一、bellman-ford
1、算法
bellmen-ford算发,这是一个可以求负权图的算法了!要说它好吧,还算好理解,要说不好吧,复杂度是O(n*m)的!有点大。思想类似Kruskal,都是以边为中心展开。
先在这里介绍一下松弛操作:
松弛操作:
对于一条边:u->v:w,dist[v] = min( dist[v] + w, dist[v]);
将从源点(S)到v点的最短路更新成:(从S到u的最短路+u到v的路径、原最短路)之间的最小值
开始演示:
首先:设dist[i] 表示 从源点到i点的路径长度;
然后,你有一个图(拿啥建图自己决定)
然后,将源点(这里设成点1)dist更新成0,其他赋值∞;
然后,我们从右向左对每一条变进行松弛操作(为啥是从右向左,怎样的一个顺序看自己,作者的目的是从左往右没啥更新的了,一次就要出结果了,所以从右往左好看点)
之后进行最多n-1轮这样的操作,就可以结束了!
2、解释点原因
a、为什么跑n-1次就可以了
因为一共n个点,所以一共需要n-1条边就可以完成这张图的更新了。每进行一次操作,至少可以更新一条边,所以最多n-1次ok了。
b、如何判负环
当跑完n-1轮上述操作后,再跑一论,如果还能出现松弛操作,就证明一定有负环!如果含有负环,就会一直跑下去,上述的操作次数就一定不止是n-1次!这样就能证明负环了
二、代码实现
#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
const int N = 1e5+7;
const int inf = 1e8;
struct edge
{
int u,v,w;
}tu[N];//存边用
int dist[N],backup[N];//dist存距离,backup存更新前的路径
int n,m;
int bellman_ford(int s)
{
memset(dist,0x3f,sizeof(dist));//初始化距离长
dist[s] = 0;
for(int i = 0; i < n; i++)
{
memcpy(backup,dist,sizeof(dist));//复制上一个距离长
for(int j = 0; j < m; j++)
dist[tu[j].v] = min(dist[tu[j].v], backup[tu[j].u] + tu[j].w);//松弛操作
}
for(int j = 0; j < m; j++)//判负环
{
if(dist[tu[j].v] > dist[tu[j].u] + tu[j].w)
return - inf;
}
return dist[n];
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=0;i<m;i++)
{
int u,v,w;scanf("%d%d%d",&u,&v,&w);
tu[i] = {u,v,w};
}
int t = bellman_ford(1);
if(t <= -inf)printf("impossible\n");
else printf("%d\n",t);
return 0;
}