标签:
Time Limit: 1000MS | Memory Limit: 30000K | |
Total Submissions: 22123 | Accepted: 7990 |
Description
Input
Output
Sample Input
3 2 1 20.0 1 2 1.00 1.00 1.00 1.00 2 3 1.10 1.00 1.10 1.00
Sample Output
YES
Source
解析
题意:
有多种汇币,汇币之间可以交换,这需要手续费,当你用100A币交换B币时,A到B的汇率是29.75,手续费是0.39,那么你可以得到(100 - 0.39) * 29.75 = 2963.3975 B币。问s币的金额经过交换最终得到的s币金额数能否增加
货币的交换是可以重复多次的,所以我们需要找出是否存在正权回路,且最后得到的s金额是增加的
怎么找正权回路呢?(正权回路:在这一回路上,顶点的权值能不断增加即能一直进行松弛)
分析:
一种货币就是一个点
一个“兑换点”就是图上两种货币之间的一个兑换方式,是双边,但A到B的汇率和手续费可能与B到A的汇率和手续费不同。
唯一值得注意的是权值,当拥有货币A的数量为V时,A到A的权值为K,即没有兑换
而A到B的权值为(V-Cab)*Rab
本题是“求最大路径”,之所以被归类为“求最小路径”是因为本题题恰恰与bellman-Ford算法的松弛条件相反,求的是能无限松弛的最大正权路径,但是依然能够利用bellman-Ford的思想去解题。
因此初始化dis(S)=V 而源点到其他点的距离(权值)初始化为无穷小(0),当s到其他某点的距离能不断变大时,说明存在最大路径;如果可以一直变大,说明存在正环。判断是否存在环路,用Bellman-Ford和spfa都可以。
spfa算法:
1 #include<stdio.h> 2 #include<string.h> 3 #include<iostream> 4 #include<queue> 5 #include<algorithm> 6 using namespace std; 7 double cost[105][105],rate[105][105]; 8 int n,vis[105]; 9 double v,dis[105]; 10 bool bellman_ford(int start){ 11 memset(dis,0,sizeof(dis)); 12 memset(vis,0,sizeof(vis)); 13 dis[start]=v; 14 queue<int>q; 15 q.push(start); 16 vis[start]=1; 17 while(!q.empty()){ 18 int x=q.front(); 19 q.pop(); 20 vis[x]=0; 21 for(int i=1;i<=n;i++){ 22 if(dis[i]<(dis[x]-cost[x][i])*rate[x][i]){ 23 dis[i]=(dis[x]-cost[x][i])*rate[x][i]; 24 if(dis[start]>v) 25 return true; 26 if(!vis[i]){ 27 q.push(i); 28 vis[i]=1; 29 } 30 } 31 } 32 } 33 return false; 34 } 35 int main(){ 36 int m,s; 37 while(scanf("%d%d%d%lf",&n,&m,&s,&v)!=EOF){ 38 memset(cost,0,sizeof(vis)); 39 memset(rate,0,sizeof(rate)); 40 41 for(int i=1;i<=n;i++){ 42 for(int j=1;j<=n;j++) 43 if(i==j) 44 rate[i][j]=1.0; 45 } 46 int x,y; 47 double rab,rba,cab,cba; 48 for(int i=1;i<=m;i++){ 49 cin>>x>>y>>rab>>cab>>rba>>cba; 50 cost[x][y]=cab; 51 cost[y][x]=cba; 52 rate[x][y]=rab; 53 rate[y][x]=rba; 54 } 55 if(bellman_ford(s)) 56 printf("YES\n"); 57 else printf("NO\n"); 58 } 59 return 0; 60 }
#include<stdio.h> #include<string.h> #include<iostream> #include<queue> #include<algorithm> using namespace std; double cost[105][105],rate[105][105]; int n,vis[105]; double v,dis[105]; bool bellman_ford(int start){ memset(dis,0,sizeof(dis)); memset(vis,0,sizeof(vis)); dis[start]=v; queue<int>q; q.push(start); vis[start]=1; while(!q.empty()){ int x=q.front(); q.pop(); vis[x]=0; for(int i=1;i<=n;i++){ if(dis[i]<(dis[x]-cost[x][i])*rate[x][i]){ dis[i]=(dis[x]-cost[x][i])*rate[x][i]; if(dis[start]>v) return true; if(!vis[i]){ q.push(i); vis[i]=1; } } } } return false; } int main(){ int m,s; while(scanf("%d%d%d%lf",&n,&m,&s,&v)!=EOF){ memset(cost,0,sizeof(vis)); memset(rate,0,sizeof(rate)); for(int i=1;i<=n;i++){ for(int j=1;j<=n;j++) if(i==j) rate[i][j]=1.0; } int x,y; double rab,rba,cab,cba; for(int i=1;i<=m;i++){ cin>>x>>y>>rab>>cab>>rba>>cba; cost[x][y]=cab; cost[y][x]=cba; rate[x][y]=rab; rate[y][x]=rba; } if(bellman_ford(s)) printf("YES\n"); else printf("NO\n"); } return 0; }
poj1860 bellman—ford队列优化 Currency Exchange
标签:
原文地址:http://www.cnblogs.com/13224ACMer/p/4633599.html