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]
更短?
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];
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;
for (int k = 0; k < v; ++k)
{
for (int i = 0; i < v; ++i)
{
for (int j = 0; j < v; ++j)
{