算是菜鸟的第一个图的运算吧。。。还好bellman算法还好,不过说服自己算法是正确的却是花费了很久的时间。。。
如果这个计算式没有这么麻烦的话,可以把权值全部变成负值来求最小的,反过来就是最大的。
#include<iostream>
#define MAX 202
using namespace std;
typedef struct Node
{
int a;
int b;
double rab,cab;
}T;
T ver[MAX];
int n,m,s,k;
double v;
double dis[MAX/2];
bool bellman_ford()
{
int flag=0;
memset(dis,0,sizeof(dis));
dis[s]=v;//记得初始化的时候对源点进行特殊的初始化
for(int i=0;i<n-1;i++)
{
flag=0;
for(int j=0;j<k;j++)//这个是对每一条边进行的
{
if(dis[ver[j].b]<(dis[ver[j].a]-ver[j].cab)*ver[j].rab)
{
dis[ver[j].b]=(dis[ver[j].a]-ver[j].cab)*ver[j].rab;
flag=1;
}
}
if(!flag)
break;
}
for(int p=0;p<k;p++)
{
if(dis[ver[p].b]<(dis[ver[p].a]-ver[p].cab)*ver[p].rab)
return true;
}
return false;
}
int main()
{
int a,b;
double rab,cab,rba,cba;
bool is;
while(cin >> n >>m >> s >>v)
{
k=0;
while(k<2*m)//注意这里的m和k的关系
{
cin >> a >> b >> rab >> cab >> rba >> cba;
ver[k].a=a;ver[k].b=b;ver[k].rab=rab;ver[k++].cab=cab;
ver[k].a=b;ver[k].b=a;ver[k].rab=rba;ver[k++].cab=cba;
}
is=bellman_ford();
if(is)
cout << "YES" << endl;
else cout << "NO" << endl;
}
return 0;
}