浅谈SPFA判负环

浅谈(SPFA)判负环

1. 关于(SPFA)：它死了

​ (SPFA(Shortest Path Faster Algorithm))是一种优秀的单源最短路算法。它可以在(O(km))的时间复杂度内求出单源最短路。

​ 但同时，(SPFA)有一个令人诟病的点，就是可以被特殊数据卡为(O(n^2))，所以我们一般会选择更稳定的堆优化(Dijkstra)。但(Dijkstra)无法处理负权边与负环，这时候我们就需要采用(SPFA)

---

2.(SPFA)判负环的原理

​ (SPFA)的本质是利用队列，重复松弛节点，来求最短路径。它更像是一个能够使节点重复入队的(BFS)。而(BFS)的搜索树的深度最大不会超过节点数。而(SPFA)同理。

​ 设(cnt[i])表示从(1)到(x)的最短路包含的边数。当执行更新(dist[y]=dist[x]+z)时，同时更新(cnt[y]=cnt[x]+1)，此时若发现(cnt[y]>=n)，则该图中有负环。

​ 另一种方式是记录每个点入队的次数，次数到达n时说明有负环。这种判定方式效率一般不如上面的高。

---

代码如下：

洛谷(P3385)

#include<bits/stdc++.h>
using namespace std;
#define maxn 60100
int head[maxn],nxt[maxn],to[maxn],tot=0;
int val[maxn];
void add(int x,int y,int z)
{
	to[++tot]=y;
	nxt[tot]=head[x];
	val[tot]=z;
	head[x]=tot;
}
void init()
{
	memset(head,0,sizeof(head));
	memset(nxt,0,sizeof(nxt));
	memset(to,0,sizeof(to));
	tot=0;
	memset(val,0,sizeof(val));
}
int dis[maxn],inq[maxn],cnt[maxn];
int n,m;
bool spfa(int s)
{
	queue<int> q;
	q.push(s);
	memset(inq,0,sizeof(inq));
	memset(cnt,0,sizeof(cnt));
	inq[s]=1;
	memset(dis,0x3f,sizeof(dis));
	dis[s]=0;
	while(!q.empty())
	{
		int x=q.front();
		q.pop();
		inq[x]=0;
		for(int i=head[x];i;i=nxt[i])
		{
			int y=to[i];
			if(dis[y]>dis[x]+val[i])
			{
				dis[y]=dis[x]+val[i];
				cnt[y]=cnt[x]+1;
				if(cnt[y]>n+1)
				return 0;
				if(!inq[y])
				{
					q.push(y);
					inq[y]=1;
				}	
				
			}
		}
	}
	return 1;
}
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)
	{
		init();
		scanf("%d%d",&n,&m);
		for(int i=1;i<=m;i++)
		{
			int x,y,z;
			scanf("%d%d%d",&x,&y,&z);
			add(x,y,z);
			if(z>=0)
			add(y,x,z);
		}
		
		if(spfa(1))
		printf("NO
");
		else 
		printf("YES
");
		
	}	
	return 0;
}