#include<bits/stdc++.h>
#define pii std::pair<int,int>
#define mkp std::make_pair
#define fir first
#define sec second
typedef long long ll;
inline void rd(){}
template<typename T,typename ...U>
inline void rd(T &x,U &...args){
	int ch=getchar();
	T f=1;x=0;
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
	x*=f;rd(args...);
}
const int N=35,M=6005,mod=998244353,V=6000;
int n,a[N],b[N],d[N],t;
int fac[M],INV[M],inv[M],Inv[M][205];
inline int KSM(int x,int n){
	int ans=1;
	while(n){
		if(n&1)ans=1ll*ans*x%mod;
		x=1ll*x*x%mod;
		n>>=1;
	}return ans;
}
inline void Init(){
	fac[0]=inv[0]=INV[0]=1;
	for(int i=1;i<=V;i++)fac[i]=1ll*fac[i-1]*i%mod;
	inv[V]=KSM(fac[V],mod-2);
	for(int i=V-1;i>=1;i--)inv[i]=1ll*inv[i+1]*(i+1)%mod;
	for(int i=1;i<=V;i++)INV[i]=1ll*inv[i]*fac[i-1]%mod;
	
}
int Y[M],A[M],B[M],mp[N][N];
inline int Det(){
	int f=1;
	for(int i=1;i<=n;i++){
		int cur=i;
		while(!mp[cur][i]&&cur<n)++cur;
		if(!mp[cur][i])return 0;
		if(cur!=i)std::swap(mp[cur],mp[i]),f=mod-f;
		for(int j=i+1;j<=n;j++){
			int div=1ll*mp[j][i]*KSM(mp[i][i],mod-2)%mod;
			for(int k=i;k<=n;k++){
				mp[j][k]=(mp[j][k]+1ll*mp[i][k]*(mod-div)%mod)%mod;
			}
		}
	}
	int ans=1;
	for(int i=1;i<=n;i++)ans=1ll*ans*mp[i][i]%mod;
	return 1ll*ans*f%mod;
}
inline int Calc(int x){
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if(b[j]<a[i]){
				mp[i][j]=0;
				continue;
			}
			mp[i][j]=Inv[x][b[j]-a[i]];
		}
	}
	return Det();
}
int map[M];
signed main(){
	rd(n);
	for(int i=1;i<=n;i++)rd(a[i]);
	for(int i=1;i<=n;i++)rd(b[i]),map[a[i]]=b[i];
	for(int i=1;i<=n;i++){
		d[i]=b[i]-a[i];
		t+=d[i];
		if(d[i]<0)return puts("0"),0;
	}
	if(t&1)return puts("0"),0;
	t/=2;
	Init();
	for(int x=1;x<=t+1;x++){
		for(int r=0;r<=200;r++){
			int now=0,v=1;
			for(int k=0;k<=r/2;k++){
				now=(now+1ll*v*inv[k]%mod*inv[r-2*k]%mod)%mod;
				v=1ll*v*(mod-x)%mod;
			}
			Inv[x][r]=now;
		}
	}
	for(int i=1;i<=t+1;i++){
		Y[i]=Calc(i);
		for(int j=1;j<i;j++)Y[i]=1ll*Y[i]*INV[i-j]%mod;
		for(int j=i+1;j<=t+1;j++)Y[i]=1ll*Y[i]*(mod-INV[j-i])%mod;
	}
	B[0]=1;
	for(int i=1;i<=t+1;i++){
		for(int j=i-1;j>=0;j--){
			A[j]=(1ll*B[j]*Y[i]+1ll*A[j]*(mod-i))%mod;
			if(j)(A[j]+=A[j-1])%=mod;
		}
		for(int j=i;j>=0;j--){
			B[j]=1ll*B[j]*(mod-i)%mod;
			if(j)(B[j]+=B[j-1])%=mod;
		}
	}
	int ans=0;
	for(int i=0;i<=t;i++){
		ans=(ans+1ll*fac[2*t-2*i]%mod*inv[t-i]%mod*A[i]%mod)%mod;
	}ans=1ll*ans*fac[t]%mod*KSM(KSM(2,t),mod-2)%mod;
	std::sort(a+1,a+n+1);
	std::sort(b+1,b+n+1);
	for(int i=1;i<=n;i++)if(map[a[i]]!=b[i])ans=0;
	printf("%d\n",ans);
	return 0;
}