[洛谷P1119]灾后重建

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

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

输入格式：
输入文件rebuild.in的第一行包含两个正整数N，M，表示了村庄的数目与公路的数量。

第二行包含N个非负整数t[0], t[1], …, t[N – 1]，表示了每个村庄重建完成的时间，数据保证了t[0] ≤ t[1] ≤ … ≤ t[N – 1]。

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

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

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

输出格式：

输出文件rebuild.out包含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
-1
5
4

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

对于30%的数据，有t[i] = 0，其中有20%的数据有t[i] = 0且N＞50；

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

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

for k=1 -> n    for i = 1 -> n        for j = 1 -> n            f[i][j] = min(f[i][j],f[i][k]+f[k][j])

PS：上述所示的K为i->j的路径中点
题解 ： 对于这类点数很少的问题，我们第一个想到的就是floyd算法，对于floyd算法每枚举的一个路径中点计算最短路之前，我们都可以算出询问时间在当前枚举点重建时间之前的所有询问的答案，这里感谢作者，把询问的时间设成不下降的。那么就直接这么做了。

#include<iostream>#include<cstring>#include<algorithm>#include<cstdio>#include<cstdlib>#include<vector>#define LiangJiaJun mainusing namespace std;int n,m,t[204],Q;bool vis[204];int mp[204][204],INF,x[50004],y[50004],t2[50004];void output(int exford){     if(mp[x[exford]][y[exford]]==INF)puts("-1");     else printf("%d\n",mp[x[exford]][y[exford]]);}int LiangJiaJun(){    scanf("%d%d",&n,&m);    for(int i=0;i<n;i++)scanf("%d",&t[i]);    memset(mp,127/3,sizeof(mp));    INF =mp[1][1];    for(int i=0;i<n;i++)mp[i][i]=0;    for(int i=1;i<=m;i++){        int u,v,w;        scanf("%d%d%d",&u,&v,&w);        mp[u][v] = mp[v][u] = w;    }    scanf("%d",&Q);    for(int s=1;s<=Q;s++)scanf("%d%d%d",&x[s],&y[s],&t2[s]);    int now=1;    for(int k=0;k<n;k++){        while(t2[now] < t[k]&&now <= Q){            if(t2[now] < t[x[now]] || t2[now] < t[y[now]])puts("-1");            else output(now);            ++now;        }        for(int i=0;i<n;i++)            for(int j=0;j<n;j++)                mp[i][j]=min(mp[i][j],mp[i][k]+mp[k][j]);    }    while(now <= Q){        if(t2[now] < t[x[now]] || t2[now] < t[y[now]])puts("-1");        else output(now);        ++now;    }    return 0;}