luoguP3379(倍增LCA)模板

https://www.luogu.org/problemnew/show/P3379
题目如上

关于LCA,利用二进制的特性解决问题,与快速幂有些相似;
须知定理:
任一实数n,有2^k+1 > n >= 2^k 那么n一定等于2^k , 2^k-1, 2^k-2, … 2^1, 2^0,中的任意几个的的和;

例如 9=2^3 + 2^0,
10=2^3 + 2^1;
11= 2^3 + 2^1 + 2^0;
那么求LCA可以不用一个一个往上跳了,利用这一规律便可以快速的找到LCA;

具体步骤如下:
1.将两个节点调整到同一高度
2.判断是否为同一节点
3.若2不成立将x,y一起往上跳,直到他们的父节点为同一个为止

f[i][j] 表示i节点向上跳2^i所到达的节点,先dfs求出f[i][0]即i的父节点,在递推求出其他;
递推公式如下
f[i][j]=f[f[i][j-1]][j-1],意思是i节点向上跳了2^j-1 后再向上跳2^j-1;
代码如下

// luogu-judger-enable-o2
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath> 
#include<vector>
#include<algorithm>
using namespace std;
vector<int> E[500010];
int deep[500010],f[500010][22];
int n,m,s;
void dfs(int now,int last){
	deep[now]=deep[last]+1;
	f[now][0]=last;
	for(int i=1;i*i<=deep[now];i++) f[now][i]=f[f[now][i-1]][i-1];
	for(int i=0;i<E[now].size();i++) if(!deep[E[now][i]]) dfs(E[now][i],now);
}
int Find_Lca(int x,int y){
	if(deep[x]<deep[y]) swap(x,y);//保证x较深
	for(int i=sqrt(deep[x]-deep[y])+1;i>=0;i--) if(deep[f[x][i]]>=deep[y]) x=f[x][i];//调整到同一高度
	if(x==y) return y;//情况2
	for(int i=sqrt(deep[x])+1;i>=0;i--) if(f[x][i]!=f[y][i]){x=f[x][i]; y=f[y][i];}
	return f[x][0];
}
int read(){
	int x=0,f=1;
	char ch=getchar();
	if(ch=='-') f=-1;
	while(ch<'0'||ch>'9') ch=getchar();
	while(ch>='0'&&ch<='9') x*=10,x+=ch-'0',ch=getchar();
	return x*f;
}
int main(){
	n=read(); m=read(); s=read();
	for(int i=1;i<=n-1;i++){
		int u=read(),v=read(); 
		E[u].push_back(v); E[v].push_back(u);
	}
	dfs(s,0);//预处理出f[i][0],和深度
	for(int i=1;i<=m;i++){
		int u=read(),v=read(); 
		printf("%d\n",Find_Lca(u,v));
	}
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值