Задание 27. Информатика. ЕГЭ. Поляков-7651
- Просмотры: 120
- Изменено: 24 ноября 2024
(В. Ланская, Р. Ягафаров) Учёный решил провести кластеризацию некоторого множества звёзд по их расположению на карте звёздного неба. 
Кластер звёзд – это набор звёзд (точек) на графике. Каждая звезда обязательно принадлежит только одному из кластеров. Центр кластера, или центроид, – это одна из звёзд на графике, сумма расстояний от которой до всех остальных звёзд кластера минимальна. Расстояние между двумя точками \(A(x_1, \, y_1)\) и \(B(x_2, \, y_2)\) вычисляется по формуле:
$$
d(A, \, B) = \sqrt{(x_1 - x_2)^2 + (y_1 - y_2)^2}
$$
Даны два входных файла (файл A и файл Б). В файле A хранятся данные о звёздах двух кластеров. В каждой строке записана информация о расположении на карте одной звезды: сначала координата \(x\), затем координата \(y\) (в условных единицах). Известно, что количество звёзд не превышает \(1000\). В файле Б хранятся данные о звёздах четырёх кластеров. Известно, что количество звёзд не превышает \(10~000\). Структура хранения информации о звездах в файле Б аналогична файлу А. Возможные данные одного из файлов иллюстрированы графиком.
Для каждого файла определите координаты центра каждого кластера, затем вычислите два числа: \(P_x\) – среднее арифметическое абсцисс центров кластеров, и \(P_y\) – среднее арифметическое ординат центров кластеров. В ответе запишите четыре числа: в первой строке сначала целую часть произведения \(P_x \times 10~000\), затем целую часть произведения \(P_y \times 10~000\) для файла А, во второй строке – аналогичные данные для файла Б.
Решение:
Точечная диаграмма, построенная для файла А в табличном процессоре, имеет вид
Видно, что кластеры можно разделить, например, прямой \(y = -4\)
Точечная диаграмма для файла Б, построенная в табличном процессоре, имеет вид
Левый кластер отделяется от остальных прямой \(y = -6\). Правые можно последовательно разделить прямыми \(y = -6\) и \(y =-2.5\).
Python (3.10+)
def which_cluster(point, task):
x, y = point
match task:
case 0:
return int(y > -4)
case 1:
if x < -6:
return 0
elif y < -6:
return 1
elif y < -2.5:
return 2
else:
return 3
def dist(point1, point2):
x1, y1 = point1
x2, y2 = point2
return ((x1 - x2)**2 + (y1 - y2)**2)**0.5
base = ''
files = ['27A.txt', '27B.txt']
clusters = [2, 4]
for t in [0, 1]:
fd = open(base + files[t])
fd.readline()
data = [[], [], [], []]
for line in fd:
x, y = map(float, line.replace(',', '.').split())
data[which_cluster([x, y], t)].append([x, y])
centers = []
for k in range(clusters[t]):
dm = 10**100
c = [0, 0]
for p in data[k]:
d = sum(dist(p, pt) for pt in data[k])
if d < dm:
dm = d
c = p
centers.append(c)
px = int(sum(c[0] for c in centers) / clusters[t] * 10000)
py = int(sum(c[1] for c in centers) / clusters[t] * 10000)
print(px, py)
Ответ:
\(-95531\) \(-47828\)
\(-38905\) \(-46034\)