#include<bits/stdc++.h>
#define DEBUG
#define mid ((l+r)>>1)
#define lson pos<<1
#define rson pos<<1|1
#define afor(x,y,z) for(int x=y;x<=z;x++)
#define sfor(x,y,z) for(int x=y;x>=z;x--)
using namespace std;
typedef double dou;
typedef long long ll;
typedef const int cint;
typedef unsigned int uint;
typedef const long long cll;
typedef unsigned long long ull;
cint N=1e5+10;
int n,m,r,mod,a[N];
int d[N<<2],b[N<<2];
int tot,head[N],to[N<<1],nxt[N<<1];
int cnt,fa[N],son[N],dep[N],siz[N],p[N],flag[N],bel[N];
void add(int x,int y) {
	to[++tot]=y;
	nxt[tot]=head[x];
	head[x]=tot;
}
void dfs1(int x,int y) {
	siz[x]=1;
	fa[x]=y;
	dep[x]=dep[y]+1;
	for(int i=head[x];i;i=nxt[i]) {
		int z=to[i];
		if(z!=y) {
			dfs1(z,x);
			siz[x]+=siz[z];
			if(siz[z]>siz[son[x]]) son[x]=z;
		}
	}
}
void dfs2(int x,int y) {
	p[x]=++cnt;
	flag[cnt]=x;
	bel[x]=y;
	if(!son[x]) return ;
	dfs2(son[x],y);
	for(int i=head[x];i;i=nxt[i]) {
		int z=to[i];
		if(z!=son[x]&&z!=fa[x]) dfs2(z,z);
	}
}
void build(int l,int r,int pos) {
	if(l==r) {
		d[pos]=a[flag[l]];
		return;
	}
	build(l,mid,lson);
	build(mid+1,r,rson);
	d[pos]=d[lson]+d[rson];
}
void pushdown(int l,int r,int pos) {
	d[lson]+=b[pos]*(mid-l+1);
	d[rson]+=b[pos]*(r-mid);
	b[lson]+=b[pos];
	b[rson]+=b[pos];
	b[pos]=0;
}
void update(int l,int r,int pos,int c,int h,int t) {
	if(h<=l&&r<=t) {
		d[pos]+=c*(r-l+1);
		b[pos]+=c;
		return;
	}
	if(b[pos]) pushdown(l,r,pos);
	if(h<=mid) update(l,mid,lson,c,h,t);
	if(t>mid) update(mid+1,r,rson,c,h,t);
	d[pos]=d[lson]+d[rson];
}
int querysum(int l,int r,int pos,int h,int t) {
	if(h<=l&&r<=t) return d[pos];
	if(b[pos]) pushdown(l,r,pos);
	int sum=0;
	if(h<=mid) sum+=querysum(l,mid,lson,h,t);
	if(t>mid) sum+=querysum(mid+1,r,rson,h,t);
	return sum;
}
void fix(int x,int y,int z) {
	while(bel[x]!=bel[y]) {
		if(dep[bel[x]]<dep[bel[y]]) swap(x,y);
		update(1,n,1,z,p[bel[x]],p[x]);
		x=fa[bel[x]];
	}
	if(dep[x]<dep[y]) swap(x,y);
	update(1,n,1,z,p[y],p[x]);
}
int getsum(int x,int y) {
	int sum=0;
	while(bel[x]!=bel[y]) {
		if(dep[bel[x]]<dep[bel[y]]) swap(x,y);
		sum+=querysum(1,n,1,p[bel[x]],p[x]);
		x=fa[bel[x]];
	}
	if(dep[x]<dep[y]) swap(x,y);
	sum+=querysum(1,n,1,p[y],p[x]);
	return sum;
}
int main() {
	int opt,x,y,z;
	scanf("%d%d%d%d",&n,&m,&r,&mod);
	afor(i,1,n) scanf("%d",&a[i]);
	dfs1(r,r);
	dfs1(r,0);
	build(1,n,1);
	afor(i,1,n-1) {
		scanf("%d%d",&x,&y);
		add(x,y);
		add(y,x);
	}
	afor(i,1,m) {
		scanf("%d",&opt);
		if(opt==1) {
			scanf("%d%d%d",&x,&y,&z);
			fix(x,y,z);
		}
		else if(opt==2) {
			scanf("%d%d",&x,&y);
			printf("%d\n",getsum(x,y)%mod);
		}
		else if(opt==3) {
			scanf("%d%d",&x,&z);
			update(1,n,1,z,p[x],p[x]+siz[x]-1);
		}
		else {
			scanf("%d",&x);
			printf("%d\n",querysum(1,n,1,p[x],p[x]+siz[x]-1)%mod);
		}
	}
	return 0;
}

为什么炸了？