首页 > 代码库 > BZOJ3732 Network
BZOJ3732 Network
这貌似是13年的noip最后一道吧?、、、
蒟蒻只会这种题呢、、、
Kruskal求出MST,然后倍增就好了
1 /************************************************************** 2 Problem: 3732 3 User: rausen 4 Language: C++ 5 Result: Accepted 6 Time:268 ms 7 Memory:4412 kb 8 ****************************************************************/ 9 10 #include <cstdio> 11 #include <algorithm> 12 13 using namespace std; 14 const int N = 16005; 15 const int M = 30005; 16 17 struct Edge { 18 int x, y, v; 19 20 inline bool operator < (const Edge &x) const { 21 return v < x.v; 22 } 23 } E[M]; 24 25 struct edge { 26 int next, to, v; 27 edge() {} 28 edge(int _n, int _t, int _v) : next(_n), to(_t), v(_v) {} 29 } e[M]; 30 31 struct tree_node { 32 int fa[16], mx[16], dep; 33 } tr[N]; 34 35 int n, m; 36 int fa[N]; 37 int tot, first[N]; 38 39 inline int read() { 40 int x = 0; 41 char ch = getchar(); 42 while (ch < ‘0‘ || ‘9‘ < ch) 43 ch = getchar(); 44 while (‘0‘ <= ch && ch <= ‘9‘) { 45 x = x * 10 + ch - ‘0‘; 46 ch = getchar(); 47 } 48 return x; 49 } 50 51 inline void Add_Edges(int x, int y, int v) { 52 e[++tot] = edge(first[x], y, v), first[x] = tot; 53 e[++tot] = edge(first[y], x, v), first[y] = tot; 54 } 55 56 int find_fa(int x) { 57 return x == fa[x] ? x : fa[x] = find_fa(fa[x]); 58 } 59 60 void Kruskal() { 61 int i, cnt, fa1, fa2; 62 sort(E + 1, E + m + 1); 63 for (i = 1; i <= n; ++i) 64 fa[i] = i; 65 for (i = 1, cnt = 0; i <= m; ++i) { 66 fa1 = find_fa(E[i].x), fa2 = find_fa(E[i].y); 67 if (fa1 != fa2) { 68 fa[fa1] = fa2, ++cnt; 69 Add_Edges(E[i].x, E[i].y, E[i].v); 70 if (cnt == n - 1) break; 71 } 72 } 73 } 74 75 void dfs(int p) { 76 int x, y; 77 for (x = 1; x < 16; ++x) { 78 tr[p].fa[x] = tr[tr[p].fa[x - 1]].fa[x - 1]; 79 tr[p].mx[x] = max(tr[p].mx[x - 1], tr[tr[p].fa[x - 1]].mx[x - 1]); 80 } 81 for (x = first[p]; x; x = e[x].next) 82 if ((y = e[x].to) != tr[p].fa[0]) { 83 tr[y].dep = tr[p].dep + 1; 84 tr[y].fa[0] = p, tr[y].mx[0] = e[x].v; 85 dfs(y); 86 } 87 } 88 89 int query(int x, int y) { 90 int res = 0, i; 91 if (tr[x].dep < tr[y].dep) swap(x, y); 92 for (i = 15; ~i; --i) 93 if (tr[tr[x].fa[i]].dep >= tr[y].dep) { 94 res = max(res, tr[x].mx[i]); 95 x = tr[x].fa[i]; 96 } 97 for (i = 15; ~i; --i) 98 if (tr[x].fa[i] != tr[y].fa[i]) { 99 res = max(res, max(tr[x].mx[i], tr[y].mx[i]));100 x = tr[x].fa[i], y = tr[y].fa[i];101 }102 if (x != y)103 res = max(res, max(tr[x].mx[0], tr[y].mx[0]));104 return res;105 }106 107 int main() {108 n = read(), m = read();109 int Q = read(), i, x, y;110 for (i = 1; i <= m; ++i)111 E[i].x = read(), E[i].y = read(), E[i].v = read();112 Kruskal();113 tr[1].dep = 1;114 dfs(1);115 while (Q--) {116 x = read(), y = read();117 printf("%d\n", query(x, y));118 }119 return 0;120 }
BZOJ3732 Network
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。