模板 树链刨分
题目描述
如题，已知一棵包含N个结点的树（连通且无环），每个节点上包含一个数值，需要支持以下操作：
操作1： 格式： 1 x y z 表示将树从x到y结点最短路径上所有节点的值都加上z
操作2： 格式： 2 x y 表示求树从x到y结点最短路径上所有节点的值之和
操作3： 格式： 3 x z 表示将以x为根节点的子树内所有节点值都加上z
操作4： 格式： 4 x 表示求以x为根节点的子树内所有节点值之和
输入输出格式
输入格式：第一行包含4个正整数N、M、R、P，分别表示树的结点个数、操作个数、根节点序号和取模数（即所有的输出结果均对此取模）。
接下来一行包含N个非负整数，分别依次表示各个节点上初始的数值。
接下来N-1行每行包含两个整数x、y，表示点x和点y之间连有一条边（保证无环且连通）
接下来M行每行包含若干个正整数，每行表示一个操作，格式如下：
操作1： 1 x y z
操作2： 2 x y
操作3： 3 x z
操作4： 4 x
输出格式：输出包含若干行，分别依次表示每个操作2或操作4所得的结果（对P取模）
输入输出样例
说明
时空限制：1s，128M
数据规模：
对于30%的数据： N≤10,M≤10 N \leq 10, M \leq 10 N≤10,M≤10
对于70%的数据： N≤103,M≤103 N \leq {10}^3, M \leq {10}^3 N≤103,M≤103
对于100%的数据： N≤105,M≤105 N \leq {10}^5, M \leq {10}^5 N≤105,M≤105
（ 其实，纯随机生成的树LCA+暴力是能过的，可是，你觉得可能是纯随机的么233 ）
(重要的事情说三遍：记得取模)
#include<iostream> #include<cstring> #include<cstdio> using namespace std; #define ls node<<1,L,(L+R)/2 #define rs node<<1|1,(L+R)/2+1,R int w[161000],head[161000],next[311000],to[311000],c[561000],lazy[561000], x,y,z,op,n,m,root,p,son[111000],size[111000],f[111000],top[111000],cnt, deep[111000],dfn[111000],hash[111000]; inline void add(int x,int y) { next[++cnt]=head[x]; to[cnt]=y; head[x]=cnt; } void pushup(int node) { c[node]=c[node<<1]+c[node<<1|1]; c[node]%=p; } void pushdown(int node,int l,int r) { if(lazy[node]) { int len=r-l+1; lazy[node<<1]+=lazy[node]; lazy[node<<1|1]+=lazy[node]; c[node<<1]+=lazy[node]*(len-(len>>1)); c[node<<1|1]+=lazy[node]*(len>>1); lazy[node]=0; } } void getson(int u,int fa) { int v; deep[u]=deep[fa]+1; size[u]=1; for(int i=head[u];i;i=next[i]) { v=to[i]; if(v!=fa) { f[v]=u; getson(v,u); size[u]+=size[v]; if(!son[u]||size[son[u]]<size[v]) { son[u]=v; } } } } void gettop(int u,int fa,int tp) { top[u]=tp; dfn[u]=++cnt; hash[cnt]=u; if(son[u]) gettop(son[u],u,tp); int v; for(int i=head[u];i;i=next[i]) { v=to[i]; if(v!=fa&&v!=son[u]) { gettop(v,u,v); } } } void add(int node,int L,int R,int l,int r,int val) { if(L>r||R<l) return; if(l<=L&&r>=R) { lazy[node]+=val; c[node]+=(R-L+1)*val; return; } pushdown(node,L,R); add(ls,l,r,val); add(rs,l,r,val); pushup(node); } int que(int node,int L,int R,int l,int r) { if(L>r||R<l) return 0; if(l<=L&&r>=R) { return c[node]%p; } pushdown(node,L,R); return (que(ls,l,r)+que(rs,l,r))%p; } void build(int node,int L,int R) { if(L==R) { c[node]=w[hash[L]]; return; } else { build(ls); build(rs); pushup(node); } } inline int get(int x,int y,int z,int o) { int f1=top[x],f2=top[y],ans=0; while(f1!=f2) { if(deep[f1]<deep[f2]) { swap(x,y); swap(f1,f2); } if(o) { add(1,1,cnt,dfn[f1],dfn[x],z); } else { ans+=que(1,1,cnt,dfn[f1],dfn[x]); ans%=p; } x=f[f1]; f1=top[x]; } if(deep[x]>deep[y]) { swap(x,y); } if(o) add(1,1,cnt,dfn[x],dfn[y],z); else ans+=que(1,1,cnt,dfn[x],dfn[y]); return ans%p; } int main() { //freopen("testdata.in","r",stdin); scanf("%d%d%d%d",&n,&m,&root,&p); for(int i=1;i<=n;i++) scanf("%d",&w[i]); for(int i=1;i<n;i++) { scanf("%d%d",&x,&y); add(x,y); add(y,x); } cnt=0; getson(root,0); gettop(root,0,root); build(1,1,cnt); while(m--) { scanf("%d",&op); if(op==1) { scanf("%d%d%d",&x,&y,&z); get(x,y,z,1); } else if(op==2) { scanf("%d%d",&x,&y); printf("%d\n",get(x,y,0,0)); } else if(op==3) { scanf("%d%d",&x,&z); add(1,1,cnt,dfn[x],dfn[x]+size[x]-1,z); } else { scanf("%d",&x); printf("%d\n",que(1,1,cnt,dfn[x],dfn[x]+size[x]-1)%p); } } }
