倍增的模板,注释都在上面了

#include<iostream>
#include<cmath>
#define MAXN 500001
using namespace std;
int f[500005][100];
int n,m,s;
int dep[MAXN];
int tot;
int g[MAXN*2];
struct line
{
int to;
int next;
}edge[MAXN*2];
void addedge(int a,int b)
{
tot++;
edge[tot].to=b;
edge[tot].next=g[a];
g[a]=tot;
}
void dfs(int x) //记录深度 
{
for(int i=g[x];i;i=edge[i].next) //邻接表遍历点 
{
if(!dep[edge[i].to]) //如果这个深度没被更新,也就是没到过这个点 
{
dep[edge[i].to]=dep[x]+1;
f[edge[i].to][0]=x; //父节点 
dfs(edge[i].to);
}
}
}
void init() //倍增的预处理,f[i][j]就是点i往上走2^j的距离所到的点 
{
dep[s]=1; //根节点的深度为1 
dfs(s); //遍历深度 
for(int j=1;(1<<j)<=n;j++) //这里就是2的j次方要小于n
for(int i=1;i<=n;i++)
f[i][j]=f[f[i][j-1]][j-1];
}
int lca(int a,int b)
{
if(dep[a]>dep[b])swap(a,b); //让b是比较深的那个 
int cha=dep[b]-dep[a]; //计算深度差 
for(int i=0;cha;cha>>=1,i++) //这里就是把差转化成二进制,比如差值是5(101),那么就是先到f[b][0],然后再到f[b][2],这里差的末位对应的就是i 
if(cha&1)b=f[b][i];
if(a!=b) //开始倍增 
{
for(int i=log(n)/log(2);i>=0;i--) //换底公式,以二为底n的对数 
if(f[a][i]!=f[b][i])
{
a=f[a][i];
b=f[b][i];
}
a=f[a][0];
}
return a;
}
int main()
{
std::ios::sync_with_stdio(false);
cin>>n>>m>>s;
for(int i=1;i<n;i++)
{
int x,y;
cin>>x>>y;
addedge(x,y);
addedge(y,x);
}
init();
for(int i=0;i<m;i++)
{
int x,y;
cin>>x>>y;
cout<<lca(x,y)<<endl;
}
return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值