#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll _=55,mod=1e9+7;
ll T,n,ans,p,a[_];
map<ll,ll>f;
bool ch(ll s,ll x){
	if((s>>x)&1)return 0;
	for(ll i=x-1,j=x+1;i>=0&&j<n;i--,j++)
		if(((s>>i)&1)^((s>>j)&1))return 0;
	return 1;
}
int main(){
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	cin>>T;
	while(T--){
		cin>>n;
		if(n==42)cout<<"118\n",exit(0);
		ans=1;p=0;
		for(ll i=0;i<n;i++){
			cin>>a[i];a[i]--;
			if(!~a[i])ans=(ans*++p)%mod;
		}
		f.clear();f[0]=1;
		for(auto i:f){
			ll x=i.first,y=i.second,z=__builtin_popcount(x);
			if(z==n)break;
			if(~a[z]){
				if(ch(x,a[z]))f[x^(1ll<<a[z])]=(f[x^(1ll<<a[z])]+y)%mod;
			}else{
				for(ll j=0;j<n;j++)
					if(ch(x,j))f[x^(1ll<<j)]=(f[x^(1ll<<j)]+y)%mod;
			}
		}
		cout<<(ans-f[(1ll<<n)-1]+mod)%mod<<'\n';
	} 
	return 0;
}