#include<bits/stdc++.h>
#define int long long
#define mod 998244353
using namespace std;
int quickpow(int a,int b){
	if(a<0) a+=mod;
	int ret=1;
	while(b){
		if(b&1) ret=ret*a%mod;
		a=a*a%mod;
		b>>=1;
	}return ret;
} 
int fac[10001],inv[10001];
void init(int n=10000){
	fac[0]=1;for(int i=1;i<=n;i++) fac[i]=fac[i-1]*i%mod;
	inv[n]=quickpow(fac[n],mod-2);
	for(int i=n-1;~i;i--) inv[i]=inv[i+1]*(i+1)%mod;
}
int C(int n,int m){
	return fac[n]*inv[m]%mod*inv[n-m]%mod;
}
int y[10001];
#define poly vector<int>
poly operator +(poly A,poly B){
	poly ret;int len=max(A.size(),B.size());
	A.resize(len),B.resize(len),ret.resize(len);
	for(int i=0;i<len;i++){
		ret[i]=A[i]+B[i];
		if(ret[i]>=mod) ret[i]-=mod;
	}
	return ret;
}
poly operator *(poly A,int b){
	for(int i=0;i<(int)A.size();i++) A[i]=A[i]*b%mod;
	return A;
}
poly add(poly A,int b){
	if(b<0) b+=mod;
	A.resize(A.size()+1);
	for(int i=(int)A.size()-1;~i;i--){
		if(i) A[i]=A[i-1]+A[i]*b%mod;
		else A[i]=A[i]*b%mod;
		if(A[i]>=mod) A[i]-=mod;
	}
	return A;
}
poly del(poly A,int b){
	if(b<0) b+=mod;
	poly ret;ret.resize(A.size()-1);
	for(int i=(int)ret.size()-1;~i;i--){
		ret[i]=A[i+1];
		A[i]-=ret[i]*b%mod;
		if(A[i]<0) A[i]+=mod;
	}
	return ret;
}
poly get(int n){
	poly ans;
	poly V={1};for(int i=0;i<=n;i++) V=add(V,-i);
	for(int i=0;i<=n;i++){
		int S=y[i];
		for(int j=0;j<=n;j++){
			if(i==j) continue;
			S=S*quickpow(i-j,mod-2)%mod;
		}
		if(S<0) S+=mod;
		ans=ans+(del(V,-i)*S);
	}
	return ans;
} 
int po(poly A,int x){
	int ret=0,X=1;
	for(int i=0;i<(int)A.size();i++){
		ret+=A[i]*X%mod;
		if(ret>=mod) ret-=mod;
		X=X*x%mod;
	}return ret;
}
int n;
int c[31][31];
int det(){
	int ret=1;
	for(int i=1;i<=n;i++){
		if(!c[i][i]){
			for(int j=i+1;j<=n;j++){
				if(c[j][i]){
					swap(c[i],c[j]);
					ret=mod-ret;
				}
			}
		}
		for(int j=i+1;j<=n;j++){
			int t=c[j][i]*quickpow(c[i][i],mod-2)%mod;
			for(int k=i;k<=n;k++){
				c[j][k]-=c[i][k]*t%mod;
				if(c[j][k]<0) c[j][k]+=mod;
			}
		}
	}
	for(int i=1;i<=n;i++) ret=ret*c[i][i]%mod;
	return ret;
}
int a[10001],b[10001];
poly G[201];
mt19937 rd(time(NULL));
signed main(){
	init();
	cin>>n;
	for(int i=1;i<=n;i++){
		cin>>a[i];
	}
	int t=0;
	for(int i=1;i<=n;i++){
		cin>>b[i];
		if(b[i]<a[i]){
			cout<<0<<'\n';
			return 0;
		}
		t+=b[i]-a[i];
	}
	for(int i=1;i<=n;i++){
		for(int j=i+1;j<=n;j++){
			if(a[i]<=a[j]&&b[j]<=b[i]){
				cout<<0;
				return 0;
			}
		}
	}
	if(t&1){
		cout<<0<<'\n';
		return 0;
	}
	t>>=1;
	int maxn=0;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			maxn=max(maxn,b[j]-a[i]);
		}
	}
	for(int d=0;d<=maxn;d++){
		for(int i=0;i*2<=d;i++){
			G[d].push_back(inv[i]*inv[d-2*i]%mod);
		}
	}
	for(int x=0;x<=200*n;x++){
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				if(b[j]-a[i]<0){
					c[i][j]=0;continue;
				}
				c[i][j]=po(G[b[j]-a[i]],x);
			}
		}
		int Det=det();
		y[x]=Det;
	}
	poly A=get(200*n);
	int ans=0;
	for(int i=0;i<=t;i++){
		if(i&1){
			ans-=C(t,i)*fac[i]%mod*fac[2*t-2*i]%mod*A[i]%mod;
		}
		else{
			ans+=C(t,i)*fac[i]%mod*fac[2*t-2*i]%mod*A[i]%mod;
		}
	}
	ans=(ans%mod+mod)%mod;
	ans=ans*quickpow(quickpow(2,mod-2),t)%mod;
	cout<<ans;
}
/*
3
1 2 3
3 4 5
*/