POJ 1860 bellman_ford

#include<iostream>
#include<fstream>
#include<cstdlib>

using namespace std;

//#define DEBUG

struct exchgpnt
{
	int a;
	int b;
	double rab;
	double cab;
};

static struct exchgpnt edges[100 * 2];
static int m, n, s;
static double v;
static double dst[101];

int bellman_ford(int edgenum)
{
	int i, j;
	memset(dst, 0, sizeof dst);
	dst[s] = v;
	for (i = 1; i < n; i++) /* relax edges */
	{
		for (j = 0; j < edgenum; j++)
		{
			int a, b;
			double rab, cab;
			a = edges[j].a; 
			b = edges[j].b;
			rab = edges[j].rab; 
			cab = edges[j].cab;
			if (dst[a] == 0.0)	continue;  /* 加上这句后运行时间减半 */
			if (dst[b] < (dst[a] - cab) * rab)
			{
				if (b == s) return 1; /* 存在正向回路 */
				dst[b] =  (dst[a] - cab) * rab;
			}
		}
	}
	for (j = 0; j < edgenum; j++)
	{
		int a, b;
		double rab, cab;
		a = edges[j].a; 
		b = edges[j].b;
		rab = edges[j].rab; 
		cab = edges[j].cab;
		if (dst[b] <  (dst[a] - cab) * rab)
		{
			return 1; /* 存在正向回路 */
		}
	}
	return 0;
}

int main()
{
#ifdef DEBUG
	fstream cin("G:\\book\\algorithms\\acm\\Debug\\dat.txt");
#endif
	while (cin >> n >> m >> s >> v)
	{
		int i, p;
		for (p = i = 0; i < m; i++)
		{
			int a, b;
			double rab, rba, cab, cba;
			//cin >> a >> b >> rab >> cab >> rba >> cba;
			/* 采用下面的输入方式时间减到0MS */
			scanf("%d%d%lf%lf%lf%lf", &a, &b, &rab, &cab, &rba, &cba);
			edges[p].a = a; edges[p].b = b;
			edges[p].cab = cab, edges[p++].rab = rab;

			edges[p].a = b; edges[p].b = a;
			edges[p].cab = cba, edges[p++].rab = rba;
		}
		if (bellman_ford(p))
			printf("YES\n");
		else
			printf("NO\n");
	}

	return 0;
}

n种货币看成n个顶点，货币之间的转换就是顶点之间的边。这样输入数据就构成了一张图，问在图中是否存在从s------>s的正权回路。

就是s货币经过若干次兑换之后再换回s货币，数量会增加。基本思想和最短路径中的bellman_ford算法相同。bellman_ford算法以边

为输入。通过对全部的边做Relax操作n-1次，如果图中不存在负权回路，就能够得到从源点s到其他节点的最短路径的值。如果存在负

权回路，可以发现。