0705 - 2D

Enunt

Gigel îşi imaginează lumea în varianta 2d, adică reprezentată în sistem de coordonate cartezian XOY. Fiecare persoană din grupul celor N prieteni ai săi este reprezentată în plan printr-un punct identificat prin abscisa şi ordonata sa. În lumea sa 2d, plouă ca în Anglia, iar picăturile de ploaie pică paralel cu axa OY, de la o înălţime infinită. Ca să îi ferească pe prietenii săi de ploaie, îşi propune să le construiască apărători pe care le va reprezenta pe hartă prin segmente de dreaptă.

Cerința

Ştiind că nu poate să deseneze pe hartă decât segmente de lungimi egale, determinaţi care este lungimea minimă a unui segment astfel încât trasând cel mult K segmente, toți cei N prieteni ai săi să fie protejați de ploaie.

Date de intrare

Fişierul de intrare 2din.txt conţine:

• pe prima linie numerele naturale N şi K;
• pe următoarele N linii se află perechi de numere reale x y, cu exact trei zecimale fiecare, reprezentând abscisa şi ordonata punctului corespunzător fiecărui prieten al lui Gigel;

Date de ieșire

Fişierul de ieşire 2dout.txt va conţine pe prima linie un singur număr real, cu trei zecimale reprezentând lungimea minimă a unui segment din cele trasate de Gigel.

Restricții și precizări

• 1 ⩽ K < N ⩽ 200000
• -1000000.000 ⩽ x,y ⩽ 1000000.000
• Lungimea minimă a unui segment pe care îl va trasa Gigel este 1.000
• Rezultatul va fi afişat cu o eroare de cel mult 0.001
• Un segment trasat între punctele (x1, y1) şi (x2, y2) protejează toate punctele de pe segment şi din semiplanul “de sub” dreapta ce trece prin aceste puncte şi care au abscisele în intervalul închis [x1, x2].

Exemplu 1

2din.txt

5 3

-5.000 1.000

-2.000 3.000

3.000 2.000

3.000 -2.000

1.000 2.000

2d.out

2.000

Exemplu 2

2din.txt

0 0

0 0

0 0

2dout.txt

Nu au fost respectate cerintele impuse

Rezolvare

#0705 - 2D
import math

def is_covered(points, length, K):
    points.sort()

    last_point = points[0]
    segments = 1

    for point in points:
        if point - last_point > length:
            segments += 1
            last_point = point

    return segments <= K

def find_min_segment_length(N, K, points):
    left, right = 0, max(points)

    while abs(left - right) > 0.001:
        mid = (left + right) / 2

        if is_covered(points, mid, K):
            right = mid
        else:
            left = mid

    return round(right, 3)

def check_restrictions(N, K, points):
    if not(1 <= K < N <= 200000):
        return False
    for point in points:
        if not (-1000000.000 <= point <= 1000000.000):
            return False
    return True

def main():
    with open("2din.txt", "r") as file_in, open("2dout.txt", "w") as file_out:
        N, K = map(int, file_in.readline().split())
        points = [float(file_in.readline().split()[0]) for _ in range(N)]

        if not check_restrictions(N, K, points):
            file_out.write("Nu au fost respectate cerintele impuse\n")
            return

        result = find_min_segment_length(N, K, points)
        file_out.write(f"{result:.3f}\n")

if __name__ == "__main__":
    main()