【Bzoj3732】Network

题意

给你N个点的无向图 (1 <= N <= 15,000)，记为：1…N。
图中有M条边 (1 <= M <= 30,000) ，第j条边的长度为： d_j ( 1 < = d_j < = 1,000,000,000).

现在有 K个询问 (1 < = K < = 20,000)。
每个询问的格式是：A B，表示询问从A点走到B点的所有路径中，最长的边最小值是多少？

---

解析

类似于那啥货车运输，嗯，懒得倍增就打了树链剖分。

---

#include <cstdio>#include <algorithm>#define Rep( i , _begin , _end ) for(int i=(_begin),i##_END=(_end);i<=(i##_END);i++)#define For( i , _begin , _end ) for(int i=(_begin),i##_END=(_end);i!=(i##_END);i++)#define Lop( i , _begin , _end ) for(int i=(_begin),i##_END=(_end);i>=(i##_END);i--)#define Dnt( i , _begin , _end ) for(int i=(_begin),i##_END=(_end);i!=(i##_END);i--)using std :: max;using std :: min;const int maxx = 15000 + 25;const int maxm = 30000 + 25;const int Inf = (unsigned)(-1) >> 1;typedef int Array[maxx];int n,m,x,y,z,num,cnt,ans,tot,t;int head[maxx],to[maxm],nxt[maxm],w[maxm];int Tx[maxx<<2];Array v,top,dpt,ftr,son,size,fky,rnk;class Edges{    public:        int x;        int y;        int v;}a[maxm<<1];bool cmp(Edges a,Edges b) {return a.v < b.v;}namespace Segment_Tree{    void modify(int i,int pos,int l,int r,int k){        if(l == r) {Tx[i] = k;return;}        int mid = (l+r) >> 1;        if(pos <= mid) modify(i<<1,pos,l,mid,k);        if(pos >  mid) modify(i<<1|1,pos,mid+1,r,k);        Tx[i] = max(Tx[i<<1],Tx[i<<1|1]);    }    int Query(int i,int x,int y,int l,int r){        if(x <= l && r <= y) return Tx[i];        int ans = -Inf;        int mid = (l+r) >> 1;        if(x <= mid) ans = max(ans,Query(i<<1,x,y,l,mid));        if(y >  mid) ans = max(ans,Query(i<<1|1,x,y,mid+1,r));        return ans;    }}namespace Kruskal{    int find(int x){        return x == fky[x]? x : fky[x] = find(fky[x]);    }}namespace Cute{    using namespace Segment_Tree;    void Ins(int x,int y,int z){        to[++num] = y;w[num] = z;nxt[num] = head[x];head[x] = num;    }    void Dfs(int x){        size[x] = 1;        for(int i=head[x];i;i=nxt[i]){            if(to[i] == ftr[x]) continue;            dpt[to[i]] = dpt[x] + 1;            ftr[to[i]] = x;            v[to[i]] = w[i];            Dfs(to[i]);size[x] += size[to[i]];            if(size[to[i]] > size[son[x]]) son[x] = to[i];        }    }    void __Dfs(int x,int brn){        rnk[x] = ++cnt;top[x] = brn;        if(son[x]) __Dfs(son[x],brn);        for(int i=head[x];i;i=nxt[i])            if(to[i] != son[x] && to[i] != ftr[x])                __Dfs(to[i],to[i]);    }    int Get(int x,int y){        int ans = -Inf;        while(top[x] != top[y]){            if(dpt[top[x]] > dpt[top[y]]) std :: swap(x,y);            ans = max(ans,Query(1,rnk[top[y]],rnk[y],1,n));            y = ftr[top[y]];        }        if(rnk[x] > rnk[y]) std :: swap(x,y);        if(x != y) ans = max(ans,Query(1,rnk[x]+1,rnk[y],1,n));        return ans;    }}using namespace Kruskal;using namespace Segment_Tree;using namespace Cute;int main(){    scanf("%d%d%d",&n,&m,&t);    Rep( i , 1 , n ) fky[i] = i;    Rep( i , 1 , m ) scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].v);    std :: sort(a+1,a+m+1,cmp);    Rep( i , 1 , m ){        int u = find(a[i].x),v = find(a[i].y);        if(u != v){            Ins(a[i].x,a[i].y,a[i].v);            Ins(a[i].y,a[i].x,a[i].v);            fky[u] = v;            tot ++;        }        if(tot == n-1) break;    }    Dfs(1);__Dfs(1,1);    Rep( i , 1 , n ) modify(1,rnk[i],1,n,v[i]);    while( t-- ){        scanf("%d%d",&x,&y);        printf("%d\n",Get(x,y));    }    return 0;}