[LNOI2014]LCA
这道题难道不是数据结构水题吗 逃）
首先想到一个i点的\(dep\)相当于根节点走到i 进过的点之和，所以这个点对答案的贡献就是从这个点走到根节点。因此我们可以转换为当前节点\(i\)到根节点所经过的所有节点权值++。
显然可以用线段树+树链剖分的数据结构进行优化 对于每一个询问 建一颗权值线段树 存储当前区间\(sum\),\(sum\)存储的就是\(sum\)
然后我们需要想到用 类似于前缀和的思想来进行优化
例如 \(l=2 ,r=4 ,z=4\) 我们可以用前缀和的思想进行优化 计算\(1到4\)的\(sum\)\(-1到1的sum\) 就可以得到答案了
想到了前缀和 想到了树剖+线段树 如果每一次查询 你都要新建两次树
那就会GG!
因此我们又想到了莫队的思想，将询问区间的\(l,r\)排序（因为询问区间为\(1到l\)，\(1到r\)，左区间都是一样的，所以我们可以将l，r一起排序），然后一个一个的加点 这样就可以不用每一次询问都新建一颗树了
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#define int long long
using namespace std;
const int mod=201314,maxn=100010;
int n,q,tot,fa[maxn],cnt_q;
int first[maxn],nxt[maxn<<1],to[maxn<<1],dep[maxn<<1],cnt_id;
int id[maxn<<1],rk[maxn<<1],size[maxn<<1],son[maxn<<1],top[maxn<<1],now;
int tag[maxn<<1],sum[maxn<<1];
struct node{
int id,pos,z,flag;
inline bool operator <(const node &x)const{
return pos<x.pos || (x.pos==pos && id<x.id);
}
}ques[maxn<<4];
struct node1{
int ans1,ans2;
}ans[maxn<<2];
void dfs1(int x,int father){
size[x]=1;fa[x]=father;
for(int i=first[x];i;i=nxt[i]){
int y=to[i];if(y==father) continue;
dep[y]=dep[x]+1;
dfs1(y,x);size[x]+=size[y];
if(size[son[x]]<size[y])
son[x]=y;
}
}
void dfs2(int x,int topo){
id[x]=++cnt_id;rk[cnt_id]=x;top[x]=topo;
if(son[x]) dfs2(son[x],topo);
for(int i=first[x];i;i=nxt[i]){
int y=to[i];if(y==fa[y] || y==son[x]) continue;
dfs2(y,y);
}
}
void add(int x,int y){nxt[++tot]=first[x];first[x]=tot;to[tot]=y;}
void tagg(int node,int l,int r,int ta)
{
sum[node]=(sum[node]+(r-l+1)*ta)%mod;
if(l<r) tag[node]=(tag[node]+ta)%mod;
}
void push_down(int node,int l,int r){
if(l<r && tag[node]){
int mid=(l+r)>>1;
tagg(node<<1,l,mid,tag[node]);
tagg(node<<1|1,mid+1,r,tag[node]);
}
tag[node]=0;
}
void push_up(int node){sum[node]=(sum[node<<1]+sum[node<<1|1])%mod;}
void modify(int node,int l,int r,int ql,int qr){
if(qr<l || ql>r) return ;
if(ql<=l && r<=qr){
tagg(node,l,r,1);
return ;
}
push_down(node,l,r);
int mid=(l+r)>>1;
modify(node<<1,l,mid,ql,qr);modify(node<<1|1,mid+1,r,ql,qr);
push_up(node);
}
void line_modify(int x,int y){
int tx=top[x],ty=top[y];
while(tx!=ty){
if(dep[tx]<dep[ty]){
x^=y^=x^=y;tx^=ty^=tx^=ty;
}
modify(1,1,n,id[tx],id[x]);
x=fa[tx];tx=top[x];
}
if(dep[x]<dep[y])
x^=y^=x^=y;
modify(1,1,n,id[y],id[x]);
}
int query(int node,int l,int r,int ql,int qr){
if(ql<=l && r<=qr) return (sum[node]%mod);
push_down(node,l,r);
int mid=(l+r)>>1,anss=0;
if(ql<=mid) anss+=query(node<<1,l,mid,ql,qr);
if(mid<qr) anss+=query(node<<1|1,mid+1,r,ql,qr);
return anss%mod;
}
int line_query(int x,int y){
int tx=top[x],ty=top[y],anss=0;
while(tx!=ty){
if(dep[tx]<dep[ty]){
x^=y^=x^=y;tx^=ty^=tx^=ty;
}
anss+=query(1,1,n,id[tx],id[x]);
anss%=mod;
x=fa[tx];tx=top[x];
}
if(dep[x]<dep[y])
x^=y^=x^=y;
anss+=query(1,1,n,id[y],id[x]);
return anss%mod;
}
signed main(){
scanf("%lld %lld",&n,&q);
for(int i=2,x;i<=n;i++){
scanf("%lld",&x);add(++x,i);
}
dep[1]=1;
dfs1(1,0);dfs2(1,1);
for(int i=1,l,r,z;i<=q;i++){
scanf("%lld %lld %lld",&l,&r,&z);r++;z++;
ques[++cnt_q]=((node){i,l,z,0});ques[++cnt_q]=((node){i,r,z,1});
}
sort(ques+1,ques+1+cnt_q);now=0;
for(int i=1;i<=cnt_q;i++){
while(now<ques[i].pos)
line_modify(1,++now);
int num=ques[i].id;
if(ques[i].flag==1) ans[num].ans1=line_query(1,ques[i].z);
else ans[num].ans2 =line_query(1,ques[i].z);
}
for(int i=1;i<=q;i++)
printf("%lld\n",(ans[i].ans1-ans[i].ans2+mod)%mod);
return 0;
}