本地测数据能过，但一交全RE？

#include<bits/stdc++.h>
using namespace std;
const int N=1e7+1;
int n,m,r,mod,num;
int w[N];
struct point
{
    int id,w,head,son,f,total,depth,topp;
}p[N];
struct line
{
    int f,t,next;
}l[N];
int tot;
void fun(int f,int t)
{
    l[++tot].f=f;
    l[tot].t=t;
    l[tot].next=p[f].head;
    p[f].head=tot;
}
void dfs1(int x,int depth)
{
    p[x].total=1;
    p[x].depth=depth;
    int maxn=-1;
    for(int i=p[x].head;i;i=l[i].next)
    {
        if(l[i].t==p[x].f) continue;
        p[l[i].t].f=x;
        dfs1(l[i].t,depth+1);
        p[x].total+=p[l[i].t].total;
        if(p[l[i].t].total>maxn)
        {
            maxn=p[l[i].t].total;
            p[x].son=l[i].t;
        }
    }
}
void dfs2(int x,int topp)
{
    p[x].id=++num;
    w[num]=p[x].w;
    p[x].topp=topp;
    if(!p[x].son) return ;
    dfs2(p[x].son,topp);
    for(int i=p[x].head;i;i=l[i].next)
    {
        if(l[i].t==p[x].f||l[i].t==p[x].son) continue;
        dfs2(l[i].t,l[i].t);
    }
}
struct tree
{
    int w,tag;
}t[N];
int ls(int x)
{
    return x<<1;
}
int rs(int x)
{
    return x<<1|1;
}
void push_up(int rt)
{
    t[rt].w=(t[rt*2].w+t[rt*2+1].w)%mod;
}
void push_down(int rt,int l,int r)
{
    t[ls(rt)].tag+=t[rt].tag;
    t[rs(rt)].tag+=t[rt].tag;
    t[ls(rt)].w+=((r+l)/2-l+1)*t[rt].tag;
    t[ls(rt)].w%=mod;
    t[rs(rt)].w+=(r-(r+l)/2)*t[rt].tag;
    t[rs(rt)].w%=mod;
    t[rt].tag=0;
}
void build(int rt,int l,int r)
{
    if(l==r)
    {
        t[rt].w=w[l]%mod;
        return ;
    }
    int mid=(l+r)/2;
    build(ls(rt),l,mid);
    build(rs(rt),mid+1,r);
    push_up(rt);
}
void update(int rt,int l,int r,int ql,int qr,int k)
{
    if(ql<=l&&qr>=r)
    {
        t[rt].tag+=k;
        t[rt].w+=k*(r-l+1);
        t[rt].w%=mod;
        return ;
    }
    if(t[rt].tag) push_down(rt,l,r);
    int mid=(r+l)/2;
    if(ql<=mid) update(ls(rt),l,mid,ql,qr,k);
    if(qr>mid) update(rs(rt),mid+1,r,ql,qr,k);
    push_up(rt);
    return ;
}
int query(int rt,int l,int r,int ql,int qr)
{
    if(ql<=l&&qr>=r) return t[rt].w;
    if(t[rt].tag) push_down(rt,l,r);
    int mid=(r+l)/2,ans=0;
    if(ql<=mid) ans+=query(ls(rt),l,mid,ql,qr),ans%=mod;
    if(qr>mid) ans+=query(rs(rt),mid+1,r,ql,qr),ans%=mod;
    push_up(rt);
    return ans;
}
int qrange(int x,int y)
{
    int ans=0;
    while(p[x].topp!=p[y].topp)
    {
        if(p[p[x].topp].depth<p[p[y].topp].depth) swap(x,y);
        ans+=query(1,1,n,p[p[x].topp].id,p[x].id);
        ans%=mod;
        x=p[p[x].topp].f;
    }
    if(p[x].depth>p[y].depth) swap(x,y);
    ans+=query(1,1,n,p[x].id,p[y].id);
    return ans%mod;
}
void uprange(int x,int y,int z)
{
    while(p[x].topp!=p[y].topp)
    {
        if(p[p[x].topp].depth<p[p[y].topp].depth) swap(x,y);
        update(1,1,n,p[p[x].topp].id,p[x].id,z);
        x=p[p[x].topp].f;
    }
    if(p[x].depth>p[y].depth) swap(x,y);
    update(1,1,n,p[x].id,p[y].id,z);
}
int qson(int x)
{
    return query(1,1,n,p[x].id,p[x].id+p[x].total-1);
}
int upson(int x,int y)
{
    update(1,1,n,p[x].id,p[x].id+p[x].total-1,y);
}
int main()
{
    //freopen(".in","r",stdin);
    cin>>n>>m>>r>>mod;
    for(int i=1;i<=n;i++)
    cin>>p[i].w;
    for(int i=1;i<n;i++)
    {
        int x,y;
        cin>>x>>y;
        fun(x,y);
        fun(y,x);
    }
    dfs1(r,1);
    dfs2(r,r);
    build(1,1,n);
    for(int i=1;i<=m;i++)
    {
        int op,x,y,z;
        cin>>op;
        switch(op)
        {
            case 1:
            {
                cin>>x>>y>>z;
                uprange(x,y,z);
                break;
            }
            case 2:
            {
                cin>>x>>y;
                cout<<qrange(x,y)<<endl;
                break;
            }
            case 3:
            {
                cin>>x>>z;
                upson(x,z);
                break;
            }
            case 4:
            {
                cin>>x;
                cout<<qson(x)<<endl;
                break;
            }
        }
    }
    return 0;
 }