Задание 27. Информатика. ЕГЭ. ЕГКР. 19.04.2025

Фрагмент звёздного неба спроецирован на плоскость с декартовой системой координат. Учёный решил провести кластеризацию полученных точек, являющихся изображениями звёзд, то есть разбить их множество на \(N\) непересекающихся непустых подмножеств (кластеров), таких что точки каждого подмножества лежат внутри прямоугольника со сторонами длиной \(H\) и \(W,\) причём эти прямоугольники между собой не пересекаются. Стороны прямоугольников не обязательно параллельны координатным осям. Гарантируется, что такое разбиение существует и единственно для заданных размеров прямоугольников.

Будем называть центром кластера точку этого кластера, сумма расстояний от которой до всех остальных точек кластера минимальна. Для каждого кластера гарантируется единственность его центра. Расстояние между двумя точками на плоскости \(A(x_1, \, y_1)\) и \(B(x_2, \, y_2)\) вычисляется по формуле: $$d(A, \, B) = \sqrt{(x_2 - x_1)^2 + (y_2 - y_1)^2}$$ В файле A хранятся данные о звёздах двух кластеров, где \(H=11, \, W=11\) для каждого кластера. В каждой строке записана информация о расположении на карте одной звезды: сначала координата \(x,\) затем координата \(y.\) Значения даны в условных единицах. Известно, что количество звёзд не превышает \(1000.\)

В файле B хранятся данные о звёздах трёх кластеров, где \(H=13, \, W=13\) для каждого кластера. Известно, что количество звёзд не превышает \(10~000.\) Структура хранения информации о звездах в файле B аналогична файлу A.

Для каждого файла определите координаты центра каждого кластера, затем вычислите два числа: \(P_x\) — среднее арифметическое абсцисс центров кластеров, и \(P_y\) — среднее арифметическое ординат центров кластеров. В ответе запишите четыре числа: в первой строке сначала абсолютное значение целой части произведения \(P_x \times 10000,\) затем абсолютное значение целой части произведения \(P_y \times 10000\) для файла A, во второй строке — аналогичные данные для файла B. Возможные данные одного из файлов иллюстрированы графиком.

Внимание! График приведён в иллюстративных целях для произвольных значений, не имеющих отношения к заданию. Для выполнения задания используйте данные из прилагаемого файла.

from math import dist

files = ['27_A.txt', '27_B.txt']

for task in (0, 1):
    data = []
    fd = open(files[task])
    fd.readline()
    for line in fd:
        x, y = map(float, line.replace(',', '.').split())
        data.append((x, y))
    clusters = []
    eps = [2, 2]
    while data:
        clusters.append([data.pop()])
        for p in clusters[-1]:
            neigh= [pt for pt in data if dist(p, pt) < eps[task]]
            clusters[-1] += neigh
            for pt in neigh:
                data.remove(pt)
    #print(len(clusters), [len(cl) for cl in clusters])
    centers = []
    for cl in clusters:
        dmin = float('inf')
        c = None
        for p in cl:
            d = sum(dist(p, pt) for pt in cl)
            if d < dmin:
                dmin = d
                c = p
        centers.append(c)
    px = sum(c[0] for c in centers) / len(centers)
    py = sum(c[1] for c in centers) / len(centers)
    print(abs(int(10_000 * px)), abs(int(10_000 * py)))