POJ 3259 Wormholes 邻接表的SPFA判断负权回路

http://poj.org/problem?id=3259

题目大意：

一个农民有农场，上面有一些虫洞和路，走虫洞可以回到 T秒前，而路就和平常的一样啦，需要花费时间走过。问该农民可不可能从某个点出发后回到该点，并且于出发时间之前?

思路：

就是让我们判断存不存在一条总权值未负的回路。

你想想，如果我们一直走这个回路，就可以无限的回到过去。（ (╯‵□′)╯︵┻━┻世上哪有时空隧道。想回到过去么？如果可以？ 你会去做什么，挽回让自己后悔的事情？ ，咳咳，继续题解。）

那就用SPFA呗，如果一个点入队次数超过n，那么就存在回路。（因为SPFA求最短路，如果存在环，那么会一直入队。。出。。，也就是说环的是没有最小值的，而正环只会走一次。）

SPFA常见的优化有SLF和LLL

SLF（Small Label First）是指在入队时如果当前点的dist值小于队首， 则插入到队首， 否则插入到队尾。

未优化的版本：

157MS

#include<cstdio>#include<string>#include<queue>using namespace std;const int INF=9999999;const int MAXN=520;const int MAXM=5200;struct edge{int to;int val;int next;}e[MAXM];int len,head[MAXN];int dis[MAXN];int n,m,w;bool SPFA(){for(int i=1;i<=n;i++)dis[i]=INF;bool vis[MAXN]={0};int cnt[MAXN]={0};int cur=1;queue<int> q;q.push(cur);vis[cur]=true;cnt[cur]=1;dis[cur]=0;while(!q.empty()){cur=q.front();q.pop();vis[cur]=false;for(int i=head[cur]  ;i!=-1; i=e[i].next){int id=e[i].to;if( dis[cur] + e[i].val < dis[ id ] ){dis[ id ] = dis[cur] + e[ i ].val;if(!vis[id]){cnt[id]++;vis[id]=true;q.push(id);if(cnt[cur]>n)return true;}}}}return false;}void add(int from,int to,int val){e[len].to=to;e[len].val=val;e[len].next= head[from];head[from]=len++;}int main(){int T;scanf("%d",&T);while(T--){memset(head,-1,sizeof(head));len=0;scanf("%d%d%d",&n,&m,&w);for(int i=0;i<m;i++){int from,to,val;scanf("%d%d%d",&from,&to,&val);add(from,to,val);add(to,from,val); //双向的}for(int i=0;i<w;i++){int from,to,val;scanf("%d%d%d",&from,&to,&val);add(from,to,-val);}if( SPFA())puts("YES");elseputs("NO");}return 0;}

采用优化的版本：

63MS

#include<cstdio>#include<string>#include<queue>using namespace std;const int INF=9999999;const int MAXN=520;const int MAXM=5200;struct edge{int to;int val;int next;}e[MAXM];int len,head[MAXN];int dis[MAXN];int n,m,w;bool SPFA(){for(int i=1;i<=n;i++)dis[i]=INF;bool vis[MAXN]={0};int cnt[MAXN]={0};int cur=1;deque<int> q;q.push_back(cur);vis[cur]=true;cnt[cur]=1;dis[cur]=0;while(!q.empty()){cur=q.front();q.pop_front();vis[cur]=false;for(int i=head[cur]  ;i!=-1; i=e[i].next){int id=e[i].to;if( dis[cur] + e[i].val < dis[ id ] ){dis[ id ] = dis[cur] + e[ i ].val;if(!vis[id]){cnt[id]++;vis[id]=true;if(cnt[cur]>n)return true;if(!q.empty() && dis[id] < dis[q.front()])q.push_front(id);elseq.push_back(id);}}}}return false;}void add(int from,int to,int val){e[len].to=to;e[len].val=val;e[len].next= head[from];head[from]=len++;}int main(){int T;scanf("%d",&T);while(T--){memset(head,-1,sizeof(head));len=0;scanf("%d%d%d",&n,&m,&w);for(int i=0;i<m;i++){int from,to,val;scanf("%d%d%d",&from,&to,&val);add(from,to,val);add(to,from,val); //双向的}for(int i=0;i<w;i++){int from,to,val;scanf("%d%d%d",&from,&to,&val);add(from,to,-val);}if( SPFA())puts("YES");elseputs("NO");}return 0;}