#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=1e5+10;
inline int read()
{
int s=0,w=1;
char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();
return s*w;
}
int n,m,rt,mod;
int r[maxn],dep[maxn],fa[maxn],top[maxn],son[maxn],siz[maxn],id[maxn],w[maxn];
int head[maxn];
struct Eage
{
int to;
int nt;
}e[maxn*2];
int cnt;
inline void add_eage(int x,int y)
{
e[++cnt].nt=head[x];
e[cnt].to=y;
head[x]=cnt;
}
inline void dfs1(int rt,int d,int pre)
{
int max_siz=0;
fa[rt]=pre;
dep[rt]=d+1;
siz[rt]=1;
for(int i=head[rt];i;i=e[i].nt)
{
int y=e[i].to;
if(fa[y])
continue;
dfs1(y,dep[rt],rt);
siz[rt]+=siz[y];
if(siz[y]>max_siz)
{
max_siz=siz[y];
son[rt]=y;
}
}
}
int dfn;
inline void dfs2(int rt,int ld)
{
id[rt]=++dfn;
top[rt]=ld;
w[dfn]=r[rt];
if(!son[rt])
return;
dfs2(son[rt],ld);
for(int i=head[rt];i;i=e[i].nt)
{
int y=e[i].to;
if(y==fa[rt]||y==son[rt])
continue;
dfs2(y,y);
}
}
int t[maxn*4],tag[maxn*4];
inline void push_up(int p)
{
t[p]=t[p<<1|1]+t[p<<1];
}
inline void push_down(int p,int l,int r)
{
tag[p<<1|1]+=tag[p];
tag[p<<1]+=tag[p];
int mid=(l+r)>>1;
t[p<<1]+=tag[p]*(mid-l+1)%mod;
t[p<<1]%=mod;
t[p<<1|1]+=tag[p]*(r-mid)%mod;
t[p<<1|1]%=mod;
tag[p]=0;
}
inline void build(int p,int l,int r)
{
if(l==r)
{
t[p]=w[l]%mod;
return;
}
int mid=(l+r)>>1;
build(p<<1,l,mid);
build(p<<1|1,mid+1,r);
push_up(p);
}
inline void update(int p,int l,int r,int ul,int ur,int x)
{
if(ul<=l&&ur>=r)
{
tag[p]+=x%mod;
t[p]+=(r-l+1)*x;
t[p]%=mod;
return;
}
if(tag[p])
push_down(p,l,r);
int mid=(l+r)>>1;
if(ul<=mid)
update(p<<1,l,mid,ul,ur,x);
if(ur>mid)
update(p<<1|1,mid+1,r,ul,ur,x);
push_up(p);
}
inline int query(int p,int l,int r,int ql,int qr)
{
int res=0;
if(l>=ql&&r<=qr)
{
return t[p]%mod;
}
if(tag[p])
push_down(p,l,r);
int mid=(l+r)>>1;
if(ql<=mid)
res+=query(p<<1,l,mid,ql,qr)%mod;
if(qr>mid)
res+=query(p<<1|1,mid+1,r,ql,qr)%mod;
return res%mod;
}
int main()
{
n=read(),m=read(),rt=read(),mod=read();
for(int i=1;i<=n;i++)
{
r[i]=read();
}
for(int i=1;i<n;i++)
{
int x,y;
x=read(),y=read();
add_eage(x,y);
add_eage(y,x);
}
dfs1(rt,0,rt);
dfs2(rt,rt);
build(1,1,n);
for(int i=1;i<=m;i++)
{
int opt,x,y,z;
opt=read();
if(opt==1)
{
x=read(),y=read(),z=read();
while(top[x]!=top[y])
{
if(dep[top[x]]<dep[top[y]])
swap(x,y);
update(1,1,n,id[top[x]],id[x],z);
x=fa[top[x]];
}
if(dep[x]<dep[y])
swap(x,y);
update(1,1,n,id[y],id[x],z);
}
else if(opt==2)
{
x=read(),y=read();
int res=0;
while(top[x]!=top[y])
{
if(dep[top[x]]<dep[top[y]])
swap(x,y);
res+=query(1,1,n,id[top[x]],id[x])%mod;
x=fa[top[x]];
}
if(dep[x]<dep[y])
swap(x,y);
res+=query(1,1,n,id[y],id[x])%mod;
printf("%d\n",res);
}
else if(opt==3)
{
x=read(),y=read();
update(1,1,n,id[x],id[x]+siz[x]-1,y);
}
else
{
x=read();
printf("%d\n",query(1,1,n,id[x],id[x]+siz[x]-1));
}
}
return 0;
}