#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
typedef unsigned long long ull;

const int mod = 998244353;
int ord[35], a[35], b[35], c[35];
int omg[1 << 15], rev[1 << 15];
int val[205][1 << 13], ans[1 << 13], fac[1 << 14], inv[1 << 14];
int w[35][35];
ull tmp[1 << 15];
int qpow(int x, int y)
{
	int ans = 1;
	while(y)
	{
		if(y & 1)
			ans = (ll)ans * x % mod;
		x = (ll)x * x % mod;
		y >>= 1; 
	}
	return ans;
}
void prep(int l)
{
	int i, j, w;
	for(i = 0; i <= l; i++)
	{
		for(j = 0; j < (1 << i); j++)
			rev[(1 << i) + j] = (rev[(1 << i) + (j >> 1)] >> 1) | ((j & 1) << (i - 1));
		omg[1 << i] = 1;
		w = qpow(3, (mod - 1) >> i);
		for(j = 1; j < (1 << i); j++)
			omg[(1 << i) + j] = (ll)omg[(1 << i) + j - 1] * w % mod;
	}
	fac[0] = 1;
	for(i = 1; i <= (1 << l); i++)
		fac[i] = (ll)fac[i - 1] * i % mod;
	inv[1 << l] = qpow(fac[1 << l], mod - 2);
	for(i = (1 << l); i >= 1; i--)
		inv[i - 1] = (ll)inv[i] * i % mod;
}
inline void ntt(int *A, int n, bool t) // 0 DFT 1 IDFT 
{
	int i, j, k, q, x;
	for(i = 0; i < (1 << n); i++)
		tmp[i] = A[i];
	for(i = 0; i < (1 << n); i++)
		if(rev[(1 << n) + i] > i)
			swap(tmp[rev[(1 << n) + i]], tmp[i]); 
	for(i = 0; i < n; i++)
	{
		for(j = 0; j < (1 << n); j += (2 << i))
			for(k = 0; k < (1 << i); k++)
			{
				x = tmp[j + k + (1 << i)] * omg[(2 << i) + k] % mod;
				tmp[j + k + (1 << i)] = tmp[j + k] + mod - x;
				tmp[j + k] = tmp[j + k] + x;
			}
		if(i == 11 || i == n - 1)
			for(j = 0; j < (1 << n); j++)
				tmp[j] %= mod;
	}
	for(i = 0; i < (1 << n); i++)
		A[i] = tmp[i];
	if(t)
	{
		for(i = 1; i < (1 << n); i++)
			if(i < (1 << n) - i)
				swap(A[i], A[(1 << n) - i]);
		q = qpow((1 << n), mod - 2);
		for(i = 0; i < (1 << n); i++)
			A[i] = (ll)A[i] * q % mod;
	}
}
bool cmp(int x, int y)
{
	return a[x] < a[y];
}
int det(int n)
{
	int i, j, k, r, q, ans = 1;
	for(i = 1; i <= n; i++)
	{
		for(j = i; j <= n; j++)
			if(w[j][i])
				break;
		if(j > n)
			return 0;
		if(j != i)
		{
			swap(w[i], w[j]);
			ans = mod - ans;
		}
		ans = (ll)ans * w[i][i] % mod;
		r = mod - qpow(w[i][i], mod - 2);
		for(j = i + 1; j <= n; j++)
		{
			q = (ll)r * w[j][i] % mod;
			for(k = i; k <= n; k++)
				w[j][k] = (w[j][k] + (ll)q * w[i][k]) % mod;
		}
	}
	return ans;
}

int main()
{
	prep(13);
	int n, m, i, j, k, l, d = 0, tmp;
	m = 0;
	scanf("%d", &n);
	for(i = 1; i <= n; i++)
	{
		scanf("%d", a + i);
		m = max(m, a[i]);
	}
	for(i = 1; i <= n; i++)
	{
		scanf("%d", b + i);
		m = max(m, b[i]);
	}
	for(i = 1; i <= n; i++)
		ord[i] = i;
	sort(ord + 1, ord + n + 1, cmp);
	for(i = 1; i <= n; i++)
		c[i] = a[ord[i]];
	for(i = 1; i <= n; i++)
		a[i] = c[i];
	for(i = 1; i <= n; i++)
		c[i] = b[ord[i]];
	for(i = 1; i <= n; i++)
		b[i] = c[i];
	for(i = 1; i <= n; i++)
	{
		if(b[i - 1] > b[i] || a[i] > b[i])
		{
			printf("0\n");
			return 0;
		}
		d += b[i] - a[i];
	}
	if(d & 1)
	{
		printf("0\n");
		return 0;
	}
	d /= 2;
	l = 0;
	while((n * m) >> l)
		l++;
	for(i = 0; i <= m; i++)
	{
		for(j = 0; j <= i / 2; j++)
		{
			val[i][j] = (ll)inv[j] * inv[i - 2 * j] % mod;
			if(j & 1)
				val[i][j] = mod - val[i][j];
		}
		ntt(val[i], l, 0);
	}
	for(i = 0; i < (1 << l); i++)
	{
		for(j = 1; j <= n; j++)
			for(k = 1; k <= n; k++)
				w[j][k] = ((b[k] >= a[j]) ? val[b[k] - a[j]][i] : 0);
		ans[i] = det(n);
	}
	ntt(ans, l, 1);
	tmp = 0;
	for(i = 0; i <= d; i++)
		tmp = (tmp + (ll)fac[2 * (d - i)] * ans[i] % mod * inv[d - i]) % mod;
	tmp = (ll)tmp * qpow(2, mod - 1 - d) % mod * fac[d] % mod;
	printf("%d\n", tmp);
	return 0;
}