#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define REP(i,n) for(ll i=0;i<(n);i++)
const long long MOD = 998244353;

vector<ll> fact(200200);

ll inverse(ll x){
    ll e = MOD-2;
    ll b = x, r = 1;
    while(e>0){
        if(e%2==1){
            r = (r*b)%MOD;
        }
        e /= 2;
        b = (b*b)%MOD;
    }
    return r;
}

ll comb(ll n, ll r){
    ll Ans = fact[n];
    Ans = (Ans*inverse(fact[r]))%MOD;
    Ans = (Ans*inverse(fact[n-r]))%MOD;
    return Ans;
}

// N,M のグリッドに K 個置く。ただし縦横ビンゴができてはならない
ll solve2(ll N, ll M, ll K){

    if(!(0<=K && K<=N*M))return 0;
    ll Ans = 0;
    vector<ll> d(N*M,0);
    REP(i,K){d[N*M-1-i]=1;}

    do{
        bool c = true;
        REP(i,N){
            bool e = true;
            REP(j,M){
                if(d[i*M+j]==0)e=false;
            }
            if(e)c=false;
        }
        REP(j,M){
            bool e = true;
            REP(i,N){
                if(d[i*M+j]==0)e=false;
            }
            if(e)c=false;
        }

        if(c)Ans++;
    }while(next_permutation(d.begin(),d.end()));

    Ans %= MOD;
    return Ans;

}

// N,M のグリッドに K 個置く。ただし縦横ビンゴができてはならない
// ll solve2(ll N, ll M, ll K){

//     ll Ans = 1;
//     if(N>M)swap(N,M);

//     if(K<N){
//         return fact[N*M];
//     }

//     ll rem = N*M-K;
//     REP(i,N){

//     }


// }

void solve(){

    ll N,M;
    cin >> N >> M;
    vector<ll> A(N*M+1,0);
    REP(i,N*M){cin >> A[i+1];}
    sort(A.begin(),A.end());

    vector<ll> R(N*M+1,0);
    
    REP(K,N*M+1){
        if(K==0)continue;
        ll num = 0,num2;
        // タテ
        num2 = solve2(N,M-1,K-N);
        num += num2;

        // ヨコ
        num2 = solve2(N-1,M,K-M);
        num += num2;

        // タテヨコ
        num2 = solve2(N-1,M-1,K-N-M+1);
        num += num2;
        
        num = (num * fact[K-1])%MOD;
        num = (num * fact[N*M-K])%MOD; 
        num = (num*N*M)%MOD;
        R[K] = num;
    }

    // cout << "NM: " << N << " " << M << endl;
    // cout << "R: ";
    // REP(i,N*M){cout << R[i+1] << " ";}cout << endl;

    ll Ans = 0;
    REP(i,N*M+1){
        ll num = (R[i]%MOD)*A[i];
        num %= MOD;
        Ans = (Ans+num)%MOD;
    }
    cout << Ans << endl;

}

signed main(void){
    //cin.tie(nullptr);
    //ios::sync_with_stdio(false);


    // for(ll N=1;N<=3;N++){
    //     for(ll M=1;M<=4;M++){

    //         vector<ll> R(N*M+1,0);

    //         cout << "N=" << N << ", M=" << M << endl;
    //         vector<ll> d(N*M,0);
    //         REP(i,N*M){d[i]=i;}
    //         do{

    //             ll Ans = 1ll<<61;
    //             REP(i,N){
    //                 REP(j,M){
    //                     bool e = false;
    //                     bool c = true;
    //                     REP(k,N){
    //                         if(d[i*M+j] < d[k*M+j])c=false;
    //                     }
    //                     if(c)e=true;
    //                     c = true;
    //                     REP(l,M){
    //                         if(d[i*M+j] < d[i*M+l])c=false;
    //                     }
    //                     if(c)e=true;

    //                     if(e){
    //                         Ans = min(Ans,d[i*M+j]);
    //                     }
    //                 }
    //             }

    //             R[Ans]++;

    //         }while(next_permutation(d.begin(),d.end()));

    //         REP(i,N*M){
    //             cout << R[i] << " ";
    //         }cout << endl << endl;

    //     }
    // }

    fact[0] = 1;
    REP(i,200100){
        fact[i+1] = (fact[i] * (i+1))%MOD;
    }

    ll T;
    cin >> T;
    REP(i,T){solve();}
    
    return 0;
    
}