洛谷P1967 货车运输

题目传送门

额自己用的kruskal+倍增LCA，AC后看大佬都是用的什么LCT (好像是这个名字)、树链剖分这些我这个蒟蒻不会的算法

思路很简单就是求一个最大生成树，然后用倍增LCA维护路径最小值，属于那种一眼题

但是我觉得对我这种蒟蒻很不友好，就是他的细节很多（OZT），现在公开处刑一下我自己

• 这个题没说一定联通，起码我是没想出用bfs维护的时候哪个点当起点，然后就用的深搜,记得起点的w是INF
• 因为一条边$(u,v)$的w我存在了v上面，所以在类似dp维护的时候应该是$w[j][i]=min(w[j][i-1],w[f[j][i-1]][i-1])$，类似st表分成了2段，我一开始就写错了
• 最后是LCA里面要先求ans在往上跳。。。

#include<bits/stdc++.h>
#define ll long long
#define INF 0x3f3f3f3f
using namespace std;
const int N=1e4+10,M=1e5+10;
int n,m,q,k;
int head[N],ver[M],nxt[M],val[M],idx;
int f[N][20],d[N],fa[N],vis[N],w[N][20];
struct node{
    int x,y,z;
    bool operator<(const node&h)const{return z>h.z;}
}e[M>>1];
inline void add(int u,int v,int w){
    ver[idx]=v,val[idx]=w,nxt[idx]=head[u],head[u]=idx++;
}
inline int Find(int x){
    return x==fa[x]?x:fa[x]=Find(fa[x]);
}
inline void kruskal(){
    int num=0;
    sort(e+1,e+m+1);
    for(int i=1;i<=m;++i){
        int x=Find(e[i].x),y=Find(e[i].y);
        int u=e[i].x,v=e[i].y,w=e[i].z;
        if(x==y) continue;
        add(u,v,w),add(v,u,w);
        fa[x]=y;
        if(++num==n-1) return;
    }
}
inline int lca(int x,int y){
    int ans=INF;
    if(Find(x)!=Find(y)) return -1;
    if(d[x]>d[y]) swap(x,y);
    for(int i=k;i>=0;i--)
        if(d[f[y][i]]>=d[x]){
            ans=min(ans,w[y][i]);
            y=f[y][i];
        }
    //printf("%d\n",ans);
    if(x==y) return ans;
    for(int i=k;i>=0;i--){
        if(f[x][i]!=f[y][i]){
            ans=min(ans,w[x][i]),ans=min(ans,w[y][i]);
            x=f[x][i],y=f[y][i];
        }
    }
    ans=min(ans,w[x][0]),ans=min(ans,w[y][0]);
    return ans;
}
void dfs(int x){
    vis[x]=1;
    for(int i=head[x];~i;i=nxt[i]){
        int y=ver[i];
        if(vis[y]) continue;
        d[y]=d[x]+1;
        f[y][0]=x;
        w[y][0]=val[i];
        dfs(y);
    }
    return ;
}
int main(){
    memset(head,-1,sizeof head);
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;++i) fa[i]=i;
    for(int i=1;i<=m;++i) scanf("%d%d%d",&e[i].x,&e[i].y,&e[i].z);
    kruskal();
    for(int i=1;i<=n;++i){
        if(!vis[i]){
            d[i]=1;
            dfs(i);
            //f[i][0]=i;
            w[i][0]=INF;
        }
    }
    k=(int)(log(n)/log(2))+1;
    for(int i=1;i<=k;++i){
        for(int j=1;j<=n;++j){
            f[j][i]=f[f[j][i-1]][i-1];
            w[j][i]=min(w[j][i-1],w[f[j][i-1]][i-1]);
           // printf("%d %d\n",f[j][i],w[j][i]);
        }
    }
    scanf("%d",&q);
    while(q--){
        int x,y;
        scanf("%d%d",&x,&y);
        printf("%d\n",lca(x,y));
    }
    return 0;
}