Как найти координаты наибольшего квадрата из единиц в матрице.Python

Сейчас решаю задачу. Решил применить формулу a[i][j] = min(a[i - 1][j - 1], a[i][j - 1], a[i - 1][j]) + 1. Теперь задаюсь вопросом, как мне найти координаты(x,y) этой квадратной матрицы?Код ниже:

coor = input().split()
N = coor[0]
W = coor[1] #Высота
H = coor[2] #Ширина
cl = []
score = 0
a = [[1 for y in range(int(W))] for x in range(int(H)+1)]
for x in range(int(N)):
    tree_coor = input().split()
    x = int(tree_coor[0])
    y = int(tree_coor[1])
    a[y][x] = 0
for i in range(len(a)):
    for j in range(len(a[i])):
        if i*j != 0:
            if a[i][j] == 1:
                a[i][j] = min(a[i - 1][j - 1], a[i][j - 1], a[i - 1][j]) + 1
            if (a[i][j] > score):
                score=a[i][j]
print(score)