#include<bits/stdc++.h>
using namespace std;
int K;
int st[50],t[50];
int jc[1010],ij[1010];
int mod;
void dfs(int x,int z,int b);
int g[50][40],f[110][40];
void ad(int &x,int y){
	x+=y;if(x>=mod)x-=mod;
}
void dfs(int x,int z,int b){
	if(x==b){
		b=0;
		for(int i=x-1;i;i--)if((st[x]>>i)&1)b=i;
		if(!b)return;
		t[b]++;
		if(z){
			int A=jc[x-1];
			for(int i=1;i<=x;i++)A/=jc[t[i]];
			g[x][z]+=A;
		}
			//计入答案 
			if(x<K+2)dfs(x+1,z,b);
			t[b]--;
		return;
	}
	int z_=z+__builtin_popcount(st[x]&st[b]);
	//for(int i=1;i<b;i++)if(a[i][x]&&a[i][b])z_++;
	if(z_<=K)dfs(x,z_,b+1);
	z_=z+__builtin_popcount((st[x]^st[b])&((1<<b)-1));
	//for(int i=1;i<b;i++)if(a[i][x]^a[i][b])z_++;
	st[x]|=(1<<b),st[b]|=(1<<x);
	if(z_<=K)dfs(x,z_,b+1);
	st[x]^=(1<<b),st[b]^=(1<<x);
}
int n,c[1010][1010],B[1010];
void sol(){
	scanf("%d%d%d",&n,&K,&mod);
	c[0][0]=B[0]=1;
	for(int i=1;i<=n;i++){
		c[i][0]=1;
		for(int j=1;j<=i;j++)
		c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;
		B[i]=0;
		for(int j=1;j<=i;j++)
		ad(B[i],1ll*B[i-j]*c[i-1][j-1]%mod);
	}
	memset(f,0,sizeof(f));
	f[0][0]=1;
	for(int i=1;i<=K*3&&i<=n;i++){
		for(int j=1;j<=i;j++)for(int z=1;z<=K;z++)if(g[j][z])
		for(int o=0;o+z<=K;o++)if(f[i-j][o])
		ad(f[i][o+z],1ll*f[i-j][o]*g[j][z]%mod*c[i-1][j-1]%mod);
	}
	int ans=0;
	for(int i=0;i<=K*3&&i<=n;i++){
		int e=0;
		for(int j=0;j<=K;j++)ad(e,f[i][j]);
		ad(ans,1ll*c[n][i]*B[n-i]%mod*e%mod);
	}
	printf("%d\n",ans);
}
int main(){
	jc[0]=1;
	for(int i=1;i<=10;i++)jc[i]=jc[i-1]*i;
	K=8,dfs(2,0,1);
	int T;scanf("%d",&T);while(T--)sol();
	return 0;
}