第六章 最短路径

Floyd-Warshall

多源最短路径问题：求任意两点之间的最短路径。

可以通过深度和广度优先搜索求出两点之间的最短路径：进行n^2遍深度或广度优先搜索，即对每两个点都进行一次深度或广度优先搜索，便可以求得任意两点之间的最短路径。

如果要让任意两点（例如a->b）之间的路程变短，只能引入第三个点（例如k），并通过这个顶点k中转即a->k->b，才可能缩短原来从a到b的路程。

甚至有时候不只通过一个点，而是经过两个点或者更多点中转会更短。

每个顶点都有可能使得另外两个顶点之间的路程变短。

当任意两点之间不允许经过第三个点时，这些点之间的最短路程就是初始路程。

假如现在只允许经过1号顶点中转，求任意两点的最短路程。只需判断e[i][1]+e[1][j]是否比e[i][j]要小即可。其中，e[i][j]表示从i到j顶点之间的路程，e[i][1]+e[1][j]表示的是从i先到1，再从1到j的路程之和。

for(i = 1; i <= n; ++i)
    for(j = 1; j <= n; ++j)
        if(e[i][j] > e[i][1] + e[1][j])
            e[i][j] = e[i][1] + e[1][j];

接下来只允许经过1和2两个顶点的情况下任意两点之间的最短路程。只需在只允许经过1时任意两点的最短路径的结果下，再判断如果经过2是否可以使得i到j之间的路程变得更短，即判断e[i][2]+e[2][j]是否比e[i][j]更短？

//经过1
for(i = 1; i <= n; ++i)
    for(j = 1; j <= n; ++j)
        if(e[i][j] > e[i][1] + e[1][j])
            e[i][j] = e[i][1] + e[1][j];

//经过2
for(i = 1; i <= n; ++i)
    for(j = 1; j <= n; ++j)
        if(e[i][j] > e[i][2] + e[2][j])
            e[i][j] = e[i][2] + e[2][j];

综上所述，基本思想为：最开始只允许经过1进行中转，接下来只允许经过1和2进行中转……允许经过1~n号所有顶点进行中转，求任意两点之间的最短路程。这是一种“动态规划”思想，从i到j只经过前k号点的最短路程。

for(k = 1; k <= n; ++k)
    for(i = 1; i <= n; ++i)
        for(j = 1; j <= n; ++j)
            if(e[i][j] > e[i][k] + e[k][j])
                e[i][j] = e[i][k] + e[k][j];

示例

#include<iostream>
#include<vector>
#include<climits>
#include<iomanip>

using namespace std;

int main(void)
{
    int v = 4;  //顶点数
    int e = 8;  //边数

    vector<vector<int>> graph(v, vector<int>(v, INT_MAX));
    for (int i = 0; i < v; ++i)
        graph[i][i] = 0;
    graph[0][1] = 2;
    graph[0][2] = 6;
    graph[0][3] = 4;
    graph[1][2] = 3;
    graph[2][0] = 7;
    graph[2][3] = 1;
    graph[3][0] = 5;
    graph[3][2] = 12;
    /*
    0  2  6  4
    ∞ 0  3  ∞
    7  ∞ 0  1
    5  ∞ 12 0
    */
    //1~k作为中转点
    for (int k = 0; k < v; ++k)
    {
        for (int i = 0; i < v; ++i)
        {
            for (int j = 0; j < v; ++j)
            {