P1119 灾后重建

导言

一道非常好的Floyd最短路练习题，从算法本质上出的题目，对于初学Floyd算法的人来说是绝佳的练习题

解题难度：普及+/提高

重要程度：提高+/省选-

算法导入

关于floyed

首先我们从Floyd算法谈起，这是一个看上去很简单的算法-事实上也的确很简单，整个算法一共只有五行，三重循环＋一个判断就能求出图中任意两点之间的最短路径。

很多人都是直接把这个算法背了下来，然后要求最短路时直接默写，没有考虑到这个算法的本质意义，而这个题目，正好考了这个算法的本质，如果只会背的话，遇见这题直接就GG。

这个算法的主要思路，就是通过其他的点进行中转来求的两点之间的最短路。因为我们知道，两点之间有多条路，如果换一条路可以缩短距离的话，就更新最短距离。而它最本质的思想，就是用其他的点进行中转，从而达到求出最短路的目的。

floyed中转

那么，如何进行中转呢？两点之间可以由一个点作为中转点更新最短路径，也可以通过多个点更新最短路径。

结合代码：

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];
//核心代码，仅仅只有5行

　　

这段代码的基本思想就是：最开始只允许经过1号顶点进行中转，接下来只允许经过1和2号顶点进行中转……允许经过1~n号所有顶点进行中转，求任意两点之间的最短路程。用一句话概括就是：从i号顶点到j号顶点只经过前k号点的最短路程。

（仔细理解这段话，它揭露了这个算法的本质并为本题提供了很好的方法）

到这里我们已经知道，Floyd算法就是一个利用其它点进行中转来求最短路的步骤。

而我们再回头看题意：

所有的边全部给出，按照时间顺序更新每一个可用的点（即修建好村庄），对于每个时间点进行两点之间询问，求对于目前建设的所有村庄来说任意两点之间的最短路

不正好就是Floyd算法中使用前k个节点更新最短路的思维吗？

落实代码

于是到了这里，我们差不多也就知道这题如何写了。

出题人还是很良心的，保证所有的数据都是用时间顺序给出的，所以我们只用读取+操作就可以了，不用再储存＋排序。

总体思路：

int main(){
    读入，存下每个村庄修复的时间
    读入所有的边并使用邻接矩阵存图
    初始化
    对于每次询问，将在当前时间前的所有点全部用Floyd更新。
    特殊判断并输出
}

读入

cin>>n>>m;
    for(int i=0;i<n;i++)
    scanf("%d",a+i);//依次输入每一个村庄建立完成时需要的时间
 /*******************初始化*************************/
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
    f[i][j]=1e9;//初始化为一个很大的值
    for(int i=0;i<n;i++)
    f[i][i]=0;//一个点到自己的距离为0
 /*******************读入边长*************************/
    int s1,s2,s3;
    for(int i=1;i<=m;i++){
        scanf("%d%d%d",&s1,&s2,&s3);
        f[s1][s2]=f[s2][s1]=s3;//无向边，存两次
    }

　　

处理

//s1,s2,s3为临时变量，表示两村庄和时间
    int q;cin>>q;//一共q个询问
    int now=0;//用于记录当前时间对应的编号，从第一个村庄开始
    for(int i=1;i<=q;i++){//处理各询问 
        scanf("%d%d%d",&s1,&s2,&s3);
        while(a[now]<=s3&&now<n){//如果目前更新的点的时间在询问点之前
            updata(now);
            now++;
        }//处理在它之前建立的村庄
        if(a[s1]>s3||a[s2]>s3)cout<<-1<<endl;//村庄未建好
        else {
            if(f[s1][s2]==1e9)cout<<-1<<endl;//两点不连通
            else cout<<f[s1][s2]<<endl;//输出路径长
        }
    }

 核心

inline void updata(int k){
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
    if(f[i][j]>f[i][k]+f[j][k])
    f[i][j]=f[j][i]=f[i][k]+f[j][k];//用这个新的更新所有前面的 
    return;
}

 最后

return 0
    //没错，完了，没了，就两个步骤加一个函数.

　　

标准的Floyd思想，内层循环不变，外层循环按照k的顺序给出。

这真的是一个很优秀的题目。同时它也提醒我们，学习算法时要弄懂各个细节，把整个算法的思想掌握，而不是仅仅把代码背下来.

链接floyed详解https://www.cnblogs.com/GumpYan/p/5540549.html

(这篇文章对Floyd算法进行了图文并茂的很详细的讲解，如果想更详细地了解此算法可以去看看）

代码

#include<iostream>
#include<cstdio>
#define N 205
using namespace std;
int n,m;
int a[N];
int f[N][N];//邻接矩阵存边
inline void updata(int k){
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++)
    if(f[i][j]>f[i][k]+f[j][k])
    f[i][j]=f[j][i]=f[i][k]+f[j][k];//用这个新的更新所有前面的 
    return;
}
int main(){
    cin>>n>>m;
    for(int i=0;i<n;i++)
    scanf("%d",a+i);//依次输入每一个村庄建立完成时需要的时间
    for(int i=0;i<n;i++)
    for(int j=0;j<n;j++){
        f[i][j]=1e9;//初始化为保证它不爆炸范围内的最大值 
    }
    for(int i=0;i<n;i++)
    f[i][i]=0;
    int s1,s2,s3;
    for(int i=1;i<=m;i++){
        scanf("%d%d%d",&s1,&s2,&s3);
        f[s1][s2]=f[s2][s1]=s3;//初始化边长 
    }
    int q;
    cin>>q;
    int now=0;
    for(int i=1;i<=q;i++){//处理各询问 
        scanf("%d%d%d",&s1,&s2,&s3);
        while(a[now]<=s3&&now<n){
            updata(now);//依次更新点，使它可以被用来更新其他的点 
            now++;
        }
        if(a[s1]>s3||a[s2]>s3)cout<<-1<<endl;
        else {
            if(f[s1][s2]==1e9)cout<<-1<<endl;
            else cout<<f[s1][s2]<<endl;
        }
    }
    return 0;
}