Vijos lxhgww的奇思妙想--求K级祖先
给出一棵树求K级祖先。O(N*logN+Q)
更详细的讲解见:https://www.cnblogs.com/cjyyb/p/9479258.html
/*
要求k级祖先,我们可以把k拆成"2^highbit(x)+tmp 形式
(highbit(x)为x在二进制位下的最高位),然后用倍增的方法把highbit(x)的部分跳了
剩下tmp的同样可以预处理掉,这样预处理就是O(n*logn)的效率,
所以对于每个询问就是O(1)回答,这样的效率就是O(n*logn+q)。
于是就考虑用长链剖分。
讲讲具体的细节操作。
数组:
f[i][j]:i的倍增得到的祖先
d[i]:i i的最深的儿子的深度(用于处理len lenlen)
dep[i]:i的节点深度
son[i]:i的长儿子
len[i]:i为长链顶点的链长
top[i]:i所在长链的顶点
hb[i]:i的highbit highbithighbit值
twi[i]:2^i
预处理:
对于highbit(x) 选择你喜欢的方法做。
对于tmp,用两个动态数组维护,存每条长链的顶点的向上len
向下len的祖先和儿子
(对于每条长链,存链中的节点显然无意义,因为你可以在询问时先跳到top)。
深搜:
一共dfs两遍
dfs1:预处理一些关于树的常规信息,如深度、长儿子
dfs2:处理链长和长链顶点
原文链接:https://blog.csdn.net/hzq_oi/article/details/88595343
*/
#include<bits/stdc++.h>
#define N 300005
using namespace std;
inline int rd()
{
int data=0,w=1;static char ch=0;
while(!isdigit(ch)&&ch!='-')ch=getchar();
if(ch=='-')w=-1,ch=getchar();
while(isdigit(ch))data=(data<<1)+(data<<3)+ch-'0',ch=getchar();
return data*w;
}
int first[N],cnt;
struct node
{int v,nxt;}
e[N<<1];
inline void add(int u,int v)
{
e[++cnt].v=v;
e[cnt].nxt=first[u];
first[u]=cnt;
}
int n,m,ans;
int f[N][20],d[N],dep[N],len[N],top[N],son[N],hb[N],twi[30];
void dfs1(int u,int fa) //根的深度为1
{
d[u]=dep[u]=dep[fa]+1;
f[u][0]=fa;
for(int register i=1;i<=19;i++)
f[u][i]=f[f[u][i-1]][i-1];
for(int register i=first[u];i;i=e[i].nxt)
{
int register v=e[i].v;
if(v==fa)continue;
dfs1(v,u);
if(d[u]<d[v]) //找出重儿子
{
d[u]=d[v];
son[u]=v;
}
}
}
void dfs2(int u,int fa)
{
len[u]=d[u]-dep[top[u]]+1;
if(!son[u])return;
top[son[u]]=top[u];
dfs2(son[u],u);
for(int register i=first[u];i;i=e[i].nxt)
if(e[i].v!=son[u]&&e[i].v!=fa)
top[e[i].v]=e[i].v,dfs2(e[i].v,u);
}
vector<int>up[N],down[N];
inline int query(int x,int k)//求k级祖先
{
if(k>=dep[x])return 0;
if(!k)return x;
x=f[x][hb[k]];//hb[k]代表k这个数字的最高位是2的多少次方
k^=twi[hb[k]];
if(!k)return x;
int register tmp=dep[x]-dep[top[x]];
if(tmp==k)
return top[x];
else
if(tmp<k)return
up[top[x]][k-tmp-1];
//x在某条轻链上,跳到top点后,还不够,不要向上跳
else
return down[top[x]][tmp-k-1];
//x在重链上,跳到top点后,跳过头了,所以还要向下移动下
}
int main()
{
n=rd();
for(int register i=1;i<n;i++)
{int register x=rd(),y=rd();add(x,y);add(y,x);}
dfs1(1,0);
top[1]=1;
dfs2(1,0); for(int i=1;i<=n;i++)
{
if(i!=top[i]) continue;
//找出每条重链的顶点
int register tmp=0,prv=i;
while(tmp<len[i]&&prv)
{
prv=f[prv][0],up[i].push_back(prv),tmp++;
}
tmp=0;prv=i;
while(tmp<len[i])
{
prv=son[prv],down[i].push_back(prv),tmp++;
//向下跳,跳到它的重儿子 }
}
twi[0]=1;
for(int register i=1;i<=20;i++)
twi[i]=twi[i-1]<<1;
for(int register i=1;i<=n;i++)
for(int register j=20;j>=0;j--)
if(twi[j]&i)
{
hb[i]=j;
//对于数字i来说,它转成2进制后最高位是2^j
//例如7的最高位是2^2,8的是2^3
break;
}
m=rd();
while(m--)
{
int register x=rd(),k=rd();
x^=ans;
k^=ans;
ans=query(x,k);
printf("%d\n",ans);
}
return 0;
}
Vijos lxhgww的奇思妙想--求K级祖先的更多相关文章
- [vijos]lxhgww的奇思妙想(长链剖分)
题意 题目链接 Sol 长链剖分 又是一个用各种花式技巧优化的暴力 它的主要思想是:对于每个节点,把深度最深的子节点当做重儿子,它们之间的边当做重边 这样就会有一些非常好的轻质 所有链长总和是\(O( ...
- 2019.01.06 vijos lxhgww的奇思妙想(长链剖分)
传送门 长链剖分模板题. 题意简述:允许O(nlogn)O(nlog_n)O(nlogn)预处理,让你支持O(1)O(1)O(1)查找任意一个点的kkk级祖先. 思路:因为要O(1)O(1)O(1) ...
- Vijos.lxhgww的奇思妙想(k级祖先 长链剖分)
题目链接 https://blog.bill.moe/long-chain-subdivision-notes/ http://www.cnblogs.com/zzqsblog/p/6700133.h ...
- 【Vijos】lxhgww的奇思妙想
题面 题解 求$k$级祖先孙子 为什么要用长链剖分啊??? 倍增并没有慢多少... 其实是我不会 长链剖分做这道题还是看这位巨佬的吧. 代码 #include<bits/stdc++.h> ...
- lxhgww的奇思妙想 长链剖分板子
https://vijos.org/d/Bashu_OIers/p/5a79a3e1d3d8a103be7e2b81 求k级祖先,预处理nlogn,查询o1 //#pragma GCC optimiz ...
- 「vijos」lxhgww的奇思妙想(长链剖分)
传送门 长链剖分的板子(又是乱搞优化暴力) 对于每一个点,我们定义它深度最深的子节点为它的重儿子(为什么不叫长儿子……),他们之间的连边为重边 然后长链剖分有几个性质 1.总链长为$O(n)$ 2.一 ...
- 【Vijos】lxhgww的奇思妙想(长链剖分)
题面 给定一棵树,每次询问一个点的\(k\)次祖先,强制在线. Vijos 题解 长链剖分. 链接暂时咕咕咕了. 现在可以戳链接看题解了 #include<iostream> #inclu ...
- 【LCA求最近公共祖先+vector构图】Distance Queries
Distance Queries 时间限制: 1 Sec 内存限制: 128 MB 题目描述 约翰的奶牛们拒绝跑他的马拉松,因为她们悠闲的生活不能承受他选择的长长的赛道.因此他决心找一条更合理的赛道 ...
- POJ 1986 Distance Queries (Tarjan算法求最近公共祖先)
题目链接 Description Farmer John's cows refused to run in his marathon since he chose a path much too lo ...
随机推荐
- python算法(一)基本知识&冒泡排序&选择排序&插入排序
本节内容: 算法基本知识 冒泡排序 选择排序 插入排序 1. 算法基本知识 1.1 什么是算法? 算法(algorithm):就是定义良好的计算过程,他取一个或一组的值为输入,并产生出一个或一组值作为 ...
- 微信支付之扫码支付、公众号支付、H5支付、小程序支付相关业务流程分析总结
前言 很久以来,一直想写一篇微信支付有关的总结文档:一方面是总结自己的一些心得,另一方面也可以帮助别人,但是因种种原因未能完全理解透彻微信支付的几大支付方式,今天有幸做一些总结上的文章,也趁此机会,将 ...
- STM32CubeIDE Debug Configurations
新建完工程并编译后,配置Debug Configurations 此时双击STM32 MCU Debugging,如下图 此时就可以下载调试生成的LED.elf文件了
- ZROI 19.08.05模拟赛
传送门 写在前面:为了保护正睿题目版权,这里不放题面,只写题解. A \(21pts:\) 随便枚举,随便爆搜就好了. \(65pts:\) 比较显然的dp,设\(f_{i,j,k}\)表示在子树\( ...
- mysql——非主键自增
今天遇到一个问题: 要创建一张表,其中我想将ip和date这两列作为一个复合主键,可以通过如下语句修改表结构: alter table tb_name add primary key (字段1,字段2 ...
- 更改pip源地址为阿里云
1.在用户名目录创建pip目录,在pip目录下创建pip.ini. 2.pip.ini中输入: [global] index-url = http://mirrors.aliyun.com/pypi/ ...
- UI编辑
UI编辑 基本部件介绍 (1)Layout(布局) (2)Space(空间) (3)Button (4)ItemView (5)ItemWidget Widget继承自View,即ListWidget ...
- 网页结构树DOM
引入 window对象 所有浏览器都支持 window 对象.它表示浏览器窗口. *如果文档包含框架(frame 或 iframe 标签),浏览器会为 HTML 文档创建一个 window 对象,并为 ...
- windows 安装 Mongodb 数据库及操作图形化软件 Robo 3T
1 下载系统对应的正确 Mongodb 和 Robo 3T 版本 2 选中 Mongodb 需要安装的路径(后续会使用路径) 3 启动 Mongodb 服务器(到安装相关的路径) 可以参考 菜鸟教程 ...
- Mongodb副本集集群搭建
一.环境准备 1.1.主机信息(机器配置要求见硬件及开发标准规范文档V1.0) 序号 主机名 IP 1 DB_01 10.202.105.52 2 DB_02 10.202.105.53 3 DB_0 ...