zoukankan      html  css  js  c++  java
  • 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
    ", query(x, y));
    118     }
    119     return 0;
    120 }
    View Code
    By Xs酱~ 转载请说明 博客地址:http://www.cnblogs.com/rausen
  • 相关阅读:
    I Hate It HDU
    线段树为什么要开4倍空间
    敌兵布阵 HDU
    线段树 模板
    线段树原理详解
    HDU5037 Frog
    HDU5187 zhx's contest
    HDU5307 He is Flying
    UOJ#34. 多项式乘法
    Bzoj4766 文艺计算姬
  • 原文地址:https://www.cnblogs.com/rausen/p/4160979.html
Copyright © 2011-2022 走看看