#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define MP make_pair
const int MAXN=3e3+5;
const int MOD=998244353;
int n,k,a[MAXN],fa[MAXN],dep[MAXN],dis[MAXN];
bool vis[MAXN];
ll s[MAXN],f[MAXN];
vector<int> edg[MAXN],blk[MAXN]; 
int find(int x){
	while(x^fa[x]) x=fa[x]=fa[fa[x]];
	return x;
}
void dfs(int u,int fa){
	for(int v:edg[u]) if(v!=fa){
		dep[v]=dep[u]+1;
		dfs(v,u);
	}
}
void dfs0(int u,int fa,int o){
	if(a[u]==o) s[dis[u]]+=f[u];
	for(int v:edg[u]) if(v!=fa)
		dis[v]=dis[u]+1,dfs0(v,u,o);
}
void dfs1(int id){
	for(int u:blk[id]) for(int v:edg[u]) if(a[v]!=id&&dep[v]>dep[u]){
		dfs1(a[v]);
		memset(s,0,sizeof(s));
		dis[u]=0;dfs0(u,0,a[v]); 
		for(int w:blk[id]){
			if(a[v]>id){
				f[w]=(s[dis[w]]+s[dis[w]+1])%MOD*f[w]%MOD;
			}else{
				f[w]=(s[dis[w]+1]+s[dis[w]+2])%MOD*f[w]%MOD;
			}
		}
	}
}
void solve(){
	memset(vis,0,sizeof(vis));
	cin>>n>>k;
	for(int i=1;i<=n;i++) edg[i].clear(),blk[i].clear(),fa[i]=i;
	for(int i=1;i<n;i++){
		int u,v;cin>>u>>v;
		edg[u].push_back(v);
		edg[v].push_back(u);
	}
	for(int i=1;i<=n;i++){
		cin>>a[i];
		vis[a[i]]=true;
		blk[a[i]].push_back(i);
	}
	for(int i=1;i<=k;i++) if(!vis[i]){
		cout<<0<<endl;return;
	}
	int cnt=n;
	for(int i=1;i<=n;i++) for(int j:edg[i])
		if(a[i]==a[j]&&find(i)!=find(j)) fa[find(i)]=find(j),cnt--; 
	if(cnt!=k){cout<<0<<endl;return;}
	for(int i=1;i<=n;i++) f[i]=1;
	dfs(1,0);dfs1(a[1]);
	ll ans=0;
	for(int i:blk[1]) ans+=f[i];
	cout<<ans%MOD<<endl;
	return;
}
int main(){
	ios::sync_with_stdio(false);
	int _;cin>>_;
	while(_--) solve();
	return 0;
}