#define _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <iomanip>
#include <cmath>
#include <algorithm>
#include <vector>
#define vd vector<double>
#define vvd vector<vector<double> >
#define ui unsigned int
using namespace std;
vvd a;
vvd b;
vd temp;
int n;
double eps = 1e-6;
vd operator /( vd a, double num)
{
vd rez;
rez.resize(a.size());
for (int i = 0; i < a.size(); ++i) {
rez[i] = a[i] /num;
}
return rez;
}
vd operator -( vd a, vd b)
{
vd rez;
rez.resize(a.size());
for (int i = 0; i < a.size(); ++i) {
rez[i] = a[i] - b[i];
}
return rez;
}
vd operator *( vd a, double num)
{
vd rez;
rez.resize(a.size());
for (int i = 0; i < a.size(); ++i) {
rez[i] = a[i] * num;
}
return rez;
}
int main ()
{
freopen("input.txt", "r", stdin);
cin >> n;
a.resize(n, vd(n));
for (int i = 0; i < n; ++i) {
for (int j = 0; j < n ; ++j) {
cin >> a[i][j];
}
}
for (int i = 0; i < n; ++i) {
for (int j = 0; j < n ; ++j) {
cout << a[i][j] <<" ";
}
cout << endl;
}
for (int N = 0; N < n; ++N) {
for (int j = N; j < n; ++j) {
for (int i = N; i < n ; ++i) {
if (fabs(a[i][j]) > eps)
{
if (fabs(a[N][j]) < eps)
{
temp = a[N];
a[N] = a[i];
a[i] = temp;
}
a[N] = a[N]/a[N][N];
for (int k = N + 1; k < n; ++k)
{
a[k] = a[k] - a[N]*a[k][N];
}
}
}
}
}
for (int i = 0; i < n; ++i) {
for (int j = 0; j < n ; ++j) {
cout << a[i][j] <<" ";
}
cout << endl;
}
for (int M = n - 1; M > 0; --M)
{
for (int s = M - 1; s >= 0; --s)
{
a[s] = a[s] - a[M] * a[s][M];
}
}
for (int i = 0; i < n; ++i) {
for (int j = 0; j < n ; ++j) {
cout << a[i][j] <<" ";
}
cout << endl;
}
}