关于Floyd最短路的一些思考

正确性

floyd算法本质是动态规划,它的正确性不那么"显然"是因为进行了压缩,状态 $f[i][j]$ 没有体现阶段,所以阶段间的转移也很难看出来.

不过,即使不理解正确性,floyd的代码也很好记

for (int k = 1; k <= n; k++)
  for (int i = 1; i <= n; i++)
    for (int j = 1; j <= n; j++)
      if (i != j && i != k && j != k)
        f[i][j] = min(f[i][j],f[i][k]+f[k][j]);

只要记住先枚举中间点$k$, 再枚举两端点,这个程序基本就不会有错了.

为什么k不能放内层for循环?
因为它代表着阶段!
压缩前动态规划的状态 $f[i][j][k]$ 表示从 $i$ 到 $j$ 路径经过点的最大标号为 $k$ 的最短路长度.

1->2 没有经过其他点,可以看作中间点最大为0
1->4->3 有一个中间点,最大中间点4
5->1->2->9 有1,2两个中间点,最大的是2

边界:$f[i][j][0]=a[i][j]$
状态转移方程:$f[i][j][k]=min(f[i][j][k-1],f[i][k][k-1]+f[k][j][k-1])$

//f[i][j][0] = a[i][j]
for (int k = 1; j <= n; k++) {
  for (int i = 1; i <= n; i++)
    for (int j = 1; j <= n; j++)
      f[i][j][k] = f[i][j][k-1];
  for (int i = 1; i <= n; i++)
    for (int j = 1; j <= n; j++)
      f[i][j][k] = min(f[i][j][k],f[i][k][k-1]+f[k][j][k-1]);
}

并行优化

理论上,只要$f[][][k-1]$比$f[][][k]$先算就可以.剩下,先枚举$i$或者先枚举$j$都行,甚至$(i,j)$乱序更新都行,内层for循环顺序不重要,应该可以并行计算,留坑