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");
else
puts("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);
else
q.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");
else
puts("NO");
}
return 0;
}