#include <bits/stdc++.h>

#define Max 100005

using namespace std;

const int mod=998244353;

int T,n,k,N1,N2,son[Max];
int f[Max][350],F[Max][350],g[Max][350];
vector<int>E[Max];

inline int dfs(int u,int fa){
	int p=1;
	son[u]=1;
	f[u][1]=1;
	for(int v:E[u]){
		if(v==fa)continue;
		int siz=dfs(v,u);
		for(int i=min(p,k+1);i>=0;i--){
			for(int j=0;j<=siz&&i+j<=k+1;j++){
				F[u][i+j]=(F[u][i+j]+1ll*f[u][i]*f[v][j])%mod;
			}
			f[u][i]=0;
		}
		for(int i=0;i<=siz+p&&i<=k+1;i++){
			f[u][i]=F[u][i];
			F[u][i]=0;
//			cout<<u<<" "<<i<<" "<<f[u][i]<<endl;
		}
		son[u]+=son[v];
		p+=siz;
	}
	f[u][0]=((f[u][0]+f[u][k]%mod)+f[u][k+1])%mod;
	f[u][k+1]=0;
//	for(int i=0;i<=k+1;i++)cout<<u<<" "<<i<<" "<<f[u][i]<<endl;
	return p;
}

inline int dfs1(int u,int fa){
	f[u][0]=1;
	g[u][0]=0;
	int p=0;
	for(int v:E[u]){
		if(v==fa)continue;
		int siz=dfs1(v,u);
		for(int i=min(p,N1+1);i>=0;i--){
			for(int j=0;j<=siz&&i+j<=N1+1;j++){
				F[u][i+j]=(F[u][i+j]+1ll*f[u][i]*f[v][j])%mod;
				g[u][i+j]=g[u][i]+g[v][j];
			}
		}
		for(int i=0;i<=siz+p&&i<=N1+1;i++){
			f[u][i]=F[u][i];
			F[u][i]=0;
		}
		p+=siz;
	}
	for(int i=p;i>=0;i--){
		int now=son[u]-(k+1)*g[u][i]-k*i;
		if(now==k+1){
			g[u][i]=g[u][i]+1;
		}else if(now==k){
			f[u][i+1]=(f[u][i+1]+f[u][i])%mod;
			g[u][i+1]=g[u][i];
		}
	}
	return p;
}

int main(){
	ios::sync_with_stdio(false);
	cin>>T;
	while(T--){
		cin>>n>>k;
		for(int i=1;i<=n;i++)E[i].clear();
		for(int i=1;i<n;i++){
			int u,v;
			cin>>u>>v;
			E[u].push_back(v);
			E[v].push_back(u);
		}
		if(k<=sqrt(n)){
			dfs(1,1);
			cout<<f[1][0]<<'\n';
		}else{
			N1=n/k;
			N2=n%k;
			N1-=N2;
			dfs1(1,1);
			cout<<f[1][N1]<<'\n';
		}
	}
	return 0;
}