#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);}

const int N=5050;
int dp[N][N][2],cnt[N],F[N],I[N];
int binom(int n,int k){return mul(F[n],mul(I[k],I[n-k]));}
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;
		scanf("%i",&n);
		int mn=n,mx=1;
		for(int i=1;i<=n;i++){
			int x;
			scanf("%i",&x);
			cnt[x]++;
			mn=min(mn,x);
			mx=max(mx,x);
		}

		bool ok=true;
		//if(mn==mx)ok=false;
		cnt[mx]--;
		if(cnt[mx]>0)ok=false;
		dp[0][0][0]=1;

		for(int i=1;i<=n;i++){
			for(int j=0;j<=n;j++){
				for(int t=0;t<2;t++){
					if(j>=cnt[i])ckadd(dp[i][j][t],mul(dp[i-1][j-cnt[i]][t],binom(j,cnt[i])));
				}
			}
			if(cnt[i]>=1){
				int add=cnt[i]-1;
				for(int j=0;j<=n;j++){
					if(j-add+1>=0)ckadd(dp[i][j][1],mul(dp[i-1][j-add+1][1],binom(j,add)));
					if(j-add+2>=0)ckadd(dp[i][j][1],mul(dp[i-1][j-add+2][0],binom(j,add)));
				}
			}
			//for(int j=0;j<=n;j++){
			//	printf("(%i %i) ",dp[i][j][0],dp[i][j][1]);
			//}
			//printf("\n");
		}

		if(!ok)printf("0\n");
		else{
			if(n==2)printf("1\n");
			else printf("%i\n",dp[n][0][1]);
		}

		for(int i=0;i<=n;i++){
			cnt[i]=0;
			for(int j=0;j<=n;j++){
				for(int t=0;t<2;t++){
					dp[i][j][t]=0;
				}
			}
		}
	}
	return 0;
}