#include <cstdio>
#include <iostream>
#include <algorithm>
#include <vector>
#include <set>
#include <map>
#include <string>
#include <cstring>
#include <stack>
#include <queue>
#include <cmath>
#include <ctime>
#include <bitset>
#include <utility>
#include <assert.h>
using namespace std;
#define rank rankk
#define mp make_pair
#define pb push_back
#define xo(a,b) ((b)&1?(a):0)
#define tm tmp

//#define LL ll
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef long long ll;
typedef pair<ll,int> pli;
typedef pair<ll,ll> pll;
const int INF=0x3f3f3f3f;
const ll INFF=0x3f3f3f3f3f3f3f3fll;
const int MAX=405;
//const ll MAXN=2e8;
//const int MAX_N=MAX;
const int MOD=1e9+7;
//const long double pi=acos(-1.0);
//const double eps=0.00000001;
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
template<typename T>inline T abs(T a) {return a>0?a:-a;}
template<class T> inline
void read(T& num) {
    bool start=false,neg=false;
    char c;
    num=0;
    while((c=getchar())!=EOF) {
        if(c=='-') start=neg=true;
        else if(c>='0' && c<='9') {
            start=true;
            num=num*10+c-'0';
        } else if(start) break;
    }
    if(neg) num=-num;
}
inline ll powMM(ll a,ll b,ll M){
    ll ret=1;
    a%=M;
//    b%=M;
    while (b){
        if (b&1) ret=ret*a%M;
        b>>=1;
        a=a*a%M;
    }
    return ret;
}
void open()
{
//    freopen("1009.in","r",stdin);
    freopen("out.txt","w",stdout);
}
inline void add(int &x,int y)//保证x、y均已经取过模
{
    x=x+y;if(x>=MOD)x-=MOD;
}
int n,n0;
int a[MAX],b[MAX],pa[MAX],pb[MAX];
int g[MAX];
int dp[405][405][135],cnt[405][405];
bitset<405>st;
//lin;
int main()
{
    scanf("%d",&n);n0=n/3;
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        pa[a[i]]=i;
    }
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&b[i]);
        pb[b[i]]=i;
    }
    g[0]=1;
    for(int i=1;i<=n0;i++)
    {
        g[i]=(ll)(3*(i-1)+1)*(3*(i-1)+2)%MOD*g[i-1]%MOD;
//        cout<<g[i]<<endl;
    }
//    lin.reset();
//    for(int i=1;i<=n;i++)
//    {
//        st.reset();
//        lin.set(a[i]);
//        st=lin;
//        for(int j=1;j<=n;j++)
//        {
//            st.set(b[j]);
//            cnt[i][j]=st.count();
//        }
//    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        {
            st.reset();
            for(int p=1;p<=i;p++)
//                st[a[p]]=1;
                st.set(a[p]);
            for(int p=1;p<=j;p++)
//                st[b[p]]=1;
                st.set(b[p]);
            cnt[i][j]=st.count();
//            cout<<i<<" "<<j<<" "<<cnt[i][j]<<endl;
        }
    dp[0][0][0]=1;
    for(int i=0;i<=n;i++)
        for(int j=0;j<=n;j++)
            for(int s=0;s<=n0;s++)
            {
                if(!i&&!j&&!s)continue;
                if(i)add(dp[i][j][s],dp[i-1][j][s]);
                if(j)add(dp[i][j][s],dp[i][j-1][s]);
                if(i&&j)add(dp[i][j][s],MOD-dp[i-1][j-1][s]);
                if(pb[a[i]]<=j||pa[b[j]]<=i)continue;
                if(i&&j&&s)
                {
                    int ge=n-3*(n0-s)-cnt[i][j];
                    if(ge>=1)
                        add(dp[i][j][s],(ll)dp[i-1][j-1][s-1]*ge%MOD);
                }
            }
    printf("%lld\n",(ll)dp[n][n][n0]*g[n0]%MOD);
    return 0;
}


Measured memory_bytes: 86272