//#pragma GCC optimize("O3")
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define rep(i,n) for (ll i=0;i<(ll)n;i++)
#define rrep(i,n) for (ll i=(n)-1;i>=(ll)0;i--)
#define loop(i,m,n) for(ll i=m;i<=(ll)n;i++)
#define rloop(i,m,n) for(ll i=m;i>=(ll)n;i--)
#define vl vector<ll>
#define vvl vector<vl>
#define vvvl vector<vvl>
#define vdbg(a) rep(ii,a.size()){cout<<a[ii]<<" ";}cout<<endl;
#define vpdbg(a) rep(ii,a.size()){cout<<"{"<<a[ii].first<<","<<a[ii].second<<"} ";}cout<<endl;
#define vvdbg(a) rep(ii,a.size()){rep(jj,a[ii].size()){cout<<a[ii][jj]<<" ";}cout<<endl;}
#define setdbg(a) for(const auto & ii:a){cout<<ii<<" ";}cout<<endl;
#define inf 4000000000000000000LL
#define mod 998244353LL
//#define mod 1000000007LL
#define eps 0.000000001
#define circlepi 3.14159265358979323846
random_device rnd;// 非決定的な乱数生成器
mt19937 mt(rnd());// メルセンヌ・ツイスタの32ビット版、引数は初期シード

//#include<boost/multiprecision/cpp_int.hpp>
//#define bbi boost::multiprecision::cpp_int
//#include<atcoder/lazysegtree>


//整数同士の累乗の計算をする。
ll power(ll A, ll B) {
	ll result = 1;
	for (ll i=0;i<B;i++){
		result *= A;
	}
	return result;
}

// nのk乗をmodで割った余りを計算
ll power_mod(ll n, ll k){
	n%=mod;
	ll ans = 1;
	while (k > 0){
		if ((k&1) ==1)ans=(ans*n)%mod;
		n=n*n%mod;
		k >>= 1;
	}
	return ans;
}

//受け取った2次元文字の外側に、文字pをコーティングする。
vector<string> pad(vector<string> &s,char p){
	ll h=s.size();
	ll w=s[0].size();
	vector<string> res(h+2,string(w+2,p));
	rep(i,h)rep(j,w)res[i+1][j+1]=s[i][j];
	return res;
}

// Union-Find
struct UnionFind {
	vector<int> par, siz;
	UnionFind(int n) : par(n, -1) , siz(n, 1) { }
	// 根を求める
	int root(int x) {
		if (par[x] == -1) return x;
		else return par[x] = root(par[x]);
	}
	// x と y が同じグループに属するかどうか (根が一致するかどうか)
	bool issame(int x, int y) {
		return root(x) == root(y);
	}
	// x を含むグループと y を含むグループとを併合する
	bool unite(int x, int y) {
		x = root(x), y = root(y);
		if (x == y) return false; 
		if (siz[x] < siz[y]) swap(x, y);
		par[y] = x;
		siz[x] += siz[y];
		return true;
	}
	// x を含むグループのサイズ
	int size(int x) {
		return siz[root(x)];
	}
};


//グリッド問題等用
vl dx={1,0,-1,0};
vl dy={0,1,0,-1};

