#include<bits/stdc++.h>
using namespace std;
const int M=2e5+5;
const int P=1e9+7;

int n;
int A[M],Cnt[M];

long long Fact[M],Inv[M],Finv[M];
void Init()
{
    Fact[0]=Fact[1]=Inv[0]=Inv[1]=Finv[0]=Finv[1]=1;
    for(int i=2;i<M;i++)
    {
        Fact[i]=Fact[i-1]*i%P;
        Inv[i]=(P-P/i)*Inv[P%i]%P;
        Finv[i]=Finv[i-1]*Inv[i]%P;
    }
}
long long Combination(int x,int y){return Fact[x]*Finv[y]%P*Finv[x-y]%P;}

#define HF Combination(Cnt[i],Cnt[i]/2)
#define BT Combination(Cnt[i],Cnt[i]/2-1)
long long dp[M][2];
void Solve()
{
    for(int i=1;i<=2*n;i++)Cnt[i]=0;
    for(int i=1;i<=2*n;i++)Cnt[A[i]]++;
    dp[0][0]=1;int c=0;
    for(int i=1;i<=2*n;i++)
    {
        if(Cnt[i]==0){dp[i][0]=dp[i-1][0],dp[i][1]=dp[i-1][1];continue;}
        if(Cnt[i]&1)
        {
            if(c&1)dp[i][0]=(dp[i-1][0]+dp[i-1][1])*HF%P;
            else dp[i][0]=dp[i][1]=dp[i-1][0]*HF%P;
        }
        else
        {
            if(c&1)dp[i][0]=(dp[i-1][0]*HF+dp[i-1][1]*BT)%P,dp[i][1]=(dp[i-1][0]*BT+dp[i-1][1]*HF)%P;
            else dp[i][0]=dp[i-1][0]*HF%P;
        }
        c+=Cnt[i];
        // printf("dp %d : %I64d,%I64d   %I64d,%I64d\n",i,dp[i][0],dp[i][1],HF,BT);
    }
    printf("%I64d\n",dp[2*n][0]);
}

int main()
{
    Init();
    int T;scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        for(int i=1;i<=2*n;i++)scanf("%d",&A[i]);
        Solve();
    }
    return 0;
}