倍增的模板,注释都在上面了
来源:互联网 发布:淘宝退货纸条怎么写 编辑:程序博客网 时间:2024/05/20 02:54
#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;
}
#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
- 倍增的模板,注释都在上面了
- App服务化, 10倍增长,你想知道的都在这里了!
- 倍增求LCA的模板
- 如何在数据库字段上面加注释
- 我在CSDN上面的Blog开张了……
- 我的blog已经放在MSN空间上面了.
- :)在CSDN上面有自己的窝了
- 今天第一次在csdn上面来发表的问题了
- 开始在CSDN上面写自己的博客了
- 倍增法求LCA的模板
- Android--在显示ImageView时图片上面和下面都出现一段空白区间的解决办法
- 【Android】在显示ImageView时图片上面和下面都出现一段空白区间的解决办法
- Xcode 修改文件上面的注释信息
- 后缀数组的倍增模板与DC3模板
- 在Eclipse中创建java类的注释模板
- 在eclipse中配置属于自己的注释模板
- 怎样在Android Studio中进行注释模板的设置
- eclipse的注释模板
- [LeetcCode]Top 100 Liked Questions(1-10)
- Faster-RCNN遇到的问题
- IOS开发中,SVN如何恢复到某一个版本(以Cornerstone为例)
- Android APP内部跳转到微信
- Redis和Memcached的区别
- 倍增的模板,注释都在上面了
- 深度学习中momentum的作用
- python设计模式(备忘录模式)
- 原版Win7注入USB3.0驱动和NVME驱动教程
- solr搜索关键字高亮显示
- React组件规范
- 6-17 Shortest Path [2](25 point(s))(dijkstra)
- KNN算法-小明逆袭的故事(1)
- java仿QQ截图带窗体检测