#include<cstdio>

#include<cctype>

#include<algorithm>

#define R register int

using namespace std;

const int maxn=;

const int MOD=1e9+;

int n;

struct Matrix{

    int m[maxn][maxn];

    void SWAP(int x,int y){   //交换两行

        for(R i=;i<=n;++i)

            swap(m[x][i],m[y][i]);

    }

    void MUL(int x,int k){    //对一行乘上k

        for(R i=;i<=n;++i)

            m[x][i]=(1LL*m[x][i]*k%MOD+MOD)%MOD;

    }

    void MULADD(int x,int y,int k){ //将第y行乘上k加到第x行上去

        for(R i=;i<=n;++i)

            m[x][i]=((m[x][i]+1LL*m[y][i]*k%MOD)%MOD+MOD)%MOD;

    }

}A,B;

int qpow(int a,int b){   //求逆元

    int ret=;

    while(b){

        if(b&) ret=1LL*ret*a%MOD;

        a=1LL*a*a%MOD;

        b>>=;

    }

    return ret;

}

void Invmatrix(){

    for(R i=;i<=n;++i){

        if(!A.m[i][i]){

            for(R j=i+;j<=n;++j)

                if(A.m[j][i]){

                    B.SWAP(i,j);

                    A.SWAP(i,j);

                    break;

                }

        }

        if(!A.m[i][i]){  //没有逆矩阵

            puts("No Solution");

            return;

        }

        int tmp=qpow(A.m[i][i],MOD-);

        B.MUL(i,tmp);

        A.MUL(i,tmp);    //系数化为1

        for(R j=i+;j<=n;++j){  //消元

            tmp=-A.m[j][i];

            B.MULADD(j,i,tmp);

            A.MULADD(j,i,tmp);

        }

    }

    for(R i=n-;i>=;--i)  //回带

        for(R j=i+;j<=n;++j){

            int tmp=-A.m[i][j];

            B.MULADD(i,j,tmp);

            A.MULADD(i,j,tmp);

        }

    for(R i=;i<=n;++i){

        for(R j=;j<=n;++j){

            printf("%d",B.m[i][j]);

            if(j!=n) printf(" ");

        }

        printf("\n");

    }

}

int main(){

    scanf("%d",&n);

    for(R i=;i<=n;++i)

        for(R j=;j<=n;++j)

            scanf("%d",&A.m[i][j]);

    for(R i=;i<=n;++i){

        for(R j=;j<=n;++j)

            B.m[i][j]=;

        B.m[i][i]=;

    }

    Invmatrix();

    return ;

}