我们熟悉的算法dijkstra可用于求单源点到其他各点的最短距离。但是有一个问题在于,图中所有边的权之都应该为正值。因为如果有负权值的边,则只要绕着这条负值边一直走,最算路径很可能越来越短。而这正是bellman ford算法的思想:如果路径可以缩短,那么就按照dijkstra的算法让路径进行N-1次(N为图中的顶点个数)缩短。如果N-1次缩短之后仍能继续缩短,那么图中必定存在负边。
算法的思想参见维基百科里关于bellman Ford算法的解释:
把每个货币币种看作图的一个顶点,把兑换过程看作顶点之间的动态的有向边(注意,是动态的有向边,因为货币数量与兑换完的数量之间是一个函数关系,而不是固定的)。那么:
poj1860这道题,则刚好把“缩短”这个概念换为“货币能增值”,如果经过某个环,让货币增值了,那么就可以让主人公的货币一直增长下去。进行N-1次“增值”之后,如果还能继续增值,则说明图中存在能无限把货币增值的环。
poj1860代码:
#include<iostream>
#include<fstream>
using namespace std;
int N,M,S;
int all;
double V;
double dis[202];
struct Edge
{
int v1;
int v2;
double rat;
double com;
}edges[404];
bool bellmanford()
{
for(int i=0;i<N-1;i++)
{
bool flag=false;
for(int j=0;j<all;j++)
{
int v1=edges[j].v1;
int v2=edges[j].v2;
if((dis[v1]-edges[j].com)*edges[j].rat>dis[v2])
{
dis[v2]=(dis[v1]-edges[j].com)*edges[j].rat;
flag=true;
}
}
if(flag==false)
break;
}
for(int j=0;j<all;j++)
{
int v1=edges[j].v1;
int v2=edges[j].v2;
if((dis[v1]-edges[j].com)*edges[j].rat>dis[v2])
return true;
}
return false;
}
int main()
{
while(cin>>N>>M>>S>>V)
{
all=0;
for(int i=0;i<M;i++)
{
cin>>edges[all].v1>>edges[all].v2>>edges[all].rat>>edges[all].com;
edges[all+1].v1=edges[all].v2;
edges[all+1].v2=edges[all].v1;
all++;
cin>>edges[all].rat>>edges[all].com;
all++;
}
for(int i=1;i<=N;i++)
{
dis[i]=0;
}
dis[S]=V;
if(bellmanford())
{
cout<<"YES"<<endl;
}
else
cout<<"NO"<<endl;
}
return 0;
}