【NOIP2013提高组 day1】货车运输

题目

---

题解

–这道题的目的其实就是求所有路径中最小值的最大值，用贪心的思想，那就要使每条边都尽量大，又要保证连通，这时候就可以找一个最大生成树了
然后就是要查询答案，直接dfs的话太慢了，我们只需要把两个点之间的联通路的最小值找出来，所以可以在求lca的时候求出答案

---

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
using namespace std;
const int MAXN=1e4+5;

int n,m;

struct hehe{
    int x;
    int y;
    int z;
}edge[50005];
int f[MAXN];

int fa[20][MAXN],next[2*MAXN],to[2*MAXN],l[2*MAXN],head[MAXN];
int tot;
int w[20][MAXN];
int deep[MAXN];
int q;

bool comp(const hehe &a,const hehe &b){
    return a.z>b.z;
}

int Find(int x){
    if(x==f[x])
        return x;
    else{
        f[x]=Find(f[x]);
        return f[x];
    }
}

void add(int x,int y,int z){
    tot++;
    next[tot]=head[x];
    to[tot]=y;
    l[tot]=z;
    head[x]=tot;
}

void dfs(int a,int b){
    deep[a]=deep[b]+1;
    for(int i=head[a];i;i=next[i]){
        if(to[i]!=b){
            fa[0][to[i]]=a;
            w[0][to[i]]=l[i];
            dfs(to[i],a);
        }
    }
}

int lca(int x,int y){
    if(deep[x]<deep[y])
        swap(x,y);
    for(int i=15;i>=0;i--){
        if(deep[x]-(1<<i)>=deep[y])
            x=fa[i][x];
    }
    if(x==y)
        return x;
    for(int i=15;i>=0;i--){
        if(fa[i][x]!=fa[i][y]){
            x=fa[i][x];
            y=fa[i][y];
        }
    }
    return fa[0][x];
}

int ans(int x,int y,int LCA){
    int ANS=0x7f7f7f7f;
    for(int i=15;i>=0;i--){
        if(deep[x]-(1<<i)>=deep[LCA]){
            ANS=min(ANS,w[i][x]);
            x=fa[i][x];
        }
    }
    for(int i=15;i>=0;i--){
        if(deep[y]-(1<<i)>=deep[LCA]){
            ANS=min(ANS,w[i][y]);
            y=fa[i][y];
        }
    }
    return ANS;
}

int main(){
//  freopen("truck.in","r",stdin);
//  freopen("truck.out","w",stdout);
    cin>>n>>m;
    for(int i=1;i<=m;i++)
        scanf("%d%d%d",&edge[i].x,&edge[i].y,&edge[i].z);
    sort(edge+1,edge+1+m,comp);
    for(int i=1;i<=n;i++)
        f[i]=i;
    for(int i=1;i<=m;i++){
        int fx=Find(edge[i].x);
        int fy=Find(edge[i].y);
        if(fx!=fy){
            f[fy]=fx;
            add(edge[i].x,edge[i].y,edge[i].z);
            add(edge[i].y,edge[i].x,edge[i].z);
        }
    }
    memset(w,0x7f7f7f7f,sizeof(w));
    dfs(1,0);
    for(int i=1;(1<<i)<=n;i++){
        for(int j=1;j<=n;j++){
            fa[i][j]=fa[i-1][fa[i-1][j]];
            w[i][j]=min(w[i-1][j],w[i-1][fa[i-1][j]]);
        }
    }
    cin>>q;
    for(int i=1;i<=q;i++){
        int x,y;
        scanf("%d%d",&x,&y);
        int fx=Find(x);
        int fy=Find(y);
        if(fx!=fy)
            printf("-1\n");
        else
            printf("%d\n",ans(x,y,lca(x,y)));
    }
    return 0;
}