Nearest Common Ancestors POJ - 1330 LCA （树上倍增模板）

这是一道模板题

树上倍增的写法和应用

#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
const int maxn=10010;
vector<int>g[maxn];
int fa[maxn][15];  //表示节点i往上跳2^j次后的节点  注意j的取值范围
int d[maxn];    //当前节点所在的深度
int n;
bool vis[maxn];
int root;   //根结点
void init()  //初始化
{
    memset(fa,0,sizeof(fa));
    memset(vis,0,sizeof(vis));
    memset(d,0,sizeof(d));
}
void dfs(int u)
{
    int i;
    vis[u]=1;
    for(i=0;i<g[u].size();i++)  //每一个与u相连的节点i
    {
        int v=g[u][i];
        if(!vis[v])
        {
            d[v]=d[u]+1;
            dfs(v);
        }
    }
}//深搜出各点的深度，存在d中 
void bz()
{
    int i,j;
    for(j=1;j<=14;j++)
        for(i=1;i<=n;i++)
        fa[i][j]=fa[fa[i][j-1]][j-1];  //i的第2^j祖先就是i的第2^(j-1)祖先的第2^(j-1)祖先
           //i的父辈们的fa数组都已经计算完毕，所以可以用来计算i
}
int lca(int u,int v)
{
    if(d[u]<d[v])
    swap(u,v);   //保证深度大的点为u，方便操作
    int dc=d[u]-d[v];   //计算深度差
    int i;
    for(i=0;i<=14;i++)  //值得注意的是，这里需要从零枚举 
    {                    //此循环用于提到深度相同
        if((1<<i)& dc)
            u=fa[u][i];
    }
    if(u==v)   //如果深度一样时，两个点相同，直接返回
        return u;
    for(i=14;i>=0;i--)  //如果从小到大的话就有可能无法提到正确位置
    {
        if(fa[u][i]!=fa[v][i])//跳2^j步不一样，就跳，否则不跳 
        {
            u=fa[u][i];
            v=fa[v][i];
        }
    }
    u=fa[u][0]; //上述过程做完，两点都在LCA下一层，所以走一步即可
    return u;
}
int main()
{
    int t;
    int x,y;
    scanf("%d",&t);
    for(int k=1;k<=t;k++)
    {
        init();
        scanf("%d",&n);
        for(int i=0;i<=n;i++)
            g[i].clear();
        for(int i=1;i<n;i++)
        {
            scanf("%d%d",&x,&y);
            g[x].push_back(y);
            fa[y][0]=x;
            if(fa[x][0]==0)    //查找根节点
                root=x;
        }
        d[root]=1;
        dfs(root);
        bz();
        int a,b;
        scanf("%d%d",&a,&b);
        printf("%d\n",lca(a,b));
    }
    return 0;
}

本文参考博客

https://blog.csdn.net/dad3zz/article/details/49386465