#include<bits/stdc++.h>
using namespace std;
const int mod=998244353,_=3005;
int c[_],cnt,d[_],f[_][_],g[_][11],h[_],k,lg[_],n,p[_],s,T,u[_],v[_],w[_],z;
vector<int>a[_],b[_],o[_];
int get(int x,int y){
	return d[x]<d[y]?x:y;
}
void dfs(int x,int fa){
	d[x]=d[fa]+1;h[x]=fa;w[x]=++cnt;g[cnt][0]=x;
	for(auto y:a[x])
		if(y!=fa)dfs(y,x);
}
int dis(int x,int y){
	if(x==y)return 0;
	if(w[y]<w[x])swap(x,y);
	int l=lg[w[y]-w[x]];
	return d[x]+d[y]-d[h[get(g[w[x]+1][l],g[w[y]-(1<<l)+1][l])]]*2;
}
void dp(int x,int fa){
	for(auto i:o[x])f[x][i]=1;
	for(auto y:b[x])
		if(y!=fa){
			dp(y,x);
			for(auto i:o[x]){
				int sum=0;
				for(auto j:o[y]){
					int dj=dis(j,p[y]),di=dis(i,h[p[y]]);
					if(dj==di||x>y&&dj==di+1||x<y&&dj+1==di)sum=(sum+f[y][j])%mod;
				}
				f[x][i]=1ll*f[x][i]*sum%mod;
			}
		}
}
int main(){
	//freopen("qwq.in","r",stdin);
	//freopen("qwq.out","w",stdout);
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	cin>>T;
	while(T--){
		cin>>n>>k;z=s=d[0]=cnt=0;lg[0]=-1;
		for(int i=1;i<=n;i++){
			a[i].clear();b[i].clear();o[i].clear();p[i]=0;
			for(int j=1;j<=n;j++)f[i][j]=0,lg[i]=lg[i>>1]+1; 
		}
		memset(f,0,sizeof(f));z=0;
		for(int i=1;i<n;i++)cin>>u[i]>>v[i],a[u[i]].push_back(v[i]),a[v[i]].push_back(u[i]);
		for(int i=1;i<=n;i++)cin>>c[i],o[c[i]].push_back(i);
		dfs(1,0);d[0]=1e9;
		for(int i=1;i<=10;i++)
			for(int j=1;j<=n-(1<<i-1);j++)g[j][i]=get(g[j][i-1],g[j+(1<<i-1)][i-1]);
		for(int i=1;i<n;i++)
			if(c[u[i]]^c[v[i]])b[c[u[i]]].push_back(c[v[i]]),b[c[v[i]]].push_back(c[u[i]]),z++;
		for(int i=1;i<=n;i++)
			if(d[p[c[i]]]>d[i])p[c[i]]=i;
		for(int i=1;i<=k;i++)
			if(!p[i])z=0;
		if(z!=k-1){
			cout<<"0\n";
			continue;
		}
		dp(c[1],0);
		//for(int i=1;i<=k;i++)
		//	for(auto j:o[i])cout<<i<<' '<<j<<' '<<f[i][j]<<'\n';
		for(auto i:o[c[1]])s=(s+f[c[1]][i])%mod;
		cout<<s<<'\n';
	}
	return 0;
}