【LNOI2014】LCA
题面
题解
考察\(dep[\mathrm{LCA}(i, x)]\)的性质,发现它是\(i\)和\(x\)的链交的长度。
那么对每个\(i\)所在的链打一个区间加标记,询问时算一下\(x\)所在的链的区间和即可。
如果有\(l \leq i \leq r\)的限制,就进行离线处理即可。
代码
好久之前的代码,有点丑见谅。
#include<bits/stdc++.h>
#define RG register
#define file(x) freopen(#x".in", "r", stdin);freopen(#x".out", "w", stdout);
#define for_edge(i, x) for(RG int i=head[x];i;i=e[i].next)
#define clear(x, y) memset(x, y, sizeof(x));
using namespace std;
template<typename T = int>
inline T read()
{
T data=0, w=1;
char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
if(ch=='-') w=-1, ch=getchar();
while(ch>='0'&&ch<='9') data=(data<<3)+(data<<1)+(ch^48), ch=getchar();
return data*w;
}
const int maxn(50010), mod(201314);
struct edge { int next, to; } e[maxn];
int head[maxn], e_num, n, q;
inline void add_edge(int from, int to) { e[++e_num]={head[from], to}; head[from]=e_num; }
int fa[maxn], heavy[maxn], size[maxn];
void dfs(int x)
{
size[x]=1;
int _max=0;
for_edge(i, x)
{
int to=e[i].to; dfs(to);
size[x]+=size[to];
if(size[to]>_max) _max=size[to], heavy[x]=to;
}
}
int pos[maxn], belong[maxn], cnt_node;
void dfs(int x, int chain)
{
pos[x]=++cnt_node;
belong[x]=chain;
int k=heavy[x];
if(!k) return;
dfs(k, chain);
for_edge(i, x)
{
int to=e[i].to;
if(to==k) continue;
dfs(to, to);
}
}
int sum[maxn<<2], lazy[maxn<<2];
#define son(i) ((root<<1)|i)
inline void pushdown(int root, int l, int r)
{
if(l==r) lazy[root]=0;
if(!lazy[root]) return;
int mid(l+r>>1);
sum[son(0)]+=(mid-l+1)*lazy[root];
sum[son(1)]+=(r-mid)*lazy[root];
lazy[son(0)]+=lazy[root];
lazy[son(1)]+=lazy[root];
lazy[root]=0;
}
inline void update(int l, int r, int root=1, int nl=1, int nr=cnt_node)
{
if(nr<l || nl>r) return;
if(l<=nl && nr<=r)
{
sum[root]+=nr-nl+1;
lazy[root]++;
return;
}
int mid(nl+nr>>1);
pushdown(root, nl, nr);
update(l, r, son(0), nl, mid); update(l, r, son(1), mid+1, nr);
sum[root]=sum[son(0)]+sum[son(1)];
}
inline void update_chain(int x) { while(x) update(pos[belong[x]], pos[x]), x=fa[belong[x]]; }
inline int query(int l, int r, int root=1, int nl=1, int nr=cnt_node)
{
if(nr<l || nl>r) return 0;
if(l<=nl && nr<=r) return sum[root];
int mid(nl+nr>>1);
pushdown(root, nl, nr);
return query(l, r, son(0), nl, mid)+query(l, r, son(1), mid+1, nr);
}
inline int query_chain(int x)
{
int ans=0;
while(x) ans+=query(pos[belong[x]], pos[x]), x=fa[belong[x]];
return ans;
}
vector<int> le[maxn], ri[maxn];
int ans[maxn], que[maxn];
int main()
{
n=read(); q=read();
for(RG int i=2;i<=n;i++) fa[i]=read()+1, add_edge(fa[i], i);
dfs(1); dfs(1, 1);
RG int l, r, z;
for(RG int i=1;i<=q;i++) l=read(), r=read(), z=read(), le[l].push_back(i), ri[r+1].push_back(i), que[i]=z+1;
for(RG int i=1;i<=n;i++)
{
update_chain(i);
for(RG auto j : ri[i]) (ans[j]+=query_chain(que[j]))%=mod;
for(RG auto j : le[i]) ans[j]=(ans[j]-query_chain(que[j])+mod)%mod;
}
for(RG int i=1;i<=q;i++) printf("%d\n", (ans[i]+mod)%mod);
return 0;
}
【LNOI2014】LCA的更多相关文章
- bzoj3626【LNOI2014】LCA
3626: [LNOI2014]LCA Time Limit: 10 Sec Memory Limit: 128 MB Submit: 1266 Solved: 448 [Submit][Stat ...
- 【BZOJ3626】LCA(树链剖分,Link-Cut Tree)
[BZOJ3626]LCA(树链剖分,Link-Cut Tree) 题面 Description 给出一个n个节点的有根树(编号为0到n-1,根节点为0).一个点的深度定义为这个节点到根的距离+1. ...
- 【LNOI2014】【BZOJ3626】NOIp2018模拟(三) LCA
Description 给出一个n个节点的有根树(编号为0到n-1,根节点为0).一个点的深度定义为这个节点到根的距离+1.设$dep[i]$表示点i的深度,$lca(i,j)$表示i与j的最近公共祖 ...
- 【P4211 LNOI2014】LCA——树链剖分 +询问离线
(7.16晚)更完先在B站颓一会儿-- --------------------------------------------------------------- (以下为luogu题面) 题目描 ...
- 【BZOJ-3626】LCA 树链剖分
3626: [LNOI2014]LCA Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1428 Solved: 526[Submit][Status ...
- 【Homework】LCA&RMQ
我校是神校,作业竟然选自POJ,难道不知道“珍爱生命 勿刷POJ”么? 所有注明模板题的我都十分傲娇地没有打,于是只打了6道题(其实模板题以前应该打过一部分但懒得找)(不过感觉我模板还是不够溜要找个时 ...
- 【模板】LCA
代码如下 #include <bits/stdc++.h> using namespace std; const int maxn=5e5+10; inline int read(){ i ...
- 【BZOJ3626】LCA(树上差分,树链剖分)
题意:给出一个n个节点的有根树(编号为0到n-1,根节点为0).一个点的深度定义为这个节点到根的距离+1.设dep[i]表示点i的深度,LCA(i,j)表示i与j的最近公共祖先.有q次询问,每次询问给 ...
- 【模板】Lca倍增法
Codevs 1036 商务旅行 #include<cstdio> #include<cmath> #include<algorithm> using namesp ...
随机推荐
- Oracle GoldenGate DDL 详细说明 使用手册(较早资料)
一. 概述 DDL 相关的参数包括:DDL.DDLERROR.DDLOPTIONS.DDLSUBST.DDLTABLE.GGSCHEMA. PURGEDDLHISTORY.PURGEMARKERHIS ...
- 解决win 和 ubuntu 虚拟机之间 无法 复制粘贴的问题,以及重装vmtool后,还是无法解决的办法
第一步:重新安装vmware-tool 我这里已经装过了,所以显示这个,不然应该是显示:安装 VMware Tools 它会打开一个文件夹,把压缩包复制到任一个地方解压,执行.pl的那个文件:sudo ...
- Python学习---模拟微信网页登录180410
WEB微信 网页登录的猜想: a. 访问页面出现二维码 b. 长轮询监听是否已经扫码并且点击确认 c. 如何进行会话保持 d. 如何获取用户列表 e. 如何发送消息(接收消息) 过程:访问微信官网[h ...
- Lua脚本语法说明(转):
Lua脚本语法说明(增加lua5.1部份特性) 转自:http://www.cnblogs.com/ly4cn/archive/2006/08/04/467550.html Lua 的语法比较简单,学 ...
- 我做SAP CRM One Order redesign的一些心得体会
框架开发和应用程序的开发完全不一样. 举个具体的最近折腾我的例子: 创建新的service order,维护header的shipping data,此时order和shipping data的mod ...
- Jenkins 配合 GitLab 实现分支的自动合并、自动创建 Tag
Jenkins 配合 GitLab 实现分支的自动合并.自动创建 Tag 背景 GitFlow工作流简介 Gitflow工作流定义了一个围绕项目发布的严格分支模型,它会相对复杂一点,但提供了用于一个健 ...
- JVM虚拟机21: 1.8中废弃永久代(PermGen)迎来元空间(Metaspace)
1.JDK8永久代的废弃 JDK8 永久代变化如下图: 1.新生代:Eden+From Survivor+To Survivor 2.老年代:OldGen 3.永久代(方法区的实现) : PermGe ...
- virtualbox+vagrant学习-3-Vagrant Share-3-SSH Sharing
SSH Sharing vagrant share通过向vagrant share提供--SSH标志,使远程SSH访问vagrant环境变得非常容易. 如果你想让同事访问你的SSH,以便对ops问题进 ...
- django 模板关闭自动转义
Django的模板中会对HTML标签和JS等语法标签进行自动转义,原因显而易见,这样是为了安全.但是有的时候我们可能不希望这些HTML元素被转义,比如我们做一个内容管理系统,后台添加的文章中是经过修饰 ...
- unittest中setUp与setUpClass执行顺序
最基础的概念 1.setUP(self)看下面的执行顺序 import unittest class TestGo(unittest.TestCase): def setUp(self): print ...