# include <bits/stdc++.h>
using namespace std;
using ll=long long;
using Poly=vector<int>;
constexpr int mod=988244353;
int a[110],b[110],p[110],c1[110],c2[110],fac[110],inv[110],f[110],g[110][110],y[110];
template<typename T>
inline void add(int &x,const T &y){x=(x+y)%mod;}
Poly operator+(const Poly &F,const Poly &G)
{
	int n=F.size(),m=G.size();
	Poly H(max(n,m));
	copy(F.begin(),F.end(),H.begin());
	for(int i=0;i<m;i++) add(H[i],G[i]);
	return H;
}
Poly operator*(const Poly &F,const int &x)
{
	Poly G=F;
	for(int &i:G) i=(ll)i*x%mod;
	return G;
}
Poly operator*(const Poly &F,const Poly &G)
{
	int n=F.size(),m=G.size();
	Poly H(n+m-1);
	for(int i=0;i<n;i++)
		for(int j=0;j<m;j++)
			add(H[i+j],(ll)F[i]*G[j]);
	return H;
}
Poly div(Poly F,const int &a)
{
	int n=F.size();
	Poly G(n-1);
	for(int i=n-1;i;i--) G[i-1]=F[i],add(F[i-1],(ll)a*F[i]);
	return G;
}
Poly Langrange(int n)
{
	Poly mul={1},F;
	for(int i=0;i<=n;i++) mul=mul*Poly{mod-i,1};
	for(int i=0;i<=n;i++)
	{
		Poly G=div(mul,i)*y[i];
		for(int j=0;j<i;j++) G=G*inv[i-j];
		for(int j=n;j>i;j--) G=G*(mod-inv[j-i]);
		F=F+G;
	}
	return F;
}
int main()
{
	int n,ans=0;
	cin>>n;
	fac[0]=fac[1]=inv[1]=1;
	for(int i=2;i<=n;i++) fac[i]=(ll)fac[i-1]*i%mod;
	for(int i=2;i<=n;i++) inv[i]=(ll)(mod-mod/i)*inv[mod%i]%mod;
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	for(int i=1;i<=n;i++) scanf("%d",&b[i]);
	for(int i=1;i<=n;i++) p[a[i]]=b[i];
	p[n+1]=n+1;
	fill(c1+1,c1+n+1,1);
	fill(c2+1,c2+n+1,1);
	for(int i=1;i<=n;i++)
		if(p[i]) c1[i]=c2[p[i]]=0;
	partial_sum(c1+1,c1+n+2,c1+1);
	partial_sum(c2+1,c2+n+2,c2+1);
	for(int x=0;x<=n;x++)
	{
		for(int i=0;i<=n;i++) g[i][0]=g[0][i]=1;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				g[i][j]=(g[i-1][j]+g[i][j-1]+(ll)(mod+x-1)*g[i-1][j-1])%mod;
		f[0]=1;
		for(int i=1;i<=n+1;i++)
		{
			if(!p[i]) continue;
			for(int j=0;j<i;j++)
				if(p[j]<p[i]) add(f[i],(ll)f[j]*g[c1[i]-c1[j]][c2[p[i]]-c2[p[j]]]);
		}
		y[x]=f[n+1];
		memset(f,0,sizeof(int)*(n+2));
	}
	Poly F=Langrange(n);
	for(int i=0,sz=F.size();i<sz && i<=c1[n];i++) add(ans,(ll)F[i]*fac[c1[n]-i]);
	cout<<ans<<endl;
	return 0;
}