#include<bits/stdc++.h>
using namespace std;
const int N=2e5+11;
int n,m,r,p,tot,tim;
int son[N],siz[N],dep[N],fa[N],v[N];
int w[N],top[N],dfn[N];
int head[N],nxt[N],ver[N];
inline int read(){
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=x*10+ch-48;ch=getchar();}
return x*f;
}
inline void write(int x) {
if(x<0){putchar('-'); x = -x;}
if(x>9) write(x / 10);
putchar(x % 10 + '0');
}
struct node{
int l,r;
int sum,lazy;
}t[N<<2];
inline void add(int u,int v){//链式前向星
ver[++tot]=v;
nxt[tot]=head[u],head[u]=tot;
return;
}
//线段树
inline void build(int x,int l,int r){//建树
t[x].l=l;t[x].r=r;
if(l==r){
t[x].sum=w[l]%p;
return;
}
int mid=(l+r)>>1;
build(x<<1,l,mid);
build(x<<1|1,mid+1,r);
t[x].sum=(t[x<<1].sum+t[x<<1|1].sum)%p;
return;
}
inline void update(int x){//下传懒标
if(t[x].lazy){
t[x<<1].sum+=t[x].lazy*(t[x<<1].r-t[x<<1].l+1);
t[x<<1|1].sum+=t[x].lazy*(t[x<<1|1].r-t[x<<1|1].l+1);
t[x<<1].sum%=p;t[x<<1|1].sum%=p;
t[x<<1].lazy+=t[x].lazy;
t[x<<1|1].lazy+=t[x].lazy;
t[x].lazy=0;
}
return;
}
inline void modify(int x,int l,int r,int k){//区修
if(l<=t[x].l&&t[x].r<=r){
t[x].sum+=k*(t[x].r-t[x].l+1);
t[x].sum%=p;
t[x].lazy+=k;
return;
}
if(t[x].lazy) update(x);
int mid=(t[x].l+t[x].r)>>1;
if(l<=mid) modify(x<<1,l,r,k);
if(r>mid) modify(x<<1|1,l,r,k);
t[x].sum=(t[x<<1].sum+t[x<<1|1].sum)%p;
return;
}
inline int query(int x,int l,int r){//区查
if(l<=t[x].l&&t[x].r<=r) return t[x].sum;
if(t[x].lazy) update(x);
int res=0;
int mid=(t[x].l+t[x].r)>>1;
if(l<=mid) res+=query(x<<1,l,r);
if(r>mid) res+=query(x<<1|1,l,r);
return res%p;
}
//树链剖分
inline void dfs1(int u,int f){
fa[u]=f;
dep[u]=dep[f]+1;
siz[u]=1;
int maxsize=-1;
for(register int i=head[u];i;i=nxt[i]){
int v=ver[i];
if(v==f) continue;
dfs1(v,u);
siz[u]+=siz[v];
if(siz[v]>maxsize){
maxsize=siz[v];
son[u]=v;
}
}
return;
}
inline void dfs2(int u,int t){
dfn[u]=++tim;
top[u]=t;
w[tim]=v[u];
if(!son[u]) return;
dfs2(son[u],t);
for(register int i=head[u];i;i=nxt[i]){
int v=ver[i];
if(v==fa[u]||v==son[u]) continue;
dfs2(v,v);
}
return;
}
inline void mson(int x,int z){//操作三
modify(1,dfn[x],dfn[x]+siz[x]-1,z);
return;
}
inline int qson(int x){//操作四
return query(1,dfn[x],dfn[x]+siz[x]-1);
}
inline void mchain(int x,int y,int z){//操作一
z%=p;
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]]) swap(x,y);
modify(1,dfn[top[x]],dfn[x],z);
x=fa[top[x]];
}
if(dep[x]>dep[y]) swap(x,y);
modify(1,dfn[x],dfn[y],z);
return;
}
inline int qchain(int x,int y){//操作二
int res=0;
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]]) swap(x,y);
res+=query(1,dfn[top[x]],dfn[x]);
x=fa[top[x]];
}
if(dep[x]>dep[y]) swap(x,y);
res+=query(1,dfn[x],dfn[y]);
return res%p;
}
int main(){
int x,y,z,op;
n=read(),m=read(),r=read(),p=read();
for(register int i=1;i<=n;i++) v[i]=read();
for(register int i=1;i<n;i++){
x=read(),y=read();
add(x,y);
add(y,x);
}
dfs1(r,r);
dfs2(r,r);
build(1,1,n);
while(m--){
op=read();
if(op==1){
x=read(),y=read(),z=read();
mchain(x,y,z);
}
if(op==2){
x=read();y=read();
write(qchain(x,y));
putchar('\n');
}
if(op==3){
x=read(),z=read();
mson(x,z);
}
if(op==4){
x=read();
write(qson(x));
putchar('\n');
}
}
return 0;
}
树链剖分板子
于 2023-02-18 15:14:18 首次发布
这篇文章展示了一段使用C++编写的代码,实现了线段树和树链剖分的数据结构,用于处理区间查询和修改操作。主要功能包括链式前向星、线段树的建立、更新和查询,以及树链剖分的深度优先搜索。这些是图论和算法竞赛中常见的高效算法。
摘要由CSDN通过智能技术生成