【题解】[LNOI2014]LCA
\(\text{Solution:}\)
这题的转化思想很巧妙……
考虑把\(dep\)给拆掉。
首先可以明确的是,每一个\(LCA\)一定在\(root\to z\)的路径上。
考虑一个\(LCA\)被选中,意味\(root\to LCA\)这个路径上所有的值都\(+1\).
于是,我们可以把询问看成将点\(p\in [l,r],add(root\to p,1)\)并从\(z\)向根的路径上求和即可。
这个区间修改,区间查询显然的树剖线段树板子。而若对于每一个询问都做一次,复杂度显然不够。
考虑将询问离线。
将一个区间询问\([l,r]\)拆成区间\([l,l-1] and [1,r]\)并用差分的思想进行处理。于是,对每一个拆开的区间按照右端点排序(因为左端点都是\(1\))进行离线处理即可,这样就避免了每一次的清空操作。
总结:对于一个\(dep\)的询问,我们可以看成:
从\(i\)向根的路径上进行区间\(+1\).
从\(z\)向根的路径上进行区间求和。
这样,我们可以巧妙地通过区间加来规避对\(dep\)的询问。
笔者一开始硬挖掘编号连续的性质往线段树维护深度和的歪路上走了很久……一度自闭
值得记录的好题。
#include<bits/stdc++.h>
using namespace std;
#define int long long
const int MAXN=5e5+10;
const int mod=201314;
int dep[MAXN],siz[MAXN],top[MAXN],head[MAXN],rt;
int tot,cnt,n,m,rk[MAXN],id[MAXN],dfstime,ls[MAXN];
int rs[MAXN],pa[MAXN],son[MAXN],qcnt,Ans[MAXN];
struct E{int nxt,to;}e[MAXN];
inline void link(int x,int y){e[++tot]=(E){head[x],y};head[x]=tot;}
inline int add(int x,int y){return (x+y)%mod;}
inline int mul(int x,int y){return 1ll*x*y%mod;}
struct SGT{int l,r,sum,tag;}tr[MAXN];
inline void pushup(int x){tr[x].sum=tr[ls[x]].sum+tr[rs[x]].sum;}
inline int LL(int x){return tr[x].r-tr[x].l+1;}
inline void pushdown(int x){
if(tr[x].tag){
int p=tr[x].tag;tr[x].tag=0;
tr[ls[x]].tag+=p;tr[rs[x]].tag+=p;
tr[ls[x]].tag%=mod;tr[rs[x]].tag%=mod;
tr[ls[x]].sum=add(tr[ls[x]].sum,mul(LL(ls[x]),p));
tr[rs[x]].sum=add(tr[rs[x]].sum,mul(LL(rs[x]),p));
}
}
void build(int l,int r,int x){
tr[x].l=l,tr[x].r=r;
if(l==r)return;
int mid=(l+r)>>1;
ls[x]=cnt++;rs[x]=cnt++;
build(l,mid,ls[x]);
build(mid+1,r,rs[x]);
pushup(x);
}
void dfs1(int x,int fa){
pa[x]=fa,siz[x]=1,dep[x]=dep[fa]+1;
for(int i=head[x];i;i=e[i].nxt){
int j=e[i].to;
if(j==fa)continue;
dfs1(j,x);siz[x]+=siz[j];
if(siz[j]>siz[son[x]])son[x]=j;
}
}
void dfs2(int x,int t){
top[x]=t,rk[id[x]=++dfstime]=x;
if(!son[x])return;dfs2(son[x],t);
for(int i=head[x];i;i=e[i].nxt){
int j=e[i].to;
if(j!=pa[x]&&j!=son[x])dfs2(j,j);
}
}
void update(int l,int r,int v,int x){
if(tr[x].l>=l&&tr[x].r<=r){
tr[x].tag=add(tr[x].tag,v);
tr[x].sum=add(tr[x].sum,mul(LL(x),v));
return;
}
int mid=(tr[x].l+tr[x].r)>>1;pushdown(x);
if(l<=mid)update(l,r,v,ls[x]);
if(mid<r)update(l,r,v,rs[x]);
pushup(x);
}
int query(int l,int r,int x){
if(tr[x].l>=l&&tr[x].r<=r)return tr[x].sum;
int T=0,mid=(tr[x].l+tr[x].r)>>1;pushdown(x);
if(l<=mid)T=add(T,query(l,r,ls[x]));
if(mid<r)T=add(T,query(l,r,rs[x]));
pushup(x);return T;
}
void updatechain(int x,int y,int v){
while(top[x]!=top[y]){
if(dep[x]<dep[y])swap(x,y);
update(id[top[x]],id[x],v,rt);
x=pa[top[x]];
}
if(id[x]>=id[y])swap(x,y);
update(id[x],id[y],v,rt);
}
int querychain(int x,int y){
int ans=0;
while(top[x]!=top[y]){
if(dep[x]<dep[y])swap(x,y);
ans=add(ans,query(id[top[x]],id[x],rt));
x=pa[top[x]];
}
if(id[x]>id[y])swap(x,y);
ans=add(ans,query(id[x],id[y],rt));
return ans;
}
struct Query{int pos,num,fg,z;}q[MAXN];
inline bool cmp(Query A,Query B){return A.pos<B.pos;}
pair<int,int>ans[MAXN];
signed main(){
scanf("%lld%lld",&n,&m);
for(int i=2;i<=n;++i){
int x;
scanf("%lld",&x);
x++;
link(x,i);link(i,x);
}
dfs1(1,0);dfs2(1,1);rt=cnt++;build(1,n,rt);
for(int i=1;i<=m;++i){
int a,b,c;
scanf("%lld%lld%lld",&a,&b,&c);
a++,b++,c++;
q[++qcnt]=(Query){a-1,i,0,c};
q[++qcnt]=(Query){b,i,1,c};
}
sort(q+1,q+qcnt+1,cmp);
int now=0;
for(int i=1;i<=qcnt;++i){
while(now<q[i].pos)updatechain(++now,1,1);
if(!q[i].fg)ans[q[i].num].first=querychain(q[i].z,1);
else ans[q[i].num].second=querychain(q[i].z,1);
}
for(int i=1;i<=m;++i)Ans[i]=(ans[i].second-ans[i].first+mod)%mod;
for(int i=1;i<=m;++i)printf("%lld\n",Ans[i]);
return 0;
}
【题解】[LNOI2014]LCA的更多相关文章
- 题解 LNOI2014 LCA
题目:传送门 这道题根本不用lca,也没有部分分... 考虑求两个点xy的lca的深度. 我们将x到树根所有点的值都加1,然后查询y到根的和,其实就是lca的深度. 所以本题离线一下上树剖乱搞就可以了 ...
- 【BZOJ3626】[LNOI2014]LCA 离线+树链剖分+线段树
[BZOJ3626][LNOI2014]LCA Description 给出一个n个节点的有根树(编号为0到n-1,根节点为0).一个点的深度定义为这个节点到根的距离+1.设dep[i]表示点i的深度 ...
- BZOJ 3626: [LNOI2014]LCA [树链剖分 离线|主席树]
3626: [LNOI2014]LCA Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2050 Solved: 817[Submit][Status ...
- bzoj 3626 [LNOI2014]LCA(离线处理+树链剖分,线段树)
3626: [LNOI2014]LCA Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1272 Solved: 451[Submit][Status ...
- bzoj 3626: [LNOI2014]LCA 离线+树链剖分
3626: [LNOI2014]LCA Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 426 Solved: 124[Submit][Status] ...
- BZOJ 3626: [LNOI2014]LCA( 树链剖分 + 离线 )
说多了都是泪啊...调了这么久.. 离线可以搞 , 树链剖分就OK了... -------------------------------------------------------------- ...
- 洛谷 P4211 [LNOI2014]LCA 解题报告
[LNOI2014]LCA 题意 给一个\(n(\le 50000)\)节点的有根树,询问\(l,r,z\),求\(\sum_{l\le i\le r}dep[lca(i,z)]\) 一直想启发式合并 ...
- P4211 [LNOI2014]LCA
P4211 [LNOI2014]LCA 链接 分析: 首先一种比较有趣的转化是,将所有点到1的路径上都+1,然后z到1的路径上的和,就是所有答案的deep的和. 对于多次询问,要么考虑有把询问离线,省 ...
- P4211 [LNOI2014]LCA LCT
P4211 [LNOI2014]LCA 链接 loj luogu 思路 多次询问\(\sum\limits_{l \leq i \leq r}dep[LCA(i,z)]\) 可以转化成l到r上的点到根 ...
- [BZOJ3626] [LNOI2014]LCA(树链剖分)
[BZOJ3626] [LNOI2014]LCA(树链剖分) 题面 给出一棵N个点的树,要求支持Q次询问,每次询问一个点z与编号为区间[l,r]内的点分别求最近公共祖先得到的最近公共祖先深度和.N, ...
随机推荐
- Android java程序员必备技能,集合与数组中遍历元素,增强for循环的使用详解及代码
Android java程序员必备技能,集合与数组中遍历元素, 增强for循环的使用详解及代码 作者:程序员小冰,CSDN博客:http://blog.csdn.net/qq_21376985 For ...
- Google谷歌在根据流量统计分析当年的2008年汶川大地震
这是一张2008年的老图,Google当时的博文说道:"当我们依照惯例整理和分析谷歌搜索引擎的流量数据时,一条从未见过的曲线出现在我们面前.当意识到发生了什么事情时,我们的眼睛湿润了.&qu ...
- RVO+CA
http://gamma.cs.unc.edu/RVO/ http://gamma.cs.unc.edu/CA/ https://arongranberg.com/astar/docs/writing ...
- 如何编写高质量的C#代码(一)
从"整洁代码"谈起 一千个读者,就有一千个哈姆雷特,代码质量也同样如此. 想必每一个对于代码有追求的开发者,对于"高质量"这个词,或多或少都有自己的一丝理解.当 ...
- 跟着兄弟连系统学习Linux-【day08】
day08-20200605 p27.软件包管理简 windows 和 linux 软件是不同的版本. Linux源码包,开源的.绝大部分都是C语言写的.源码包安装速度比较慢.需要先编译后再安装.脚本 ...
- 读取文本文件中的中文打印到Eclipse控制台为何显示问号
原因:未将文本文件存为utf-8编码格式而是ascii编码格式.
- 支持MySQL数据库的agumaster版本
下载地址:https://files.cnblogs.com/files/xiandedanteng/agumaster20200501.zip
- IT人35岁危机:到底是因为爱还是责任?
互联网蚕食世界,未来属于IT人. 这是属于互联网的时代,每个人都是网络的弄潮儿,由于网络越来越被需要,互联网IT行业被推上了较高的位置,这也导致IT行业的就业环境火热,越来越多的人被IT行业广阔的发展 ...
- Vue企业级优雅实战04-组件开发01-SVG图标组件
(后续的文章 公众号会提前一周更新,欢迎关注文末的微信公众号:程序员搞艺术) 预览本文的实现效果: # gitee git clone git@gitee.com:cloudyly/dscloudy- ...
- 学习 | css3基本动画之demo篇
移动端使用的框架是zepto,但是zepto的内置对象没有传统的animate这个方法,效果都是需要css3来实现的,zepto也不支持fadeIn和fadeOut等一些基本的动画,基于这一现状,我自 ...