求助！关于建树部分的疑问

在调试的时候发现一个很神奇的问题，

void build(int p,int l,int r){
    t[p].l=l;t[p].r=r;
    t[p].len=r-l+1;t[p].flag=0;
    if(l==r){t[p].v=a[b[l]];return;}
    mid=(l+r)>>1;cout<<l<<' '<<r<<' '<<mid<<endl;
    build(p*2,l,mid);
    build(p*2+1,mid+1,r);
    t[p].v=(t[p*2].v+t[p*2+1].v)%mod;
}

这一段代码不正确，会出现1 3，1 2,2 3这种重复计算的情况，

void build(int p,int l,int r){
    t[p].l=l;t[p].r=r;
    int mid=(l+r)/2;t[p].flag=0;
    t[p].len=r-l+1;
//  cout<<l<<' '<<r<<' '<<mid<<endl;
    if(l==r){t[p].v=a[b[l]];return ;}
    build(p*2,l,mid);
    build(p*2+1,mid+1,r);
    t[p].v=(t[p*2].v+t[p*2+1].v)%mod;
}

而这一段就是正确的。蒟蒻的我百思不得其解，求助各位大佬帮助

下面是ac的代码

#include<bits/stdc++.h>
using namespace std;
#define N 500860
#define ll long long
#define in read() 
ll a[N],b[N],mod;
int fa[N],n,m,dfn[N],cnt,x,y,z,v,l,r,mid,root,op,deep[N],head[N],tp[N],size[N],son[N],tt;
struct node{
    int next,to;
}e[N*2];
struct tree{
    int l,r,len;
    ll v,flag;
}t[N*4];
void down(int p){
    if(t[p].flag==0)return;
    ll v=t[p].flag;
    t[p*2].v+=v*t[p*2].len;    t[p*2].v%=mod;
    t[p*2+1].v+=v*t[p*2+1].len;t[p*2+1].v%=mod;
    t[p*2].flag+=v;  t[p*2+1].flag+=v;
    t[p*2].flag%=mod;t[p*2+1].flag%=mod;
    t[p].flag=0;
} 
void add(int x,int y){
    e[++cnt].next=head[x];
    e[cnt].to=y;
    head[x]=cnt;
}
ll read(){ 
    ll x=0,f=1;char c=getchar();
    while(c<'0'||c>'9'){
        if(c=='-')f=-1;c=getchar();
    } 
    while(c<='9'&&c>='0'){
        x=(x<<3)+(x<<1)+c-'0';
        c=getchar();
    }
    return x*f;
}
void dfs1(int x){
    size[x]=1;
    int big=0;
    for(int i=head[x];i;i=e[i].next){
        int y=e[i].to;
        if(y==fa[x])continue;
        fa[y]=x;deep[y]=deep[x]+1;
        dfs1(y);
        size[x]+=size[y];
        if(size[y]>big)big=size[y],son[x]=y;
    }
}
void dfs2(int x,int tou){
    dfn[x]=++tt;
    b[tt]=x;
    tp[x]=tou;
    if(son[x]) dfs2(son[x],tou);
    for(int i=head[x];i;i=e[i].next){
        int y=e[i].to;
        if(y!=fa[x]&&y!=son[x]){
            dfs2(y,y);
        }
    }
}
//void build(int p,int l,int r){
//  t[p].l=l;t[p].r=r;
//  t[p].len=r-l+1;t[p].flag=0;
//  if(l==r){t[p].v=a[b[l]];return;}
//  mid=(l+r)>>1;cout<<l<<' '<<r<<' '<<mid<<endl;

//  build(p*2,l,mid);
//  build(p*2+1,mid+1,r);
//  t[p].v=(t[p*2].v+t[p*2+1].v)%mod;
//}
void build(int p,int l,int r){
    t[p].l=l;t[p].r=r;
    int mid=(l+r)/2;t[p].flag=0;
    t[p].len=r-l+1;
//  cout<<l<<' '<<r<<' '<<mid<<endl;
    if(l==r){t[p].v=a[b[l]];return ;}
    build(p*2,l,mid);
    build(p*2+1,mid+1,r);
    t[p].v=(t[p*2].v+t[p*2+1].v)%mod;
}
void add(int p,int l,int r,int v){
    if(l>t[p].r||r<t[p].l)return;
    else if(l<=t[p].l&&r>=t[p].r) {
        t[p].v+=t[p].len*v;t[p].v%=mod;
        t[p].flag+=v; t[p].flag%=mod;
    }
    else{
        down(p);
        add(p*2,l,r,v);add(p*2+1,l,r,v);
        t[p].v=(t[p*2].v+t[p*2+1].v)%mod;
    }
}
ll query(int p,int l,int r){
    if(l>t[p].r||r<t[p].l)return 0;
    else if(l<=t[p].l&&r>=t[p].r) return t[p].v;
    else{
        down(p);
        return query(p*2,l,r)+query(p*2+1,l,r);
    }
}
void op1(int x,int y,int z){
    while(tp[x]!=tp[y]){
        if(deep[tp[x]]>deep[tp[y]])swap(x,y);
        add(1,dfn[tp[y]],dfn[y],z);
        y=fa[tp[y]];
    }
    if(deep[x]>deep[y])swap(x,y);
    add(1,dfn[x],dfn[y],z);
}
void op2(int x,int y){
    ll ans=0;
    while(tp[x]!=tp[y]){
        if(deep[tp[x]]>deep[tp[y]])swap(x,y);
        ans+=query(1,dfn[tp[y]],dfn[y]);ans%=mod;
        y=fa[tp[y]];
    }
    if(deep[x]>deep[y])swap(x,y);
    ans+=query(1,dfn[x],dfn[y]);
    cout<<ans%mod<<endl;
}
void op3(int x,int z){
    add(1,dfn[x],dfn[x]+size[x]-1,z);
}
void op4(int x){
    cout<<query(1,dfn[x],dfn[x]+size[x]-1)%mod<<endl;
}
int main(){
    n=in,m=in,root=in,mod=in;
    for(int i=1;i<=n;i++) a[i]=in;
    for(int i=1;i<n;i++){
        x=in,y=in;
        add(x,y);add(y,x);
    }
    dfs1(root);
    dfs2(root,root);
    build(1,1,n);
    while(m--){
        op=in;
        if(op==1){ x=in,y=in,z=in;op1(x,y,z);}
        else if(op==2){ x=in,y=in;op2(x,y);}
        else if(op==3){ x=in,y=in;op3(x,y);}
        else{ x=in;op4(x);}
    }
    return 0; 
}