You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

63 lines
2.6 KiB

2 years ago
#include <bits/stdc++.h>
using namespace std;
const int INF = 0x3f3f3f3f;
const int N = 1003;
2 years ago
int g[N][N]; // 邻接矩阵
int n; // n个点
int w[N]; // 额外费用
int path[N][N]; // 路径
2 years ago
void floyd() {
for (int k = 1; k <= n; k++)
for (int i = 1; i <= n; i++)
2 years ago
if (g[i][k] != INF) // floyd优化
2 years ago
for (int j = 1; j <= n; j++) {
2 years ago
if (g[i][j] > g[i][k] + g[k][j] + w[k]) { // w[k]:点权
g[i][j] = g[i][k] + g[k][j] + w[k]; // k的加入使得i->j的路径变短
2 years ago
path[i][j] = path[i][k];
}
2 years ago
if (g[i][j] == g[i][k] + g[k][j] + w[k]) { // 如果存在多条最短路径也就是除了k还有其它k1,k2使得i->j距离一样小
if (path[i][j] > path[i][k]) path[i][j] = path[i][k]; // 字典序,谁更小就留下谁
2 years ago
}
}
}
int main() {
while (cin >> n && n) {
for (int i = 1; i <= n; i++) {
for (int j = 1; j <= n; j++) {
2 years ago
path[i][j] = j; // 路径初始化,记录整条路径上离i节点最近的最短路径上的下一个点只有i->j时下一个点可不就是j
2 years ago
cin >> g[i][j];
if (g[i][j] == -1) g[i][j] = INF;
2 years ago
}
}
2 years ago
for (int i = 1; i <= n; i++) cin >> w[i]; // 读入点权
// 多源最短路径
2 years ago
floyd();
2 years ago
// 处理询问
2 years ago
int x, y;
2 years ago
while (cin >> x >> y) {
if (x == -1 && y == -1) break;
printf("From %d to %d :\n", x, y);
printf("Path: %d", x);
2 years ago
int u = x, v = y;
2 years ago
// 理解路径思路:
// (1) 从起点x出发,用循环打印路径,最后一个打印的肯定是y
// (2) 从起点x出发,第二个点应该是离x最近的并且是最短路径上的那个点,这个点就是path[x][y]!
// path[x][y]从起点x出发到终点y有多条最短路径我们选择字典序最小的那条最短路径然后path[x][y]就是从x出发离x最近的这条最短路径上的点。
2 years ago
while (x != y) {
2 years ago
printf("-->%d", path[x][y]); // 输出距离x最近的那个点
x = path[x][y]; // 更换x概念向y逼近让循环跑起来
2 years ago
}
2 years ago
puts("");
if (g[u][v] < INF)
printf("Total cost : %d\n", g[u][v]);
2 years ago
else
2 years ago
puts("-1");
puts("");
2 years ago
}
}
return 0;
}