Bellman-Ford算法

Bellman-Ford

1.基本思路

for n-1:
   for 所有边(a, b, w)边 //随便存，如定义一个结构体
      dist[b] = min[dist[b], dist[a] + w]

循环后一定满足：对于边ab，有dist[b] <= dist[a] + w(a,b)
(三角不等式 松弛操作)
定理
最外层的第k次迭代得到的距离dist[b]为从起点到顶点b所经过的弧数不超过k的最短路的路长。
注：
该定理展示了dist的实际意义。可以看到dist随着迭代的增加是非增的，且对于无负圈(圈的总权重小于0)的图，所有最短有向路中弧数一定不超过n-1（n-1条弧的路是Hamilton-路），所以Bellman-Ford算法一定在第n-1次迭代后收敛(最短有向路的路长确定)。
若迭代第n次时，dist的值发生变动，则说明存在一条长度为n的最短路到达顶点n，该路有n+1个顶点，则由抽屉原理知有两个点重合，即存在圈，且为负圈，继续迭代也将无法收敛，因此该算法也可判断负圈的存在。
适用范围
Bellman-Ford是一种采用逐次逼近的方法进行求解的算法，可以解决权重有负数但是没有负有向圈的最短路问题，或者解决有边数限制的最短路问题。

2.例题

有边数限制的最短路
给定一个 n 个点 m 条边的有向图，图中可能存在重边和自环， 边权可能为负数。
请你求出从 1 号点到 n 号点的最多经过 k 条边的最短距离，如果无法从 1 号点走到 n 号点，输出 impossible。
注意：图中可能存在负权回路 。
输入格式
第一行包含三个整数 n,m,k。
接下来 m 行，每行包含三个整数 x,y,z，表示存在一条从点 x 到点 y 的有向边，边长为 z。
输出格式
输出一个整数，表示从 1 号点到 n 号点的最多经过 k 条边的最短距离。
如果不存在满足条件的路径，则输出 impossible。
数据范围
1≤n,k≤500,
1≤m≤10000,
任意边长的绝对值不超过 10000。
输入样例：

3 3 1
1 2 1
2 3 1
1 3 3

输出样例：

3

思路
存在负权边，考虑用Bellman-Ford算法。
用结构体存储边，在每次循环里遍历所有条边来更新dist即可，为保证不串联更新dist(串联更新则不满足题中对路长的限制)，则用backup数组备份上一次的dist。
时间复杂度
$O(km)$
代码

#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

const int N = 510, M = 10010;

int n, m, k;
int dist[N], backup[N];

struct Edge
{
    int a, b, w;
}edges[M];

void bellman_ford()
{
    memset(dist, 0x3f, sizeof dist);
    dist[1] = 0;
    
    for(int i = 0; i < k; i ++) // 弧数不超过k的最短路
    {
        memcpy(backup, dist, sizeof dist); // 备份上次的dist
        for(int j = 0; j < m; j ++) // 遍历每条边
        {
            int a = edges[j].a, b = edges[j].b, w = edges[j].w;
            dist[b] = min(dist[b], backup[a] + w); 
        }
    }
    // 若最终不到达，且终点的内邻居与其有负权边，则无穷会更新，即减小
}

int main()
{
    scanf("%d%d%d", &n, &m, &k);
    
    for(int i = 0; i < m; i ++)
    {
        int a, b, w;
        scanf("%d%d%d", &a, &b, &w);
        edges[i] = {a, b, w};
    }
    
    bellman_ford();
    
    if(dist[n] > 0x3f3f3f3f / 2) puts("impossible");
    else printf("%d\n", dist[n]);
    
    return 0;
}

tips
为何无穷常取0x3f3f3f3f?
参考关于0x3f3f3f3f(0x四个3f)
主要是满足无穷+常数=无穷、无穷+无穷=无穷、防止溢出变负数，在最短路松弛操作(min)时减少出错；
以及0x3f3f3f3f的十进制是1061109567，也就是10^9级别的，基本满足无穷大；
memset(数组,0x3f,sizeof(数组))置无穷大很方便。
（负无穷大：0xcfcfcfcf）。