#include <stdio.h>
#include <assert.h>
#pragma GCC target "avx2"
#pragma GCC optimize "Ofast"
#include <immintrin.h>

typedef __m256i m256;
typedef unsigned uint;
typedef long long ll;
typedef unsigned long long ull;
#define ms(a, x) memset(a, x, sizeof(a))
#define pb push_back
#define fs first
#define sn second
#define ALL(v) begin(v), end(v)
#define SZ(v) ((int) (v).size())
#define lbv(v, x) (lower_bound(ALL(v), x) - (v).begin())
#define ubv(v, x) (upper_bound(ALL(v), x) - (v).begin())
const int INF = 0x3f3f3f3f;
const ll LLINF = 0x3f3f3f3f3f3f3f3f;
#define FR(i, n) for(int i = 0; i < (n); i++)
#define FOR(i, a, b) for(int i = (a); i < (b); i++)
#define FORR(i, a, b) for(int i = (a); i >= (b); i--)
#define dbg(x) {cerr << #x << ' ' << x << endl;}
#define dbgArr(arr, n) {cerr << #arr; FR(_i, n) cerr << ' ' << (arr)[_i]; cerr << endl;}

#define vload(ptr) _mm256_loadu_si256((m256*) (ptr))
#define vstore(ptr, v) _mm256_storeu_si256((m256*) (ptr), (v))
#define vand(a, b) _mm256_and_si256(a, b)
#define vor(a, b) _mm256_or_si256(a, b)
#define vxor(a, b) _mm256_xor_si256(a, b)
#define vblendPs(a, b, c) (m256) _mm256_blendv_ps((__m256) (a), (__m256) (b), (__m256) (c))

const int S12 = 9740686;;
const int C24_11 = 2496144;
const int C24_12 = 2704156;
const int MASK = ((1 << 24) - 1);
const int MN = 24;
int n, mask;
int adjMat[MN], iAdjMat[MN];
bool ans[MN][MN];
int dpMem[C24_11+C24_12][8], (*dpMem1)[8] = dpMem, (*dpMem2)[8] = dpMem+C24_11;
int popPtr[MN+1], popSz[MN+1];
char mp[3<<MN], popMem[3*S12];

#define atMp(i) (*(int*) (mp + (i)*3))
#define atPopMem(i) (*(int*) (popMem + (i)*3))

void solvePop(int pop, int oldDp[][8], int newDp[][8]){
    m256 zeroV = _mm256_setzero_si256();

    m256 vAdj[MN];
    FR(i, MN) vAdj[i] = _mm256_set1_epi32(adjMat[i]);

    int inv = pop > n/2 ? ((1<<n)-1) : 0;

    FR(i2, popSz[pop]){
        int s2 = (atPopMem(popPtr[pop]+i2) ^ inv) & MASK;
        m256 res = zeroV;
        m256 s2v = _mm256_set1_epi32(s2);

        #define proc(n2){ \
            int s1 = s2 & ~(1 << n2); \
            int i1 = atMp(s1) & MASK; \
            m256 maskV = _mm256_slli_epi32(vand(vload(oldDp[i1]), s2v), 31-n2); \
            res = vor(res, vblendPs(zeroV, vAdj[n2], maskV)); \
        }

        proc(0);
        proc(1);
        proc(2);
        proc(3);
        proc(4);
        proc(5);
        proc(6);
        proc(7);
        proc(8);
        proc(9);
        proc(10);
        proc(11);
        proc(12);
        proc(13);
        proc(14);
        proc(15);
        proc(16);
        proc(17);
        proc(18);
        proc(19);
        proc(20);
        proc(21);
        proc(22);
        proc(23);

        vstore(newDp+i2, res);
    }
}

void solve(int *ord){
    FR(i, n) adjMat[i] = 0;
    FR(i, n) FR(j, n) adjMat[i] |= ((iAdjMat[ord[i]] >> ord[j]) & 1) << j;

    int (*oldMem)[8] = dpMem2, (*newMem)[8] = dpMem1;

    FR(n1, 8) oldMem[0][n1] = adjMat[n1];

    FR(pop, n-2){
        solvePop(pop+1, oldMem, newMem);
        int (*tmp)[8] = oldMem;
        oldMem = newMem, newMem = tmp;
    }

    FR(i, 8) FR(j, n) if(i != j){
        ans[ord[i]][ord[j]] = (oldMem[atMp(mask^(1<<i)^(1<<j)) & MASK][i] >> j) & 1;
    }
}

int choose[MN+1][MN+1];

void init(){
    FR(a, MN+1){
        choose[a][0] = 1;
        FOR(b, 1, a+1) choose[a][b] = choose[a-1][b-1] + choose[a-1][b];
    }
    FR(i, n+1) popSz[i] = choose[n][i];
    FR(i, n/2) popPtr[i+1] = popPtr[i] + popSz[i];
    FOR(i, n/2+1, n+1) popPtr[i] = popPtr[n-i];

    int ind[MN+1] = {};

    FR(s, 1<<n){
        int pop = __builtin_popcount(s);
        int i = ind[pop]++;
        if(pop > n/2) i = popSz[pop]-1 - i;
        else atPopMem(popPtr[pop]+i) |= s;
        atMp(s) |= i;
    }
}

int main(){
    scanf("%d\n", &n);
    FR(i, n){
        FR(j, n) iAdjMat[i] |= (getchar()-'0') << j;
        getchar();
    }

    //n = 24;
    //FR(i, n) FR(j, i) if(rng() & 1) iAdjMat[i] |= 1 << j, iAdjMat[j] |= 1 << i;

    mask = (1<<n)-1;
    init();

    int ord[MN];
    FR(i, MN) ord[i] = i;
    solve(ord);
    FR(i, n) ord[i] = (i+8) % n;
    solve(ord);
    FR(i, n) ord[i] = (i+16) % n;
    solve(ord);

    FR(i, n){
        FR(j, n) putchar('0'+ans[i][j]);
        putchar('\n');
    }
}