#include<bits/stdc++.h>
#define ci const int
#define ll long long
using namespace std;
ci mod=998244353;
ll qk(ll x,int y=mod-2){
	ll ans=1;
	while(y){
		if(y&1)ans=ans*x%mod;
		x=x*x%mod,y>>=1;
	}
	return ans;
}
inline void add(int &x,ci v){
	x+=v,x-=x<mod?0:mod;
}
inline void sub(int &x,ci v){
	x-=v,x+=x<0?mod:0;
}
ci N=35,M=N*200;
int n,a[N],b[N],t;
int mat[N][N][105];
int fac[M],dfac[M],inv[M],ans[M];
int co[M],cur[M],val[M];
void upd(ci i,ci j,ci d){
	for(int x=0;x*2<=d;++x)
		mat[i][j][x]=(ll)dfac[x]*dfac[d-2*x]%mod;
}
int A[N][N];
int det(){
	int ans=1;
	for(int i=1;i<=n;++i){
		for(int j=i;j<=n;++j)
			if(A[j][i]){
				swap(A[i],A[j]);
				if(i!=j)ans=mod-ans;
				break;
			}
		if(!A[i][i])return 0;
		ci inv=qk(A[i][i]);
		ans=(ll)ans*A[i][i]%mod;
		for(int j=i;j<=n;++j)A[i][j]=(ll)A[i][j]*inv%mod;
		for(int j=i+1;j<=n;++j)
			for(int k=n;k>=i;--k)
				sub(A[j][k],(ll)A[i][k]*A[j][i]%mod);
	}
	return ans;
}
int main(){
	fac[0]=fac[1]=dfac[0]=dfac[1]=inv[1]=1;
	for(int i=2;i<M;++i)
		fac[i]=(ll)fac[i-1]*i%mod,
		inv[i]=(ll)(mod-mod/i)*inv[mod%i]%mod,
		dfac[i]=(ll)dfac[i-1]*inv[i]%mod;
	scanf("%d",&n);
	for(int i=1;i<=n;++i)scanf("%d",&a[i]);
	for(int i=1;i<=n;++i){
		scanf("%d",&b[i]),t+=b[i]-a[i];
		if(a[i]>b[i])return 0*puts("0");
	}
	for(int i=1;i<=n;++i)
		for(int j=i+1;j<=n;++j)
			if((a[i]<a[j])!=(b[i]<b[j]))
				return 0*puts("0");
	if(t&1)return 0*puts("0");
	t>>=1;
	for(int i=1;i<=n;++i)
		for(int j=1;j<=n;++j){
			if(a[i]>b[j])continue;
			upd(i,j,b[j]-a[i]);
		}
	for(int i=0;i<=t;++i){
		for(int x=1;x<=n;++x)
			for(int y=1;y<=n;++y){
				A[x][y]=0;
				for(int p=0,tmp=1;p*2<=b[y]-a[x];++p)
					add(A[x][y],(ll)tmp*mat[x][y][p]%mod),
					tmp=(ll)tmp*i%mod;
			}
		ans[i]=det();
	}
	cur[0]=1;
	for(int j=0;j<=t;++j)
		for(int k=t+1;~k;--k)
			cur[k]=((k?cur[k-1]:0)-(ll)cur[k]*j%mod+mod)%mod;
	for(int i=0;i<=t;++i){
		int tmp=ans[i];
		for(int j=0;j<=t;++j)
			if(i!=j)
				tmp=tmp*qk((i-j+mod)%mod)%mod;
		for(int j=0;j<=t+1;++j)val[j]=cur[j];
		for(int j=t+1;j;--j){
			add(co[j-1],(ll)val[j]*tmp%mod);
			add(val[j-1],(ll)val[j]*i%mod);
		}
//		for(int j=0;j<=t;++j)add(co[j],(ll)cur[j]*tmp%mod);
	}
	int ans=0;
	for(int i=0;i<=t;++i){
		ci tmp=(ll)co[i]*((ll)fac[t]*dfac[i]%mod*dfac[t-i]%mod)%mod*fac[i]%mod*fac[2*(t-i)]%mod;
		if(i&1)sub(ans,tmp);
		else add(ans,tmp);
	}
	printf("%lld",ans*qk((mod+1)/2,t)%mod);
	return 0;
}