# include <bits/stdc++.h>
using namespace std;
const int mod=998244353,g=3,_2=(mod+1)/2;
typedef long long ll;
typedef unsigned long long ull;
typedef vector<int> Poly;
ll fac[6010],inv[6010];
int n,m,a[40],b[40],ta[40],tb[40],id[40],rev[6010],w[6010],maxn;
Poly F[40][40],G;
ll C(int n,int m){return fac[n]*inv[n-m]%mod*inv[m]%mod;}
ll power(ll a,int b)
{
	ll ans=1;
	for(;b;b>>=1,a=a*a%mod)
		if(b&1) ans=ans*a%mod;
	return ans;
}
void init(int len)
{
	int l=0;maxn=1;
	while(maxn<=len) maxn<<=1,l++;
	for(int i=0;i<maxn;i++) rev[i]=(rev[i>>1]>>1)|((i&1)<<(l-1));
	for(int i=1;i<maxn;i<<=1)
	{
		ll wm=power(g,(mod-1)/(i<<1));w[i]=1;
		for(int j=1;j<i;j++) w[i+j]=w[i+j-1]*wm%mod;
	}
}
void NTT(Poly &p,bool flag)
{
	static ull a[6010];
	p.resize(maxn);
	for(int i=0;i<maxn;i++) a[i]=p[rev[i]];
	for(int i=1;i<maxn;i<<=1)
		for(int j=0;j<maxn;j+=(i<<1))
			for(int k=j;k<j+i;k++)
			{
				ull x=a[k],y=a[k+i]*w[i+k-j]%mod;
				a[k]=x+y;a[k+i]=x+mod-y;
			}
	if(flag) for(int i=0;i<maxn;i++) p[i]=a[i]%mod;
	else
	{
		int inv=power(maxn,mod-2);
		reverse(a+1,a+maxn);
		for(int i=0;i<maxn;i++) p[i]=a[i]%mod*inv%mod;
	}
}
int Det(int x)
{
	static int a[40][40];
	int ans=1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			a[i][j]=F[i][j][x];
	for(int i=1;i<=n;i++)
	{
		for(int j=i;j<=n;j++)
			if(a[j][i])
			{
				for(int k=i;k<=n;k++) swap(a[i][k],a[j][k]);
				if(i!=j) ans=mod-ans;
				break;
			}
		if(!a[i][i]) return 0;
		ans=(ll)ans*a[i][i]%mod;
		ll t=power(a[i][i],mod-2);
		for(int j=i;j<=n;j++) a[i][j]=a[i][j]*t%mod;
		for(int j=i+1;j<=n;j++)
		{
			ll t=mod-a[j][i];
			for(int k=i;k<=n;k++) a[j][k]=(a[j][k]+t*a[i][k])%mod;
		}
	}
	return ans;
}
int main()
{
	int len=0,ans=0;
	cin>>n;
	for(int i=1;i<=n;i++) scanf("%d",&a[i]),id[i]=i,len-=a[i];
	for(int i=1;i<=n;i++) scanf("%d",&b[i]),m=max(m,b[i]),len+=b[i];
	fac[0]=fac[1]=inv[0]=inv[1]=1;
	for(int i=2;i<=n*m;i++) fac[i]=fac[i-1]*i%mod;
	for(int i=2;i<=n*m;i++) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
	for(int i=2;i<=n*m;i++) inv[i]=inv[i-1]*inv[i]%mod;
	sort(id+1,id+n+1,[](int i,int j){return a[i]<a[j];});
	for(int i=1;i<=n;i++) ta[i]=a[id[i]],tb[i]=b[id[i]];
	swap(a,ta);swap(b,tb);
	for(int i=1;i<n;i++)
		if(b[i]>b[i+1]) return puts("0"),0;
	for(int i=1;i<=n;i++)
		if(a[i]>b[i]) return puts("0"),0;
	if(len&1) return puts("0"),0;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		{
			if(a[i]>b[j]) continue;
			int l=b[j]-a[i];
			F[i][j].resize(l/2+1);
			for(int k=0;k<=l/2;k++) F[i][j][k]=inv[k]*inv[l-2*k]%mod;
		}
	init(n*m/2);G.resize(maxn);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			NTT(F[i][j],1);
	for(int i=0;i<maxn;i++) G[i]=Det(i);
	NTT(G,0);
	for(int i=0;i<=len/2;i++) ans=(ans+(i&1?mod-inv[len/2-i]:inv[len/2-i])*fac[len-2*i]%mod*G[i])%mod;
	cout<<power(_2,len/2)*ans%mod*fac[len/2]%mod<<endl;
	return 0;
}