suoi31 最近公共祖先2 (倍增lca)

根为r时x、y的公共祖先，就是lca(x,r),lca(x,y),lca(r,y)中深度最大的那一个，不要再在倍增的时候判来判去还判不对了...

#include<bits/stdc++.h>
#define pa pair<int,int>
#define CLR(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
const int maxn=3e5+10;

inline ll rd(){
    ll x=0;char c=getchar();int neg=1;
    while(c<'0'||c>'9'){if(c=='-') neg=-1;c=getchar();}
    while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*neg;
}

int eg[maxn*2][2],egh[maxn],ect;
int dfn[maxn][2],tot;
int dep[maxn],fa[maxn][20];
int N,M;

inline void adeg(int a,int b){
    eg[++ect][0]=b;eg[ect][1]=egh[a];egh[a]=ect;
}

void dfs(int x){
    for(int i=0;fa[x][i]&&fa[fa[x][i]][i];i++){
        fa[x][i+1]=fa[fa[x][i]][i];
    }
    for(int i=egh[x];i;i=eg[i][1]){
        int b=eg[i][0];if(b==fa[x][0]) continue;
        dep[b]=dep[x]+1,fa[b][0]=x;
        dfs(b);
    }
}

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

int main(){
    //freopen("","r",stdin);
    int i;
    N=rd(),M=rd();
    for(i=1;i<N;i++){
        int a=rd(),b=rd();
        adeg(a,b);adeg(b,a);
    }
    dep[1]=1;dfs(1);
    for(i=1;i<=M;i++){
        int r=rd(),x=rd(),y=rd();
        int a=lca(x,y),b=lca(x,r),c=lca(y,r);
        int mm=max(dep[a],max(dep[b],dep[c]));
        if(dep[a]==mm) printf("%d
",a);
        else if(dep[b]==mm) printf("%d
",b);
        else printf("%d
",c);
    }
    return 0;
}