#include<cstdio>

 #include<cstdlib>

 #include<algorithm>

 #include<cstring>

 using namespace std;

 int a[][];

 int ans[][];

 int c(int x,int y){return (x*+y+);}

 int d[]={,-,,,};

 void gauss(int m,int n){

     int line=;

     for(int k=;k<=m;k++){

         int i=line;

         while(i<=m){if(a[i][k])break;i++;}

         if(i>m)continue;

         if(i!=line){swap(a[i],a[line]);}

         for(i=;i<=m;i++){

             if(i!=line&&a[i][k]){

                 for(int j=k;j<=n;j++){

                     a[i][j]^=a[line][j];

                 }

             }

         }

         line++;

     }

 //    for(int i=1;i<=30;i++){

 //        for(int j=1;j<=30;j++){

 //            printf("%d ",a[i][j]);

 //        }

 //        printf("\n");

 //    }

     for(int i=;i<;i++){

         for(int j=;j<;j++){

             ans[i][j]=a[c(i,j)][n];

         }

     }

 }

 void solve(){

     memset(a,,sizeof(a));

     memset(ans,,sizeof(ans));

     int x;

     for(int i=;i<;i++){

         for(int j=;j<;j++){

             scanf("%d",&x);

             a[c(i,j)][]=x;

             a[c(i,j)][c(i,j)]=;

             for(int k=;k<;k++){

                 int x=i+d[k],y=j+d[k+];

                 if(x>=&&x<&&y>=&&y<)a[c(i,j)][c(x,y)]=;

             }

         }

     }

     gauss(,);

     for(int i=;i<;i++){

         for(int j=;j<;j++){

             printf("%d ",ans[i][j]);

         }

         printf("%d\n",ans[i][]);

     }

 }

 int main()

 {

 //    freopen("data.in","r",stdin);

     int T;

     scanf("%d",&T);

     for(int i=;i<=T;i++){

         printf("PUZZLE #%d\n",i);

         solve();

     }

     return ;

 }