链接: POJ - 1860 Currency Exchange
题意:题目是说给出几种货币的兑换比率和兑换手续费,现有其中一种货币 s ,其数量为 v ,问能不能经过一些兑换之后换到更多货币s,即兑换后数量大于 v。货币兑换的计算为:(当前货币 – 手续费)X 兑换率。注意 A 到 B 和 B 到 A 间的兑换比率和手续费是不同的。
可以看一下输入数据,第一行为 4 个数,n , m , s , v,n 为货币类型数量,m 为下面有几种货币间的兑换渠道,接下来有 m 行,每一行有 6 个数,x , y , r1, c1, r2 , c2 , 代表货币 x 和货币 y 之间的兑换,货币 x 换成 y 需要收 c1 手续费,兑换率为 r1,货币 y 换成 x 需要收 c2 手续费,兑换率为 r2.
解法:可以把题目抽象为一个图,每一种货币就是一个顶点,而货币间的兑换就是一条带权有向边,权值可正可负。题目所求就是在这个有向图上从顶点 s 出发,初始距离为 v, 经过一些路径和顶点后回到顶点 s 时,距离大于 v。那么,只要我们能在 s 所能到达的路径中,找到一个正环,经过这个环可以让手中所持货币数量不断增加,那么,就一定能够通过不断地走这个环获得更多货币从而使回到 s 时,货币数量大于 v.
所以,我们可以很容易联系到负权图的最短路算法 Bellman-flod 或者 SPFA。但这两个算法是找负环的,需要做一定的修改,变成求正环。我们以 SPFA 为例,修改的地方有两处,一个是初始化,各个顶点的距离初始值改为 0,顶点 s 的距离为 v, 另一个是根据当前顶点的边来更新距离值的时候,应该用 > 来判断,即兑换后的货币大于原来货币时进行更新。
还有一个需要注意的是精度问题。货币间的兑换率和手续费是浮点数,所以浮点数的比较要特殊处理,可以用 le-6 或 le-7 来作为精度就可以了。如定义 eps = le-7 ,那么 a > b 应该写成 a – b > eps。 POJ 不考虑精度就可以过,Ural 必须处理精度问题。
#include <iostream>
#include <cstring>
#include <string>
#include <queue>
#include <vector>
#include <cstdio>
using namespace std;
#define maxn 1003
#define INF 0x3f3f3f3f
#define eps 1e-7
struct Edge
{
int x,y;
double r,c; // x 到 y 的换算率,c 为手续费
Edge(int xx,int yy,double rr,double cc)
{
x = xx;
y = yy;
r = rr;
c = cc;
}
};
vector<Edge> ways[maxn]; //存储边
double dis[maxn]; //顶点的距离
int cnt[maxn]; //顶点进队次数
bool ever[maxn]; //顶点是否在队列中
int n,m,s;
double v;
void init()
{
memset(ways,0,sizeof(ways));
memset(cnt,0,sizeof(cnt));
memset(ever,0,sizeof(ever));
for(int i = 1;i <= n;++i)
{
dis[i] = 0;
}
for(int i = 0;i < m;++i)
{
int t1,t2;
double r1,c1,r2,c2;
cin>>t1>>t2;
cin>>r1>>c1>>r2>>c2;
ways[t1].push_back(Edge(t1,t2,r1,c1));
ways[t2].push_back(Edge(t2,t1,r2,c2));
}
}
bool SPFA()
{
ever[s] = 1;
cnt[s]++;
queue<int> q;
q.push(s);
dis[s] = v;
while(!q.empty())
{
int cur = q.front();
q.pop();
int len = ways[cur].size();
ever[cur] = 0;
for(int i = 0;i < len;++i)
{
int y = ways[cur][i].y; //y 为边的下一个顶点
//根据汇率换算公式进行计算,判断换算后的数量是否大于原有数量,此处注意精度
if((dis[cur] - ways[cur][i].c)*ways[cur][i].r - dis[y] > eps)
{
dis[y] = (dis[cur] - ways[cur][i].c)*ways[cur][i].r;
if(!ever[y])
{
q.push(y);
ever[y] = 1;
cnt[y]++;
//进队次数大于 n-1 则判断有正环
if(cnt[y] >= n) return true;
}
}
}
}
return false;
}
int main()
{
cin>>n>>m>>s>>v;
init();
bool flag = SPFA();
if(flag)
printf("YES");
else
printf("NO");
return 0;
}