Floyd最短路径

无负权值最短路径问题的Floyd算法。

二维数组d表示有向图边信息，d[i][j]表示从结点i到结点j的最短路径长度，若i==j，则d[i][j]=0，若i，j之间无直接可达路径，则d[i][j]取值无穷大（为计算方便，取值1000000）

二维数组p表示路径走向信息，若结点i直接可达结点j，则p[i][j]=j，否则p[i][j]=0

计算如下有向图中任意两点的最短路径长度和走向：

define H 1000000

d的初始值：d = {
{0,7,5,12,H,H},
{H,0,H,3,H,H},
{H,H,0,6,H,15},
{12,H,6,0,H,8},
{6,H,13,H,0,H},
{H,4,15,H,9,0}}

p的初始值：p = {
{1,2,3,4,0,0},
{0,2,0,4,0,0},
{0,0,3,4,0,6},
{1,0,3,4,0,6},
{1,0,3,0,5,0},
{0,2,3,0,5,6}};

Floyd算法步骤：

1、边信息矩阵d，路径信息矩阵p
2、i = 1,j = 1,k = 1
3、u = d[i][k]+d[k][j],若u < d[i][j],d[i][j] = u,p[i][j] = p[i][k];否则转到4
4、若j < n，则j = j + 1,转到3；否则转到5
5、若i < n, 则i = i + 1,转到3；否则转到6
6、若k < n, 则k = k + 1,i = 1,j = 1,转到3，否则end

其实就是一个三重循环。代码如下：

#include<iostream>
using namespace std;
#define H 100000
void fun(int d[6][6],int p[6][6])
{
	int i,j,k;
	for(k=0;k<6;k++)
	{
		for(i=0;i<6;i++)
		{
			for(j=0;j<6;j++)
			{
				if(d[i][k]+d[k][j]<d[i][j])
				{
					d[i][j]=d[i][k]+d[k][j];
					p[i][j]=p[i][k];
				}
			}
		}
	}
}
int main()
{
	int d[6][6] = {
		{0,7,5,12,H,H},
		{H,0,H,3,H,H},
		{H,H,0,6,H,15},
		{12,H,6,0,H,8},
		{6,H,13,H,0,H},
		{H,4,15,H,9,0}};
	int p[6][6] = {
		{1,2,3,4,0,0},
		{0,2,0,4,0,0},
		{0,0,3,4,0,6},
		{1,0,3,4,0,6},
		{1,0,3,0,5,0},
		{0,2,3,0,5,6}};
	fun(d,p);
	int i,j;
	for(i=0;i<6;i++)
	{
		for(j=0;j<6;j++)
			cout<<d[i][j]<<' ';
		cout<<endl;
	}
	cout<<endl;
	for(i=0;i<6;i++)
	{
		for(j=0;j<6;j++)
			cout<<p[i][j]<<' ';
		cout<<endl;
	}
	return 0;
}

结果截图：

此时边信息矩阵d中d[i][j]即表示从结点i到结点j的最大路径长度。而路径信息矩阵p记录路径走向。

比如查找结点1到结点5的最短路径走向：

第1行第5列，值为2，表示经过结点2；第2行第5列，值为4，表示经过结点4；第4行第5列，值为6，表示经过结点6，第6行第5列，值为5结束。

得到路径走向：1->2->4->6->5

Floyd算法，可以得到所以结点之间的最短路径长度及走向（若结点之间可达）。