#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;

using ll=long long;
using ldub=long double;
using lldub=__float128;
using str=string;

using mint=modint;

template<class T=ll>
using tup2=tuple<T,T>;
template<class T=ll>
using tup3=tuple<T,T,T>;
template<class T=ll>
using tup4=tuple<T,T,T,T>;
template<class T=ll>
using tup5=tuple<T,T,T,T,T>;
template<class T=ll>
using tup6=tuple<T,T,T,T,T,T>;
template<class T=ll>
using tup7=tuple<T,T,T,T,T,T,T>;

template<class T=ll>
using vec=vector<T>;
template<class T=ll>
using vec2=vector<vec<T>>;
template<class T=ll>
using vec3=vector<vec2<T>>;
template<class T=ll>
using vec4=vector<vec3<T>>;
template<class T=ll>
using vec5=vector<vec4<T>>;
template<class T=ll>
using vec6=vector<vec5<T>>;

template<class T>
using que=queue<T>;
template<class T>
using Pque=priority_queue<T>;
template<class T>
using pque=priority_queue<T, vector<T>, greater<T>>;

struct Edge{
  ll from,to,w=1,num=-1;
};
using gvec=vector<Edge>;
using gvec2=vector<gvec>;

template<class T>
bool chmax(T &a,T b){if(a<b){a=b;return 1;}return 0;}
template<class T>
bool chmin(T &a,T b){if(b<a){a=b;return 1;}return 0;}

#define pb push_back
#define pf push_front
#define pob pop_back
#define pof pop_front
#define ins insert
#define tup make_tuple
#define low lower_bound
#define pops(a) __builtin_popcountll(a)
#define ctz(a) __builtin_ctzll(a)
#define all(v) v.begin(),v.end()
#define rall(v) v.rbegin(),v.rend()
#define nperm(v) next_permutation(all(v))
#define Yes cout << "Yes" << endl
#define No cout << "No" << endl
#define YN(f) cout << ((f)?"Yes":"No") << endl

const int IINF=2e9;
const ll INF=4e18;

#ifndef LOCAL
#define debug(x)
#endif




// =================================================




map<vec<ll>,vec2<ll>> M;

void solve(){
  ll N;
  cin >> N;
  vec X(N),Y(N),Z(2*N-1);
  for(ll i=0;i<N;i++) cin >> X[i];
  for(ll j=0;j<N;j++) cin >> Y[j];
  for(ll t=0;t<2*N-1;t++) cin >> Z[t];
  vec2 A(N,vec(N));
  for(ll k=0;k<30;k++){
    ll f=0,g=0;
    for(ll i=0;i<N;i++) if(((X[i]>>k)&1)==1) f=1;
    for(ll j=0;j<N;j++) if(((Y[j]>>k)&1)==0) g=1;
    if(f&&g){
      cout << -1 << endl;
      return;
    }
    vec2<ll> F(N,vec<ll>(N,-1));
    if(f){
      for(ll i=0;i<N;i++){
        for(ll j=0;j<N;j++) if(((X[i]>>k)&1)==1) F[i][j]=1;
      }
      vec C(N);
      for(ll i=0;i<N;i++) if(((X[i]>>k)&1)==1) C[i]=N;
      for(ll t=0;t<2*N-1;t++){
        ll s=0;
        for(ll i=0;i<=t;i++){
          ll j=t-i;
          if(i>=N||j>=N) continue;
          if(F[i][j]!=-1) s^=1;
        }
        for(ll i=0;i<=t;i++){
          ll j=t-i;
          if(i>=N||j>=N) continue;
          if(F[i][j]!=-1) continue;
          if((((Z[t])>>k)&1)==s||C[i]==N-1) F[i][j]=0;
          else F[i][j]=1,s^=1,C[i]++;
        }
      }
    }
    if(g){
      for(ll i=0;i<N;i++){
        for(ll j=0;j<N;j++) if(((Y[j]>>k)&1)==0) F[i][j]=0;
      }
      vec C(N);
      for(ll j=0;j<N;j++) if(((Y[j]>>k)&1)==0) C[j]=N;
      for(ll t=0;t<2*N-1;t++){
        ll s=0;
        for(ll j=0;j<=t;j++){
          ll i=t-j;
          if(i>=N||j>=N) continue;
          if(F[i][j]==-1) s^=1;
        }
        for(ll j=0;j<=t;j++){
          ll i=t-j;
          if(i>=N||j>=N) continue;
          if(F[i][j]!=-1) continue;
          if((((Z[t])>>k)&1)==s||C[j]==N-1) F[i][j]=1;
          else F[i][j]=0,s^=1,C[j]++;
        }
      }
    }
    if(!f&&!g){
      if(N<=3){
        vec E(2*N-1);
        for(ll t=0;t<2*N-1;t++) E[t]=(Z[t]>>k)&1;
        if(!M.count(E)){
          cout << -1 << endl;
          return;
        }
        F=M[E];
      }else{
        for(ll i=0;i<N;i++){
          for(ll j=0;j<N;j++){
            if(0<=i-j&&i-j<=1) continue;
            F[i][j]=(i+j)&1;
          }
        }
        for(ll t=0;t<2*N-1;t++){
          ll s=0;
          for(ll i=0;i<=t;i++){
            ll j=t-i;
            if(i>=N||j>=N) continue;
            if(F[i][j]!=-1) s^=1;
          }
          if(((Z[t]>>k)&1)==s) F[t-t/2][t/2]=0;
          else F[t-t/2][t/2]=1;
        }
      }
    }
    for(ll i=0;i<N;i++){
      for(ll j=0;j<N;j++) A[i][j]+=F[i][j]<<k;
    }
  }
  {
    vec B(N);
    for(ll i=0;i<N;i++){
      B[i]=A[i][0];
      for(ll j=1;j<N;j++) B[i]&=A[i][j];
    }
    if(B!=X){
      cout << -1 << endl;
      return;
    }
  }
  {
    vec B(N);
    for(ll j=0;j<N;j++){
      B[j]=A[0][j];
      for(ll i=1;i<N;i++) B[j]|=A[i][j];
    }
    if(B!=Y){
      cout << -1 << endl;
      return;
    }
  }
  {
    vec B(2*N-1);
    for(ll i=0;i<N;i++){
      for(ll j=0;j<N;j++) B[i+j]^=A[i][j];
    }
    if(B!=Z){
      cout << -1 << endl;
      return;
    }
  }
  for(ll i=0;i<N;i++){
    for(ll j=0;j<N;j++) cout << A[i][j] << ' ';
    cout << endl;
  }
}

int main(){
  ios::sync_with_stdio(false);
  cin.tie(nullptr);

  mint::set_mod(998244353);
  // mint::set_mod(1000000007);

  for(ll n=1;n<=3;n++){
    for(ll b=0;b<(1<<(n*n));b++){
      vec2 F(n,vec(n));
      for(ll k=0;k<(n*n);k++) F[k/n][k%n]=(b>>k)&1;
      ll f=1;
      for(ll i=0;i<n;i++){
        ll g=0;
        for(ll j=0;j<n;j++) if(F[i][j]==0) g=1;
        if(!g) f=0;
      }
      for(ll j=0;j<n;j++){
        ll g=0;
        for(ll i=0;i<n;i++) if(F[i][j]==1) g=1;
        if(!g) f=0;
      }
      if(!f) continue;
      vec E(2*n-1);
      for(ll i=0;i<n;i++) for(ll j=0;j<n;j++) E[i+j]^=F[i][j];
      M[E]=F;
    }
  }

  ll T=1;
  cin >> T;
  while(T--) solve();
  return 0;
}