#include<bits/stdc++.h>
using namespace std;
using ll=long long;
const int N=3e6+10,V=9e6+10,mod=998244353,i2=(mod+1)/2;
int T,n,m,k,ans,a[N];
int fac[N],inv[N];
int cnt,vis[V],pri[V],mn[V];
int nex[N],mx[V],now[V];
int p[N],q[N],ta[N],tb[N],fa[N*2];
int find(int x){
	return fa[x]==x?x:fa[x]=find(fa[x]);
}
void init(int n=9e6,int m=3e6){
	for(int i=2;i<=n;i++){
		if(!vis[i])pri[++cnt]=i,mn[i]=i;
		for(int j=1;j<=cnt&&i*pri[j]<=n;j++){
			vis[i*pri[j]]=1,mn[i*pri[j]]=pri[j];
			if(i%pri[j]==0)break;
		}
	}
	for(int i=fac[0]=1;i<=m;i++)fac[i]=1ll*fac[i-1]*i%mod;
	inv[1]=1;
	for(int i=2;i<=m;i++)inv[i]=1ll*inv[mod%i]*(mod-mod/i)%mod;
	memset(vis,0,sizeof vis);
}
#define id(i,j) (((i)-1)*m+(j))
int trs(int n,int i){
	if(i&1)return (i+1)/2;
	return i/2+(n+1)/2;
}
void get(){
	scanf("%d%d",&n,&m),k=n*m,ans=1;
	for(int i=1;i<=k;i++)scanf("%d",&a[i]);
	for(int i=1;i<=n;i++)p[i]=trs(n,i);
	for(int i=1;i<=m;i++)q[i]=trs(m,i);
	vector<vector<int> >P,Q;
	for(int i=1;i<=n;i++)if(p[i]==i){
		int lcm=1;
		vector<int>ers;
		for(int j=1;j<=m;j++)if(!vis[j]){
			vector<int>col{0};
			for(int u=j;!vis[u];u=q[u])vis[u]=1,col.push_back(a[id(i,u)]);
			int cnt=col.size()-1;
			for(int i=2,j=0;i<=cnt;i++){
				for(;j&&col[j+1]!=col[i];j=nex[j]);
				nex[i]=j+=col[j+1]==col[i];
			}
			int len,las;
			if(cnt%(cnt-nex[cnt])||nex[cnt]*2<cnt)len=cnt;
			else len=cnt-nex[cnt];
			las=len;
			ers.push_back(len);
			for(;len>1;len/=mn[len]){
				if(++now[mn[len]]>mx[mn[len]]){
					lcm=1ll*lcm*mn[len]%mod;
					++mx[mn[len]];
				}
			}
			for(len=las;len>1;len/=mn[len])now[mn[len]]=0;
		}
		ans=1ll*ans*lcm%mod;
		fill(vis,vis+1+m,0);
		for(int len:ers)for(;len>1;len/=mn[len])mx[mn[len]]=0;
	}
	for(int j=1;j<=m;j++)if(q[j]==j){
		int lcm=1;
		vector<int>ers;
		for(int i=1;i<=n;i++)if(!vis[i]){
			vector<int>col{0};
			for(int u=i;!vis[u];u=p[u])vis[u]=1,col.push_back(a[id(u,j)]);
			int cnt=col.size()-1;
			for(int i=2,j=0;i<=cnt;i++){
				for(;j&&col[j+1]!=col[i];j=nex[j]);
				nex[i]=j+=col[j+1]==col[i];
			}
			int len,las;
			if(cnt%(cnt-nex[cnt])||nex[cnt]*2<cnt)len=cnt;
			else len=cnt-nex[cnt];
			ers.push_back(len);
			las=len;
			for(;len>1;len/=mn[len]){
				if(++now[mn[len]]>mx[mn[len]]){
					lcm=1ll*lcm*mn[len]%mod;
					mx[mn[len]]++;
				}
			}
			for(len=las;len>1;len/=mn[len])now[mn[len]]=0;
		}
		ans=1ll*ans*lcm%mod;
		fill(vis,vis+1+n,0);
		for(int len:ers)for(;len>1;len/=mn[len])mx[mn[len]]=0;
	}
	for(int i=1;i<=n;i++)if(!vis[i]){
		P.push_back({});
		for(int u=i;!vis[u];u=p[u])vis[u]=1,P.back().push_back(u);
	}
	fill(vis,vis+1+n,0);
	for(int i=1;i<=m;i++)if(!vis[i]){
		Q.push_back({});
		for(int u=i;!vis[u];u=q[u])vis[u]=1,Q.back().push_back(u);
	}
	fill(vis,vis+1+m,0);
	iota(fa,fa+1+n+m,0);
	for(auto &a:P)for(auto &b:Q){
		if(a.size()==1||b.size()==1)continue;
		bool flag=0;
		for(int i:a)for(int j:b){
			flag|=vis[::a[id(i,j)]];
			ans=1ll*ans*inv[++vis[::a[id(i,j)]]]%mod;
		}
		int mx=0;
		for(int i:a)for(int j:b){
			mx=max(mx,vis[::a[id(i,j)]]);
		}
		if(mx==a.size()*b.size()){
			ans=1ll*ans*fac[a.size()*b.size()]%mod;
		}
		else if(flag){
			ans=1ll*ans*fac[a.size()*b.size()]%mod;
		}else if(a.size()%2==1&&b.size()%2==1){
			ans=1ll*ans*fac[a.size()*b.size()]%mod*i2%mod;
		}else if(a.size()%2==1){
			if(!ta[a.back()])ans=ans*2%mod;
			ans=1ll*ans*fac[a.size()*b.size()]%mod*i2%mod;
			ta[a.back()]=1;
		}else if(b.size()%2==1){
			if(!tb[b.back()])ans=ans*2%mod;
			ans=1ll*ans*fac[a.size()*b.size()]%mod*i2%mod;
			tb[b.back()]=1;
		}else{
			ans=1ll*ans*fac[a.size()*b.size()]%mod*i2%mod;
			int fx=find(a.back()),fy=find(b.back()+n);
			if(fx!=fy)ans=ans*2%mod,fa[fx]=fy;
		}
		for(int i:a)for(int j:b)vis[::a[id(i,j)]]=0;
	}
	fill(ta,ta+1+n,0),fill(tb+1,tb+1+m,0);
	cout<<ans<<endl;
}
int main(){
	init();
	for(scanf("%d",&T);T--;)get();
	return 0;
}