import sys
import numpy
from math import inf
from array import array
S = 16384

mapa = []
input_file = open("01.in", "rb")
for y in range(S):
    # Musíme dát pozor na endianitu (data jsou
    # little-endian) a tak načtení svěříme numpy
    numpy_radek = numpy.fromfile(
        input_file,
        dtype="<u2",  # unsigned 2bajtová čísla
        count=S       # v little-endian(<)
    )
    # Pro rychlejší přístup k prvkům převedeme
    # na interní Pythoní pole
    radek = array('H')  # minimálně 2b číslo
    radek.fromlist(numpy_radek.tolist())
    mapa.append(radek)

print("loaded")
koupene = []
yk = 0
while yk <= S:
    new_yk = 2*S
    xk = 0
    while xk <= S:
        print(f"\rkvadrant {xk} {yk} ", end="", flush=True)
        best = (inf, (-1, -1))
        for x in range(xk, min(xk+501, S)):
            for y in range(yk, min(yk+501, S)):
                if mapa[y][x] != 0:
                    score = mapa[y][x] - 3.5*(x - xk) - 0.08*(y-yk)  # 3 0.1
                    if score < best[0]:
                        best = (score, (y, x))
        if best[0] != inf:
            new_yk = min(new_yk, best[1][0] + 501)
            koupene.append(best[1])
            xk = best[1][1] + 501
        else:
            new_yk = yk + 501
            xk += 501
    yk = new_yk

print("\nVýpis", flush=True)
f = open("PPMrizka.out", "w")
f.write(f"{len(koupene)}\n")
for dum in koupene:
    f.write(f"{dum[0]} {dum[1]}\n")
f.close()