//bellman-ford 算法
//判断有无环
//对于每条路进行更新,最大的环为n-1,所以最多进行n-1次更新(其实就是以s为原点的一个最短路的树)
//那么对于第n次,还能进行更新那么就能存在环
//模板:
dis[s] = 0;
MEM(dis,inf);
for(int i=1;i<n;i++){
int flag = 0; // 优化
for(int j=1;j<=m;j++){
if(dis[edge[j].v] > dis[edge[j].u] + edge[j].cost){
dis[edge[j].v] = dis[edge[j].u] + edge[j].cost
flag = 1;
}
}
if(!flag){
break;//如果不能再更新了,那么就直接跳出,优化
}
}
//判断是否有环
for(int j=1;j<=m;j++){
if(dis[edge[j].v] > dis[edge[j].u] + edge[j].cost){
return true;//还能更新,有环
}
}
return false;
bellman-ford 判断是针对图,整个图可以有多个联通分块
这道题要找正环
同理可得,把每次更新的条件换一下就行
代码如下
<pre name="code" class="cpp">#include <cstdio>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <climits>
#include <string>
#include <vector>
#include <cmath>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <sstream>
#include <cctype>
using namespace std;
typedef long long ll;
typedef pair<int ,int> pii;
#define MEM(a,b) memset(a,b,sizeof a)
#define CLR(a) memset(a,0,sizeof a);
const int inf = 0x3f3f3f3f;
const int MOD = 1e9 + 7;
struct Node
{
int v;
int u;
double rate;
double com;
};
double dis[200];
Node edge[400];
double get(double v,double rate, double com){
return (v-com)*rate;
}
int main()
{
#ifdef LOCAL
freopen("in.txt", "r", stdin);
// freopen("out.txt","w",stdout);
#endif
int n,m,s;
double v;
while(cin >> n >>m >> s >> v){
int cnt = 0;
for(int i=1;i<=m;i++){
int u,v;
cin >> u >> v;
Node node;node.v = v,node.u = u;
cin >> node.rate >> node.com;
edge[++cnt] = node;
node.v = u,node.u = v;
cin >> node.rate >> node.com;
edge[++cnt] = node;
}
CLR(dis);
dis[s] = v;
int ans = 0;
for(int i=1;i<n;i++){
int flag = 0;
for(int j=1;j<=m*2;j++){
double ans = get(dis[edge[j].u],edge[j].rate,edge[j].com);
if(dis[edge[j].v] < ans){
dis[edge[j].v] = ans;
flag = 1;
}
}
if(!flag)break;
}
for(int j=1;j<=m*2;j++){
if(dis[edge[j].v] < get(dis[edge[j].u],edge[j].rate,edge[j].com))ans = 1;
}
if(ans){
printf("YES\n");
}
else printf("NO\n");
}
return 0;
}