crout_parcial

n = int(raw_input())
A = []

for i in range(n):
    A.append([])
    for j in range (n):
        A[i].append(float(raw_input()))

a = []
L = []
U = []

for i in range(n):
    a.append([])
    for j in range(n):
        a[i].append(A[i][j])

for i in range(n):
    L.append([])
    U.append([])
    for j in range(n):
        L[i].append(0.0)
        U[i].append(0.0)


print'\n\n\nMatriz A:\n\n'

for i in range(n):
    for j in range(n):
        print a[i][j],
    print'\n'

print '\n\n'


cont = 0

for k in range (n):
    pos = k
    maxi = -1000000000000.0
    for i in range(k,n):
        res = 0.0
        for p in range(k):
            res += L[i][p]*U[p][k]
            cont+=2
        L[i][k] = a[i][k]-res
        #a[i][k]=L[i][k]
        cont+=1

    for i in range(k,n):
        if maxi < abs(a[i][k]):
            pos = i

    a[pos], a[k] = a[k], a[pos]

    for i in range (k,n):
        res = 0.0
        for p in range(k):
            res += L[k][p]*U[p][i]
            cont+=2
        U[k][i] = (a[k][i]-res)/L[k][k]
        cont+=2


for i in range(n):
    for j in range(n):
        print L[i][j],
    print'\n'

print '\n\n'

for i in range(n):
    for j in range(n):
        print U[i][j],
    print '\n'
print '\n\n'

print cont