【模板】最近公共祖先（LCA）

题目描述

如题，给定一棵有根多叉树，请求出指定两个点直接最近的公共祖先。

输入输出格式

输入格式：

第一行包含三个正整数N、M、S，分别表示树的结点个数、询问的个数和树根结点的序号。

接下来N-1行每行包含两个正整数x、y，表示x结点和y结点之间有一条直接连接的边（数据保证可以构成树）。

接下来M行每行包含两个正整数a、b，表示询问a结点和b结点的最近公共祖先。

输出格式：

输出包含M行，每行包含一个正整数，依次为每一个询问的结果。

思路

推荐用倍增法求LCA。
先构建倍增数组(f[i][j]=f[f[i-1][j-1]][j-1])
根据倍增数组，如果跳到的点相同，则不跳，否则跳。
最后输出father[x]。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=500007,ln=20;
struct E
{
    int to,next;
}e[maxn*2];
int p[maxn][ln+1],d[maxn],list[maxn],cnt=1,n,q,u;
void add(int u,int v)
{
    e[++cnt].to=v; e[cnt].next=list[u]; list[u]=cnt;
}
void dfs(int r,int fa)
{
    d[r]=d[fa]+1;
    p[r][0]=fa;
    for(int i=1; (1<<i)<=d[r]; i++) p[r][i]=p[p[r][i-1]][i-1];
    for(int i=list[r]; i!=-1; i=e[i].next)
    {
        int u=e[i].to;
        if(u!=fa) dfs(u,r);
    }
}
int gcd(int x,int y)
{
    if(d[x]>d[y]) swap(x,y);
    for(int i=ln; i>=0; i--)
    {
        if(d[x]<=d[y]-(1<<i)) y=p[y][i];
    }
    if(x==y) return x;
    for(int i=ln; i>=0; i--)
    {
        if(p[x][i]==p[y][i]) continue;else
        {
            x=p[x][i]; y=p[y][i];
        }
    }
    return p[x][0];
}
int main()
{
    memset(list,-1,sizeof(list));
    scanf("%d%d%d",&n,&q,&u);
    for(int i=1; i<=n-1; i++) 
    {
        int x,y;
        scanf("%d%d",&x,&y);
        add(x,y); add(y,x);
    }
    dfs(u,0);
    while(q--)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        printf("%d\n",gcd(x,y));
    }
}