Фрагмент звёздного неба спроецирован на плоскость с декартовой системой координат. Учёный решил провести кластеризацию полученных точек, являющихся изображениями звёзд, то есть разбить их множество на \(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=6,\) \(W=4{,}5\) для каждого кластера. В каждой строке записана информация о расположении на карте одной звезды: сначала координата \(x,\) затем координата \(y.\) Значения даны в условных единицах. Известно, что количество звёзд не превышает \(1000.\)
В файле B хранятся данные о звёздах трёх кластеров, где \(H=6,\) \(W = 5\) для каждого кластера. Известно, что количество звёзд не превышает \(1000.\) Структура хранения информации о звездах в файле B аналогична файлу A.
Известно, что в файле B имеются координаты ровно трёх «лишних» точек, являющихся аномалиями, возникшими в результате помех при передаче данных. Эти три точки не относятся ни к одному из кластеров, их учитывать не нужно.
Для файла A определите координаты центра каждого кластера, затем найдите два числа: \(P_x\) — сумму абсцисс центров кластеров, и \(P_y\) — сумму ординат центров кластеров. Для файла B найдите два числа: \(Q_1\) — минимальное расстояние между точками, принадлежащими двум различным кластерам, и \(Q_2\) — максимальное расстояние между точками, принадлежащими двум различным кластерам.
В ответе запишите четыре числа: в первой строке — сначала абсолютную величину целой части произведения \(P_x \times 10~000,\) затем абсолютную величину целой части произведения \(P_y \times 10~000;\) во второй строке — сначала целую часть произведения \(Q_1 \times 10~000,\) затем целую часть произведения \(Q_2 \times 10~000.\)
Возможные данные одного из файлов иллюстрированы графиком. Внимание! График приведён в иллюстративных целях для произвольных значений, не имеющих отношения к заданию. Для выполнения задания используйте данные из прилагаемого файла.

Решение:
Python
from math import dist
base = ''
files = ['27A.txt', '27B.txt']
eps = 1
for t in [0, 1]:
data = [tuple(map(float, line.replace(',', '.').split())) for line in open(base + files[t])]
clusters = []
while data:
clusters.append([data.pop()])
for p in clusters[-1]:
neigh = [pt for pt in data if dist(p, pt) < eps]
clusters[-1] += neigh
for pt in neigh:
data.remove(pt)
clusters = [cl for cl in clusters if len(cl) > 1]
#print(len(clusters), [len(cl) for cl in clusters])
if t == 0:
centers = []
for cl in clusters:
dmin = 10 ** 100
c = None
for p in cl:
d = sum(dist(p, pt) for pt in cl)
if d < dmin:
c = p
dmin = d
centers.append(c)
px = sum(x for x, y in centers)
py = sum(y for x, y in centers)
print(int(abs(px * 10_000)), int(abs(py * 10_000)))
else:
q1 = 10**100
q2 = 0
for i in range(2):
for j in range(i+1, 3):
q1 = min(q1, min(dist(p1, p2) for p1 in clusters[i] for p2 in clusters[j]))
q2 = max(q2, max(dist(p1, p2) for p1 in clusters[i] for p2 in clusters[j]))
print(int(q1 * 10_000), int(q2 * 10_000))
Ответ:
\(92256 \,\, 258611\)
\(33863 \,\, 170816\)
Эффективно готовьтесь к ЕГЭ по информатике с новым тренажёром, эмулирующем работу станции КЕГЭ, которая используется на реальном экзамене