#include<bits/stdc++.h>
using namespace std;
#define MAXN 200005
#define ll long long
#define P 1000000007
char c[MAXN];
int T,ans,n,m,i,j,k,f[MAXN],g[MAXN],h[MAXN],a[MAXN],s[2][MAXN],u[MAXN],v[MAXN],x[MAXN];
int Pow(int x,int y)
{
    int s=1;
    for(;y;y>>=1,x=(ll)x*x%P)if(y&1)s=(ll)s*x%P;
    return s;
}
int main()
{
    scanf("%d",&T);
    for(i=f[0]=1;i<MAXN;i++)f[i]=(ll)f[i-1]*i%P;
    for(g[MAXN-1]=Pow(f[MAXN-1],P-2),i=MAXN-2;~i;i--)g[i]=(ll)g[i+1]*(i+1)%P;
    for(h[1]=1,i=2;i<MAXN;i++)h[i]=(ll)(P-P/i)*h[P%i]%P;
    while(T--)
    {
        scanf("%d",&n);
        fill(s[0],s[0]+(n<<1)+1,0);
        fill(s[1],s[1]+(n<<1)+1,0);
        for(i=0,ans=1;i<n<<1;i++)
        {
            scanf("%d",a+i);
            s[i&1][a[i]]++;
        }
        for(i=1;i<=n<<1;i++)ans=(ll)ans*f[s[0][i]+s[1][i]]%P*g[s[0][i]]%P*g[s[1][i]]%P;
        for(i=m=0;i<n<<1;i+=2)if(a[i]!=a[i|1])
        {
            u[++m]=a[i];
            v[m]=a[i|1];
        }
        for(i=x[0]=1,j=0;i<=m;i++)
        {
            if(u[i]!=v[i-1])j=(ll)x[i-1]*s[0][u[i]]%P*h[s[1][u[i]]+1]%P;
            else j=(j+(ll)x[i-2]*s[0][u[i]]%P*h[s[1][u[i]]+1])%P;
            x[i]=(x[i-1]+(ll)j*s[1][v[i]]%P*h[s[0][v[i]]+1])%P;
        }
        printf("%d\n",(ll)ans*x[m]%P);
    }
    return 0;
}