[最大生成树+LCA]NOIP 2013——货车运输

[最大生成树+LCA]NOIP 2013——货车运输

题目梗概

给出一张有权无向图。

每次询问x->y的所有路径，最小值边权最大的路径。

解题思路

构造一棵最大生成树（可能是多棵）。

每次在树上求一条路径的最小值就可以了。

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=10005,maxm=50005;
struct jz{
    int x,y,w;
    bool operator<(const jz &b)const{return w>b.w;}
}a[maxm];
int n,m,Q,fa[maxn],f[maxn][15],g[maxn][15],dep[maxn];
int tot,lnk[maxn],son[2*maxn],nxt[2*maxn],w[2*maxn];
bool vis[maxn];
int get(int x){
    if (fa[x]==x) return x;
    fa[x]=get(fa[x]);
    return fa[x];
}
void add(int x,int y,int z){nxt[++tot]=lnk[x];lnk[x]=tot;son[tot]=y;w[tot]=z;}
void DFS(int x){
    vis[x]=1;
    for (int j=lnk[x];j;j=nxt[j]) if (!vis[son[j]])
    g[son[j]][0]=w[j],f[son[j]][0]=x,dep[son[j]]=dep[x]+1,DFS(son[j]);
}
int query(int x,int y){
    if (dep[x]<dep[y]) swap(x,y);int ans=g[0][0];
    for (int j=log2(n);j>=0;j--) if (dep[f[x][j]]>=dep[y]) ans=min(ans,g[x][j]),x=f[x][j];
    if (x==y) return ans;
    for (int j=log2(n);j>=0;j--) if (f[x][j]!=f[y][j]) ans=min(ans,min(g[x][j],g[y][j])),x=f[x][j],y=f[y][j];
    return min(ans,min(g[x][0],g[y][0]));
}
int main(){
    freopen("exam.in","r",stdin);
    freopen("exam.out","w",stdout);
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++) scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].w);
    sort(a+1,a+1+m);
    for (int i=1;i<=n;i++) fa[i]=i;
    for (int i=1;i<=m;i++){
        int x=get(a[i].x),y=get(a[i].y);
        if (x==y) continue;fa[x]=y;
        add(a[i].x,a[i].y,a[i].w);
        add(a[i].y,a[i].x,a[i].w);
    }
    memset(g,63,sizeof(g));
    for (int i=1;i<=n;i++) if (!vis[i]) DFS(i);
    for (int j=1;j<=log2(n);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-1],g[f[i][j-1]][j-1]);
    scanf("%d",&Q);
    while(Q--){
        int x,y;scanf("%d%d",&x,&y);
        if (get(x)!=get(y)){printf("-1\n");continue;}
        printf("%d\n",query(x,y));
    }
    return 0;
}