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+;

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

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

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

 void dfs(int x){
     for(int i=;fa[x][i]&&fa[fa[x][i]][i];i++){
         fa[x][i+]=fa[fa[x][i]][i];
     }
     for(int i=egh[x];i;i=eg[i][]){
         int b=eg[i][];if(b==fa[x][]) continue;
         dep[b]=dep[x]+,fa[b][]=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]+);i>=;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>=;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][];
 }

 int main(){
     //freopen("","r",stdin);
     int i;
     N=rd(),M=rd();
     for(i=;i<N;i++){
         int a=rd(),b=rd();
         adeg(a,b);adeg(b,a);
     }
     dep[]=;dfs();
     for(i=;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\n",a);
         else if(dep[b]==mm) printf("%d\n",b);
         else printf("%d\n",c);
     }
     return ;
 }