poj 3259 最短路判负环 spfa算法和Bellman_ford算法

题目：http://poj.org/problem?id=3259

题意：

有个人在一个图上旅行，当它从一点出发后，可能会经过一些虫洞边，然后时光倒流，回到原点，看到自己，看了好久样例一才看清题意，他必须要看到“自己”，也就是要有个负环回路，所以这就是个求负环回路的问题

分析：

spfa算法和Bellman_ford算法都可以

spfs代码：

#include <iostream>
#include<cstdio>
#include<cstring>
#include<map>
#include<queue>
#include<algorithm>
using namespace std;
const int N=3000;
#define MP(x,y) make_pair(x,y)
typedef pair<int,int> pii;
vector<pii>e
;
bool vis[555];
int dis[555],cnt[555];
int n,m,w;
void init()
{
scanf("%d%d%d",&n,&m,&w);
for(int i=0;i<=n;i++)e[i].clear();
int u,v,t;
for(int i=0;i<m;i++){
scanf("%d%d%d",&u,&v,&t);
e[u].push_back(MP(v,t));
e[v].push_back(MP(u,t));
}
for(int i=0;i<w;i++){
scanf("%d%d%d",&u,&v,&t);
e[u].push_back(MP(v,-t));
}
}
bool spfa(int s)
{
memset(vis,0,sizeof(vis));
memset(cnt,0,sizeof(cnt));
memset(dis,0x3f,sizeof(dis));
queue<int>q;
dis[s]=0;
q.push(s);
vis[s]=1;
cnt[s]++;
while(!q.empty()){
int u=q.front();
q.pop();
vis[u]=0;
for(int i=0;i<e[u].size();i++){
int v=e[u][i].first,d=e[u][i].second;
if(dis[v]>dis[u]+d){
dis[v]=dis[u]+d;
if(!vis[v]){
q.push(v);
vis[v]=1;
cnt[v]++;
if(cnt[v]>N-1)return 1;
}
}
}
}
return 0;
}
int main()
{
int T;
scanf("%d",&T);
while(T--){
init();
if(spfa(1)||dis[1]<0)printf("YES\n");
else printf("NO\n");
}
return 0;
}
Bellman_ford算法：

推荐一篇很好的介绍Bellman_ford算法的文章http://www.wutianqi.com/?p=1912

#include <iostream>
#include<cstdio>
#include<cstring>
#include<map>
#include<queue>
#include<algorithm>
using namespace std;
const int N=5300;
typedef struct edge
{
int u,v,w;
};
edge e
;
int dis[555];
int n,m,w,cnt;
void init()
{
scanf("%d%d%d",&n,&m,&w);
int u,v,t;
cnt=0;
for(int i=0;i<m;i++){
scanf("%d%d%d",&u,&v,&t);
e[cnt].u=u;e[cnt].v=v;e[cnt++].w=t;
e[cnt].u=v;e[cnt].v=u;e[cnt++].w=t;
}
for(int i=0;i<w;i++){
scanf("%d%d%d",&u,&v,&t);
e[cnt].u=u;e[cnt].v=v;e[cnt++].w=-t;
}
}
bool Bellman_Ford()
{
memset(dis,0x3f,sizeof(dis));
for(int i=1;i<=n;i++){
bool flag=1;
for(int j=0;j<cnt;j++){
if(dis[e[j].v]>dis[e[j].u]+e[j].w)
dis[e[j].v]=dis[e[j].u]+e[j].w,flag=0;
}
if(flag)break;
}
for(int j=0;j<cnt;j++){
if(dis[e[j].v]>dis[e[j].u]+e[j].w)return 1;
}
return 0;
}
int main()
{
int T;
scanf("%d",&T);
while(T--){
init();
if(Bellman_Ford())printf("YES\n");
else printf("NO\n");
}
return 0;
}