题目戳我

\(\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的更多相关文章

  1. 题解 LNOI2014 LCA

    题目:传送门 这道题根本不用lca,也没有部分分... 考虑求两个点xy的lca的深度. 我们将x到树根所有点的值都加1,然后查询y到根的和,其实就是lca的深度. 所以本题离线一下上树剖乱搞就可以了 ...

  2. 【BZOJ3626】[LNOI2014]LCA 离线+树链剖分+线段树

    [BZOJ3626][LNOI2014]LCA Description 给出一个n个节点的有根树(编号为0到n-1,根节点为0).一个点的深度定义为这个节点到根的距离+1.设dep[i]表示点i的深度 ...

  3. BZOJ 3626: [LNOI2014]LCA [树链剖分 离线|主席树]

    3626: [LNOI2014]LCA Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 2050  Solved: 817[Submit][Status ...

  4. bzoj 3626 [LNOI2014]LCA(离线处理+树链剖分,线段树)

    3626: [LNOI2014]LCA Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1272  Solved: 451[Submit][Status ...

  5. bzoj 3626: [LNOI2014]LCA 离线+树链剖分

    3626: [LNOI2014]LCA Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 426  Solved: 124[Submit][Status] ...

  6. BZOJ 3626: [LNOI2014]LCA( 树链剖分 + 离线 )

    说多了都是泪啊...调了这么久.. 离线可以搞 , 树链剖分就OK了... -------------------------------------------------------------- ...

  7. 洛谷 P4211 [LNOI2014]LCA 解题报告

    [LNOI2014]LCA 题意 给一个\(n(\le 50000)\)节点的有根树,询问\(l,r,z\),求\(\sum_{l\le i\le r}dep[lca(i,z)]\) 一直想启发式合并 ...

  8. P4211 [LNOI2014]LCA

    P4211 [LNOI2014]LCA 链接 分析: 首先一种比较有趣的转化是,将所有点到1的路径上都+1,然后z到1的路径上的和,就是所有答案的deep的和. 对于多次询问,要么考虑有把询问离线,省 ...

  9. P4211 [LNOI2014]LCA LCT

    P4211 [LNOI2014]LCA 链接 loj luogu 思路 多次询问\(\sum\limits_{l \leq i \leq r}dep[LCA(i,z)]\) 可以转化成l到r上的点到根 ...

  10. [BZOJ3626] [LNOI2014]LCA(树链剖分)

    [BZOJ3626] [LNOI2014]LCA(树链剖分) 题面 给出一棵N个点的树,要求支持Q次询问,每次询问一个点z与编号为区间[l,r]内的点分别求最近公共祖先得到的最近公共祖先深度和.N, ...

随机推荐

  1. Linux文件描述符与重定向

    文件描述符可以理解为linux跟踪打开文件,而分配的一个数字,这个数字有点类似c语言操作文件时候的句柄,通过句柄就可以实现文件的读写操作. 当Linux启动的时候会默认打开三个文件描述符,分别是: 标 ...

  2. LWPR

    Scriptable Render Pipeline https://docs.unity3d.com/Manual/ScriptableRenderPipeline.html     Unity轻量 ...

  3. vue <compent> 结合is实现动态切换子组件

    这里就很自然的产生一个新问题:如何给相应的子组件传props?

  4. LiteOS间歇计算技术:IOT终端真正感受“电量自由”

    摘要:LiteOS间歇计算为什么能使物联网终端实现长续航? 物联网设备.场景复杂多样,小到智能穿戴的耳机,大到大型基建设备,更有我们陌生而又熟悉场景,例如深海探测.森林监控.野生动物跟踪等等能量采集场 ...

  5. JAVA JDK 环境变量配置 入门详解 - 精简归纳

    JAVA JDK 环境变量配置 入门详解 - 精简归纳 JERRY_Z. ~ 2020 / 9 / 13 转载请注明出处!️ 目录 JAVA JDK 环境变量配置 入门详解 - 精简归纳 一.为什么j ...

  6. App测试理论简介

    一.App测试常见关注点 1.App的功能测试 功能测试都是我们首要测试的,只有功能实现了才算符合上线发布的最低标准.我们需要检测产品功能是否已实现.产品功能是否符合设计要求.产品功能是否有重复.产品 ...

  7. [LeetCode]647. 回文子串(DP)

    ###题目 给定一个字符串,你的任务是计算这个字符串中有多少个回文子串. 具有不同开始位置或结束位置的子串,即使是由相同的字符组成,也会被计为是不同的子串. 示例 1: 输入: "abc&q ...

  8. 深入理解python语言

    2008年,安卓操作系统诞生:PC时代向移动时代转换 互联网,视窗 2017/5/27柯洁最终0:3AlphaGo 计算机技术的演进过程 不同编程语言的设计初心和适用对象 C语言核心解决的是性能问题, ...

  9. JMeter5.0在windows(含插件安装)

    一.jmeter下载 前提:已经安装jdk8+ jmeter下载地址:http://jmeter.apache.org/download_jmeter.cgi 有Binaries和Source版本 前 ...

  10. SVN检出maven项目

    (一)直接单击项目,右键选择configure,选择convert to maven project (二)删除project explorer中的项目,并重新从工作区间导入maven项目.