P1119 灾后重建 Floyd算法本质

题目背景

B 地区在地震过后，所有村庄都造成了一定的损毁，而这场地震却没对公路造成什么影响。但是在村庄重建好之前，所有与未重建完成的村庄的公路均无法通车。换句话说，只有连接着两个重建完成的村庄的公路才能通车，只能到达重建完成的村庄。

题目描述

给出 B 地区的村庄数 N，村庄编号从 0 到 N-1，和所有 MM条公路的长度，公路是双向的。并给出第 i个村庄重建完成的时间 t​，你可以认为是同时开始重建并在第 t 天重建完成，并且在当天即可通车。若 t​ 为 00 则说明地震未对此地区造成损坏，一开始就可以通车。之后有 Q 个询问 (x,y,t)，对于每个询问你要回答在第 t 天，从村庄 x 到村庄 y 的最短路径长度为多少。如果无法找到从 x 村庄到 y 村庄的路径，经过若干个已重建完成的村庄，或者村庄 x 或村庄 y 在第 t 天仍未重建完成，则需要返回 -1。

输入格式

第一行包含两个正整数N,M，表示了村庄的数目与公路的数量。

第二行包含N个非负整数t0, t1,…, tN-1，表示了每个村庄重建完成的时间，数据保证了t0 ≤ t1 ≤ … ≤ tN-1。

接下来M行，每行3个非负整数i, j, w，w为不超过10000的正整数，表示了有一条连接村庄i与村庄j的道路，长度为w，保证i≠j，且对于任意一对村庄只会存在一条道路。

接下来一行也就是M+3行包含一个正整数Q，表示Q个询问。

接下来Q行，每行33个非负整数x, y, t，询问在第t天，从村庄x到村庄y的最短路径长度为多少，数据保证了t是不下降的。

输出格式

共Q行，对每一个询问(x, y, t)输出对应的答案，即在第t天，从村庄x到村庄y的最短路径长度为多少。如果在第t天无法找到从x村庄到y村庄的路径，经过若干个已重建完成的村庄，或者村庄x或村庄y在第t天仍未修复完成，则输出−1。

输入输出样例

输入 

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

输出 

-1
-1
5
4

说明/提示

对于30%的数据，有N≤50N≤50；

对于30%的数据，有i​=0，其中有%20%的数据有t_i = 0ti​=0且N>50N>50；

对于50%的数据，有Q≤100Q≤100；

对于100%的数据，有N≤200N≤200，M≤N×(N−1)/2，Q≤50000Q≤50000，所有输入数据涉及整数均不超过100000100000。

思路

我们很多时候都只是背下代码，而忘了去思考一个算法的思想本质是什么，而Floyd算法的本质是什么，我们先看它的代码

//互相影响，互相更新
for (int k = 1; k <= G.n; k++)
    for (int i = 1; i <= G.n; i++)
        for (int j = 1; j <= G.n; j++)
        {
            if ((*D)[i][j] > (*D)[i][k] + (*D)[k][j]) //如果经过下标为k的顶点比原两点路径更短
            {
                (*D)[i][j] = (*D)[i][k] + (*D)[k][j]; //将当前两点权值设置为更小的那个
                (*P)[i][j] = (*P)[i][k];              //路径设置为经过下标为k的顶点
            }
        }

关于这个代码，我们可以把经过下标为k的顶点作为中转点，k从1开始，表示若把1作为中转点，任意两点间的距离会不会更短，会的话就更新；然后k=2，表示1和2都可以作为中转点，再次判断任意两点间的距离会不会更短，依次类推最后k=n，表示n个点都可以作为中转点，最后一次更新后的任意两点间的距离就是最优最短路径，换句话说这个图是一点点生成的，每添加一个中转点就重新更新一次最短路径

总结出来就是：不断的动态更新i号顶点到j号顶点只能经过前k个点的距离，到最后的时候就是最短路径

而对题意我们可以把它理解为：题目按照时间顺序给出村庄的修好时间（题目给的数据已经是排好序的），我们需要查询在指定天数时，村庄A和村庄B的最短路径。那么我们只能把能在指定天数前修好的村庄添加到图中作为中转点，其他的不添加，那么村庄的修好顺序就是k的顺序。

我们直接来看代码

AC CODE

#include <bits/stdc++.h>
using namespace std;
#define inf 1e7 + 5
int n, m, vis[205], a[205][205], shortpath[205][205], vex[205];

void floyd(int k)
{
    for (int i = 0; i < n; i++)
        for (int j = 0; j < n; j++)
        {
            if (shortpath[i][j] > shortpath[i][k] + shortpath[k][j])
                shortpath[i][j] = shortpath[j][i] = shortpath[i][k] + shortpath[k][j];
        }
}
int main()
{
    cin >> n >> m;
    //初始化
    for (int i = 0; i < n; i++)
        cin >> vex[i];
    for (int i = 0; i < n; i++)
        for (int j = 0; j < n; j++)
        {
            if (i == j)
                a[i][j] = 0;
            else
                a[i][j] = inf;
        }

    for (int k = 0; k < m; k++)
    {
        int i, j, w;
        cin >> i >> j >> w;
        a[i][j] = a[j][i] = w;
    }
    for (int i = 0; i < n; i++)
        for (int j = 0; j < n; j++)
            shortpath[i][j] = a[i][j];

    //Floyd算法 核心code
    int q;
    cin >> q;   
    int num = 0;
    while (q--)
    {
        int x, y, t;
        cin >> x >> y >> t;
        while (vex[num] <= t && num < n)   //把修好天数小于t的添加到图中作为中转点，并且不能超出范围
        {
            floyd(num);
            num++;
        }
        if (vex[x] > t || vex[y] > t) //如果x或y本身都没有修完，也就是都没有加入floyd算法，就表示没有路径
            cout << -1 << endl;
        else
        {
            if (shortpath[x][y] == inf) //x和y虽然都修好了，但是不连通无法到达，也表示没有路径
                cout << -1 << endl;
            else
                cout << shortpath[x][y] << endl;
        }
    }
    return 0;
}

一定要好好思考算法的本质！！！