#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;

const int MAXN = 4e5;

int N, M, P[MAXN+10];

const int MOD = 998244353;

struct mint
{
    int x;
    mint() : x(0) {}
    mint(int _x) : x(_x) {}
    mint operator + (int ot) const { return x+ot>=MOD ? x+ot-MOD : x+ot; }
    mint operator - (int ot) const { return x<ot ? x+MOD-ot : x-ot; }
    mint operator * (int ot) const { return 1ll*x*ot%MOD; }
    mint operator += (int ot) { return *this = *this + ot; }
    mint operator -= (int ot) { return *this = *this - ot; }
    mint operator *= (int ot) { return *this = *this * ot; }
    operator int() const { return x; }
};
mint mpow(mint a, int b)
{
    mint ret=1;
    while(b)
    {
        if(b&1) ret=ret*a;
        a=a*a; b>>=1;
    }
    return ret;
}

const mint G = 3;

void dft(vector<mint> &F, bool inv)
{
    int n=F.size();
    int b=__lg(n);
    vector<int> rev(n);
    for(int i=0; i<n; i++)
    {
        rev[i]=(rev[i/2]>>1) | ((i&1)<<(b-1));
        if(i<rev[i]) swap(F[i], F[rev[i]]);
    }
    vector<mint> w(n);
    w[0]=w[1]=1;
    for(int k=2; k<n; k<<=1)
    {
        mint z[2]={1, mpow(G, inv ? MOD-1-(MOD-1)/k/2 : (MOD-1)/k/2)};
        for(int i=k; i<(k<<1); i++) w[i]=w[i>>1]*z[i&1];
    }
    for(int d=1; d<n; d<<=1)
    {
        for(int i=0; i<n; i+=d+d)
        {
            for(int j=0; j<d; j++)
            {
                mint b=F[i|j|d]*w[j|d];
                F[i|j|d]=F[i|j]-b;
                F[i|j]+=b;
            }
        }
    }
    if(inv)
    {
        mint val=mpow(n, MOD-2);
        for(int i=0; i<n; i++) F[i]*=val;
    }
}

vector<mint> multiply(vector<mint> F, vector<mint> G)
{
    int n=1;
    int sz=F.size()+G.size()-1;
    while(n<sz) n<<=1;
    F.resize(n); G.resize(n);
    dft(F, false); dft(G, false);
    for(int i=0; i<n; i++) F[i]*=G[i];
    dft(F, true);
    F.resize(sz);
    return F;
}

mint fac[MAXN+10], ifac[MAXN+10];
mint ncr(int n, int r)
{
    if(0<=r && r<=n) return fac[n]*ifac[r]*ifac[n-r];
    return 0;
}

int main2()
{
    cin >> N >> M;
    for(int i=1; i<=N*M; i++) cin >> P[i];
    sort(P+1, P+N*M+1);

    vector<mint> A(N*M+1), B(N*M+1), C(N*M+1);
    for(int i=0; i<=N; i++) for(int j=0; j<=M; j++) A[i*j]+=ncr(N, i)*ncr(M, j)*((N+M+i+j)%2 ? mint(MOD-1) : mint(1))*fac[i*j];
    for(int i=0; i<=N*M; i++) B[i]=ifac[N*M-i];

    auto V = multiply(A, B);
    for(int i=0; i<=N*M; i++) C[N*M-i]=V[i+N*M]*fac[N*M-i];

    // for(int i=0; i<=N*M; i++) cout << A[i] << " ";
    // cout << " : A\n";
    // for(int i=0; i<=N*M; i++) cout << B[i] << " ";
    // cout << " : B\n";
    // for(int i=0; i<=N*M; i++) cout << C[i] << " ";
    // cout << " : C\n";

    mint ans=0;
    for(int i=0; i<N*M; i++) ans+=C[i]*mint(P[i+1]-P[i]);

    cout << ans << "\n";

    return 0;
}

void reset()
{

}

int main()
{
    // ios_base::sync_with_stdio(false); cin.tie(NULL);

    fac[0]=1;
    for(int i=1; i<=MAXN; i++) fac[i]=fac[i-1]*i;
    ifac[MAXN]=mpow(fac[MAXN], MOD-2);
    for(int i=MAXN; i>=1; i--) ifac[i-1]=ifac[i]*i;

    int TC;
    cin >> TC;
    while(TC--)
    {
        main2();
        reset();
    }
}