Bellman-Ford算法可以完美的解决带有负权边的图,用一句话概括就是:对所有边进行n-1次”松弛“操作。他的核心代码只有四行如下:
for(int k=1;k<=n-1;k++)
for(int i=1;i<=m;i++)
if(dis[v[i]]>dis[u[i]]+w[i])
dis[v[i]]>dis[u[i]]+w[i];
上边的代码外循环共循环了n - 1 次(n为顶点的个数),内循环共循环了m次(m代表边的个数)即枚举每一条边, dis 数组是的作用和dijkstra 算法一样,也是用来记录源点到其余各个顶点的最短路径,u,v,w 三个数组用来记录边的信息。例如第i条边存储在u[i]、v[i]、w[i]中,表示从顶点u[i]到顶点v[i]这条边(u[i] --> v[i])权值为w[i]。
两层for循环的意思是:看看能否通过u[i]—>v[i] (权值为w[i])这条边,使的1号顶点的距离变短。即1号顶点到u[i]号顶点的距离(dis[u[i]]) 加上 u[i] —> v[i]这条边(权值为w[i])的值是否比原来1号顶点到v[i]号距离(dis[v[i]])要小。
注意:最短路径肯定是一个不包含回路的简单路径。回路分为正权回路(即回路权值之和为正)和负权回路(即回路权值之和为负)。我们分别来讨论一下为什么这两种回路都不可能有。如果最短路径中包含正权回路,那么去掉这个回路,一定可以得到更短的路径。如果最短路径中包含负权回路,那么肯定没有最短路径,因为每多走一次?负权回路就可以得到更短的路径。??因此,最短路径肯定是-一个不包含回路的简单路径,即最多包含n-1条边,所以进行n-1轮松弛就可以了。
//Bellman-Ford算法
#include<iostream>
using namespace std;
int main()
{
int u[50],v[50],w[50];
int dis[50],book[50];
int n,m;
//n为顶点数,m为边数
cin>>n>>m;
for(int i=1;i<=m;i++)
{
cin>>u[i]>>v[i]>>w[i];
}
//dis数组初始化
for(int i=1;i<=n;i++)
{
dis[i]=99999999;
}
dis[1]=0;
//核心代码
//只因为在含有n个顶点的图中,任意两点之间的最短路径最多包含n-1个,
//所以需要进行n-1轮的循环
//最短路径一定不包含回路的简单路径
for(int k=1;k<=n-1;k++)
{
int check=0; //标记在本轮松弛中数组dis是否会发生更新,作用是进行优化
for(int i=1;i<=m;i++) //枚举每一条边,进行一次松弛
{
if(dis[v[i]]>dis[u[i]]+w[i])
{
dis[v[i]]=dis[u[i]]+w[i];
check=1; //数组更新
}
}
if(check=0) // 松弛一次,检查dis是否更新,没有更新则提前退出
break;
}
int flag=0;
for(int i=1;i<=m;i++)
{
if(dis[v[i]]>dis[u[i]]+w[i])
flag=1;
}
if(flag==1)
{
cout<<"此图含有负权回路";
}
else
{
for(int i=1;i<=n;i++)
{
cout<<dis[i]<<" ";
}
}
return 0;
}
若含有负权回路,则不含有最短路径。(即在进行n-1轮松弛后,仍然可以继续松弛,那么此图必然存在负权回路)
不含负权回路,输出结果