#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define int long long
const int Mod=998244353;
const int Maxn=5010;
int qpow(int a,int b)
{
	int res=1;
	while(b)
	{
		if(b&1)
		{
			res=res*a%Mod;
		}
		a=a*a%Mod;
		b>>=1;
	}
	return res;
}
int a[Maxn],b[Maxn];
int n;
int C[Maxn][Maxn];

int work()
{
	int j,k,r=1,fh=0,l;
	for(int i=1;i<=n;i++)
	{
		k=i;
		for(j=i;j<=n;j++)
		{
			if(C[j][i])
			{
				k=j;
				break;
			}
		}
		if(C[k][i]==0)
		{
			return 0;
		}
		for(++j;j<=n;j++)
		{
			if(C[j][i] && C[j][i]<C[k][i])
			{
				k=j;
			}
		}
		if(i!=k)
		{
			swap(C[k],C[i]);
			fh^=1;
		}
		for(j=i+1;j<=n;j++)
		{
			if(C[j][i]>C[i][i])
			{
				swap(C[i],C[j]);
				fh^=1;
			}
			while(C[j][i])
			{
				l=C[i][i]/C[j][i];
				for(k=i;k<=n;k++)
				{
					C[i][k]=(C[i][k]+(ll)(Mod-l)*C[j][k])%Mod;
				}
				swap(C[j],C[i]);
				fh^=1;
			}
		}
		r=(ll)r*C[i][i]%Mod;
	}
	if(fh) return (Mod-r+Mod)%Mod;
	else return (r+Mod)%Mod;
}
void solve()
{
	scanf("%lld",&n);
	for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
	for(int i=1;i<=n;i++) scanf("%lld",&b[i]);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			C[i][j]=a[i]^b[j];
	if(n>61) puts("0");
	else printf("%lld\n",work());
}
signed main()
{
	int T;
	scanf("%lld",&T);
	while(T--)
	{
		solve();
	}
	return 0;
}