#include <bits/stdc++.h>
using namespace std;
#define ll long long
const int mod=998244353;
int add(int x,int y){x+=y;return x>=mod?x-mod:x;}
int sub(int x,int y){x-=y;return x<0?x+mod:x;}
int mul(int x,int y){return (ll)x*y%mod;}
void ckadd(int &x,int y){x=add(x,y);}
void cksub(int &x,int y){x=sub(x,y);}

const int N=400050;
int F[N],I[N];
int binom(int n,int k){
    return mul(F[n],mul(I[k],I[n-k]));
}
int c[N];
vector<int> vls;
void precalc(int n,int m){
	for(int k=0;k<=n*m;k++)c[k]=0;
	vls.clear();
    for(int i=0;i<n;i++){
        for(int j=0;j<m;j++){
        	int now=mul(binom(n,i),binom(m,j));
        	int sm=(n-i)*(m-j);
        	if((i+j)%2==0)ckadd(c[sm],now);
        	else cksub(c[sm],now);
        }
    }
	for(int k=0;k<=n*m;k++)if(c[k]!=0)vls.push_back(k);
}
int SolveInclusionExclusion(int cnt){
	int ans=0;
	for(int j=(int)vls.size()-1;j>=0;j--){
		int sm=vls[j];
		if(sm<cnt)break;
		ckadd(ans,mul(binom(sm,cnt),c[sm]));
	}
    return ans;
}

int Solve(int n,int m,int cnt){
    vector<vector<int>> dp(n+1,vector<int>(m+1,0));
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            if(i*j>=cnt){
                dp[i][j]=binom(i*j,cnt);
                for(int a=0;a<=i;a++){
                    for(int b=0;b<=j;b++){
                        if(a==0 && b==0)continue;
                        if(a==i && b!=j)continue;
                        if(b==j && a!=i)continue;
                        cksub(dp[i][j],mul(dp[i-a][j-b],mul(binom(i,a),binom(j,b))));
                    }
                }
            }
        }
    }
    return dp[n][m];
    /*if(cnt>n*m)return 0;
    if(n>m)swap(n,m);
    int ans=binom(n*m,cnt);
    for(int i=1;i<=n;i++){
        int now=mul(Solve(n-i,m,cnt),binom(n,i));
        if(i%2==1)cksub(ans,now);
        else ckadd(ans,now);
    }
    return ans;*/
}

int main(){
    F[0]=1;
    for(int i=1;i<N;i++)F[i]=mul(F[i-1],i);
    I[0]=I[1]=1;
    for(int i=2;i<N;i++)I[i]=mul(mod-mod/i,I[mod%i]);
    for(int i=1;i<N;i++)I[i]=mul(I[i-1],I[i]);
    int t;
    scanf("%i",&t);
    while(t--){
        int n,m;
        scanf("%i %i",&n,&m);
        precalc(n,m);
        vector<int> a(n*m);
        for(int i=0;i<n*m;i++){
            scanf("%i",&a[i]);
        }
        sort(a.begin(),a.end());
        int ans=0;
        for(int i=0;i<n*m;i++){
            int diff=(i==0?a[i]:a[i]-a[i-1]);
            if(diff==0)continue;
            int ways=mul(mul(F[n*m-i],F[i]),SolveInclusionExclusion(n*m-i));
            ckadd(ans,mul(diff,ways));
            //printf("bingo>=%i ways=%i\n",a[i],ways);
        }
        printf("%i\n",ans);
    }
    return 0;
}