洛谷P2633 Count on a tree(主席树,倍增LCA,树上差分)
题目大意
就是给你一棵树,每个点都有点权,每次任意询问两点间路径上点权第k小的值(强制在线)。
思路分析
第k小。。。。。。又是主席树了。但这次变成树了,无法直接维护前缀和。
又是树上差分的小套路——每一个点到根的前缀和还是很好维护对吧。
询问\(u,v\)的时候,我们可以知道\(size[root,u]\)和\(size[root,v]\)的和。
但我们需要的只是一条路径,\(lca(u,v)\)以上的全不要,\(lca(u,v)\)也只要算一次。
于是用\(size[root,u]+size[root,v]-size[root,lca(u,v)]-size[root,father(lca(u,v))]\),也就是询问的时候四个点一起跳。
求LCA最方便的是倍增法(不会的百度一下),还有每个点对应的线段树从其父亲的线段树继承而来(根节点从\(0\)号空线段树继承而来),这两个操作我们在一次dfs建树时就可以一并处理完了。
话说我好久没打过倍增LCA了,老是写挂。。。。。。
闲话
另外,本蒟蒻听说倍增的二维数组把长度小的那一维度(即表示\(2^j\)的维度)开在前面会跑的快一些。
于是本蒟蒻亲自验证了一下,开在前面1112ms,开在后面992ms(没开O2,开了以后两个差不多)。。。。。。
这又是什么鬼?!
难道有评测机的不稳定因素?(洛谷一直很稳定啊!)
或者是有其他原因?欢迎各位大佬指教。
下面贴代码(开在前面):
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define R register int
#define G c=getchar()
#define in(Z) G;\
while(c<'-')G;\
Z=c&15;G;\
while(c>'-')Z*=10,Z+=c&15,G;
const int N=100009,M=4000009;
int P,he[N],ne[N<<1],to[N<<1],f[18][N],d[N],b[N];
int n,SZ,a[N],g[N],rt[N],lc[M],rc[M],s[M];
void build(R&u,R l,R r){//建初始空线段树
u=++P;
if(l!=r){
R m=(l+r)>>1;
build(lc[u],l,m);
build(rc[u],m+1,r);
}
}
inline void insert(R*u,R v,R k){//更新
R l=1,r=SZ,m;
while(l!=r){
s[*u=++P]=s[v]+1;
m=(l+r)>>1;
if(k<=m)r=m,rc[*u]=rc[v],u=&lc[*u],v=lc[v];
else l=m+1,lc[*u]=lc[v],u=&rc[*u],v=rc[v];
}
s[*u=++P]=s[v]+1;
}
inline int ask(R u,R v,R w,R x,R k){//询问,四个点一起搞
R l=1,r=SZ,m,q;
while(l!=r){
m=(l+r)>>1;
q=s[lc[u]]+s[lc[v]]-s[lc[w]]-s[lc[x]];
if(k<=q)r=m,u=lc[u],v=lc[v],w=lc[w],x=lc[x];
else l=m+1,u=rc[u],v=rc[v],w=rc[w],x=rc[x],k-=q;
}
return g[l];
}
void dfs(R u,R fa){//建好树,顺便预处理LCA
d[u]=d[fa]+1;
f[0][u]=fa;
for(R&i=b[u];f[i+1][u]=f[i][f[i][u]];++i);
insert(&rt[u],rt[fa],lower_bound(g+1,g+SZ+1,a[u])-g);
for(R i=he[u];i;i=ne[i])
if(to[i]!=fa)dfs(to[i],u);
}
inline int getlca(R u,R v){//查LCA
R i;
if(d[u]<d[v])i=u,u=v,v=i;
for(i=b[u];d[u]>d[v]&&i>=0;--i)
if(d[f[i][u]]>=d[v])u=f[i][u];
for(i=b[u];i>=0;--i)
if(f[i][u]!=f[i][v])u=f[i][u],v=f[i][v];
return u==v?u:f[0][u];//这里老是忘记判
}
int main(){
register char c;
R p=1,m,i,u,v,k,lca,lans=0;
in(n);in(m);
for(i=1;i<=n;++i){in(a[i]);}
memcpy(g,a,(n+1)<<2);//搞出来离散化
sort(g+1,g+n+1);
SZ=unique(g+1,g+n+1)-g-1;
build(rt[0],1,SZ);
for(i=1;i<n;++i)
{
in(u);in(v);
to[++p]=v;ne[p]=he[u];he[u]=p;
to[++p]=u;ne[p]=he[v];he[v]=p;//建边
}
dfs(1,0);
while(m--){
in(u);in(v);in(k);
lca=getlca(u^=lans,v);
printf("%d\n",lans=ask(rt[u],rt[v],rt[lca],rt[f[0][lca]],k));
}
return 0;
}
洛谷P2633 Count on a tree(主席树,倍增LCA,树上差分)的更多相关文章
- 洛谷P2633 Count on a tree(主席树上树)
题目描述 给定一棵N个节点的树,每个点有一个权值,对于M个询问(u,v,k),你需要回答u xor lastans和v这两个节点间第K小的点权.其中lastans是上一个询问的答案,初始为0,即第一个 ...
- 洛谷 P2633 Count on a tree 主席树
在一棵树上,我们要求点 $(u,v)$ 之间路径的第$k$大数. 对于点 $i$ ,建立 $i$ 到根节点的一棵前缀主席树. 简单容斥后不难得出结果为$sumv[u]+sumv[v]−sumv[l ...
- 洛谷P2633 Count on a tree 主席树
传送门:主席树 解题报告: 传送门! umm这题我还麻油开始做 所以 先瞎扯一波我的想法,如果错了我就当反面教材解释这种典型错误,对了我就不管了QwQ 就直接dfs,在dfs的过程中建树 然后就直接查 ...
- 洛谷 P2633 Count on a tree
P2633 Count on a tree 题目描述 给定一棵N个节点的树,每个点有一个权值,对于M个询问(u,v,k),你需要回答u xor lastans和v这两个节点间第K小的点权.其中last ...
- 【BZOJ-2588】Count on a tree 主席树 + 倍增
2588: Spoj 10628. Count on a tree Time Limit: 12 Sec Memory Limit: 128 MBSubmit: 3749 Solved: 873[ ...
- 洛谷P2633 Count on a tree(主席树,倍增LCA)
洛谷题目传送门 题目大意 就是给你一棵树,每个点都有点权,每次任意询问两点间路径上点权第k小的值(强制在线). 思路分析 第k小......又是主席树了.但这次变成树了,无法直接维护前缀和. 又是树上 ...
- ☆ [洛谷P2633] Count on a tree 「树上主席树」
题目类型:主席树+\(LCA\) 传送门:>Here< 题意:给出一棵树.每个节点有点权.问某一条路径上排名第\(K\)小的点权是多少 解题思路 类似区间第\(K\)小,但放在了树上. 考 ...
- 洛谷P2633 Count on a tree
题目描述 给定一棵N个节点的树,每个点有一个权值,对于M个询问(u,v,k),你需要回答u xor lastans和v这两个节点间第K小的点权.其中lastans是上一个询问的答案,初始为0,即第一个 ...
- 洛谷 P2633 Count on a tree 题解
题面 对于每个点建立一颗主席树: 然后按照树上差分的思想统计主席树的前缀和: lca+主席树+前向星存表就可以了: #include <bits/stdc++.h> #define inc ...
随机推荐
- 安装Docker的三种方式
本人在安装docker时之前一直采用的是系统自带的docker软件包安装,导致下载的docker不是最新版本,因此会有很多docker命令无法使用,例如network等等,现将安装docker的方式总 ...
- java web classpath
在使用ssh等框架开发web程序时配置文件(xml和properties)存放的路径一般为src下,当部署程序时则必须存在于classes路径下,具体如下 src不是classpath, WEB-IN ...
- Vivado中xilinx_courdic IP核(求exp指数函数)使用
由于Verilog/Vhdl没有计算exp指数函数的库函数,所以在开发过程中可利用cordic IP核做exp函数即e^x值: 但前提要保证输入范围在(-pi/4—pi/4) 在cordic核中e^x ...
- Spring Data JPA、MyBatis还有Hibernate有什么区别
原文:https://www.imooc.com/article/19754?block_id=tuijian_wz Spring Data JPA.MyBatis还有Hibernate有什么区别 2 ...
- 大数据入门第二十五天——elasticsearch入门
一.概述 推荐路神的ES权威指南翻译:https://es.xiaoleilu.com/010_Intro/00_README.html 官网:https://www.elastic.co/cn/pr ...
- 20155233 《网络对抗技术》EXP3 免杀原理与实践
正确使用msf编码器,msfvenom生成如jar之类的其他文件,veil-evasion,自己利用shellcode编程等免杀工具或技巧 使用msf编码器生成jar包 输入命令msfvenom -p ...
- cli 开发记录
最近要开发一个 cli,主要作用是方便同事生成前端项目,做了一天半,基本参考的是 vue-cli. cli 要实现的功能: 用 cnpm install zt-cli -g 全局安装,这个就要将你做的 ...
- eclipse中设置项目的编码方式
1.windows->Preferences...打开"首选项"对话框,左侧导航树,导航到general->Workspace,右侧Text file encoding ...
- 二叉树 c++
树 非空树 有一个(root)根节点r 其余节点可分为m个互不相交的有限集(子树)T1....Tm 具有n个节点的树,具有(n-1)条连接(指针域),需要构成结构体,尽可能减少空间域的浪费,使用儿子兄 ...
- 微软职位内部推荐-Software Engineer-Sharepoint
微软近期Open的职位: SharePoint is a multi-billion dollar enterprise business that has grown from an on-prem ...