救救孩子，初学过不了样例

#include <bits/stdc++.h>
#define lc k*2
#define rc k*2+1
using namespace std;
int n, m, r, p, cnt=0;
const int MN=1e5+115;
struct EDGE{
    int to, nex, beg, w;
}edge[MN];
void add_node(int x, int y){
    cnt++;
    edge[cnt].to=y;
    edge[cnt].nex=edge[x].beg;
    edge[x].beg=cnt;
    return;
}
struct STN{
    int l, r, val, tag; 
}tr[MN*4];
int wt[MN], w[MN];
int res=0;
int son[MN], id[MN], father[MN], tot, dep[MN], size[MN], top[MN];
void pushup(int k){
    tr[k].val=(tr[lc].val+tr[rc].val)%p;
    return;
}
void build(int k, int l, int r){
    tr[k].l=l; tr[k].r=r;
    if(l==r){
        tr[k].val=wt[l];
        if(tr[k].val>p){tr[k].val%=p;}
        return;
    }
    int mid=(l+r)>>1;
    build(lc,l,mid);build(rc,mid+1,r);
    pushup(k);
    return;
} 
void pushdown(int k){
    if(tr[k].tag){
        tr[lc].tag+=tr[k].tag;
        tr[rc].tag+=tr[k].tag;
        tr[lc].val+=tr[k].tag*(tr[lc].r-tr[lc].l+1);
        tr[rc].val+=tr[k].tag*(tr[rc].r-tr[rc].l+1);
        tr[lc].val%=p;
        tr[rc].val%=p;
        tr[k].tag=0;
    }

}
void update(int k, int l, int r, int v){
    if(tr[k].l>=l&&tr[k].r<=r){
        tr[k].val+=v*(tr[k].r-tr[k].l+1);
        tr[k].tag+=v;
    }else{
        pushdown(k);
        int mid=(tr[k].l+tr[k].r)>>1;
        if(mid>=l){
            update(lc,l,r,v);
        }
        if(mid<r){
            update(rc,l,r,v);
        }
        pushup(k);
    }

}
void query(int k, int l, int r){
    if(tr[k].l>=l&&tr[k].r<=r){
        res+=tr[k].val;
        res%=p;
        return;
    }
    pushdown(k);
    int mid=(tr[k].l+tr[k].r)>>1;
    if(mid>=l){
        query(lc,l,r);
    }
    if(mid<r){
        query(rc,l,r);
    }
}
void dfs1(int x, int f, int depth){
    dep[x]=depth;father[x]=f;size[x]=1;
    int bs=-1;
    for(int i=edge[x].beg;i;i=edge[i].nex){
        int y=edge[i].to;
        if(y==f){
            continue;
        }
        dfs1(y,x,depth+1);
        size[x]+=size[y];
        if(size[y]>bs){
            son[x]=y;
            bs=size[y];
        }
    }
    return;
}
void dfs2(int x, int topf){//我讨厌树链剖分呜呜呜 
    id[x]=++tot;wt[tot]=w[x];top[x]=topf;
    if(!son[x]){
        return;
    }
    dfs2(son[x],topf);
    for(int i=edge[x].beg;i;i=edge[i].nex){
        int y=edge[i].to;
        if(y==father[x]||y==son[x]){
            continue;
        }
        dfs2(y,y);
    }
    return;
}
//难绷的来了 
int Range(int x, int y){
    int ans=0;
    while(top[x]!=top[y]){//两点不同链 
        if(dep[top[x]]<dep[top[y]]){
            swap(x,y);
        }
        res=0;
        query(1,id[top[x]],id[x]);ans+=res;ans%=p;
        x=father[top[x]];//x跳到顶上 

    }
    //马上处就停了 
    if(dep[x]>dep[y]){
        swap(x,y);
    }
    res=0;
    query(1,id[x],id[y]);
    ans+=res;
    return ans%p;
}
void Change(int x, int y, int k){
    k%=p;
    while(top[x]!=top[y]){
        if(dep[top[x]]<dep[top[y]]){
            swap(x,y);
        }
        update(1,id[top[x]],id[x],k);
        x=father[top[x]];
    }
    if(dep[x]>dep[y]){
        swap(x,y);
    }
    update(1,id[x],id[y],k);
    return;
}
int Son(int x){
    res=0;
    query(1,id[x],id[x]+size[x]-1);
    return res;
}
void cSon(int x, int k){
    update(1,id[x],id[x]+size[x]-1,k);
    return;
}
int main(){
    scanf("%d%d%d%d",&n,&m,&r,&p);
    for(int i=1; i<=n; i++){scanf("%d",&w[i]);}
    for(int i=1; i<=n; i++){
        int a, b;
        scanf("%d%d",&a,&b);
        add_node(a,b);
        add_node(b,a);
    }
    dfs1(r,0,1);
    dfs2(r,r);
    build(1,1,n);
    for(int i=1; i<=m; i++){
        int op;
        scanf("%d",&op);
        if(op==1){
            int a, b, c;
            scanf("%d%d%d",&a,&b,&c);
            Change(a,b,c);
        }else if(op==2){
            int a, b;
            scanf("%d%d",&a,&b);
            printf("%d\n",Range(a,b));
        }else if(op==3){
            int a, b;
            scanf("%d%d",&a,&b);
            cSon(a,b);
        }else{
            int a;
            scanf("%d",&a);
            printf("%d\n",Son(a));
        }
    }
    return 0;
}