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

#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;
 }