def main():
    grid = [[0] * (N + 2)]
    apnd = grid.append
    for n in A:
        apnd([0] + n + [0])
    apnd([0] * (N + 2))

    for l in L:
        case1 = search(l, grid)
        _, row1, col1, row2, col2 = case1
        for i in range(row1, row2 + 1):
            line = grid[i][col1:col2]
            grid[i][col1:col2] = list(map(lambda x: int(not x), line))
        print(row1, col1, row2, col2 - 1)

def search(l, grid):
    maxb = -1
    for i in range(1, N + 1):
        for j in range(1, N + 2 - l):
            black = grid[i][j:j+l].count(1)
            if l == 1 and black == 0 and grid[i][j-1] == grid[i][j+1] == 1:
                maxb = 0
                return (maxb, i, j, i, j)
            elif black == l and grid[i][j-1] == grid[i][j+l] == 0:
                maxb = black
                return (maxb, i, j, i, j + l)
            elif maxb < black:
                maxb = black
                ret = (maxb, i, j, i, j + l)
    return ret

N, K = map(int, input().split())
L = tuple(map(int, input().split()))
A = [list(map(int, list(input()))) for _ in [0] * N]
main()