#include <bits/stdc++.h>
using namespace std;
#define rep(i,n) for(int i=0;i<n;i++)
#define per(i,n) for(int i=n-1;i>=0;i--)
#define rng(i,c,n) for(int i=c;i<n;i++)
#define fi first
#define se second
#define vec vector
#define pb push_back
#define sz(a) (int)a.size()
#define all(a) a.begin(),a.end()
typedef long long ll;
typedef vector<int> vi;
typedef pair<int,int> pii;
void print(){cout<<'\n';}
template<class h,class...t>
void print(const h&v,const t&...u){cout<<v<<' ',print(u...);}

// snuke's mod int
template <ll mod>
struct modint{
	ll x;
	modint(ll x=0):x((x%mod+mod)%mod){}
	modint operator-()const{return modint(-x);}
	modint& operator+=(const modint a){if((x+=a.x)>=mod) x-=mod; return *this;}
	modint& operator-=(const modint a){if((x+=mod-a.x)>=mod) x-=mod; return *this;}
	modint& operator*=(const modint a){(x*=a.x)%=mod; return *this;}
	modint& operator/=(const modint& a){(x*=a.inv().x)%=mod; return *this;}
	modint operator+(const modint& a)const{return modint(*this)+=a;}
	modint operator-(const modint& a)const{return modint(*this)-=a;}
	modint operator*(const modint& a)const{return modint(*this)*=a;}
	modint operator/(const modint& a)const{return modint(*this)/=a;}
	modint pow(ll n)const{
		modint res=1,x(*this);
		while(n){
			if(n&1)res*=x;
			x*=x;
			n>>=1;
		}
		return res;
	}
	modint inv()const{return pow(mod-2);}
	bool operator<(const modint& a)const{return x<a.x;}
	bool operator==(const modint& a)const{return x==a.x;}
	bool operator!=(const modint& a)const{return x!=a.x;}
};

using mint=modint<1000000007>;

const int _n=6e6;

mint p4[_n];

mint fact[_n],invfact[_n];

void init_fact(){
	fact[0]=1;
	rng(i,1,_n){
		fact[i]=fact[i-1]*mint(i);
	}
	invfact[_n-1]=fact[_n-1].inv();
	per(i,_n-1){
		invfact[i]=invfact[i+1]*mint(i+1);
	}
	p4[0]=1;
	rng(i,1,_n){
		p4[i]=p4[i-1]*mint(4);
	}
}

mint cnk(ll n,ll k){
	if(k<0 or k>n) return 0;
	return fact[n]*invfact[k]*invfact[n-k];
}

void print(vi a){
	for(auto x:a){
		cout<<x<<" ";
	}
	cout<<"\n";
}

// int total=0;

const mint i4=mint(1)/mint(4);

void slv(){
	int n;
	cin>>n;
	// total+=n;

	const int n2=n*2;

	vi a(n2);
	rep(i,n2){
		cin>>a[i];
	}

	if(*max_element(all(a))==2){
		rep(i,sz(a)){
			if(a[i]==0){
				cout<<"0\n";
				return;
			}
			a[i]=(a[i]==2?0:1);
		}
		rotate(a.begin(),a.begin()+1,a.end());
	}

	if(*min_element(all(a))==-1){
		cout<<"0\n";
		return;
	}

	const mint ip2=mint(1)/mint(2).pow(n2);

	if(*max_element(all(a))==0){
		mint ans=n2;
		ans*=mint(n2-1);
		// print(ans.x);
		ans*=fact[n2*2-2];
		ans*=ip2;
		cout<<ans.x<<"\n";
		return;
	}

	rep(i,n2){
		a.pb(a[i]);
	}

	const int si=sz(a);

	// for(auto x:a)cout<<x<<" ";
	// print();

	vi ps(si);
	vi ps_bad(si);
	rep(i,si){
		if(i) ps_bad[i]+=ps_bad[i-1];
		if(i%2){ // me
			if(i+1<si and a[i]==1 and a[i+1]==0){
				ps_bad[i]++;
			}
		}else{ // other
			if(i+1<si and a[i]==0 and a[i+1]==1){
				ps_bad[i]++;
			}
		}
		if(i) ps[i]+=ps[i-1];
		ps[i]+=(i+1<si and a[i]==0 and a[i+1]==1);
	}
	vec<vi> nxt(si+1,vi(2,si));
	per(i,si){
		nxt[i][0]=nxt[i+1][0];
		nxt[i][1]=nxt[i+1][1];
		if(a[i]==0){
			nxt[i][0]=i;
		}else{
			nxt[i][1]=i;
		}
	}

	// rep(i,si) cout<<ps[i]<<" ";
	// print();

	auto get=[&](int l,int r)->mint{
		if(a[r]==0){
			return 0;
		}
		int cnt_bad=ps_bad[r-1]-(!l?0:ps_bad[l-1]);
		if(cnt_bad>0){
			return 0;
		}
		
		mint ans=0;
		int groups=ps[r-1]-(!l?0:ps[l-1]);
		int one=groups;
		int k=n2*2;
		k-=2*groups;
		mint wys=p4[groups];
		
		if(a[l]==0){
			one--;
			groups--;
			k+=2;
			wys*=i4;
			int i=nxt[l][1];
			i--;
			i-=l;
			// print(i,one,k,wys.x);
			one++;
			k-=2;
			wys*=4;
			mint old_wys=wys;
			int old_k=k;
			int len=n2-i;
			// len/=2;
			len-=one;
			wys*=mint(len);
			k-=1;
			wys*=fact[k];
			// print(wys.x);
			ans+=wys;

			wys=old_wys;
			k=old_k;
			len=i;
			wys*=mint(len);
			k--;
			len=n2;
			len-=one+1;
			// print("e",wys.x,len,i);
			wys*=mint(len);
			k--;
			wys*=fact[k];
			// print("e",wys.x,k,len,i);
			// print(wys.x);
			ans+=wys;
		}else{
			int j=nxt[l][0];
			j=min(j,r+1);
			j--;
			// auto it=zeros.lower_bound(l);
			// if(it==zeros.end()){
			// 	j=r;
			// }else{
			// 	j=(*it)-1;
			// }
			j-=l;
			mint old_wys=wys;
			int old_k=k;
			wys*=2;
			k-=1;
			int len=n2-j-1-one;
			wys*=mint(len);
			k-=1;
			wys*=fact[k];
			ans+=wys;
			
			k=old_k;
			wys=old_wys*mint(2);
			k-=1;
			len=j+1;
			wys*=mint(len);
			k-=1;
			one++;
			
			len=n2-one;
			wys*=mint(len);
			k-=1;
			wys*=fact[k];
			ans+=wys;
		}
		return ans;
	};

	mint ans=0;
	for(int i=1;i<n2;i+=2){
		ans+=get(i,i+n2-1);
	}
	if(count(all(a),1)==sz(a)){
		ans*=2;
	}
	ans*=ip2;
	cout<<ans.x<<"\n";
}

signed main(){
	ios::sync_with_stdio(0),cin.tie(0);
	init_fact();

	int t;
	t=1;
	cin>>t;

	rep(cs,t){
		slv();
	}
}