Фрагмент звёздного неба спроецирован на плоскость с декартовой системой координат. Учёный решил провести кластеризацию полученных точек, являющихся изображениями звёзд, то есть разбить их множество на N непересекающихся непустых подмножеств (кластеров), таких что точки каждого подмножества лежат внутри прямоугольника со сторонами длиной H и W, причём эти прямоугольники между собой не пересекаются. Стороны прямоугольников не обязательно параллельны координатным осям.
Гарантируется, что такое разбиение существует и единственно для заданных размеров прямоугольников.
Будем называть антицентром кластера точку этого кластера, сумма расстояний от которой до всех остальных точек кластера максимальна. Для каждого кластера гарантируется единственность его антицентра. Расстояние между двумя точками А (х1, y1) и B (х2, y2) на плоскости вычисляется по формуле:
В файле А хранятся координаты точек двух кластеров, где Н = 8, W = 4 для каждого кластера. В каждой строке записана информация о расположении на карте одной звезды: сначала координата х, затем координата у. Значения даны в условных единицах. Известно, что количество точек не превышает 1000.
В файле Б хранятся координаты точек трёх кластеров, где Н = 6, W = 7 для каждого кластера. Известно, что количество точек не превышает 10 000. Структура хранения информации о звёздах в файле Б аналогична структуре в файле А.
Известно, что в файле Б имеются координаты ровно трёх «лишних» точек, представляющих аномалии, которые возникли в результате помех при передаче данных. Эти три точки не относятся ни к одному из кластеров, их учитывать не нужно.
Для файла А определите координаты антицентра каждого кластера, затем вычислите два числа: P1 — сумма абсциссы и ординаты антицентра кластера с наименьшим количеством точек, и P2 — сумма абсциссы и ординаты антицентра кластера с наибольшим количеством точек. Гарантируется, что во всех кластерах количество точек различно.
Для файла Б определите координаты антицентра каждого кластера, затем вычислите два числа: Qx — абсциссу наиболее отдалённого антицентра кластера от начала координат, и Qy — ординату ближайшего антицентра кластера к началу координат.
В ответе запишите четыре числа: в первой строке — сначала целую часть абсолютной величины произведения P1 × 10 000, затем целую часть абсолютной величины произведения P2 × 10 000; во второй строке — сначала целую часть произведения Qx × 10 000, затем целую часть произведения Qy × 10 000.
Возможные данные одного из файлов проиллюстрированы графиком.
Ответ:
Построим диаграммы для файла А и Б. Для этого воспользуемся табличным редактором.
Диаграмма для файла А:
Диаграмма для файла Б:
Приведём решение на языке Python для файла А.
from math import *
f = open('27.txt')
claster = [[] for i in range(2)]
for s in f:
x,y= map(float, s.split())
if y < 90:
claster[0].append((x,y))
else:
claster[1].append((x,y))
def centroid(claster):
rad = []
for i in claster:
rad.append((sum(dist(i,j) for j in claster),i))
return max(rad)[1]
len_claster = sorted(claster, key = len )
centr = [centroid(i) for i in len_claster]
P1 = sum(centr[0]) *10000
P2 = sum(centr[1]) *10000
print(int(abs(P1)), int(abs(P2)))
Приведём решение на языке Python для файла В.
from math import *
f = open('27.txt')
claster = [[] for i in range(3)]
for s in f:
x,y= map(float, s.split())
if x > 16 and 32 > y > 24:
claster[0].append((x,y))
elif x > 16 and 41 > y > 33:
claster[1].append((x,y))
elif x > 16 and 50 > y > 42:
claster[2].append((x,y))
def centroid(claster):
rad = []
for i in claster:
rad.append((sum(dist(i,j) for j in claster),i))
return max(rad)[1]
centr = [centroid(i) for i in claster]
Qx = max([(dist(i, (0,0)), i) for i in centr])[1][0]
Qy = min([(dist(i, (0,0)), i) for i in centr])[1][1]
print(int(abs(Qx*10000)), int(abs(Qy*10000)))
Ответ: 1126711 1517181 213883 264132.