【Noip2013】货车运输

3287 货车运输

 

2013年NOIP全国联赛提高组

 时间限制: 1 s

 空间限制: 128000 KB

 题目等级 : 钻石 Diamond

题解

题目描述 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。

算法的话应该很明显了：生成树加树上倍增，也没必要去纠结最大最小，生成树瓶颈树。

反正就是按边从大到小排序，之后搜生成树，生成完毕就break。记得每个点Dfs一次，

之后就直接倍增跑LCA。把每次sum（a，b）换成max（a，b）就行。

#include <cstdio>  
#include <algorithm>  
  
using namespace std;  
  
const int maxx = 50000 + 100;  
const int Inf = 100000 + 500;

template<class T> T chkmin(const T &a,const T &b){return a<b? a:b;}

int head[maxx],depth[maxx],father[maxx];  
int grand[maxx][20+2],dis[maxx][20+2];  
bool done[maxx];  
int n,m,x,y,z,num,root,cnt,t;  

inline int Read(){  
    int x=0,f=1;char c=getchar();  
    while(c>'9'||c<'0') {if(c=='-') f=-1;c=getchar();}  
    while(c>='0'&&c<='9') {x=x*10+c-'0'; c=getchar();}  
    return x*f;  
}

struct Edge{  
    int next;  
    int to;  
    int value;  
}Edges[maxx<<1];  

struct Road{
    int x;
    int y;
    int v;
}Roads[maxx<<1];

int cmp(Road a,Road b){
    return a.v>b.v;
}

int find(int x){
    if(father[x]!=x)
    father[x]=find(father[x]);
    return father[x];
}

void unionn(int x,int y){
    int r1=find(x),r2=find(y);
    if(r1!=r2)
    father[r1]=r2;
}
  
void Add(int x,int y,int z){  
    Edges[++num].to=y;  
    Edges[num].next=head[x];  
    Edges[num].value=z;  
    head[x]=num;  
}  
  
void Dfs(int x){  
    done[x]=true;
    for(int i=1;i<=20;i++){  
        if((1<<i) > depth[x]) break;  
        grand[x][i]=grand[grand[x][i-1]][i-1];  
        dis[x][i]=chkmin(dis[x][i-1],dis[grand[x][i-1]][i-1]);  
    }  
    for(int i=head[x];i;i=Edges[i].next){  
        int now=Edges[i].to;  
        if(done[now]) continue;  
        depth[now]=depth[x]+1;  
        grand[now][0]=x;  
        dis[now][0]=Edges[i].value;  
        Dfs(now);
    }  
}  
  
inline int Lca(int x,int y){  
    int Ans=Inf;
    if(depth[x]>depth[y]) swap(x,y);  
    int d=depth[y]-depth[x];  
    for(int i=0;i<=20;i++){  
        if((1<<i) & d) {  
           Ans = chkmin(Ans,dis[y][i]);  
           y=grand[y][i];  
        }
    }  
    for(int i=20;i>=0;i--){  
        if(grand[x][i]!=grand[y][i]){  
            Ans = chkmin(Ans,dis[x][i]);  
            Ans = chkmin(Ans,dis[y][i]);  
            x=grand[x][i];  
            y=grand[y][i];  
        }
    }
    if(x==y) return Ans;  
    else{
        Ans = chkmin(Ans,dis[x][0]);  
        Ans = chkmin(Ans,dis[y][0]); 
        return Ans;
    }  
}  
  
int main(){  
    n=Read();m=Read();
    for(int i=1;i<=n;i++)
        father[i]=i;
    for(int i=1;i<=m;i++)
        Roads[i].x=Read(),Roads[i].y=Read(),Roads[i].v=Read();
    sort(Roads+1,Roads+m+1,cmp);
    for(int i=1;i<=m;i++){
        if(find(Roads[i].x)!=find(Roads[i].y)){
            Add(Roads[i].x,Roads[i].y,Roads[i].v);
            Add(Roads[i].y,Roads[i].x,Roads[i].v);
            unionn(Roads[i].x,Roads[i].y);
            cnt++;   
        }
        if(cnt==n-1)
            break;
    }
    for(int i=1;i<=n;i++)
        Dfs(i);  
    t = Read();
    while( t-- ){
        x=Read();y=Read();
        if(find(x)!=find(y))
            printf("-1\n");
        else
            printf("%d\n",Lca(x,y));  
    }
    return 0;
}