#include<bits/stdc++.h>
using namespace std;
const int p=998244353;
int bb[210]={};
int qow(int a,int n){
	int ans=1;
	for(;n;n>>=1){
		if(n&1) ans=((long long)ans*a)%p;
		a=((long long)a*a)%p; 
	}
	return ans;
}
inline int inverse(int a){
	return qow(a,p-2);
}
int aa[50][50];
int det(int m){
	int ans=1;
	for(int i=1;i<=m;i++){
		for(int j=i+1;j<=m;j++){
			if(aa[j][i]){
				for(int k=1;k<=m;k++) swap(aa[i][k],aa[j][k]);
                ans=p-ans;
				break;
			}
		}
		int u=inverse(aa[i][i]);
		for(int j=i+1;j<=m;j++){
			int uu=((long long)u*aa[j][i])%p;
			for(int k=1;k<=m;k++) aa[j][k]=(aa[j][k]+(long long)(p-uu)*aa[i][k])%p; 
		}
	}
	for(int i=1;i<=m;i++) ans=((long long)ans*aa[i][i])%p;
    return ans;
}
int fac[10000];
int infac[10000];
int anss[10000]={};
int t[10000];
int u[10000]={}; 
int main(){
	fac[0]=1;
	for(int i=1;i<10000;i++) fac[i]=((long long)fac[i-1]*i)%p;
	for(int i=0;i<10000;i++) infac[i]=inverse(fac[i]);
	int n,a[50],b[50];
	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]);
	for(int i=1;i<=n;i++) bb[a[i]]=b[i];
	for(int i=1,tot=0;i<=200;i++){
		if(bb[i]){
			a[++tot]=i;
			b[tot]=bb[i];
		}
	}
	for(int i=2;i<=n;i++) if(b[i]<=b[i-1]) {
		printf("0");
		return 0;
	}
	int s=0;
	for(int i=1;i<=n;i++) s+=b[i]-a[i];
	if((s<0)||(s&1)){
		printf("0");
		return 0;
	} 
	s/=2;
	u[0]=1;
	for(int i=0;i<=s;i++){
		u[i+1]=1;
		for(int j=i;j>=1;j--) u[j]=(u[j-1]+(long long)(p-i)*u[j])%p;
		u[0]=((long long)(p-i)*u[0])%p;
	}
	for(int i=0;i<=s;i++){
		for(int j=1;j<=n;j++){
			for(int k=1;k<=n;k++){
				aa[j][k]=0;
				if(b[k]>=a[j]) for(int o=(b[k]-a[j])/2;o>=0;o--) aa[j][k]=((long long)i*aa[j][k]+(long long)infac[o]*infac[b[k]-a[j]-2*o])%p;
                //if(i==0) cout<<aa[j][k]<<endl;
			}
		}
		int w=det(n);
        //if(i==1) cout<<(w*144ll)%p<<endl;
		w=((((long long)infac[i]*infac[s-i])%p)*w)%p;
		if((s-i)&1) w=(p-w)%p;
		t[s]=1;
		for(int j=s-1;j>=0;j--) t[j]=(u[j+1]+(long long)i*t[j+1])%p;
		for(int j=0;j<=s;j++) anss[j]=(anss[j]+(long long)t[j]*w)%p;
	}
	int ans=0;
	for(int j=0,jj=1;j<=s;j++,jj=p-jj) ans=(ans+(((((long long)infac[s-j]*fac[2*s-2*j])%p)*jj)%p)*anss[j])%p;
	ans=((long long)ans*qow(2,p-1-s))%p;
    printf("%lld",((long long)ans*fac[s])%p);
	return 0;
}