28pts AC 5,8,11其余WA 玄关

#include <iostream>
using namespace std;
typedef long long ll;
const ll N=100005;
ll n,m,r,mod;
ll a[N];
ll siz[N],fa[N],dep[N],seg[N],rev[N],top[N],son[N];
ll hd[N],nt[2*N],v[2*N],cnt=1;
ll tr[N*4],mk[N*4];
inline void add(ll x,ll y){
    v[++cnt]=y;
    nt[cnt]=hd[x];
    hd[x]=cnt;
}
inline void dfs1(ll u,ll f){
    siz[u]=1;
    fa[u]=f;
    dep[u]=dep[f]+1;
    for(ll i=hd[u];i;i=nt[i]){
        ll y=v[i];
        if(y!=f){
            dfs1(y,u);
            siz[u]+=siz[y];
            if(siz[y]>siz[son[u]]){
                son[u]=y;
            }
        }
    }
}
inline void dfs2(ll u,ll f){
    if(son[u]){
        seg[son[u]]=++seg[0];
        top[son[u]]=top[u];
        rev[seg[0]]=son[u];
        dfs2(son[u],u);
    }
    for(ll i=hd[u];i;i=nt[i]){
        ll y=v[i];
        if(!top[y]){
            seg[y]=++seg[0];
            rev[seg[0]]=y;
            top[y]=y;
            dfs2(y,u);
        }
    }
}
inline void build(ll l,ll r,ll p){
    if(l==r){
        tr[p]=a[rev[l]];
    }else{
        ll mid=(l+r)/2;
        build(l,mid,p*2);
        build(mid+1,r,p*2+1);
        tr[p]=tr[p*2]+tr[p*2+1];
    }
}
inline void push_down(ll p,ll len){
    mk[p*2]+=mk[p];
    mk[p*2+1]+=mk[p];
    tr[p*2]+=mk[p]*(len-len/2);
    tr[p*2+1]+=mk[p]*(len/2);
    mk[p]=0;
}
inline void update(ll l,ll r,ll d,ll p,ll cl,ll cr){
    if(cl>r||cr<l){
        return;
    }else if(cl>=l&&cr<=r){
        tr[p]+=(cr-cl+1)*d;
        if(cr>cl){
            mk[p]+=d;
        }
    }else{
        ll mid=(cl+cr)/2;
        push_down(p,cr-cl+1);
        update(l,r,d,p*2,cl,mid);
        update(l,r,d,p*2+1,mid+1,cr);
        tr[p]=tr[p*2]+tr[p*2+1];
    }
}
inline ll query(ll l,ll r,ll p,ll cl,ll cr){
    if(cl>r||cr<l){
        return 0;
    }else if(cl>=l&&cr<=r){
        return tr[p];
    }else{
        ll mid=(cl+cr)/2;
        push_down(p,cr-cl+1);
        return (ll)(query(l,r,p*2,cl,mid)+query(l,r,p*2+1,mid+1,cr));
    }
}
inline void op1(ll x,ll y,ll z){
    while(top[x]!=top[y]){
        if(dep[top[x]]<dep[top[y]]){
            swap(x,y);
        }
        update(seg[top[x]],seg[x],z,1,1,n);
        x=fa[top[x]]; 
    }
    if(dep[x]>dep[y]){
        swap(x,y);
    } 
    update(seg[x],seg[y],z,1,1,n);
}
inline ll op2(ll x,ll y){
    ll ans=0;
    while(top[x]!=top[y]){
        if(dep[top[x]]<dep[top[y]]){
            swap(x,y);
        }
        ans=(ans+query(seg[top[x]],seg[x],1,1,n));
        x=fa[top[x]]; 
    }
    if(dep[x]>dep[y]){
        swap(x,y);
    } 
    ans=(ans+query(seg[x],seg[y],1,1,n));
    return ans;
}
inline void op3(ll x,ll y){
    update(seg[x],seg[x]+siz[x]-1,y,1,1,n);
}
inline ll op4(ll x){
    return query(seg[x],seg[x]+siz[x]-1,1,1,n);
}
int main(){
    std::ios::sync_with_stdio(0);
    std::cin.tie(0);
    std::cout.tie(0);
    cin>>n>>m>>r>>mod;
    for(ll i=1;i<=n;i++){
        cin>>a[i];
    }
    for(ll i=1;i<n;i++){
        ll x,y;
        cin>>x>>y;
        add(x,y);
        add(y,x);
    }
    dfs1(r,0);
    seg[r]=++seg[0];
    top[r]=r;
    rev[r]=1;
    dfs2(r,0);
    build(1,n,1);
    for(ll i=1;i<=m;i++){
        ll k;
        cin>>k;
        if(k==1){
            ll x,y,z;
            cin>>x>>y>>z;
            op1(x,y,z);
        }else if(k==2){
            ll x,y;
            cin>>x>>y;
            cout<<op2(x,y)%mod<<'\n';
        }else if(k==3){
            ll x,y;
            cin>>x>>y;
            op3(x,y);
        }else{
            ll x;
            cin>>x;
            cout<<op4(x)%mod<<'\n';
        }
    } 
}