sideman 做好了回到 \text{Gliese}Gliese 星球的硬件准备,但是 \text{sideman}sideman 的导航系统还没有完全设计好。为了方便起见,我们可以认为宇宙是一张有 NN 个顶点和 MM 条边的带权无向图,顶点表示各个星系,两个星系之间有边就表示两个星系之间可以直航,而边权则是航行的危险程度。
\text{sideman}sideman 现在想把危险程度降到最小,具体地来说,就是对于若干个询问 (A, B)(A,B),\text{sideman}sideman 想知道从顶点 AA 航行到顶点 BB 所经过的最危险的边的危险程度值最小可能是多少。作为 \text{sideman}sideman 的同学,你们要帮助 \text{sideman}sideman 返回家园,兼享受安全美妙的宇宙航行。所以这个任务就交给你了。
输入输出格式
输入格式:
第一行包含两个正整数 NN 和 MM,表示点数和边数。
之后 MM 行,每行三个整数 AA,BB 和 LL,表示顶点 AA 和 BB 之间有一条边长为 LL 的边。顶点从 11 开始标号。
下面一行包含一个正整数 QQ,表示询问的数目。
之后 QQ 行,每行两个整数 AA 和 BB,表示询问 AA 和 BB 之间最危险的边危险程度的可能最小值。
输出格式:
对于每个询问, 在单独的一行内输出结果。如果两个顶点之间不可达, 输出 \text{impossible}impossible。
#include <bits/stdc++.h> using namespace std; const int maxn = 1e6 + 5; int head[maxn], top[maxn], father[maxn]; int son[maxn], size[maxn], depth[maxn]; int fa[maxn], val[maxn], L[maxn], R[maxn], vis[maxn]; int n, m, q, u, v, tot, cnt, Index; struct Node { int u, v, w; bool operator < (const Node &x) const { return w < x.w; } }a[maxn]; struct Edge { int u, v, next; }edge[maxn]; int find(int x) { return x == fa[x] ? x : fa[x] = find(fa[x]); } void init() { memset(head, -1, sizeof(head)); tot = 1; } void add(int u, int v) { edge[++tot].u = u; edge[tot].v = v; edge[tot].next = head[u]; head[u] = tot; } void dfs1(int u, int fa) { size[u] = 1; son[u] = 0; vis[u] = 1; father[u] = fa; depth[u] = depth[fa] + 1; int maxson = -1; for (int i = head[u]; i != -1; i = edge[i].next) { int to = edge[i].v; if (to == fa) { continue; } dfs1(to, u); size[u] += size[to]; if (size[to] > maxson) { maxson = size[to]; son[u] = to; } } } void dfs2(int u, int topf) { top[u] = topf; L[u] = R[u] = ++Index; if (son[u]) { dfs2(son[u], topf); } for (int i = head[u]; i != -1; i = edge[i].next) { int to = edge[i].v; if (top[to]) { continue; } dfs2(to, to); } R[u] = Index; } int LCA(int x, int y) { while (top[x] != top[y]) { if (depth[top[x]] < depth[top[y]]) { swap(x, y); } x = father[top[x]]; } if (depth[x] > depth[y]) { return y; } return x; } void kruskal() { for (int i = 1; i <= n; i++) { fa[i] = i; } sort(a + 1, a + m + 1); for (int i = 1; i <= m; i++) { int fu = find(a[i].u); int fv = find(a[i].v); if (fu != fv) { ++cnt; fa[cnt] = fa[fu] = fa[fv] = cnt; val[cnt] = a[i].w; add(cnt, fu); add(fu, cnt); add(cnt, fv); add(fv, cnt); } } for(int i=1;i<=cnt;++i) { if(!vis[i]){ int f=find(i); dfs1(f,0); dfs2(f,f); } } } int main() { ios::sync_with_stdio(false); cin.tie(0); cout.tie(0); init(); cin >> n >> m; cnt = n; for (int i = 1; i <= m; i++) { cin >> a[i].u >> a[i].v >> a[i].w; } kruskal(); cin >> q; while (q--) { cin >> u >> v; if (find(u) != find(v)) { cout << "impossible" << endl; continue; } cout << val[LCA(u, v)] << endl; } return 0; }