#include<bits/stdc++.h>
#define int long long
using namespace std;
const int MAXN=55,mod=1e9+7;
int T,n,a[MAXN];
bool vis[MAXN];
unordered_map<int,int> dp[MAXN];
bool check(int mask,int val){
	for(int i=max(1ll,val*2-n);i<=min(n,val*2-1);i++){
		if(mask>>i&1){
			int d=val*2-i;
			if(!(mask>>d&1)) return 0;
		}
	}
	return 1;
}
inline int dfs(int i,int mask){
	if(i==0) return 1;
	if(dp[i].find(mask)!=dp[i].end()) return dp[i][mask];
	int ans=0;
	if(a[i]){
		if(check(mask^(1ll<<a[i]),a[i])) ans=dfs(i-1,mask^(1ll<<a[i]));
	}
	else{
		for(int v=1;v<=n;v++)
			if(!vis[v]&&(mask>>v&1)&&check(mask^(1ll<<v),v))
				ans=(ans+dfs(i-1,mask^(1ll<<v)))%mod;
	}
	return dp[i][mask]=ans;
}
signed main(){
	cin>>T;
	while(T--){
		cin>>n;
		for(int i=1;i<=n;i++)
			dp[i].clear(),vis[i]=0;
		for(int i=1;i<=n;i++)
			cin>>a[i],vis[a[i]]=1;
			
		int ans=1,cnt=0;
		for(int i=1;i<=n;i++)
			if(!vis[i]) cnt++;
		for(int i=1;i<=cnt;i++) ans=ans*i%mod;
		cout<<(ans-dfs(n,(1ll<<n+1)-1)+mod)%mod<<'\n';
	}
	return 0;
}