【洛谷】 P1119 灾后重建 最短路floyd算法

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

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

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

第二行包含NN个非负整数t_0, t_1,…, t_{N-1}t
0
​
,t
1
​
,…,t
N−1
​
，表示了每个村庄重建完成的时间，数据保证了t_0 ≤ t_1 ≤ … ≤ t_{N-1}t
0
​
≤t
1
​
≤…≤t
N−1
​
。

接下来MM行，每行33个非负整数i, j, wi,j,w，ww为不超过1000010000的正整数，表示了有一条连接村庄ii与村庄jj的道路，长度为ww，保证i≠ji

​
=j，且对于任意一对村庄只会存在一条道路。

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

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

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

输入输出样例
输入 #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
-1
5
4
说明/提示
对于30%30%的数据，有N≤50N≤50；

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

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

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

题意就不多说明了。

思路（floyd算法）：

标签有最小生成树，没想到怎么用最小生成树写。
但是最短路的提示很明显了，数量级又这么小（n<=200），大好的机会这不用floyd用啥？
但是注意到有t的限制，简单来说
每次询问，询问的两个端点必须在询问时间t之前就已经重建好了。有的人觉得限制就只是这样。但是我们在用floyd做题时，中转站k必须也是重建好的，也就是说每个k要在询问时间t之前就建好了。
看似有点麻烦，但是题目给了简化条件——原Ti单调不减，询问时间也单调不减。这有什么作用呢？我可以在每次询问时，遍历k且T[k] <= t（询问时间），然后保持k在每次询问中传递下去，因为询问时间也是单调的，所以我的k不用每次都从1开始，前面询问时的k作中转已经帮我们计算过了。
大体就是这样，然后主体就是floyd的板子了。

#include <iostream>
#include <vector>
#include <cstdio>
#include <map>
#include <string>
#include <cmath>
#include <cstring>
#include <stack>
#include <queue>
#include <algorithm>
#define maxn 500
using namespace std;
int static inf = 0x3f3f3f3f;
int t[maxn];
int dp[maxn][maxn];
int n, m, q;

void init()
{
    for(int i=1; i<=n; i++)
        for(int j=1; j<=n; j++)
            if(i==j) dp[i][j] = 0;
            else dp[i][j] = inf;
}

int main()
{
    scanf("%d%d",&n,&m);
    init();
    for(int i = 1; i<=n; i++)
     {
        scanf("%d",&t[i]);
     }
    for(int i=1; i<=m; i++)
    {
        int x,y,dis;
        scanf("%d%d%d",&x,&y,&dis);
        x++,y++;
        dp[x][y] = dp[y][x] = dis;
    }

    scanf("%d",&q);
    int k = 1;
    for(int z=1; z<=q; z++)
    {
        int x,y,ti;
        scanf("%d%d%d",&x,&y,&ti);
        x++,y++;
        int tmp = max(x,y);
        if(t[tmp]<=ti)
        {
                for(; t[k]<=ti&&k<=n; k++)
                    for(int i=1; i<=n; i++)
                        for(int j=1; j<=n; j++)
                        {
                            if(dp[i][k]!=inf&&dp[k][j]!=inf&&dp[i][k] + dp[k][j] < dp[i][j])
                            {
                                dp[i][j] = dp[i][k] + dp[k][j];
                            }
                        }
            if(dp[x][y]!=inf)
            printf("%d\n",dp[x][y]);
            else
            printf("-1\n");
        }
        else
            printf("-1\n");
    }
    return 0;
}