E - Currency Exchange POJ - 1860

bellman-ford队列优化 反向运用bellman-ford 求正权回路

题解：（摘自kuangbin的博客）贴一发超链接

点击打开链接

有多种汇币，汇币之间可以交换，这需要手续费，当你用100A币
交换B币时，A到B的汇率是29.75，手续费是0.39，那么你可以得到
(100 - 0.39) * 29.75 = 2963.3975 B币。问s币的金额经过交换最终
得到的s币金额数能否增加
货币的交换是可以重复多次的，所以我们需要找出是否存在
正权回路，且最后得到的s金额是增加的
怎么找正权回路呢？（正权回路：在这一回路上，顶点的权值能不断增加即能一直进行松弛）

关键在于反向利用Bellman-Ford算法


单源最短路径算法，因为题目可能存在负边，所以用Bellman Ford算法，
原始Bellman Ford可以用来求负环，这题需要改进一下用来求正环

一种货币就是图上的一个点
一个“兑换点”就是图上两种货币之间的一个兑换环，相当于“兑换方式”M的个数，是双边
唯一值得注意的是权值，当拥有货币A的数量为V时，A到A的权值为K，即没有兑换
而A到B的权值为(V-Cab)*Rab
本题是“求最大路径”，之所以被归类为“求最小路径”是因为本题题恰恰
与bellman-Ford算法的松弛条件相反，求的是能无限松弛的最大正权路径，
但是依然能够利用bellman-Ford的思想去解题。
因此初始化d(S)=V   而源点到其他店的距离（权值）初始化为无穷小（0），
当s到其他某点的距离能不断变大时，说明存在最大路径

code:

#include <iostream>   
#include <algorithm>   
#include <cstdio>   
#include <cstring>  
#include<cmath>  
#include<cstdlib>  
#include<map>  
#include<string>  
#include<vector>  
#include<queue>  
#include<functional>  
#pragma warning(disable:4996)  
//最大边的最小权值  
using namespace std;
typedef long long ll;
const int INF = 999999999;
typedef pair<int, int> P;//第一个值表示最大权值  
struct edge {
	int to;
	double cost,r, c;
};
int mn;
double d[1001];
vector<edge>G[1005];
int inq[1005];
int cnt[1005];
int n, m, s;
double v;
bool bellmanford(int s)
{
	queue<int>que;
	memset(inq, 0, sizeof(inq));
	memset(cnt, 0, sizeof(cnt));
	for (int i = 1; i <= n; i++)
		d[i] = 0;

	d[s] = v;
	inq[s] = 1;
	que.push(s);
	while (!que.empty())
	{
		int u = que.front();
		que.pop();
		inq[u] = 0;
		for (int i = 0; i < G[u].size(); i++)
		{
			edge e = G[u][i];
			if (d[u] > 0 && d[e.to] < (d[u] - e.c)*e.r)
			{
				d[e.to] = (d[u] - e.c)*e.r;
				if (!inq[e.to])
				{
					que.push(e.to);
					inq[e.to] = 1;
					if (++cnt[e.to] > n)
						return true;
				}
			}
		}

	}

	return false;
}
int main()
{
	cin >> n >>  m >> s >> v;
	for (int i = 1; i <= n; i++)G[i].clear();
	for (int i = 1; i <= m; i++)
	{
		int x, y;
		double rab, cab, rba, cba;
		cin >> x >> y >> rab >> cab >> rba >> cba;
		edge e;
		e.to = y, e.r = rab,e.c=cab;
		G[x].push_back(e);
		e.to = x, e.r = rba, e.c = cba;
		G[y].push_back(e);

	}
	if (bellmanford(s))
		printf("YES\n");
	else
		printf("NO\n");
}