void solve(){
	ll n;
	cin>>n;
	vl x(n),y(n),z(n+n-1);
	//横ANDがX,縦ORがY,左下→右上向き斜めXORがZ
	rep(i,n)cin>>x[i];
	rep(i,n)cin>>y[i];
	rep(i,n+n-1)cin>>z[i];

	vvl ans(n,vl(n,0));

	rep(b,30){
		vl xx=x,yy=y,zz=z;
		vl xcnt=x,ycnt=y,zcnt=z;
		rep(i,n){
			xx[i]=(xx[i]>>b)%2;
			yy[i]=(yy[i]>>b)%2;
			xcnt[i]=n;
			ycnt[i]=n;
		}
		rep(i,n+n-1){
			zz[i]=(zz[i]>>b)%2;
			zcnt[i]=min(i+1,n+n-i-1);
		}

		vvl tmp(n,vl(n,-1));
		
		ll type=2;
		rep(i,n){
			if(xx[i]==1){
				if(type==1){
					cout<<-1<<endl;
					return;
				}
				type=0;
				rep(j,n){
					tmp[i][j]=1;
					xcnt[i]=0;
					ycnt[j]=0;
					zcnt[i+j]--;
					zz[i+j]^=1;
				}
			}
			if(yy[i]==0){
				if(type==0){
					cout<<-1<<endl;
					return;
				}
				type=1;
				rep(j,n){
					tmp[j][i]=0;
					xcnt[j]=0;
					ycnt[i]=0;
					zcnt[i+j]--;
				}
			}
		}
		if(type==2){
			if(n<=2){
				//3以下は全探索
				rep(i,1LL<<(n*n)){
					vvl testing(n,vl(n));
					rep(j,n*n){
						if(i&(1LL<<j))testing[j/n][j%n]=1;
						else testing[j/n][j%n]=0;
					}
					vl ansx(n,1),ansy(n,0),ansz(n+n-1,0);
					rep(i,n){
						rep(j,n){
							ansx[i]&=testing[i][j];
							ansy[j]|=testing[i][j];
							ansz[i+j]^=testing[i][j];
						}
					}
					if(ansx==xx&&ansy==yy&&ansz==zz){
						tmp=testing;
						break;
					}
				}
				if(tmp[0][0]==-1){
					cout<<-1<<endl;
					return;
				}
				
				rep(i,n)rep(j,n)ans[i][j]+=tmp[i][j]*(1LL<<b);
				continue;
			}
			if(n==3){
				//3以下は全探索
				rep(i,1LL<<((n-1)*(n-1))){
					vvl testing(n,vl(n));
					rep(j,(n-1)*(n-1)){
						if(i&(1LL<<j))testing[j/(n-1)][j%(n-1)+1]=1;
						else testing[j/(n-1)][j%(n-1)+1]=0;
					}
					
					testing[0][0]=zz[0];
					testing[1][0]=zz[1]^testing[0][1];
					testing[2][0]=zz[2]^testing[1][1]^testing[0][2];
					testing[2][1]=zz[3]^testing[1][2];
					testing[2][2]=zz[4];

					vl ansx(n,1),ansy(n,0);
					rep(i,n){
						rep(j,n){
							ansx[i]&=testing[i][j];
							ansy[j]|=testing[i][j];
						}
					}
					if(ansx==xx&&ansy==yy){
						tmp=testing;
						break;
					}
				}
				if(tmp[0][0]==-1){
					cout<<-1<<endl;
					return;
				}
				
				rep(i,n)rep(j,n)ans[i][j]+=tmp[i][j]*(1LL<<b);
				continue;
			}
			//n>=4のみ解く。
			tmp[0][1]=0;
			tmp[1][1]=0;
			tmp[n-2][1]=1;
			tmp[n-2][0]=1;
			zz[n-2]^=1;
			zz[n-1]^=1;
			zcnt[1]--;
			zcnt[2]--;
			zcnt[n-2]--;
			zcnt[n-1]--;
			loop(i,2,n-1){
				tmp[1][i]=1;
				zz[1+i]^=1;
				zcnt[1+i]--;
				tmp[i][n-2]=0;
				zcnt[n-2+i]--;
			}

			//XORを調整する。
			rep(i,n)rep(j,n){
				if(tmp[i][j]==-1){
					tmp[i][j]=zz[i+j];
					zz[i+j]=0;
				}
			}

		}else{
			if(type==0){
				rep(i,n)rep(j,n){
					if(tmp[i][j]==-1&&zcnt[i+j]>=2){
						tmp[i][j]=0;
						xcnt[i]--;
						zcnt[i+j]--;
					}
				}
			}else{
				rep(j,n)rep(i,n){
					if(tmp[i][j]==-1&&zcnt[i+j]>=2){
						tmp[i][j]=1;
						zz[i+j]^=1;
						ycnt[j]--;
						zcnt[i+j]--;
					}
				}
			}
			//zcnt埋めが終わった奴をチェック
			rep(i,n+n-1){
				if(zcnt[i]==0){
					if(zz[i]==1){
						cout<<-1<<endl;
						return;
					}
				}
			}
			//zcntが1な場所を埋める
			rep(i,n)rep(j,n){
				if(zcnt[i+j]==1&&tmp[i][j]==-1){
					tmp[i][j]=zz[i+j];
					zcnt[i+j]=0;
					xcnt[i]=max(0LL,xcnt[i]-1);
					ycnt[j]=max(0LL,ycnt[j]-1);
				}
			}
		}

		vl ansx(n,1),ansy(n,0);
		rep(i,n){
			rep(j,n){
				ansx[i]&=tmp[i][j];
				ansy[j]|=tmp[i][j];
			}
		}
		ll hokan=0;
		rep(i,n){
			if(tmp[i][n-i-1]==type){
				hokan=i;
				break;
			}
		}
		rep(i,n){
			if(ansx[i]!=xx[i]){
				swap(tmp[hokan][n-hokan-1],tmp[i][n-i-1]);
			}
			if(ansy[i]!=yy[i]){
				swap(tmp[hokan][n-hokan-1],tmp[n-i-1][i]);
			}
		}

		rep(i,n)rep(j,n)ans[i][j]+=tmp[i][j]*(1LL<<b);
	}
	vl ansx(n,(1LL<<30)-1),ansy(n,0);
	rep(i,n){
		rep(j,n){
			ansx[i]&=ans[i][j];
			ansy[j]|=ans[i][j];
		}
	}
	if(ansx!=x||ansy!=y){
		cout<<-1<<endl;
		return;
	}
	vvdbg(ans);
	return;
}

//メイン
int main(){
	ll t;
	cin>>t;
	rep(i,t)solve();
}