#include<bits/stdc++.h>
using namespace std;
const int Max_N=30,Mv=2e2,M=6e3;
const long long Mod=998244353;
int n,m,a[Max_N+5],b[Max_N+5],c[Max_N+5];
long long p[M+5],Inv[M+5],A[Max_N+5][Max_N+5],Ml[Mv+5],ans[M+5],f[M+5],g[M+5],res[M+5];
long long Mul(int x,int y){
	long long sum=1;
	for(;y;y>>=1,x=1ll*x*x%Mod)
	 if(y&1) sum=sum*x%Mod;
	return sum;
}
int main(){
	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++)
	 for(int j=i+1;j<=n;j++)
	  if((a[i]<=a[j])!=(b[i]<=b[j])){
	  	puts("0"); return 0;
	  }
	for(int i=1;i<=n;i++) c[i]=b[i]-a[i],m+=c[i];
	if(m&1){ puts("0"); return 0; }
	p[0]=Inv[0]=Inv[1]=1;
	for(int i=1;i<=m;i++) p[i]=p[i-1]*i%Mod;
	for(int i=2;i<=m;i++) Inv[i]=Inv[Mod%i]*(Mod-Mod/i)%Mod;
	for(int i=1;i<=m;i++) Inv[i]=Inv[i]*Inv[i-1]%Mod;
	for(int e=0;e<=m;e++){
		memset(A,0,sizeof(A));
		Ml[0]=1; for(int i=1;i<=Mv;i++) Ml[i]=Ml[i-1]*e%Mod;
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){
				int C=b[j]-a[i];
				for(int D=0;D<=C/2;D++){
					int E=C-2*D;	
					A[i][j]+=((D&1)?-1ll:1ll)*Inv[D]*Inv[E]%Mod*Ml[E]%Mod;
					A[i][j]=(A[i][j]+Mod)%Mod;
				}
			}
		}
		for(int i=1;i<=n;i++){
			for(int j=i+1;j<=n&&A[i][i]==0;j++)
			 if(A[j][i])
			  for(int k=1;k<=n;k++)
			   A[i][k]=(A[i][k]-A[j][k]+Mod)%Mod;
			if(!A[i][i]) break;
			long long INV=Mul(A[i][i],Mod-2);
			for(int j=i+1;j<=n;j++){
				long long d=A[j][i]*INV%Mod;
				for(int k=1;k<=n;k++)
				 A[j][k]-=A[i][k]*d%Mod,
				 A[j][k]=(A[j][k]+Mod)%Mod;
			}
		}
		ans[e]=1; for(int i=1;i<=n;i++) ans[e]=ans[e]*A[i][i]%Mod;
	}
	f[0]=1;
	for(int i=0;i<=m;i++){
		for(int j=0;j<=i;j++) res[j+1]=f[j];
		for(int j=0;j<=i;j++) res[j]-=1ll*f[j]*i%Mod,res[j]=(res[j]+Mod)%Mod;
		for(int j=0;j<=i+1;j++) f[j]=res[j],res[j]=0;
	}
	for(int i=0;i<=m;i++){
		for(int j=0;j<=m;j++)
		 if(i^j) ans[i]=ans[i]*Mul((i-j+Mod)%Mod,Mod-2)%Mod;
		for(int j=0;j<=m+1;j++) res[j]=f[j];
		for(int j=m+1;j>=1;j--)
		 res[j-1]=(res[j-1]+1ll*i*res[j])%Mod;
		for(int j=0;j<=m;j++) res[j]=res[j+1];
		for(int j=0;j<=m;j++) g[j]+=ans[i]*res[j]%Mod,g[j]%=Mod; 
	}
	long long Ans=0;
	for(int i=0;i<=m;i++)
	 Ans=(Ans+g[i]*p[i]%Mod*Inv[i>>1]%Mod)%Mod;
	Ans=Ans*p[m/2]%Mod*Mul(Mul(2,Mod-2),m/2)%Mod;
	printf("%lld\n",Ans);
	return 0;
}