LCA(Least Common Ancestors),即最近公共祖先,是指在有根树中,找出某两个结点u和v最近的公共祖先。 ———来自百度百科
如:
LCA(4,6)=4
LCA(10,11)=9
LCA(3,11)=1
暴力
求两个点的公共祖先,必定在有根树下进行,那么就是先从根结点开始遍历,先找到每个点的儿子,和更新每一个点 的 fa[i] ,更新每个点的深度,在这些预处理过后就可以开始LCA 了
就拿图中的LCA(3,8)为例:
我们LCA函数要求的就是3和8的最小公共祖先,
一开始3和8并不在同一层,所以他们可以要在同一层之后,然后开始向前一起找,才能找到LCA
1.就是将 dep[x] 和dep[y]更新成一样
2.然后一个个向上枚举,每次向上枚举一个父节点,如果一样那么就找到了最近公共祖先
inline void dfs(int x,int fa)//dfs建立每个结点之间的关系
{
dep[x]=dep[fa]+1;//这一个的深度,就是父节点的深度加1
father[x]=fa;//这个结点的父节点就是上次遍历的结点
for(int i=head[x];i;i=edge[i].ne)
{
int to=edge[i].to;
if(to!=fa)dfs(to,x);
}
}
inline int LCA(int x,int y)
{
if(x==y)return x;
if(dep[x]<dep[y])swap(x,y);//把x更新为小的那个结点
while(dep[x]>dep[y])x=father[x];//将其更新为同一层
if(x==y)return x;
while(x!=y)x=father[x],y=father[y];//不断向上找父亲
return x;
}
倍增法求LAC
我们发现,每次一个个地往上跳跃,找到父亲的时间复杂度至少为O(n),如果结合题目中的Q次询问,那么就可能会超时,那么我们可以想到如何简化网上找公共祖先的过程
我们不妨设
f[i][j]=k 为第 i 个 结点 第 2^j 的父亲结点的编号为k
f[i][j]=f[f[i][j-1][j-1] 可以从此结点的2^j-1 的父亲在向上 2^j-1就可以得到
那么j的范围是多少呢? 和i 的深度有关,第k深度的结点 i 有k个祖先,也就是 log2 k= j; 2^j=k;
所以可以预处理出每个数的log值
方便之后的运算
for(int i=1;i<=n;i++)
lg[i]=lg[i-1]+(1<<(lg[i-1]+1)==i);//处理每个数的log值
inline void dfs(int x,int fa)
{
dep[x]=dep[fa]+1;
father[x][0]=fa;
for(int i=1;i<=lg[dep[x]];i++)father[x][i]=father[father[x][i-1]][i-1];//倍增发更新祖先
for(int i=head[x];i;i=edge[i].ne)
{
int to=edge[i].to;
if(to!=fa)dfs(to,x);
}
}
在做好前期准备之后,就可以跳跃地求公共祖先了:
1.先是同暴力的一样调到同一深度, 设他们的层数相差 k层,那么最多跳跃 k个祖先,所以 x = f [ x ] [ l g [ d [ x ] - d [ y ] ]
2.在调到一层时,从最多开始跳跃,如果一样说明这是公共祖先,但不一定是最小的,然后在调小,如果不是了,那么两个点都向上一层,最后不断地调变小就可以找到了(可以用样例手推一下)
复杂度O(logn)
inline int LCA(int x,int y)
{
if(x==y)return x;
if(dep[x]<dep[y])swap(x,y);
while(dep[x]>dep[y])x=father[x][lg[dep[x]-dep[y]]];
if(x==y)return x;
for(int i=lg[dep[x]];i>=0;i--)
if(father[x][i]!=father[y][i])
x=father[x][i],y=father[y][i];
return father[x][0];
}