树pao(雾)

今天难得睡醒了，大概睡了13个小时，打算今晚把树剖学了。

嘿嘿嘿雀魂真好玩，这篇文章咕了

有学上了，找到了水平远高于我的队友，好开心的说，，，

卡空间感觉治不了了。。。

板子题是洛谷P3384

实在不知道怎么优化了，感觉应该是哪里写错了，但是是MLE又不是WA，明天再说惹

参考博客：

https://www.cnblogs.com/pks-t/p/9194322.html#_label0

https://www.cnblogs.com/George1994/p/7821357.html

下面那篇要好懂一些，码风也不错。所以为什么还要粘上面的呢？

其实树剖真不难。。就是码量大一些，如果你会 线段树，dfs序，lca。那么应该一看就懂了。

树剖其实就是把一个树剖成很多条链，然后在链上搞来搞去。

先来了解一下概念

• 重结点：子树结点数目最多的结点；
• 轻节点：父亲节点中除了重结点以外的结点；
• 重边：父亲结点和重结点连成的边；
• 轻边：父亲节点和轻节点连成的边；
• 重链：由多条重边连接而成的路径；
• 轻链：由多条轻边连接而成的路径；

然后有几个结论，不会证真的不会证，·········  记不得了，可以看上面blog

反正就是告诉你我们这样搞不会t。。。

那么我们怎么维护呢

比方说我们要  求 x 到 y 的路径长。

当x与y不在一条重链上时，我们可以 往上跳，跳到 top[x]这个位置，然后累加这一段答案，因为x与top[x]在一条重链所以很简单。

这样跳呀跳呀就跳到一条重链上去了。 然后再算一下就阔以了。

为什么不会T呢？复杂度是什么呢？这些我都不会。

同理 更新操作也是一样的。

如果对子树呢? 显然子树内的序号都是连续的。所以直接对 tid[x],tid[x]+siz[x]-1   这段区间搞来搞去就行了。

484很简单qwq

嗯敲代码的时候可就爽死你了

哦还没有说怎么算top，hson这些东西，可以先dfs一下求siz和hson。再dfs一下求top。具体看代码。

然后粘一个适合我自己码风的板子。vector的改天再说。因为自己平时喜欢vector。。。

要板子的话可以去看上面blog，真的讲得很好，每一步的代码什么的。

 #include <bits/stdc++.h>
 typedef long long ll;
 using namespace std;
 inline int read() {
     int X=,w=; char c=getchar();
     while (c<''||c>'') { if (c=='-') w=-; c=getchar(); }
     while (c>=''&&c<='') X=(X<<)+(X<<)+c-'',c=getchar();
     return X*w;
 }
 const int N = 1e5+;
 struct Edge{
     int to,next;
 }edge[N];int cnt=;int head[N];
 void addEdge(int u,int v){
     edge[++cnt].to=v;
     edge[cnt].next=head[u];
     head[u]=cnt;
 }
 struct Node{
     int lz,val;
 }tree[N];
 int siz[N],top[N],hson[N],dep[N],fa[N],tid[N],ran[N];
 int n,m,r,p,val[N],tot=;
 void dfs1(int u,int par,int d){
     dep[u]=d;
     fa[u]=par;
     siz[u]=;
     for(int i=head[u];i;i=edge[i].next){
         int v = edge[i].to;
         if(v!=fa[u]){
             dfs1(v,u,d+);
             siz[u]+=siz[v];
             if(hson[u]==-||siz[v]>siz[hson[u]]){
                 hson[u]=v;
             }
         }
     }
 }
 void dfs2(int u,int t){
     top[u]=t;
     tid[u]=++tot;
     ran[tot]=u;
     if(!hson[u])return;
     dfs2(hson[u],t);
     for(int i=head[u];i;i=edge[i].next){
         int v = edge[i].to;
         if(v!=hson[u]&&v!=fa[u]){
             dfs2(v,v);
         }
     }
 }
 void push_up(int x){
     tree[x].val=(0ll+tree[x<<].val+tree[x<<|].val)%p;
 }
 void push_down(int x,int l,int r){
     if(!tree[x].lz)return;
     int mid=l+r>>;
     tree[x<<].val=(0ll+tree[x<<].val+(mid-l+)*tree[x].lz)%p;
     tree[x<<|].val=(0ll+tree[x<<|].val+(r-mid)*tree[x].lz)%p;
     tree[x<<].lz=(0ll+tree[x<<].lz+tree[x].lz)%p;
     tree[x<<|].lz=(0ll+tree[x<<|].lz+tree[x].lz)%p;
     tree[x].lz=;
 }
 void build(int k,int l,int r){
     if(l==r){
         tree[k].val=val[ran[l]]%p;
     } else{
         int mid=l+r>>;
         build(k<<,l,mid);
         build(k<<|,mid+,r);
         push_up(k);
     }
 }
 int query(int l,int r,int ql,int qr,int rt){
     int res = ;
     if(ql<=l&&qr>=r){
         res=(0ll+res+tree[rt].val)%p;
         return res;
     }
     int mid = (l+r)>>;
     push_down(rt,l,r);
     if(ql<=mid) res=(0ll+res+query(l,mid,ql,qr,rt<<))%p;
     if(qr>mid) res=(0ll+res+query(mid+,r,ql,qr,rt<<|))%p;
     return res;
 }
 void upd(int l,int r,int ul,int ur,int rt,int k){
     if(ul<=l&&ur>=r) {
         tree[rt].val = (tree[rt].val + k * (r - l + ))%p;
         tree[rt].lz+=k;
         return;
     }
     int mid=l+r>>;
     push_down(rt,l,r);
     if(ul<=mid)upd(l,mid,ul,ur,rt<<,k);
     if(ur>mid)upd(mid+,r,ul,ur,rt<<|,k);
     push_up(rt);
 }
 int query_p(int u,int v){//x
     int ans=;
     while (top[u]!=top[v]){
         if(dep[top[u]]<dep[top[v]])swap(u,v);
         ans=(0ll+ans%p+query(,n,tid[top[u]],tid[u],))%p;
         u=fa[top[u]];
     }
     if(dep[u]>dep[v])swap(u,v);
     ans=(0ll+ans+query(,n,tid[u],tid[v],))%p;
     return ans;
 }
 void upd_p(int u,int v,int z){
     z%=p;
     while (top[u]!=top[v]){
         if(dep[top[u]]<dep[top[v]])swap(u,v);
         upd(,n,tid[top[u]],tid[u],,z);
         u=fa[top[u]];
     }
     if(dep[u]>dep[v])swap(u,v);
     upd(,n,tid[u],tid[v],,z);
 }
 int main(){
     n=read();m=read();r=read();p=read();
     for(int i=;i<=n;i++)val[i]=read();
     int op,x,y,z;
     for(int i=;i<n;i++){
         x=read();y=read();
         addEdge(x,y);
         addEdge(y,x);
     }
     dfs1(r,,);
     dfs2(r,r);
     build(,,n);
     /**while (m--){
         op=read();
         if(op==1){
             x=read();y=read();z=read();
             upd_p(x,y,z);
         } else if(op==2){
             x=read();y=read();
             printf("%d\n",query_p(x,y));
         } else if(op==3){
             x=read();z=read();
             upd(1,n,tid[x],tid[x]+siz[x]-1,1,z);
         } else{
             x=read();
             printf("%d\n",query(1,n,tid[x],tid[x]+siz[x]-1,1));
         }
     }*/
 }