倍增法求lca

 f[i][j]表示从i这个节点出发,向上走2^j步到达的点(超过了最大深度就返回0)

那么显然有f[i][j]=f[f[i][j-1]][j-1]

因为相当于是从i点先向上走2^(j-1)步,再走2^(j-1)步。等价于一共走2^j步。

 

然后求lca就是:先把深度大的那个点往上跳,使两个点的深度相同。

对应的代码:

for(int i=18;i>=0;i--) if(dep[f[x][i]]>=dep[y]) x=f[x][i];

 可以把两点之间的深度差理解为一个二进制数。

举个例子:一个点深度为3,另一个点深度为16。那么深度为16的点向上跳13步它们的深度就相同了。

13转换成二进制为1101。那么这个跳的过程就是先跳(1000)步,再跳(100)步,

然后跳(10)步,发现深度超了,就不会跳。

最后跳(1)步,它们深度就相同了。

 

之后特判一下这时是否相同。

如果不相同:它们就一起往上跳。

如果发现  f[x][i]==f[y][i]  ,那么说明它们跳到了lca或者是超过了lca,就不跳。(当前深度≤lca深度)

但是如果  f[x][i]!=f[y][i]  ,那么它们一定还都在lca下面,就继续跳。(当前深度>lca深度)

for(int i=18;i>=0;i--) if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];

这里跳完之后,可以发现,它们一定刚好跳到了lca的下面一个。(可以把深度想成二进制思考。)

那么它们的父亲就是lca。

void dfs(int u,int f){
	dep[u]=dep[f]+1,f[u][0]=f;
	for(int i=Head[u];i;i=Next[i])
		dfs(V[i],u);
}
void pre(){
	dfs(1,0);
	for(int i=1;(1<<i)<=N;++i)
		for(int u=1;u<=N;++u)
			f[u][i]=f[f[u][i-1]][i-1];
}
int lca(int x,int y){
	if(dep[x]<dep[y]) swap(x,y);
	for(int i=18;i>=0;i--) if(dep[f[x][i]]>=dep[y]) x=f[x][i];
	if(x==y) return x;
	for(int i=18;i>=0;i--) if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];
	return f[x][0];
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值