定义
Floyd-Warshall算法(Floyd-Warshall algorithm)是解决任意两点间的最短路径的一种算法,可以正确处理有向图或负权的最短路径问题,同时也被用于计算有向图的传递闭包。Floyd-Warshall算法的时间复杂度为O(N3),空间复杂度为O(N2)。
算法思想原理
Floyd算法是一个经典的动态规划算法。用通俗的语言来描述的话,首先我们的目标是寻找从点i到点j的最短路径。从动态规划的角度看问题,我们需要为这个目标重新做一个诠释(这个诠释正是动态规划最富创造力的精华所在)
从任意节点i到任意节点j的最短路径不外乎2种可能,1是直接从i到j,2是从i经过若干个节点k到j。所以,我们假设Dis(i,j)为节点u到节点v的最短路径的距离,对于每一个节点k,我们检查Dis(i,k) + Dis(k,j) < Dis(i,j)是否成立,如果成立,证明从i到k再到j的路径比i直接到j的路径短,我们便设置Dis(i,j) = Dis(i,k) + Dis(k,j),这样一来,当我们遍历完所有节点k,Dis(i,j)中记录的便是i到j的最短路径的距离。
算法描述
- 从任意一条单边路径开始。所有两点之间的距离是边的权,如果两点之间没有边相连,则权为无穷大。
- 对于每一对顶点 u 和 v,看看是否存在一个顶点 w 使得从 u 到 w 再到 v 比己知的路径更短。如果是更新它。
- Floyd算法过程矩阵的计算
用arr[i][j]来表示i->j的距离,若没有连接即为无穷。
从i到j时:i可以到k,然后从k到j,即arr[i][k]+arr[k][j]<arr[i][j]时,更新arr[i][j];
用path[i][j]标记i到j走过的路径。
#include <iostream>
#include <limits.h>
using namespace std;
int n;
const int MAXN = 1000;
int arr[MAXN][MAXN];
int path[MAXN][MAXN];
int floyd() //Floyd算法
{
for (int i = 1; i <= n; i++)//经过i到达其他地方
{
//更新arr[j][k]的值
for (int j = 1; j <= n; j++)
{
for (int k = 1; k <= n; k++)
{
//从j可以到i,然后从i到k,即arr[j][i]+arr[i][k]<arr[j][k]时,更新arr[j][k];
if (arr[j][k] > arr[j][i] + arr[i][k])
{
arr[j][k] = arr[j][i] + arr[i][k];
//标记经过的路径
path[j][k] = i;
}
}
}
}
return 0;
}
void Ppath( int i, int j)
{
int k;
k = path[i][j];
if (k == -1)
{
return;
}
Ppath(i, k);
printf("%d->", k);
Ppath(k, j);
}
int main()
{
memset(path, -1, sizeof(path));
cout << "请输入矩阵的阶数:";
cin >> n;
cout << "请输入矩阵:(若无法从i到j,请输入-1):" << endl;
for (int i = 1; i <= n; i++)
{
for (int j = 1; j <= n; j++)
{
cin >> arr[i][j];
if (arr[i][j] == -1)
{
arr[i][j] = 0x3ffff;//表示最大值
}
}
}
//Floyd算法
floyd();
//输出
cout << "距离矩阵的最终结果:" << endl;
for (int i = 1; i <= n; i++)
{
for (int j = 1; j <= n; j++)
{
printf("%6d", arr[i][j]);
}
cout << endl;
}
cout << "路径矩阵的最终结果:" << endl;
for (int i = 1; i <= n; i++)
{
for (int j = 1; j <= n; j++)
{
printf("%6d", path[i][j]);
}
cout << endl;
}
cout << endl << "请输入起始地址和最终地址:";
int x, y;
cin >> x >> y;
printf("%d->", x);
Ppath(x, y);
printf("%d\n", y);
system("pause");
return 0;
}
运行结果如下: