333
правки
Изменения
→Алгоритм
==Алгоритм ==
Сперва приведем алгоритм, корректность которого позднее будет доказана.
Рассмотрим метод <tex> MinDisc </tex>, который будет принимать на вход множество точек, а возвращать будет искомую охватывающую окружность. Идея метода в том, что он генерирует случайную перестановку из точек, которые поступили на вход и в получившимся порядке по одной добавляет их в текущую охватывающую окружность. Как будет доказано позднее, если на каком-то шаге алгоритма текущая точка оказалась вне текущей минимальной окружности, то утверждается, что она лежит на границе новой окружности, которая будет содержать в себе ее и все предыдущие точки. Опишем приведенный выше алгоритм псевдоком:
if pi <tex> \in </tex> Di−1 then Di = Di−1 // point is inside. nothing to do here
else Di = MinDiscWithPoint ({p1,...,pi−1}, pi) // point is lying on the boundary of minimal circle
}
return Dn
Как вы заметили, мы использовали выше метод <tex> MinDiscWithPoint </tex> - это метод, который по множеству точек и еще одной точке возвращают искомую окружность. При этом выделенная точка (которая передается в качестве второго аргумента) ''обязательно'' лежит ''на'' окружности. Эта функция очень похожа на предыдущую, поэтому приведу лишь ее псевдокод:
MinDiscWithPoint(S, q)
Let P - random permutation of S
Let D1 - be the smallest enclosing disc for {p1, q}.
for i = 2 to n {
if pi <tex> \in </tex> Di−1 then Di = Di−1
else Di = MinDiscWith2Points ({p1,...,pi−1}, pi, q)
}
return Dn
