codeup 问题 B: 算法7-16：弗洛伊德最短路径算法

                                 问题 B: 算法7-16：弗洛伊德最短路径算法

                                                               时间限制: 1 Sec  内存限制: 32 MB
                                                                          提交: 107  解决: 66
                                                              [提交][状态][讨论版][命题人:外部导入]

题目描述

在带权有向图G中，求G中的任意一对顶点间的最短路径问题，也是十分常见的一种问题。解决这个问题的一个方法是执行n次迪杰斯特拉算法，这样就可以求出每一对顶点间的最短路径，执行的时间复杂度为O(n3)。而另一种算法是由弗洛伊德提出的，时间复杂度同样是O(n3)，但算法的形式简单很多。可以将弗洛伊德算法描述如下：

在本题中，读入一个有向图的带权邻接矩阵（即数组表示），建立有向图并按照以上描述中的算法求出每一对顶点间的最短路径长度。

 

输入

输入的第一行包含1个正整数n，表示图中共有n个顶点。其中n不超过50。

以后的n行中每行有n个用空格隔开的整数。对于第i行的第j个整数，如果大于0，则表示第i个顶点有指向第j个顶点的有向边，且权值为对应的整数值；如果这个整数为0，则表示没有i指向j的有向边。当i和j相等的时候，保证对应的整数为0。

输出

共有n行，每行有n个整数，表示源点至每一个顶点的最短路径长度。如果不存在从源点至相应顶点的路径，输出-1。对于某个顶点到其本身的最短路径长度，输出0。

请在每个整数后输出一个空格，并请注意行尾输出换行。

样例输入

4
0 3 0 1
0 0 4 0
2 0 0 0
0 0 1 0

样例输出

0 3 2 1 
6 0 4 7 
2 5 0 3 
3 6 1 0 

提示

 

在本题中，需要按照题目描述中的算法完成弗洛伊德算法，并在计算最短路径的过程中将每个顶点是否可达记录下来，直到求出每一对顶点的最短路径之后，算法才能够结束。

相对于迪杰斯特拉算法，弗洛伊德算法的形式更为简单。通过一个三重循环，弗洛伊德算法可以方便的求出每一对顶点间的最短距离。

另外需要注意的是，为了更方便的表示顶点间的不可达状态，可以使用一个十分大的值作为标记。而在题目描述中的算法示例使用了另外一个三维数组对其进行表示，这使原本的O(n3)时间复杂度增长到了O(n4)，这也是需要自行修改的部分。

解题思路：此题是全源最短路径问题，用之前的dijkstra算法解决的话要执行n次，虽然到最后复杂度也还是O(n^3)，但是人家弗洛伊德算法简单啊，不仅思路简单，代码还简单。。。。正因为弗洛伊德算法思想及其简单，相当于暴力一遍，所以有个很大的限制就是顶点个数n被限制到200以内。

弗洛伊德算法的思想：如果存在一个顶点k，使得以k作为中介点时可以使顶点 i 到顶点 j 的距离缩短，则使用 k 作为中介点，即 if(d[i][k]+d[k][j]<d[i][j]) d[i][j]=d[i][k]+d[k][j]；

弗洛伊德算法的具体实现：

#define maxv 50
#define inf 100000000
int d[maxv][maxv]//d 除了i->i自身的距离初始化为0，其余全部初始化为inf,表示不可达

void Floyd()
{
    for(int k=0;k<n;k++)//枚举顶点k
    {
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)//以k作为中介点，看是否能优化顶点 i 到 j 的距离
            {
                if(d[i][k]!=inf&&d[k][j]!=inf&&d[i][k]+d[k][j]<d[i][j])
                    d[i][j]=d[i][k]+d[k][j];
            }
        }
    }
}

AC代码：

#include <iostream>
#define maxv 51
#define inf 100000000
using namespace std;
int d[maxv][maxv];
int n,tem;

void Floyd()//3个for循环是弗洛伊德算法的核心部分，这个算法的思想非常简单，相当于暴力了一遍
{
    for(int k=0;k<n;k++)
    {
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                if(d[i][k]!=inf&&d[k][j]!=inf&&d[i][k]+d[k][j]<d[i][j])
                    d[i][j]=d[i][k]+d[k][j];
            }
        }
    }
}

int main()
{
    fill(d[0],d[0]+maxv*maxv,inf);

    cin>>n;
    for(int i=0;i<n;i++)
        d[i][i]=0;

    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
        {
            cin>>tem;
            if(tem==0) continue;
            else
            {
                d[i][j]=tem;
            }
        }
    }
    Floyd();

    for(int i=0;i<n;i++)//请注意此题的输出格式要求，这可以算是一个坑吧
    {
        for(int j=0;j<n;j++)
        {
            if(i==j&&j==n-1) cout<<0<<endl;
            else if(i==j&&j!=n-1) cout<<0<<" ";
            else if(d[i][j]==inf&&j==n-1) cout<<-1<<endl;
            else if(d[i][j]==inf&&j!=n-1) cout<<-1<<" ";
            else if(d[i][j]!=inf&&j!=n-1) cout<<d[i][j]<<" ";
            else if(d[i][j]!=inf&&j==n-1) cout<<d[i][j]<<endl;
        }
    }
    return 0;
}