#include<bits/stdc++.h>
using namespace std;
const int N=5005,M=998244353;
int T,n,cnt[N],C[N][N],dp[N][N];
void ck(int &x){(x>=M)&&(x-=M);}
void slv(){
    scanf("%d",&n);
    for(int i=1;i<=n;++i)cnt[i]=0;
    for(int i=1;i<=n;++i){
        int x;
        scanf("%d",&x);
        cnt[x]++;
    }
    for(int i=0;i<=n;++i)for(int j=0;j<=n;++j)dp[i][j]=0;
    int m=0;
    dp[0][0]=1;
    for(int i=n;i;i--)if(cnt[i]){
        for(int j=m;j>=m-j;--j){
            int p=max(j-1,0)+((j<<1)==n)-(m-j);
            if(p<0)continue;
            // printf("___%d %d  %d\n",m-j,j,p);
            dp[m-j+cnt[i]][j]=(dp[m-j+cnt[i]][j]+1ll*dp[m-j][j]*C[p][cnt[i]])%M;
            dp[m-j+cnt[i]-1][j+1]=(dp[m-j+cnt[i]-1][j+1]+1ll*dp[m-j][j]*C[p][cnt[i]-1])%M;
        }
        m+=cnt[i];
    }
    // for(int i=0;i<=4;++i){
    //     for(int j=0;j<=4;++j)printf("%d ",dp[i][j]);
    //     puts("");
    // }
    printf("%d\n",dp[n/2][n/2]);
}
int main(){
    scanf("%d",&T);
    int n=5000;
    for(int i=0;i<=n;++i){
        C[i][0]=1;
        for(int j=1;j<=i;++j)ck(C[i][j]=C[i-1][j]+C[i-1][j-1]);
    }
    while(T--)slv();
}