Ruudukko

import random

# n = 50
# rows = []
# for i2 in range(n):
#     rows.append(random.sample(range(1, n*n+1), n*n))
# n = 3
# rows = [
#     [2,1,3], 
#     [1,1,1], 
#     [9,2,7]
# ]

# for row in rows:
#     for element in row:
#         print(element, end=" ")
#     print()

n = int(input())
rows = []
for i in range(n):
    rows.append([int(x) for x in input().split(" ")])

unique_numbers = []
for i2 in range(n):
    for j1 in range(n):
        if(rows[i2][j1] not in unique_numbers):
            unique_numbers.append(rows[i2][j1])
unique_numbers = sorted(unique_numbers)

def count_global(lambda_function):
    c = 0
    for i in range(n):
        for j in range(n):
            if lambda_function(rows[i][j]):
                c += 1
    return c

def get_numbers_on_row_and_column(x, y, include_self):
    numbers = []
    for i in range(n):
        if include_self:
            numbers.append(rows[y][i])
        else:
            if i != x:
                numbers.append(rows[y][i])
    for i in range(n):
        if i != y:
            numbers.append(rows[i][x])
    return numbers

def get_coordinates_of(lambda_function):
    coords = []
    for i in range(n):
        for j in range(n):
            if lambda_function(rows[i][j]):
                coords.append((i, j))
    return coords

def count_elements_on_same_row_and_column(x, y, lambda_function):
    c = 0
    for i in range(n):
        if lambda_function(rows[y][i]):
            c += 1
        if lambda_function(rows[i][x]):
            c += 1
    return c

cache = {}

for num in unique_numbers:

    coords = get_coordinates_of(lambda x: x == num)
    if num == unique_numbers[0]:
        for c in coords:
            cache[c] = 1
    else:
        for c in coords:
            lines = []
            i1 = c[0]
            for j1 in range(n):
                if rows[i1][j1] < num:
                    lines.append((i1, j1))
            
            j2 = c[1]
            for i2 in range(n):
                if rows[i2][j2] < num:
                    lines.append((i2, j2))



            cache[c] = sum([cache[x] for x in lines]) + 1

print(sum(cache.values()) % (10**9 + 7))