LCA—倍增法求解

LCA（Least Common Ancestors）

即最近公共祖先，是指在有根树中，找出某两个结点u和v最近的公共祖先
常见解法一般有三种，这里讲解一种在线算法—倍增

---

首先我们定义fa[u][j]表示结点u的第2^j祖先
那么要怎么求出全部的fa数组呢

不难发现fa[u][0]就是u的父亲结点
这些父亲结点我们可以直接初始化

对于其他结点则有
$fa[u][j]=fa[fa[u][j-1]][j-1]$
即u的第 2 ^ (j-1) 祖先的第 2 ^ (j-1) 祖先 就是u的第2^j祖先（有点像快速幂那样分解）

---

预处理各节点深度+初始fa[u][0]

void dfs(int u,int pa)
{
    dep[u]=dep[pa]+1;
    fa[u][0]=pa;
    for(int i=head[u];i;i=E[i].nxt)
    {
        int v=E[i].v;
        if(v!=pa) dfs(v,u);
    }
}

预处理fa数组

for(int i=1;(1<<i)<=n;i++)//由于远的祖先要由近的得到，所以幂次一定在循环外层
for(int u=1;u<=n;u++)
fa[u][i]=fa[ fa[u][i-1] ][i-1];

---

预处理之后怎么求解LCA(u，v)呢
我么先假定dep[u]>dep[v]
则两点深度差 d=dep[u]-dep[v]

首先我们要把u升到与v同样的深度

for(int i=0;(1<<i)<=d;++i)
if((1<<i)&d) u=fa[u][i];

深度差d二进制下第k位是否位1，是则向上跳2^k个祖先，和二进制数转换十进制数是一样的道理
(注意二进制最右边一位是第0号位)

if( (1<<i) & d )

这一段用来检查d的二进制下第i位是否为1

---

抬升到相同高度后就可以开始查询LCA了

for(int i=32;i>=0;--i)
if(fa[u][i]!=fa[v][i])
{
	u=fa[u][i];
	v=fa[v][i];
}

即从u和v最远的祖先开始
如果u的第2^i祖先等于 v的第2^i祖先，就不移动
否则u和v同时上移2^i个深度
最后u的父亲一定是u和v的LCA

---

最后贴上一份完整代码 模板题传送门

#include<iostream>
#include<cstdio>
#include<vector>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;

int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

const int maxn=1000010;
int n,m,s;
struct node{int v,nxt;}E[maxn];
int head[maxn],tot;
int fa[maxn][35],dep[maxn];

void add(int u,int v)
{
	E[++tot].nxt=head[u];
	E[tot].v=v;
	head[u]=tot;
}

void dfs(int u,int pa)
{
    dep[u]=dep[pa]+1;
    fa[u][0]=pa;
    for(int i=head[u];i;i=E[i].nxt)
    {
        int v=E[i].v;
        if(v!=pa) dfs(v,u);
    }
}

int LCA(int u,int v)
{
	if(dep[u]<dep[v]) swap(u,v);
	int d=dep[u]-dep[v];
	
	for(int i=0;(1<<i)<=d;++i)
	if((1<<i)&d) u=fa[u][i];
	
	if(u==v) return v;//v一开始就是u的祖先，故LCA是v
	for(int i=32;i>=0;--i)
	if(fa[u][i]!=fa[v][i])
	{
		u=fa[u][i];
		v=fa[v][i];
	}
	
	return fa[u][0];
}

int main()
{
	n=read(); m=read(); s=read();
	for(int i=1;i<n;++i)
	{
		int u=read(),v=read();
		add(u,v); add(v,u);
	}
	
	dfs(s,-1);
	for(int i=1;(1<<i)<=n;++i)
    for(int u=1;u<=n;++u)
    fa[u][i]=fa[ fa[u][i-1] ][i-1];
	
	while(m--)
	{
		int u=read(),v=read();
		printf("%d\n",LCA(u,v));
	}
	return 0;
}