Dijkstra算法虽然好,但是它不能解决带有负权边(边的权值为负数)的图。接下来
介绍一个无论是思想上还是代码实现上都堪称完美的最短路算法:Bellman-Ford。
Bellman-Ford算法非常简单,核心代码只有4行,并且可以完美地解决带有负权边的图.
思路 : 一张有向图,有n个点,m条边,用dis[]数组保存源点到各点的最短距离,可以通过对边进行n-1次的遍历,当其满足dis[v]>dis[u]+w的时候,就对其进行松弛更新,重复n-1次以后就能得到答案,如果n-1次以后还能继续更新,则可以判断图中出现了负权环,思路非常简短。
举例验算:
1.例图和边输入的顺序如下,并对dis数组进行初始化
另外实现这个算法的数据结构是边集数组.因为我们需要存储起始点,终点,权值.从而对每条边进行松弛.
参考代码
#include<iostream>
#include<queue>
using namespace std;
const int maxn = 1000;
const int INF = 0x3f3f3f3f;
int u[maxn],v[maxn],w[maxn],dis[maxn],flag;//dis:从起始点到当前点的最短路径
int n,m,x;//n:顶点数 m:边数 x:起始点
void Ford(int s){
for(int i = 1; i <= n; i++){//dis数组进行初始化
dis[i] = INF;
}
dis[s] = 0;
for(int i = 0; i < n-1;i++){//
flag = 0;//flag:用于标记此次是否所有边是否进行松弛了.
for(int j = 0; j < m; j++){
if(dis[v[j]]>=dis[u[j]]+w[j]){
dis[v[j]] = dis[u[j]]+w[j];
flag = 1;
}
}
if(!flag){//如果 这次没有进行松弛,说明所有的点 都已松弛完毕.
break;
}
}
}
int main()
{
cin>>n>>m>>x;
for(int i = 0; i<m; i++) {
cin>>u[i]>>v[i]>>w[i];
}
Ford(x);
for(int i = 1; i <= n; i++){
cout<<dis[i]<<"\t";
}
cout<<endl;
return 0;
}
运行结果:
Ford算法还可以去检测一个图是否有负权回路,如果在进行n-1轮松弛后,仍然存在:
if(dis[v[i]] > dis[u[i]] + w[i])
dis[v[i]] = dis[u[i]] + w[i];
则说明存在负环回路.
SPFA算法
思路: SPFA算法就是用队列优化过的Bellman-Ford算法,初始时将源点加入队列。每次选出队首结点,对其的所有出边进行松弛更新,更新成功的点加入队列,同一个结点可能被多次更新,但是同一个结点只能在同时在队列中出现一个,重复这个操作直到队列为空.
参考代码
#include<iostream>
#include<queue>
using namespace std;
const int maxn = 100;
const int INF = 0x3f3f3f3f;
int n,m,s,c,num;
int head[maxn],dis[maxn],vis[maxn];//vis:标记节点是否在队列中
struct node{
int next,to,w;
}e[maxn*maxn];
void add(int u,int v,int cc){
e[++num].next = head[u];
e[num].to = v;
e[num].w = cc;
head[u] = num;
}
void spfa(int u){
for(int i = 1; i <= n; i++){//dis进行初始化.
dis[i] = INF;
}
queue<int> q;
q.push(u);
vis[u] = 1;
dis[u] = 0;
while(!q.empty()){
int x = q.front();
q.pop();//怕最后忘记弹出,直接在这里处理吧..
vis[x] = 0;
for(int i = head[x]; i; i=e[i].next){
if(dis[e[i].to] > dis[x]+e[i].w){//如果该边可以松弛
dis[e[i].to] = dis[x]+e[i].w;//dis更新
if(!vis[e[i].to]) {
q.push(e[i].to);
vis[e[i].to] = 1;
}
}
}
}
}
int main()
{
int u,v,w;
cin>>n>>m>>s;
for(int i = 0; i < m; i++){
cin>>u>>v>>w;
add(u,v,w);
}
spfa(s);
for(int i = 1; i<=n;i++){
cout<<dis[i]<<"\t";
}
cout<<endl;
return 0;
}
/*
5 7 1
1 2 2
1 5 10
2 3 3
2 5 7
3 4 4
4 5 5
5 3 6
*/
运行结果:
最短路径算法分析
对于判断正环或者负环的问题,我们可以使用Ford,SPFA,或者DFS(判断能否回到曾经搜过的点).