#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,j=x;i>=0&&j<n;i--,j++)
		if((s>>i&1)^(s>>j&1))return 0;
	return 1;
}
int main(){
	cin>>T;
	while(T--){
		cin>>n;f.clear();f[p=0]=ans=1;
		for(ll i=0;i<n;i++){
			cin>>a[i];a[i]--;
			if(!~a[i])ans=(ans*++p)%mod;
		}
		for(auto i:f){
			ll x=i.first,y=i.second,z=__builtin_popcountll(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;
}