Участник:Zhu.a-v/Алгоритм Киркпатрика

Материал из Алговики
Перейти к навигации Перейти к поиску

Автор статьи: Журавская Александра

Свойства и структура алгоритма

Алгоритм Киркпатрика решает задачу построения выпуклой оболочки набора точек методом "разделяй и властвуй" [1].

Общее описание алгоритма

Выпуклой оболочкой множества [math]\displaystyle{ X }[/math] называется наименьшее выпуклое множество, содержащее множество [math]\displaystyle{ X }[/math].

Справедлива следующая интерпретация выпуклой оболочки: если вбить в деревянную доску несколько гвоздей, забивая их только до середины ножки, затем накинуть вокруг них петлю веревки и затянуть, веревка опишет выпуклую оболочку для множества точек-гвоздей.

Математическое описание алгоритма

Входные данные: множество точек [math]\displaystyle{ X = {(x_1,y_1),(x_2,y_2),...,(x_n,y_n)} }[/math].

Ясно, что выпуклая оболочка множества точек на плоскости [math]\displaystyle{ X }[/math] --- это выпуклый многоугольник, вершины которого содержатся в [math]\displaystyle{ X }[/math], и все остальные элементы [math]\displaystyle{ X }[/math] лежат внутри этого многоугольника. Значит, задача поиска выпуклой оболочки [math]\displaystyle{ conv X }[/math] сводится к выбору такого набора вершин [math]\displaystyle{ (x_{i_1},y_{i_1}),(x_{i_2},y_{i_2}),...,(x_{i_n},y_{i_n}) }[/math] , которые образуют выпуклый многоугольник, содержащий все точки множества [math]\displaystyle{ X }[/math].

Выходные данные: множество точек [math]\displaystyle{ conv X }[/math] --- вершины выпуклой оболочки для множества [math]\displaystyle{ X }[/math].


Применим рекурсивный подход "разделяй и властвуй" для решения поставленной задачи:

1. Если множество состоит из 1 элемента, его выпуклая оболочка совпадает с самим множеством. Иначе шаг 2.

2. Проведем на плоскости прямую, разделяющую множество на два примерно одинаковых подмножества [math]\displaystyle{ X_l }[/math], [math]\displaystyle{ X_r }[/math]. Выберем для этого следующую прямую: [math]\displaystyle{ x = \frac{\sum_{(a,b) \in T} a}{|T|} }[/math] --- это вертикальная линия, медиана абцисс точек множества, для которого на текущем шаге строится выпуклая оболочка. Для каждого из полученных подмножеств найдем выпуклую оболочку. Перейдем на шаг 3.

3. Построим [math]\displaystyle{ conv(X) }[/math] по [math]\displaystyle{ conv(X_l) }[/math] и [math]\displaystyle{ conv(X_r) }[/math], построенным на предыдущем шаге. Поскольку множества [math]\displaystyle{ X_l }[/math] и [math]\displaystyle{ X_r }[/math] разделены вертикальной линией, их выпуклые оболочки представляют собой непересекающиеся выпуклые многоугольники. Следовательно, для построения общей выпуклой оболочки этих многоугольников достаточно найти их общие касательные - верхнюю и нижнюю, а затем отбросить правую цепочку в многоугольнике [math]\displaystyle{ conv(X_l) }[/math] и левую цепочку в многоугольнике [math]\displaystyle{ conv(X_r) }[/math]. Пусть левый и правый многоугольники представлены списками вершин [math]\displaystyle{ P_l }[/math] и [math]\displaystyle{ P_r }[/math]. Поиск общей касательной начнем с мостика [math]\displaystyle{ (u,v) }[/math], соединяющего самую правую вершину [math]\displaystyle{ u }[/math] в [math]\displaystyle{ P_l }[/math] с самой левой вершиной [math]\displaystyle{ v }[/math] в [math]\displaystyle{ P_r }[/math]. Далее будем последовательно перемещать эти вершины вдоль границ выпуклых многоугольников против часовой стрелки в [math]\displaystyle{ P_l }[/math] и по часовой стрелке в [math]\displaystyle{ P_r }[/math] до тех пор, пока существуют вершины, лежащие слева от вектора [math]\displaystyle{ (u,v) }[/math]. Таким образом мы найдем общую верхнюю касательную. Аналогично, вычислим нижнюю касательную. Но перемещение теперь осуществляется по часовой стрелке в [math]\displaystyle{ P_l }[/math] и против часовой стрелки в [math]\displaystyle{ P_r }[/math]. Поиск нижней касательной завершается, когда мостик окажется в положении, при котором не существует вершин, лежащих справа от вектора [math]\displaystyle{ (u,v) }[/math]. [2]

Вычислительное ядро алгоритма

Макроструктура алгоритма

Схема реализации последовательного алгоритма

Последовательная сложность алгоритма

Информационный граф

Ресурс параллелизма алгоритма

Входные и выходные данные алгоритма

Свойства алгоритма

Программная реализация алгоритма

Особенности реализации последовательного алгоритма

Локальность данных и вычислений

Возможные способы и особенности параллельной реализации алгоритма

Масштабируемость алгоритма и его реализации

Динамические характеристики и эффективность реализации алгоритма

Выводы для классов архитектур

Существующие реализации алгоритма

Литература

<references \>

  1. ↑ Kirkpatrick, David G.; Seidel, Raimund (1986). "The ultimate planar convex hull algorithm". SIAM Journal on Computing. 15 (1): 287–299.
  2. ↑ Местецкий Л.М. Лекции по вычислительной геометрии. Учебное пособие. ВМК МГУ, 2013.