[NOIP2013] 货车运输 最大生成树 LCA

货车运输

题目描述 Description
A 国有 n 座城市，编号从 1 到 n，城市之间有 m 条双向道路。每一条道路对车辆都有重量限制，简称限重。现在有 q 辆货车在运输货物，司机们想知道每辆车在不超过车辆限重的情况下，最多能运多重的货物。

输入描述 Input Description
第一行有两个用一个空格隔开的整数 n，m，表示 A 国有 n 座城市和 m 条道路。
接下来 m 行每行 3 个整数 x、y、z，每两个整数之间用一个空格隔开，表示从 x 号城市到 y 号城市有一条限重为 z 的道路。注意：x 不等于 y，两座城市之间可能有多条道路。
接下来一行有一个整数 q，表示有 q 辆货车需要运货。
接下来 q 行，每行两个整数 x、y，之间用一个空格隔开，表示一辆货车需要从 x 城市运输货物到 y 城市，注意：x 不等于 y。

输出描述 Output Description
输出共有 q 行，每行一个整数，表示对于每一辆货车，它的最大载重是多少。如果货车不能到达目的地，输出-1。

样例输入 Sample Input
4 3
1 2 4
2 3 3
3 1 1
3
1 3
1 4
1 3

样例输出 Sample Output
3
-1
3

数据范围及提示 Data Size & Hint
对于 30%的数据，0 < n < 1,000，0 < m < 10,000，0 < q < 1,000；
对于 60%的数据，0 < n < 1,000，0 < m < 50,000，0 < q < 1,000；
对于 100%的数据，0 < n < 10,000，0 < m < 50,000，0 < q < 30,000，0 ≤ z ≤ 100,000。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int INF = 100000000 + 7;
const int M = 100000 + 5;
const int N = 10000 + 5;
const int P = 16;
int n,m,q,cnt,tot,dep[N],last[N],fa[N],anc[N][17],d[N][17];
bool vis[N];
struct edge{
    int u,v,w;
}a[M];
struct Edge{
    int to,v,next;
}e[M];
void insert( int u, int v, int w ){
    e[++cnt].to = v; e[cnt].v = w; e[cnt].next = last[u]; last[u] = cnt;
    e[++cnt].to = u; e[cnt].v = w; e[cnt].next = last[v]; last[v] = cnt;
}
bool cmp( edge a, edge b ){
    return a.w > b.w;
}
int find( int x ){
    return fa[x]==x ? x : fa[x]=find(fa[x]);
}
void dfs( int x ){
    vis[x] = 1;
    for( int i = 1; i <= P; i++ ){
        if( dep[x] < (1<<i) ) break;
        anc[x][i] = anc[anc[x][i-1]][i-1];
        d[x][i] = min( d[x][i-1], d[anc[x][i-1]][i-1]);
    } for( int i = last[x]; i; i = e[i].next ){
        if( vis[e[i].to] ) continue;
        anc[e[i].to][0] = x;
        d[e[i].to][0] = e[i].v;
        dep[e[i].to] = dep[x]+1;
        dfs( e[i].to );
    }
}
int lca( int x, int y ){
    if( dep[x] < dep[y] ) swap( x, y );
    int derta = dep[x] - dep[y];
    for( int i = 0; i <= P; i++ ) if( (1<<i) & derta ) x = anc[x][i];
    for( int i = P; i >= 0 ; i-- ) if( anc[x][i] != anc[y][i] ){ x = anc[x][i]; y = anc[y][i]; }
    if( x == y ) return x;
    x = anc[x][0]; return x;
}
int query( int x, int y ){
    int minn = INF;
    int derta = dep[x] - dep[y];
    for( int i = 0; i <= P; i++ )
        if( derta & (1<<i) ){
            minn = min(minn,d[x][i]);
            x = anc[x][i];
        }
    return minn;
}
int main(){
    freopen("truck.in","r",stdin);
    freopen("truck.out","w",stdout);
    scanf("%d%d", &n, &m);
    for( int i = 1; i <= n; i++ ) fa[i] = i;
    for( int i = 1; i <= m; i++ ) scanf("%d%d%d", &a[i].u, &a[i].v, &a[i].w);
    std::sort( a+1, a+m+1, cmp );
    for( int i = 1; i <= m; i++ ){
        int u = a[i].u; int v = a[i].v;
        int x = find(a[i].u); int y = find(a[i].v);
        if( x != y ){
            fa[x] = y;
            insert( u, v, a[i].w );
            tot++; if( tot == n-1 ) break;
        }
    }
    for( int i = 1; i <= n; i++ ) if( !vis[i] ) dfs(i);
    scanf("%d", &q);
    for( int i = 1; i <= q; i++ ){
        int x, y; scanf("%d%d", &x, &y);
        if( find(x) != find(y) ){ printf("-1\n"); continue; }
        int anci = lca( x, y );
        int ans = min( query( x, anci ), query( y, anci ) );
        printf("%d\n", ans);
    }
    return 0;
}