最短路径算法——清晰简单的弗洛伊德算法(Floyd)

弗洛伊德算法(Floyd)

\qquad 上一篇文章介绍了迪杰斯特拉算法(Dijkstra)。具体请看：最短路径算法——简单明了的迪杰斯特拉算法(Dijkstra)。Dijkstra适用于非负权图，并且一次只能从网络中找源点到任何一个节点的最短路径，而Floyd算法的应用更加广泛，可以求网络中任意两点之间的最短路径，而且弗洛伊德算法适用于负权图，这篇文章就用图和表的形式来介绍一下弗洛伊德算法！

1 思想（原理）

\qquad Floyd算法可以给出网络中任意两个节点之间的最短路径，因此它是比Dijkstra更一般的算法。Floyd算法的思想是将$n$个节点的网络表示为$n$行$n$列的矩阵，而矩阵中的元素$(i,j)$表示从节点$i$到节点$j$的距离$d_{ij}$，如果两点直接没有边相连，则相应的元素就是无穷$(\infty )$.

2 步骤

<1> 第0步：定义初始距离矩阵$D_0$、节点序列矩阵$S_0$，如下表。对角线上用”—“表示不需要从自身到自身。

\qquad 这里的节点序列矩阵相当于路线表，如下表，$S_{ij}=j$表示，从节点$i$到节点$j$只需经过节点$j$即可。

令$k=1$.

<2> 一般的第k步：令第$k$行为枢轴行，第$k$列为枢轴列。对于矩阵$D_{k-1}$（上一步完成后的矩阵）中对的每一个元素做三重操作。

如果满足条件：
$$d_{ik}+d_{kj}<d_{ij}(i\ne k,j\ne k,i\ne j)$$
则进行下面的操作：

<\a> 用$d_{ik}+d_{kj}$代替矩阵$D_{k-1}$中的元素$d_{ij}$，从而得到矩阵$D_k$.
<\b> 用$k$代替矩阵$S_{k-1}$中的元素$s_{ij}$，从而得到矩阵$S_k$.
<\c> 令$k=k+1$，如果$k=n+1$，停止，否则重复<2>.

3 栗子

\qquad 直接看方法步骤会感觉太抽象，这里用一个例子进行步骤的演示.
\qquad 对下图中的网络，求任意两个节点之间的最短路径，图中弧上给出了相应节点间的距离。弧(3,5)是有向的，其他边都是双边。

迭代0：矩阵$D_0$和$S_0$代表初始的网络。可以看到矩阵$D_0$除了$d_{53}=\infty$外(因为弧(3,5)是单向弧)，$D_0$是对称的。

迭代1：令$k=1$. $D_0$矩阵中的黄色阴影表示的第1行和第1列为枢轴行和枢轴列。根据三重操作发现可以改进的元素是$d_{23}$和$d_{32}$，即

(1) $d_{21}+d_{13}=3+10<\infty$，则在$d_{23}$中用$13$代替$\infty$，并令$s_{23}=1$.
(2) $d_{31}+d_{12}=10+3<\infty$，则在$d_{32}$中用$13$代替$\infty$，并令$s_{32}=1$.

此时得到$D_1$和$S_1$，得到下表

迭代2：令$k=2$. $D_0$矩阵中的黄色阴影表示的第2行和第2列为枢轴行和枢轴列。根据三重操作发现可以改进的元素是$d_{14}$和$d_{41}$，即

(1) $d_{21}+d_{42}=3+5=8<\infty$，则在$d_{41}$中用$8$代替$\infty$，并令$s_{41}=2$.
(2) $d_{12}+d_{24}=3+5=8<\infty$，则在$d_{14}$中用$8$代替$\infty$，并令$s_{14}=2$.

此时得到$D_2$和$S_2$，得到下表

迭代3：令$k=3$. $D_0$矩阵中的黄色阴影表示的第3行和第3列为枢轴行和枢轴列。根据三重操作发现可以改进的元素是$d_{15}$和$d_{25}$，即

(1) $d_{13}+d_{35}=10+15=25<\infty$，则在$d_{15}$中用$25$代替$\infty$，并令$s_{15}=3$.
(2) $d_{23}+d_{35}=13+15=28<\infty$，则在$d_{25}$中用$28$代替$\infty$，并令$s_{25}=3$.

此时得到$D_3$和$S_3$，得到下表

迭代4：令$k=4$. $D_0$矩阵中的黄色阴影表示的第4行和第4列为枢轴行和枢轴列。根据三重操作发现可以改进的元素是$d_{15}$、$d_{23}$、$d_{25}$、$d_{32}$、$d_{35}$、$d_{51}$、$d_{52}$、$d_{53}$，即

(1) $d_{14}+d_{45}=8+4=12<25$，则在$d_{15}$中用$12$代替$25$，并令$s_{15}=4$.
(2) $d_{24}+d_{43}=5+6=11<13$，则在$d_{23}$中用$11$代替$13$，并令$s_{23}=4$.
(3) $d_{24}+d_{45}=5+4=9<28$，则在$d_{25}$中用$9$代替$28$，并令$s_{25}=4$.
(4) $d_{34}+d_{42}=6+5=11<13$，则在$d_{25}$中用$11$代替$13$，并令$s_{32}=4$.
(5) $d_{34}+d_{45}=6+4=10<15$，则在$d_{35}$中用$10$代替$15$，并令$s_{35}=4$.
(6) $d_{41}+d_{54}=8+4=12<\infty$，则在$d_{51}$中用$12$代替$\infty$，并令$s_{51}=4$.
(7) $d_{42}+d_{54}=5+4=9<\infty$，则在$d_{52}$中用$9$代替$\infty$，并令$s_{52}=4$.
(8) $d_{43}+d_{54}=6+4=10<\infty$，则在$d_{53}$中用$10$代替$\infty$，并令$s_{53}=4$.

此时得到$D_4$和$S_4$，得到下表

迭代5：令$k=5$. $D_0$矩阵中的黄色阴影表示的第5行和第5列为枢轴行和枢轴列。根据三重操作发现没有可以改进的元素了。

因此，最后得到的矩阵为：

\qquad 这两个矩阵包含了网络中任意两个节点最短路径的所有信息。如从矩阵$D$中可以看出节点1到节点5的最短路径长度为12.从矩阵$S$中发现，节点1到节点5的中间节点是节点4，即节点1→节点4→节点5，再看节点1→节点4中间是节点2，即节点1需要通过节点2到达节点4，即节点1→节点2→节点4；而节点4可以直接到节点5，中间没有节点，因此可以得到节点1到节点5的最短路径是节点1→节点2→节点4→节点5.