NOIP 2013 货车运输 TYVJ 3115

---

题目大意：
中文题，不翻译了哈>_<

---

分析：
什么暴力做法30、60分算法打击都会咱就不探讨了，我们就说一说怎么拿满分？？
大家都能看出来这道题是要先求一棵最大生成树，然后每两个点之间的路径是唯一的，所以我们只需要求出xy的LCA并且求出路径上的最小值即可，重点是怎么求捏？？
我们只需要在维护f[i][j]这个记录祖先的数组的同时维护一个与ij同样意义的g[i][j]数组来记录i到它的2的j次方祖先的路径上的最小值就好啦

---

代码如下：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#define inf 0x3f3f3f3f
#define Min(a,b,c) min(a,min(b,c))
using namespace std;
const int maxm=50000+5,maxn=10000+5;
int n,m,fa[maxn],f[maxn][30],g[maxn][30],dep[maxn],cnt,to[maxm*2],nxt[maxm*2],hd[maxn],w[maxm*2];
struct edge{
    int x,y,w; 
    friend bool operator < (edge a,edge b){
        return a.w>b.w;
    }
}s[maxm];
inline int read(){
    char ch=getchar();
    int f=1,x=0;
    while(!(ch>='0'&&ch<='9')){
        if(ch=='-')
            f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
        x=x*10+ch-'0',ch=getchar();
    return f*x;
}
inline int find(int x){
    return fa[x]==x ? x : fa[x]=find(fa[x]);
}
inline void add(int x,int y,int z){
    to[cnt]=y;
    w[cnt]=z;
    nxt[cnt]=hd[x];
    hd[x]=cnt++;
}
void kruskal(){
    for(int i=1;i<=m;i++){
        int fx=find(s[i].x),fy=find(s[i].y);
        if(fx!=fy)
            fa[fy]=fx,add(s[i].x,s[i].y,s[i].w),add(s[i].y,s[i].x,s[i].w);
    }
}
inline void dfs(int root,int fa){
    for(int i=hd[root];i!=-1;i=nxt[i])
        if(to[i]!=fa)
            f[to[i]][0]=root,g[to[i]][0]=w[i],dep[to[i]]=dep[root]+1,dfs(to[i],root);
}
inline void init(){
    for(int j=1;j<=25;j++)
        for(int i=1;i<=n;i++)
            f[i][j]=f[f[i][j-1]][j-1],g[i][j]=Min(g[i][j],g[i][j-1],g[f[i][j-1]][j-1]);
}
inline int LCA(int a,int b){
    int ans=inf;
    if(dep[a]<dep[b])
        swap(a,b);
    int d=dep[a]-dep[b];
    for(int i=0;i<=25;i++)
        if((1<<i)&d)
            ans=min(ans,g[a][i]),a=f[a][i];
    if(a==b)
        return ans;
    for(int i=25;i>=0;i--)
        if(f[a][i]!=f[b][i])
            ans=Min(ans,g[a][i],g[b][i]),a=f[a][i],b=f[b][i];
    ans=Min(ans,g[a][0],g[b][0]);
    return ans;
}
signed main(void){
    n=read(),m=read();
    memset(hd,-1,sizeof(hd)),cnt=0;
    for(int i=1;i<=m;i++)
        s[i].x=read(),s[i].y=read(),s[i].w=read();
    sort(s+1,s+m+1);
    for(int i=1;i<=n;i++)
        f[i][0]=i,fa[i]=i;
    memset(g,inf,sizeof(g));
    kruskal();
    int q=read();
    for(int i=1;i<=n;i++)
        if(f[i][0]==i)
            dep[i]=0,dfs(i,-1);
    init();
    while(q--){
        int x=read(),y=read();
        int fx=find(x),fy=find(y);
        if(fx!=fy)
            cout<<-1<<endl;
        else
            cout<<LCA(x,y)<<endl;
    }
    return 0;
}

---

by >o< neighthorn