Uva 11354 Bond（最小生成树+LCA）

题目地址：https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2339

思路：

1.寻找s到t路径上最小瓶颈路。则该路径必在最小生成树上：最小生成树上为权值尽可能短的n-1条边，若该路径上最小的最大权值边不在最小生成树上，则此权值会比最小生成树上的权值大（假设存在一条路径u,p[1],......p[n],v，所求边不在最小生成树上，即其权值比最小生成树上权值小。又由于按照贪心原则，权值小的边优先加入生成树中，该边应位于MST中，矛盾）。

2.求出MST后，按倍增的方法求出s到t的最大权值即可。

#include<cstdio>
#include<vector>
#include<cstring>
#include<iostream>
#include<algorithm>
#define debug
using namespace std;
const int INF=1e9;
const int maxn=5e4+50;
struct Node1
{
    int st,ed,w;
    Node1(int st=0,int ed=0,int w=0):st(st),ed(ed),w(w) {}
};
struct Node2
{
    int v,w;
    Node2(int v=0,int w=0):v(v),w(w) {}
};
int n,m,q;
vector<Node1> edge;
vector<Node2> g[maxn];
int fa[maxn],anc[maxn][50];
int v[maxn],d[maxn],mm[maxn][50];
int Find(int root)
{
    int x=root;
    while(x!=fa[x]) x=fa[x];
    return x;
}
void Union(int r1,int r2)
{
    fa[r1]=r2;
}
void dfs(int x,int dep)
{
    v[x]=1,d[x]=dep;
    for(int i=0; i<g[x].size(); i++)
    {
        int nt=g[x][i].v;
        if(!v[nt])
        {
            anc[nt][0]=x;
            mm[nt][0]=g[x][i].w;
            int k=0;
            while(anc[anc[nt][k]][k]!=0)
            {
                anc[nt][k+1]=anc[anc[nt][k]][k];
                mm[nt][k+1]=max(mm[nt][k],mm[anc[nt][k]][k]);
                k++;
            }
            dfs(nt,dep+1);
        }
    }
}
int cmp(Node1 a,Node1 b)
{
    return a.w<b.w;
}
void build()
{
    sort(edge.begin(),edge.end(),cmp);
    int k=0;
    for(int i=0; i<edge.size(); i++)
    {
        int f1=Find(edge[i].st);
        int f2=Find(edge[i].ed);
        if(f1!=f2)
        {
            Union(f1,f2);
            g[edge[i].st].push_back(Node2(edge[i].ed,edge[i].w));
            g[edge[i].ed].push_back(Node2(edge[i].st,edge[i].w));
            k++;
        }
        if(k==n-1) break;
    }
}
void mindistanc()
{
    for(int i=1; i<=n; i++)
        if(!v[i]) dfs(i,0);
    scanf("%d",&q);
    for(int i=1; i<=q; i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        if(d[x]<d[y]) swap(x,y);
        int l=d[x]-d[y];
        int k=0,ans=-INF;
        while(l!=0)
        {
            if(l&1)
            {
                ans=max(ans,mm[x][k]);
                x=anc[x][k];
            }
            l>>=1;
            k++;
        }
        k=0;
        while(x!=y)
        {
            if((anc[x][k]!=anc[y][k])||(!k))
            {
                ans=max(ans,max(mm[x][k],mm[y][k]));
                x=anc[x][k],y=anc[y][k];
                k++;
            }
            else k--;
        }
        printf("%d\n",ans);
    }
}
void init()
{
    edge.clear();
    memset(v,0,sizeof(v));
    memset(d,0,sizeof(d));
    memset(anc,0,sizeof(anc));
    memset(mm,0,sizeof(mm));
    for(int i=0; i<=n; i++)
    {
        fa[i]=i,g[i].clear();
    }
}
int main()
{
#ifdef debu
    freopen("in.txt","r",stdin);
#endif // debug
    int cas=0;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        init();
        cas++;
        if(cas!=1) printf("\n");
        for(int i=0; i<m; i++)
        {
            int x,y,w;
            scanf("%d%d%d",&x,&y,&w);
            edge.push_back(Node1(x,y,w));
        }
        build();
        mindistanc();
    }
    return 0;
}