TIOJ-1163

發表於 更新於 Disqus：

有一天，時弦無意間聽到了黑暗騎士與憂鬱人士的對話：
「最近家裡附近的道路又在整修了，我等了三天才有通路連到這裡，真是科科。」
「哭哭啦，你比我好多了。從我家到這裡的路蓋了五天才蓋好咧。」
時弦聽了以後覺得很有趣，於是就把他們兩人的對話對照了地圖：
map
地圖上每個點都是一個地方，而每一條邊上面所寫的數字代表第幾天以後該道路才會通。
例如要有一條從A走到D的道路需要等到第三天才能通行。
「如果我隨便指定兩個點，要如何很快地知道這兩個地點要等到第幾天才有連接的道路呢？」
時弦想了很久，決定請TDY好威幫忙。哪知道TDY好威一下子就解決了這個問題，卻不告訴時弦。
他說：「你應該自己想出來的，這個題目根本就是秒殺嘛…」
四處求助無門的時弦，決定把解決這個問題的重責大任交給你！

$V\le 30000, E\le 50000, Q\le 50000$，詢問兩點不連通輸出-1。

觀察一下可以發現，我們只需要盡可能用權重小的讓圖連通，因此就可以做MST。
用完MST後會變一座森林，對每棵樹做LCA即可知道他們路徑權重最大值。
連通可以直接用kruskal時的dsu來判，真的很好用。
實作量稍大，不過這也是必備的技能(吧
複雜度$\mathcal{O}((V+E)log(V+E))$。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
/*
--------------              |   /
      |                     |  /
      |                     | /
      |             *       |/          |    |         ------            *
      |                     |           |    |        /      \
      |             |       |\          |    |       |       |\          |
   \  |             |       | \         |    |       |       | \         |
    \ |             |       |  \        |    |        \     /   \        |
     V              |       |   \        \__/|         -----     \       |
*/
#include <bits/stdc++.h>
using namespace std;

#define EmiliaMyWife ios::sync_with_stdio(0); cin.tie(NULL);
/*-----------------------------------------------------------------------------------------------------*/

const int N = 3e4 + 25;
struct ed {
	int a, b, c;
	bool operator<(ed x) { return c < x.c; }
};
int pa[N], sz[N], dep[N], lca[15][N], mx[15][N];
bool vis[N];
vector<vector<pair<int, int>>> edge;

int fnd(int x) { return pa[x] == x ? pa[x] : pa[x] = fnd(pa[x]); }
bool uni(int a, int b) {
	if((a = fnd(a)) == (b = fnd(b)))
		return false;
	if(sz[a] > sz[b])
		swap(a, b);
	return sz[b] += sz[a], pa[a] = b, true;
}
void dfs(int v) {
	vis[v] = 1;
	for(auto [u, c] : edge[v]) {
		if(!vis[u]) {
			dep[u] = dep[v] + 1;
			lca[0][u] = v;
			mx[0][u] = c;
			dfs(u);
		}
	}
}
int solve(int a, int b) {
	if(fnd(a) != fnd(b))
		return -1;
	int res = 0;
	if(dep[a] < dep[b])
		swap(a, b);
	for(int i = 14; ~i; i--)
		if(dep[lca[i][a]] >= dep[b])
			res = max(res, mx[i][a]), a = lca[i][a];
	if(a == b)
		return res;
	for(int i = 14; ~i; i--) {
		if(lca[i][a] != lca[i][b]) {
			res = max(res, mx[i][a]); a = lca[i][a];
			res = max(res, mx[i][b]); b = lca[i][b];
		}
	}
	return max({res, mx[0][a], mx[0][b]});
}

signed main() {
	EmiliaMyWife

	int n, m;
	cin >> n >> m;
	for(int i = 1; i <= n; i++)
		pa[i] = i, sz[i] = 1;
	edge.resize(n + 1);
	vector<ed> edg(m);
	vector<bool> has(m);
	for(int i = 0; i < m; i++)
		cin >> edg[i].a >> edg[i].b >> edg[i].c;
	sort(edg.begin(), edg.end());
	for(int i = 0; i < m; i++) {
		if(uni(edg[i].a, edg[i].b))
			has[i] = 1;
	}
	for(int i = 0; i < m; i++) {
		if(has[i]) {
			edge[edg[i].a].push_back({edg[i].b, edg[i].c});
			edge[edg[i].b].push_back({edg[i].a, edg[i].c});
		}
	}
	for(int i = 1; i <= n; i++) {
		if(!vis[i])
			dfs(i);
	}
	for(int i = 1; i < 15; i++)
		for(int j = 1; j <= n; j++)
			lca[i][j] = lca[i - 1][lca[i - 1][j]], mx[i][j] = max(mx[i - 1][j], mx[i - 1][lca[i - 1][j]]);
	int q;
	cin >> q;
	while(q--) {
		int a, b;
		cin >> a >> b;
		cout << solve(a, b) << '\n';
	}

	return 0;
}