为什么开了O2反而错了

#include<bits/stdc++.h>
using namespace std;
const int N=2e5+11;
int n,m,r,p,tot,tim;
int son[N],siz[N],dep[N],fa[N],v[N];
int w[N],top[N],dfn[N];
int head[N],nxt[N],ver[N];
struct node{
    int l,r;
    int sum,lazy;
}t[N<<2];
void add(int u,int v){//链式前向星 
    ver[++tot]=v;
    nxt[tot]=head[u],head[u]=tot;
}
//线段树
void build(int x,int l,int r){//建树 
    t[x].l=l;t[x].r=r;
    if(l==r){
        t[x].sum=w[l]%p;
        return;
    }
    int mid=(l+r)>>1;
    build(x<<1,l,mid);
    build(x<<1|1,mid+1,r);
    t[x].sum=(t[x<<1].sum+t[x<<1|1].sum)%p;
}
void update(int x){//下传懒标 
    if(t[x].lazy){ 
        t[x<<1].sum+=t[x].lazy*(t[x<<1].r-t[x<<1].l+1);
        t[x<<1|1].sum+=t[x].lazy*(t[x<<1|1].r-t[x<<1|1].l+1);
        t[x<<1].sum%=p;t[x<<1|1].sum%=p;
        t[x<<1].lazy+=t[x].lazy;
        t[x<<1|1].lazy+=t[x].lazy;
        t[x].lazy=0;
    }
}
void modify(int x,int l,int r,int k){//区修 
    if(l<=t[x].l&&t[x].r<=r){
        t[x].sum+=k*(t[x].r-t[x].l+1);
        t[x].sum%=p;
        t[x].lazy+=k;
        return;
    }
    if(t[x].lazy) update(x);
    int mid=(t[x].l+t[x].r)>>1;
    if(l<=mid) modify(x<<1,l,r,k);
    if(r>mid) modify(x<<1|1,l,r,k);
    t[x].sum=(t[x<<1].sum+t[x<<1|1].sum)%p;
}
int query(int x,int l,int r){//区查 
    if(l<=t[x].l&&t[x].r<=r) return t[x].sum;
    if(t[x].lazy) update(x);
    int res=0;
    int mid=(t[x].l+t[x].r)>>1;
    if(l<=mid) res+=query(x<<1,l,r);
    if(r>mid) res+=query(x<<1|1,l,r);
    return res%p;
}
//树链剖分
void dfs1(int u,int f){
    fa[u]=f;
    dep[u]=dep[f]+1;
    siz[u]=1;
    int maxsize=-1;
    for(int i=head[u];i;i=nxt[i]){
        int v=ver[i];
        if(v==f) continue;
        dfs1(v,u);
        siz[u]+=siz[v];
        if(siz[v]>maxsize){
            maxsize=siz[v];
            son[u]=v;
        }
    }
}
void dfs2(int u,int t){
    dfn[u]=++tim;
    top[u]=t;
    w[tim]=v[u];
    if(!son[u]) return;
    dfs2(son[u],t);
    for(int i=head[u];i;i=nxt[i]){
        int v=ver[i];
        if(v==fa[u]||v==son[u]) continue;
        dfs2(v,v);
    }
}
void mson(int x,int z){//操作三 
    modify(1,dfn[x],dfn[x]+siz[x]-1,z);
}
int qson(int x){//操作四 
    query(1,dfn[x],dfn[x]+siz[x]-1);
}
void mchain(int x,int y,int z){//操作一 
    z%=p;
    while(top[x]!=top[y]){
        if(dep[top[x]]<dep[top[y]]) swap(x,y);
        modify(1,dfn[top[x]],dfn[x],z);
        x=fa[top[x]];
    }
    if(dep[x]>dep[y]) swap(x,y);
    modify(1,dfn[x],dfn[y],z);
}
int qchain(int x,int y){//操作二 
    int res=0;
    while(top[x]!=top[y]){
        if(dep[top[x]]<dep[top[y]]) swap(x,y);
        res+=query(1,dfn[top[x]],dfn[x]);
        x=fa[top[x]];
    }
    if(dep[x]>dep[y]) swap(x,y);
    res+=query(1,dfn[x],dfn[y]);
    return res%p;
}
int main(){
    int x,y,z,op;
    scanf("%d%d%d%d",&n,&m,&r,&p);
    for(int i=1;i<=n;i++) scanf("%d",&v[i]);
    for(int i=1;i<n;i++){
        scanf("%d%d",&x,&y);
        add(x,y);
        add(y,x);
    }
    dfs1(r,r);
    dfs2(r,r);
    build(1,1,n);
    while(m--){
        scanf("%d",&op);
        if(op==1){
            scanf("%d%d%d",&x,&y,&z);
            mchain(x,y,z);
        }
        if(op==2){
            scanf("%d%d",&x,&y);
            printf("%d\n",qchain(x,y));
        }
        if(op==3){
            scanf("%d%d",&x,&z);
            mson(x,z);
        }
        if(op==4){
            scanf("%d",&x);
            printf("%d\n",qson(x));
        }
    }
    return 0;
}

不开O2是对的，开了O2，就RE和TLE

对的

开了O2的

求大佬解答