Каналы на Марсе

def tra(mat):
    tm = list()
    for i in range(len(mat)):
        tm.append(list(len(mat[i]) * '0'))
    for i in range(len(mat)):
        for j in range(len(mat)):
            tm[j][i] = mat[i][j]
    return tm
def scm(mat, lt):
    fug = list()
    for i in range(len(mat)):
        fg = list()
        arg = 1
        a = mat[i][0]
        for i1 in range(1, len(mat[i])):
            if mat[i][i1] == a:
                arg += 1
            else:
                fg.append([arg, a])
                arg = 1
                a = mat[i][i1]
        fg.append([arg, a])
        fg.sort()
        if len(fg) > 1:
            gh = fg[-1]
            for i in range(len(fg) - 2, -1, -1):
                if fg[i][0] == gh[0]:
                    if fg[i][1] == gh[1]:
                        pass
                    else:
                        fug.append([fg[i][0], "?"])
                        break
                else:
                    fug.append(fg[-1])
                    break
        else:
            fug.append(fg[0])
    if lt == 1:
        return fug
    elif lt == 0:
        fug.sort()
        gh = fug[-1]
        fg = list()
        f = False
        for i in range(len(fug) - 2, -1, -1):
            if fug[i][0] == gh[0]:
                if fug[i][1] == gh[1]:
                    pass
                else:
                    return [fug[i][0], "?"]
                    f = True
                    break
            else:
                return fug[-1]
                f = True
                break
        if not f:
            return fug[-1]
N = int(input())
POLE = []
ALL_0 = []
ALL_1 = []
ALL_2 = []
OTVET = []
for i in range(N):
    POLE.append(list(input())) # заполняем список POLE, который будет хранить все поле посевов, где в каждой ячейке хранится 
if N == 1:
    print(POLE[0][0])
    print(1)
else:            
    test = POLE
    max_col = len(test[0])
    max_row = len(test)
    cols = [[] for _ in range(max_col)]
    rows = [[] for _ in range(max_row)]
    fdiag = [[] for _ in range(max_row + max_col - 1)]
    bdiag = [[] for _ in range(len(fdiag))]
    min_bdiag = -max_row + 1
 
    for x in range(max_col):
        for y in range(max_row):
            cols[x].append(test[y][x])
            rows[y].append(test[y][x])
            fdiag[x+y].append(test[y][x])
            bdiag[x-y-min_bdiag].append(test[y][x])
    ALL_0.append(fdiag)
    ALL_0.append(bdiag)
    ALL_1 = [c for d in ALL_0 for c in d]
    for i in range (N):
        ALL_1.append(POLE[i])
    ALL_2 = (tra(POLE))
    for i in range (N):
        ALL_1.append(ALL_2[i])
    #print(ALL_1)
    #print(ALL_1)
    OTVET.append(scm(ALL_1,0))
    print (OTVET[0][1])
    print (OTVET[0][0])