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
        hline = grid[row1][col1:col2]
        grid[row1][col1:col2] = list(map(lambda x: int(not x), hline))
        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 == black == 1:
                if grid[i][j-1] == grid[i][j+1] == grid[i-1][j] == grid[i+1][j] == 0:
                    maxb = black
                    ret = [maxb, i, j, i, j + l]
                    return ret
            if 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()