# include <bits/stdc++.h>
using namespace std;
const int mod=64123;
typedef long long ll;
ll inv[2010];
int n,K,id[2010],ans[2010],d[2010],a[2010],f[2010],g[2010],Lg[2010][2010],sz[2010],fa[2010],son[2010],top[2010],dwn[2010],dfn[2010],idfn[2010],tot;
bool b[2010];
vector<int> G[2010];
ll power(ll a,int b)
{
	ll ans=1;
	for(;b;b>>=1,a=a*a%mod)
		if(b&1) ans=ans*a%mod;
	return ans;
}
struct Num
{
	int v,c;
	Num(int _v=1,int _c=0):v(_v),c(_c){}
	void operator*=(int x)
	{
		if(x) v=(ll)v*x%mod;
		else c++;
	}
	void operator/=(int x)
	{
		if(x) v=(ll)v*power(x,mod-2)%mod;
		else c--;
	}
	operator int(){return c?0:v;}
}Lf[2010][2010];
struct Mat
{
	int a,b,c,d;
	Mat(int _a=0,int _b=0,int _c=0,int _d=0):a(_a),b(_b),c(_c),d(_d){}
	Mat operator+(const Mat &x)const{return Mat((a+x.a)%mod,(b+x.b)%mod,(c+x.c)%mod,(d+x.d)%mod);}
	Mat operator*(const Mat &x)const{return Mat((ll)a*x.a%mod,((ll)a*x.b+b)%mod,((ll)c*x.a+x.c)%mod,((ll)c*x.b+d+x.d)%mod);}
};
typedef vector<Mat> Info;
Info val[5010],M[2010];
void work(int u)
{
	if(!b[u]) for(int i=0;i<=n;i++) M[u][i]=Mat(Lf[u][i],1,Lf[u][i],Lg[u][i]);
	else for(int i=0;i<=n;i++) M[u][i]=Mat((ll)Lf[u][i]*i%mod,1,(ll)Lf[u][i]*i%mod,Lg[u][i]);
}
Info operator*(const Info &a,const Info &b)
{
	Info ans(n+1);
	for(int i=0;i<=n;i++) ans[i]=a[i]*b[i];
	return ans;
}
void dfs1(int u)
{
	sz[u]=1;
	for(int v:G[u])
	{
		if(v==fa[u]) continue;
		fa[v]=u;dfs1(v);
		sz[u]+=sz[v];
		if(sz[v]>sz[son[u]]) son[u]=v;
	}
}
void dfs2(int u)
{
	dfn[u]=++tot;idfn[tot]=u;
	if(u==son[fa[u]]) top[u]=top[fa[u]];
	else top[u]=u;
	if(son[u]) dfs2(son[u]);
	for(int v:G[u])
		if(v!=fa[u] && v!=son[u]) dfs2(v);
	if(son[u]) dwn[u]=dwn[son[u]];
	else dwn[u]=u;
}
void dfs3(int u)
{
	f[u]=1;
	for(int v:G[u])
	{
		if(v==fa[u]) continue;
		dfs3(v);
		f[u]=(ll)f[u]*f[v]%mod;g[u]=(g[u]+g[v])%mod;
	}
	g[u]=(g[u]+f[u])%mod;f[u]++;
}
void pushup(int rt){for(int i=0;i<=n;i++) val[rt][i]=val[rt*2][i]*val[rt*2+1][i];}
void build(int rt,int l,int r)
{
	if(l==r) return val[rt]=M[idfn[l]],void();
	int mid=(l+r)/2;
	build(rt*2,l,mid);build(rt*2+1,mid+1,r);
	val[rt].resize(n+1);pushup(rt);
}
void update(int rt,int l,int r,int x)
{
	if(l==r) return val[rt]=M[idfn[l]],void();
	int mid=(l+r)/2;
	if(x<=mid) update(rt*2,l,mid,x);
	else update(rt*2+1,mid+1,r,x);
	pushup(rt);
}
Info query(int rt,int l,int r,int x,int y)
{
	if(x<=l && r<=y) return val[rt];
	int mid=(l+r)/2;
	if(y<=mid) return query(rt*2,l,mid,x,y);
	else if(x>mid) return query(rt*2+1,mid+1,r,x,y);
	else return query(rt*2,l,mid,x,y)*query(rt*2+1,mid+1,r,x,y);
}
void calc(int u)
{
	Info t=query(1,1,n,dfn[u],dfn[dwn[u]]);
	for(int i=0;i<=n;i++) f[i]=(t[i].a+t[i].b)%mod,g[i]=(t[i].c+t[i].d)%mod;
}
void update(int u)
{
	b[u]=1;work(u);
	for(;u;u=fa[top[u]])
		if(top[u]!=1)
		{
			int v=top[u];calc(v);
			for(int i=0;i<=n;i++) Lf[fa[v]][i]/=f[i],Lg[fa[v]][i]=(Lg[fa[v]][i]-g[i]+mod)%mod;
			update(1,1,n,dfn[u]);calc(v);
			for(int i=0;i<=n;i++) Lf[fa[v]][i]*=f[i],Lg[fa[v]][i]=(Lg[fa[v]][i]+g[i])%mod;
			work(fa[v]);
		}
		else update(1,1,n,dfn[u]);
}
int Lagrange()
{
	static int F[2010],mul[2010];
	mul[0]=1;
	for(int i=0;i<=n;i++)
	{
		for(int j=i+1;j;j--) mul[j]=(mul[j-1]+(ll)(mod-i)*mul[j])%mod;
		mul[0]=(ll)(mod-i)*mul[0]%mod;
	}
	for(int i=0;i<=n;i++)
	{
		int cur=mul[n+1],w=inv[i]*((n-i)&1?mod-inv[n-i]:inv[n-i])%mod;
		for(int j=n;j>=0;j--) F[j]=(F[j]+(ll)w*cur%mod*ans[i])%mod,cur=((ll)cur*i+mul[j])%mod;
	}
	int aans=0;
	for(int i=K;i<=n;i++) aans=(aans+F[i])%mod;
	return aans;
}
int main()
{
	int u,v;
	scanf("%d%d%*d",&n,&K);
	inv[0]=inv[1]=1;
	for(int i=2;i<=n;i++) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
	for(int i=2;i<=n;i++) inv[i]=inv[i-1]*inv[i]%mod;
	for(int i=1;i<=n;i++) scanf("%d",&d[i]),a[i]=d[i],id[i]=i;
	sort(a+1,a+n+1);sort(id+1,id+n+1,[](int i,int j){return d[i]>d[j];});
	int len=unique(a+1,a+n+1)-a-1;
	for(int i=1;i<n;i++)
	{
		scanf("%d%d",&u,&v);
		G[u].push_back(v);G[v].push_back(u);
	}
	dfs1(1);dfs2(1);dfs3(1);
	for(int i=1;i<=n;i++)
	{
		Num lf;int lg=0;
		for(int j:G[i])
			if(j!=fa[i] && j!=son[i]) lf*=f[j],lg=(lg+g[j])%mod;
		fill(Lf[i],Lf[i]+n+1,lf);fill(Lg[i],Lg[i]+n+1,lg);
		M[i].resize(n+1);work(i);
	}
	build(1,1,n);
	for(int i=len,j=1;i>=1;i--)
	{
		for(;j<=n && d[id[j]]>=a[i];j++) update(id[j]);
		calc(1);
		for(int k=0;k<=n;k++) ans[k]=(ans[k]+(ll)(a[i]-a[i-1])*g[k])%mod;
	}
	cout<<Lagrange()<<endl;
	return 0;
}