bellman-ford算法 时间复杂度O(nm)

在一般情况下，spfa算法都优于bf算法，但遇到最短路的边数有限制的题时，只能用bf算法

bf算法和dijkstra很像

#include<iostream>
#include<queue>
#include<cstring>
#include<algorithm>using namespace std;const int N = 510,M = 10010;int n,m,k;
int dist[N],backup[N]; //backup备份数组struct Edge{int a,b,w;
}Edge[M]; //存所有边int bellman_ford(){memset(dist,0x3f,sizeof dist);dist[1] = 0;for(int i = 0;i < k;i ++ ){memcpy(backup,dist,sizeof dist); //备份dist，不会出现串联情况for(int j = 0;j < m;j ++ ){int a = Edge[j].a,b = Edge[j].b,w = Edge[j].w;dist[b] = min(dist[b],backup[a] + w);}}if(dist[n] > 0x3f3f3f3f / 2) return 0;else return dist[n];
}int main(){cin >> n >> m >> k;for(int i = 0;i < m;i ++ ){int a,b,w;cin >> a >> b >> w;Edge[i] = {a,b,w};}int t = bellman_ford();if(!t) cout << "impossible" << endl;else cout << t << endl;return 0;
}

spfa算法 时间复杂度一般O(m), 最坏O(nm)

基本上单源最短路都可以用spfa来解决

spfa的核心优化思路是：拿我更新过的点来更新别人。一个点如果没有被更新过的话，拿它来更新别人一定是没有效果的，只有该点变小了，该点后面的点才会变小

spfa代码和堆优化dijkstra特别像

spfa算法求最短路

#include<iostream>
#include<queue>
#include<cstring>
#include<algorithm>using namespace std;typedef pair<int,int> PII;const int N = 150010;int n,m;
int h[N],w[N],e[N],ne[N],idx;
int dist[N];
bool vis[N];int add(int a,int b,int c){e[idx] = b,w[idx] = c,ne[idx] = h[a],h[a] = idx,idx ++ ;
}int spfa(){memset(dist,0x3f,sizeof dist);dist[1] = 0;queue<int> q; //队列里存的是变小的aq.push(1);vis[1] = true;while(q.size()){int t = q.front();q.pop();vis[t] = false;for(int i = h[t];i != -1;i = ne[i]){int j = e[i];if(dist[j] > dist[t] + w[i]){dist[j] = dist[t] + w[i];if(!vis[j]){q.push(j);vis[j] = true;}}}}if (dist[n] == 0x3f3f3f3f) return 0;return dist[n];
}int main(){cin >> n >> m;memset(h,-1,sizeof h);while(m -- ){int x,y,z;cin >> x >> y >> z;add(x,y,z);}int t = spfa();if(!t) cout << "impossible" << endl;else cout << t << endl;return 0;
}

spfa算法求负环

spfa算法可以求出负环用的是抽屉原理，即把多于n+1个的物体放到n个抽屉里，则至少有一个抽屉里的东西不少于两件。

代码在spfa求最短路的模板上稍加改动即可

#include<iostream>
#include<queue>
#include<cstring>
#include<algorithm>using namespace std;typedef pair<int,int> PII;const int N = 150010;int n,m;
int h[N],w[N],e[N],ne[N],idx;
int dist[N],cnt[N];
bool vis[N];int add(int a,int b,int c){e[idx] = b,w[idx] = c,ne[idx] = h[a],h[a] = idx,idx ++ ;
}bool spfa(){queue<int> q;for(int i = 1;i <= n;i ++ ){ //由于存在的负环1号点可能走不到，所以要把每一个点都推进队列vis[i] = true;q.push(i);}vis[1] = true;while(q.size()){int t = q.front();q.pop();vis[t] = false;for(int i = h[t];i != -1;i = ne[i]){int j = e[i];if(dist[j] > dist[t] + w[i]){dist[j] = dist[t] + w[i];cnt[j] = cnt[t] + 1; //最重要的一步，如果j被更新了最短路，那么意味着j点的cnt是前一个点t+1条边达到的if(cnt[j] >= n) return true;if(!vis[j]){q.push(j);vis[j] = true;}}}}return false;
}int main(){cin >> n >> m;memset(h,-1,sizeof h);while(m -- ){int x,y,z;cin >> x >> y >> z;add(x,y,z);}if(spfa()) cout << "Yes" << endl;else cout << "No" << endl;return 0;
}