#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=3030,M=1e6+6,mod=998244353;
int n,m,p[N],vis[N],cnt[N];
LL f[N],ans,fac[N];
LL ksm(LL x,int y)
{
	LL ret=1;
	while(y)
	{
		if(y&1)
		{
			ret=ret*x%mod;
		}
		x=x*x%mod;
		y>>=1;
	}
	return ret;
}
int gcd(int x,int y)
{
	if(y==0)
	{
		return x;
	}
	return gcd(y,x%y);
}
void work()
{
	ans=1;
	scanf("%d%d",&n,&m);
	fac[0]=1;
	for(int i=1;i<=n;i++)
	{
		fac[i]=fac[i-1]*i%mod;
		vis[i]=cnt[i]=0;
		scanf("%d",&p[i]);
	}
	for(int i=1;i<=n;i++)
	{
		if(vis[i]==0)
		{
			int sz=0;
			vis[i]=1;
			for(int j=p[i];j;j=p[j])
			{
				vis[j]=1;
				sz++;
				if(j==i)
				{
					break;
				}
			}
			cnt[sz]++;
		}
	}
	for(int i=1;i<=n;i++)
	{
		f[0]=1;
		for(int j=1;j<=cnt[i];j++)
		{
			f[j]=0;
		}
	//	cout<<i<<endl;
		for(int j=1;j<=cnt[i];j++)
		{
			for(int k=1;k<=j;k++)
			{
				int len=k*i;
				if(len/gcd(len,m)==i)
				{
				//	cout<<k<<" legal"<<endl;
					f[j]+=f[j-k]*fac[k-1]%mod*ksm(i,k-1)%mod;
					f[j]%=mod;
				}
			}
		//	cout<<j<<' '<<f[j]<<endl;
		}
		ans*=f[cnt[i]];
		ans%=mod;
	}
	printf("%lld\n",ans);
}
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)
	{
		work();
	}
	return 0;
}