Time Limit: 1 second
Memory Limit: 128 MB
【问题描述】
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。
【数据规模】
对于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。
Sample Input1
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
Sample Output1
-1 -1 5 4
【题解】
这道题要利用floyd的原理来做。 我们floyd的3层循环如下。 for (int k=1;k<= n;k++) for (int i = 1;i <= n;i++) for (int j=1;j <= n;j++) if (w[i][j] > w[i][k]+w[k][j]) f[i][j] = w[i][k]+w[k][j]; 这里的k层循环(第一层)枚举的是经过哪一些点作为中间点来缩短i->j的距离。 可以理解为用了前k个点作为中间点来尝试更新任意两点之间的距离。 这点可以为我们所利用。 看一下我们的询问。 x,y,time 如果t[x] 或者t[y]>time则肯定是输出-1的。 对于其他的 我们可以在k层循环中的i,j循环完毕之后加上下面这些东西。 即 for (k=1->n) { for (i=1->n) for (j =1->n) ... 在这个位置加上我们下面所说的东西; } k层循环仍是枚举n个点。 如果t[k]<=a[now].time且t[k+1] <=a[now].time 则k可以继续枚举。表示我们可以利用k和k+1来作为中间点更新任意两点之间的距离。 如果遇到t[k]<=a[now].time且t[k+1]>a[now].time。 则表示我们最多只能用k来作为中间点更新任意两点之间的距离了。 这时我们只能尝试在利用前k个点之后输出w[x][y]了。 不能再用k+1这个点了。因为k+1这个点在a[now].time时还没有修建好。 遇到这样的k之后。now++.(a[now].time是随着now的增加递增的)。 如果now递增后t[k+1]<=a[now].time了。则可以继续利用k+1来作为中间点更新任意两点之间的距离。 怎么样?理解了吧。 然后我们把k层循环的下界改为0. 因为可能有在0时刻的询问。 这个问题在脑子里想的时候有点烦。但是写下来其实还好。 然后要先将询问存下来。 村庄从0--n-1有点麻烦。直接改成1--n,即输入的时候x,y都递增1。
【代码】
#include <cstdio> #include <cstring> struct question //用结构体把询问存下来。 { int x, y, time; }; int n, m, t[201] = { 0 }, w[201][201], q; //t数组是各个节点修建好的时间。 question a[50001] = { 0 }; void input_data() { memset(w, 127 / 3, sizeof(w));//一开始w数组赋值为一个很大的数字。 scanf("%d%d", &n, &m); for (int i = 1; i <= n; i++) //输入各个节点修建好的时刻。 scanf("%d", &t[i]); for (int i = 1; i <= m; i++) //输入边权信息。 { int x, y, z; scanf("%d%d%d", &x, &y, &z); x++; y++; w[x][y] = w[y][x] = z; } scanf("%d", &q); for (int i = 1; i <= q; i++) //输入q个询问。 { scanf("%d%d%d", &a[i].x, &a[i].y, &a[i].time); a[i].x++; a[i].y++; } } void get_ans() { int now = 1; t[n + 1] = t[n] + 10000; //这是防止上溢。 for (int k = 0; k <= n; k++) //k从0开始枚举 { for (int i = 1; i <= n; i++) //以k作为中间节点尝试更新任意两点之间的距离。 for (int j = 1; j <= n; j++) if (w[i][j] > w[i][k] + w[k][j]) w[i][j] = w[i][k] + w[k][j]; while (now <= q && t[k] <= a[now].time && t[k + 1] > a[now].time) {//如果询问还没结束。且这个节点在所询问的时间内。且k+1这个节点修建的时间超过询问的时间 if (t[a[now].x] > a[now].time || t[a[now].y] > a[now].time) printf("-1 "); else //输出依靠前k个节点作为中间节点更新出的任意两点之间的距离 { if (w[a[now].x][a[now].y] >= w[0][0]) printf("-1 "); else printf("%d ", w[a[now].x][a[now].y]); } now++; //看一下下一个询问是否符合要求。 } if (now > q) //如果询问都输出了则结束。 break; } } int main() { input_data(); get_ans(); return 0; }