倍增法求lca模板

/*
倍增法求lca
*/
const int maxlogv=16;
int f[maxlogv][MAXN];
int dep[MAXN];
//先dfs
void lca_init(){
    for(int k=0;k+1<maxlogv;++k){
        for(int v=1;v<=n;++v){
            if(f[k][v]==0) f[k+1][v]=0;
            else f[k+1][v]=f[k][f[k][v]];
        }
    }
}
void dfs(int u,int fa)
{
    f[0][u]=fa;dep[u]=dep[fa]+1;
    for(int k=head[u];k!=-1;k=e[k].next){
        int v=e[k].v;
        if(v==fa) continue;
        dfs(v,u);
    }
}
 
int lca(int u,int v)
{
    if(dep[u]>dep[v]) swap(u,v);
    for(int k=0;k<maxlogv;++k){
        if( (dep[v]-dep[u])>>k&1){
            v=f[k][v];
        }
    }
    if(u==v) return u;
    for(int k=maxlogv-1;k>=0;--k){
        if(f[k][u]!=f[k][v]){
            u=f[k][u];
            v=f[k][v];
        }
    }
    return f[0][u];
}