救命！！！

哪里错了呀，万能的谷民们救救我！！！

#include<bits/stdc++.h>
#define N 100005
using namespace std;
long long a[8*N],lazy[8*N],n,m,x,y,k,h[8*N],w[8*N],fa[8*N],dep[8*N],siz[8*N],son[8*N],top[8*N],id[8*N],p;
int cnt;
struct Node{
    int ne,to,val;
}e[8*N];
void addl(int u,int v)
{
    e[++cnt].ne=h[u];
    e[cnt].to=v;
    h[u]=cnt;
}
void push_up(int num)
{
    a[num]=a[num*2]+a[num*2+1];
}
void push_down(int num,int l,int r)
{
    if(lazy[num])
    {
    lazy[num*2]+=lazy[num];
    lazy[num*2+1]+=lazy[num];
    int mid=(l+r)>>1;
    a[num*2]+=(mid-l+1)*lazy[num];
    a[num*2+1]+=(r-mid)*lazy[num];
    lazy[num]=0;
    }
}
void build(int l,int r,int num)
{
    if(l==r) 
    {
        w[l]%=p;
        a[num]=w[l];
    }
    else
    {
        int mid=(l+r)>>1;
        build(l,mid,2*num);
        build(mid+1,r,2*num+1);
        push_up(num);
    }
}
void add(int x,int y,int k,int l,int r,int num)
{
    if(x<=l&&y>=r)
    {
        lazy[num]+=k;
        a[num]+=(r-l+1)*k;
    }
    else
    {
        push_down(num,l,r);
        int mid=(l+r)>>1;
        if(x<=mid) add(x,y,k,l,mid,num*2);
        if(y>mid) add(x,y,k,mid+1,r,num*2+1);
        push_up(num);
    }
}
long long find(int x,int y,int l,int r,int num)
{
    if(x<=l && y>=r) return a[num];
    push_down(num,l,r);
    {
        int mid=(l+r)>>1;
        long long ans=0;
        if(x<=mid) ans+=find(x,y,l,mid,num*2);
        if(y>mid) ans+=find(x,y,mid+1,r,num*2+1);
        return ans;
    }
}
void dfs1(int u,int f,int deep)
{
    dep[u]=deep;
    fa[u]=f;
    siz[u]=1;
    int maxson=-1;
    for(int i=h[u];i;i=e[i].ne)
    {
        int v=e[i].to;
        if(v==f) continue;
        dfs1(v,u,deep+1);
        siz[u]+=siz[v];
        if(siz[v]>maxson)
        {
            son[u]=v;
            maxson=siz[v];
        }
    }
}
void dfs2(int u,int tto)
{
    id[u]=++cnt;
    w[cnt]=w[u];
    top[u]=tto;
    if(!son[u])
    return ;
    dfs2(son[u],tto);
    for(int i=h[u];i;i=e[i].ne)
    {
        int v=e[i].to;
        if(v==fa[u]||v==son[u])
        continue;
        dfs2(v,v); 
    }
}
int main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);
    int r;
    cin>>n>>m>>r>>p;
    for(int i=1;i<=n;i++)
    cin>>w[i];
    for(int i=1;i<n;i++)
    {
        cin>>x>>y;
        addl(x,y);
        addl(y,x);
    }
    cnt=0;
    dfs1(r,0,1);
    dfs2(r,r);
    build(1,n,1);
    for(int i=1;i<=m;i++)
    {
        int opt;
        cin>>opt;
        if(opt==1)
        {
            cin>>x>>y>>k;
            k%=p;
            while(top[x]!=top[y])
            {
                if(dep[top[x]]<dep[top[y]]) swap(x,y);
                add(id[top[x]],id[x],k,1,n,1);
                x=fa[top[x]];
            }
            if(dep[x]>dep[y]) swap(x,y);
            add(id[x],id[y],k,1,n,1);
        }
        if(opt==2)
        {
            long long ans=0;
            cin>>x>>y;
            while(top[x]!=top[y])
            {
                if(dep[top[x]]<dep[top[y]]) swap(x,y);
                ans+=find(id[top[x]],id[x],1,n,1);
                ans%=p;
                x=fa[top[x]];
            }
            if(dep[x]>dep[y]) swap(x,y);
            ans+=find(id[x],id[y],1,n,1);
            cout<<ans%p<<"\n";
        }
        if(opt==3)
        {
            cin>>x>>y;
            add(id[x],id[x]+siz[x]-1,y,1,n,1);
        }
        if(opt==4)
        {
            cin>>x;
            cout<<find(id[x],id[x]+siz[x]-1,1,n,1)<<"\n";
        }
    }
    return 0;
 }