Dijkstra算法没办法处理带有负权的图,所以这时候就需要Bellman-Ford算法了,在假设途中没有负权回路(回路的权值和为负,即回路中负权的值大于其他几遍的权值和)可以采用Bellman-Ford算法,但是该算法的时间复杂度为O(V*E),效率较低
1.初始化:将除源点外的所有顶点的最短距离估计值 d[v] ——>+∞, d[s]——>0;
2.迭代求解:反复对边集E中的每条边进行松弛操作,使得顶点集V中的每个顶点v的最短距离估计值逐步逼近其最短距离;(运行|v|-1次)
3.检验负权回路:判断边集E中的每一条边的两个端点是否收敛。如果存在未收敛的顶点,则算法返回false,表明问题无解;否则算法返回true,并且从源点可达的顶点v的最短距离保存在 d[v]中。
#include <iostream>
#include <stdio.h>
using namespace std;
const int maxnum=100; //最大边数
const int maxint=9999; //源点和某点不可达时的距离
//有向边的结构体
typedef struct Edge
{
int Start; //有向边边的起始点
int End; //有向边的终点
int Weight;//边的权重
}Edge;
Edge edge[maxnum];//有向边的数组
int dist[maxnum]; //距离数组
//节点的数目、边的数目、源节点的下标
int nodenum,edgenum,source;
//初始化函数
void Init()
{
cin>>nodenum>>edgenum>>source;
for(int i=1;i<=nodenum;i++)
dist[i]=maxint;
dist[source]=0;
for(int i=1;i<=edgenum;i++)
{
cin >> edge[i].Start >> edge[i].End >> edge[i].Weight;
if(source==edge[i].Start)
{
dist[edge[i].End]=edge[i].Weight;
}
}
}
//松弛函数
void Relax(int Start,int End,int Weight)
{
if(dist[End]>dist[Start]+Weight)
dist[End]=dist[Start]+Weight;
}
//贝尔曼福特函数
bool Bellman_Ford()
{
for(int i=1;i<=nodenum-1;i++)
{
for(int j=1;j<=edgenum;j++)
{
Relax(edge[j].Start,edge[j].End,edge[j].Weight);
}
}
bool flag=1;
for(int i=1;i<=edgenum;i++)
{
//判断是否存在负回路
if(dist[edge[i].End]>dist[edge[i].Start]+edge[i].Weight)
{
flag=0;
break;
}
}
return flag;
}
int main()
{
freopen("out.txt","r",stdin);//打开txt
Init();
if(Bellman_Ford())
{
for(int i=1;i<nodenum;i++)
cout<<dist[i]<<endl;//打印源节点到每个节点的距离
}
return 0;
}