求最短路径的基本方法——SPFA(贝尔曼-福特算法的优化算法)

前篇

关于Floyd算法与Dijkstra算法—>点击这里

SPFA算法

SPFA算法介绍

SPFA算法是求解单源最短路径问题的一种算法，由理查德·贝尔曼（Richard Bellman） 和 莱斯特·福特 创立的。有时候这种算法也被称为 Moore-Bellman-Ford 算法，因为 Edward F. Moore 也为这个算法的发展做出了贡献。它的原理是对图进行V-1次松弛操作，得到所有可能的最短路径。其优于迪杰斯特拉算法的方面是边的权值可以为负数、实现简单，缺点是时间复杂度过高，高达 O(VE) 其中V是节点数量，E是边数量。但算法可以进行若干种优化，提高了效率。

浅谈SPFA算法

我们用数组dis记录每个结点的最短路径估计值，用二维数组来存储图E。我们采取的方法是动态逼近法：设立一个队列用来保存待优化的结点，优化时每次取出队首结点u，并且用u点当前的最短路径估计值对离开u点所指向的结点v进行松弛操作，如果v点的最短路径估计值有所调整，且v点不在当前的队列中，就将v点放入队尾。这样不断从队列中取出结点来进行松弛操作，直至队列空为止
注意：如果题目明确说明存在负环，那么就必须用SPFA算法来解决，这时候我们必须另开一个数组用于判断i号点到j号点的最短路径的步数是否<n（n为节点数），否则每次经过负环，他的权值都会减少，这样会导致无限循环，无法算出结果

经典例题

题目来源

点击这里

思路

洛谷的模板题，只需要判断是否存在负环即可
由于是多实例，首先必须将$dis(存放权值数组)，vis(标记数组)，cnt（存放步数数组）以及存放图的数组$进行初始化
首先将顶点1存入队列，将顶点1标记为真，每次取出队列的队首，在将队首标记为false
用队首元素的最短路径估计值对离开u点所指向的结点v进行松弛操作，若当前节点大于松弛操作之后的权值，就将v点放入队尾，并且让v点在u点步数的基础上＋1，判断当前步数是否大于等于n(是否存在负环)，若满足则直接返回，不满足则继续循环

code

const int INF=INT_MAX;
const int N=1e4+5;
int dis[N],vis[N],cnt[N];
int n,m;
struct edge{
	int v,w;
};
vector<edge> e[N];
void Init(){
	for(int i=1;i<=n;++i){
		dis[i]=INF;
		vis[i]=0;
		cnt[i]=0;
		e[i].clear();
	} 
	
}
bool spfa(){
	queue<int> q;
	dis[1]=0;
	vis[1]=true;
	q.push(1);
	while(!q.empty()){
		int u=q.front();
		q.pop();
		vis[u]=false;
		for(int i=0;i<e[u].size();++i){
			int v=e[u][i].v,w=e[u][i].w;
			if(dis[v]>dis[u]+w){
				dis[v]=dis[u]+w;
				cnt[v]=cnt[u]+1;
				if(cnt[v]>=n){
					return true;
				}
				if(!vis[v]){
					q.push(v);
					vis[v]=true;
				}
			}
		}
	}
	return false;
}
void solve(){
	cin >> n >> m;
	Init();
	for(int i=1;i<=m;++i){
		int u,v,w;
		cin >> u >> v >> w;
		e[u].push_back({v,w});
		if(w>=0) e[v].push_back({u,w});
    }
	if(spfa()) cout << "YES" << endl;
	else cout << "NO" << endl;
	return ;
}