stringtranslate.com

Триангуляция минимального веса

Три возможных триангуляции одного и того же многоугольника. Центральная триангуляция имеет меньший вес (сумму периметров).

В вычислительной геометрии и информатике задача триангуляции с минимальным весом — это задача нахождения триангуляции с минимальной общей длиной ребра . [1] То есть входной многоугольник или выпуклая оболочка входного набора точек должны быть разделены на треугольники, которые сходятся ребром к ребру и вершиной к вершине таким образом, чтобы минимизировать сумму периметров треугольников. Задача является NP-трудной для входных наборов точек, но может быть аппроксимирована с любой желаемой степенью точности. Для входных многоугольников она может быть решена точно за полиномиальное время. Триангуляция с минимальным весом также иногда называется оптимальной триангуляцией .

История

Проблема триангуляции минимального веса множества точек была поставлена ​​Дюппе и Готтшалком (1970), которые предложили ее применение для построения триангулированных нерегулярных сетевых моделей контуров суши и использовали жадную эвристику для ее аппроксимации. Шамос и Хои (1975) предположили, что триангуляция минимального веса всегда совпадает с триангуляцией Делоне , но это было быстро опровергнуто Ллойдом (1977), и действительно Киркпатрик (1980) показал, что веса двух триангуляций могут отличаться на линейный множитель. [2]

Задача триангуляции минимального веса стала печально известной, когда Гари и Джонсон (1979) включили ее в список открытых задач в своей книге о NP-полноте , и многие последующие авторы опубликовали частичные результаты по ней. Наконец, Мульцер и Роте (2008) показали, что она NP-трудна, а Реми и Стегер (2009) показали, что точные приближения к ней могут быть эффективно построены.

Сложность

Вес триангуляции множества точек на евклидовой плоскости определяется как сумма длин ее ребер. Ее вариант решения — задача определения, существует ли триангуляция с весом, меньшим заданного веса; было доказано, что она NP-трудна Mulzer & Rote (2008). Их доказательство основано на редукции от PLANAR-1-IN-3-SAT, особого случая проблемы выполнимости булевых уравнений , в которой 3-CNF, граф которой является планарным, принимается, когда у нее есть истинностное назначение, которое удовлетворяет ровно одному литералу в каждом предложении . Доказательство использует сложные гаджеты и включает компьютерную помощь для проверки правильного поведения этих гаджетов.

Неизвестно, является ли задача решения триангуляции с минимальным весом NP-полной , поскольку это зависит от известной открытой проблемы, можно ли вычислить сумму радикалов за полиномиальное время. Однако Мульцер и Роте отмечают, что задача является NP-полной, если веса ребер округлены до целых значений.

Хотя это NP-трудно, триангуляция с минимальным весом может быть построена за субэкспоненциальное время с помощью алгоритма динамического программирования , который рассматривает все возможные простые разделители циклов точек внутри триангуляции, рекурсивно находит оптимальную триангуляцию на каждой стороне цикла и выбирает разделитель циклов, приводящий к наименьшему общему весу. Общее время для этого метода составляет . [3]

Приближение

Несколько авторов доказали результаты, связывающие триангуляцию минимального веса с другими триангуляциями с точки зрения коэффициента аппроксимации , худшего отношения общей длины ребра альтернативной триангуляции к общей длине триангуляции минимального веса. В этом ключе известно, что триангуляция Делоне имеет коэффициент аппроксимации , [4] и что жадная триангуляция (триангуляция, образованная путем добавления ребер в порядке от самого короткого к самому длинному, на каждом шаге включая ребро всякий раз, когда оно не пересекает более раннее ребро) имеет коэффициент аппроксимации . [5] Тем не менее, для случайно распределенных наборов точек как триангуляция Делоне, так и жадная триангуляция находятся в пределах логарифмического множителя минимального веса. [6]

Результат о сложности Мульцера и Роте также подразумевает NP-сложность нахождения приближенного решения с относительной ошибкой аппроксимации не более O(1/n 2 ). Таким образом, полностью полиномиальная схема аппроксимации для триангуляции с минимальным весом маловероятна. Однако возможна квазиполиномиальная схема аппроксимации: для любой константы ε 0 решение с коэффициентом аппроксимации 1 + ε может быть найдено за квазиполиномиальное время exp(O((log  n ) 9 ). [7]

Эвристика

Из-за сложности нахождения точных решений триангуляции минимального веса многие авторы изучали эвристики, которые в некоторых случаях могут найти решение, хотя нельзя доказать, что они работают во всех случаях. В частности, большая часть этих исследований была сосредоточена на проблеме нахождения наборов ребер, которые гарантированно принадлежат триангуляции минимального веса. Если подграф триангуляции минимального веса, найденный таким образом, оказывается связным, то оставшееся нетриангулированное пространство можно рассматривать как образующее простой многоугольник, и всю триангуляцию можно найти с помощью динамического программирования для нахождения оптимальной триангуляции этого многоугольного пространства. Тот же подход динамического программирования можно распространить и на случай, когда подграф имеет ограниченное число связных компонентов [8] , и тот же подход нахождения связного графа и последующего применения динамического программирования для заполнения полигональных пробелов, окружающих ребра графа, также использовался как часть эвристик для нахождения триангуляций с низким весом, но не с минимальным весом. [9]

Граф, образованный соединением двух точек, когда они являются ближайшими соседями друг друга, обязательно является подграфом триангуляции с минимальным весом. [10] Однако этот граф взаимного ближайшего соседа является соответствием и, следовательно, никогда не является связным. Смежное направление исследований находит большие подграфы триангуляции с минимальным весом с помощью основанных на окружностях β -скелетов , геометрических графов, образованных включением ребра между двумя точками u и v , когда не существует третьей точки w, образующей угол uwv, больший некоторого параметра θ. Было показано, что для достаточно малых значений θ граф, образованный таким образом, является подграфом триангуляции с минимальным весом. [11] Однако выбор θ, необходимый для того, чтобы гарантировать, что это меньше угла {{{1}}}, для которого β -скелет всегда связен.

Более сложная техника, называемая LMT-скелетом, была предложена Дикерсоном и Монтегю (1996). Она формируется итеративным процессом, в котором поддерживаются два набора ребер, набор ребер, известных как принадлежащие триангуляции с минимальным весом, и набор ребер, которые являются кандидатами на принадлежность к ней. Первоначально набор известных ребер инициализируется выпуклой оболочкой ввода, и все оставшиеся пары вершин образуют ребра-кандидаты. Затем, в каждой итерации процесса построения, ребра-кандидаты удаляются всякий раз, когда нет пары треугольников, образованных оставшимися ребрами, образующими четырехугольник, для которого ребро-кандидат является кратчайшей диагональю, и ребра-кандидаты перемещаются в набор известных ребер, когда нет другого ребра-кандидата, которое их пересекает. LMT-скелет определяется как набор известных ребер, полученных после того, как этот процесс прекращает вносить какие-либо изменения. Он гарантированно является подграфом триангуляции минимального веса, может быть эффективно построен, и в экспериментах на наборах до 200 точек он часто оказывался связанным. [12] Однако было показано, что в среднем для больших наборов точек он имеет линейное число связных компонент. [13]

Другие эвристики, которые были применены к проблеме триангуляции минимального веса, включают генетические алгоритмы [14], ветвь и границу [ 15] и алгоритмы оптимизации муравьиной колонии [16] .

Вариации

Триангуляция полигона минимального веса может быть построена за кубическое время с использованием подхода динамического программирования , о котором независимо сообщили Гилберт (1979) и Клинчек (1980). В этом методе вершины нумеруются последовательно вокруг границы полигона, и для каждой диагонали от вершины i до вершины j , которая лежит внутри полигона, оптимальная триангуляция вычисляется путем рассмотрения всех возможных треугольников ijk внутри полигона, сложения весов оптимальных триангуляций ниже диагоналей ik и jk и выбора вершины k , которая приводит к минимальному общему весу. То есть, если MWT( i , j ) обозначает вес триангуляции минимального веса полигона ниже ребра ij , то общий алгоритм выполняет следующие шаги:

После завершения этой итерации MWT(1, n ) сохранит общий вес триангуляции минимального веса. Сама триангуляция может быть получена путем рекурсивного поиска по этому массиву, начиная с MWT(1, n ), на каждом шаге выбирая треугольник ijk, который приводит к минимальному значению для MWT( i , j ) и рекурсивного поиска MWT( i , k ) и MWT( j , k ).

Аналогичные методы динамического программирования могут быть также адаптированы к входным наборам точек, где все, кроме постоянного числа точек, лежат на выпуклой оболочке входа, [17] и к наборам точек, которые лежат на постоянном числе вложенных выпуклых многоугольников или на постоянном числе линий, никакие две из которых не пересекаются внутри выпуклой оболочки. [18]

Также возможно сформулировать версию задач триангуляции набора точек или полигона, в которой разрешено добавлять точки Штейнера , дополнительные вершины, чтобы уменьшить общую длину ребра результирующих триангуляций. В некоторых случаях результирующая триангуляция может быть короче триангуляции минимального веса на линейный множитель. Можно аппроксимировать триангуляцию Штейнера минимального веса набора точек с точностью до постоянного множителя оптимального, но постоянный множитель в аппроксимации велик. [19]

Связанные проблемы, которые также были изучены, включают построение псевдотриангуляций минимального веса [20] и характеристику графов триангуляций минимального веса. [21]

Примечания

  1. ^ Обзоры проблемы см. в Xu (1998), Levcopoulos (2008) и De Loera, Rambau & Santos (2010).
  2. См. также Манахер и Зобрист (1979).
  3. ^ Лингас (1998).
  4. ^ Киркпатрик (1980). Более слабая граница была дана Манахером и Зобристом (1979).
  5. ^ Семейство примеров, доказывающих, что отношение аппроксимации равно, было дано Левкопулосом (1987), а соответствующая верхняя граница — Левкопулосом и Крзнаричем (1998). Как и в случае отношения аппроксимации для триангуляции Делоне, более слабая граница была также дана Манахером и Зобристом (1979).
  6. ^ Лингас (1986).
  7. ^ Remy & Steger (2009). Более ранние результаты по полиномиальным алгоритмам аппроксимации см. в Plaisted & Hong (1987) (логарифмическое факторное приближение) и Levcopoulos & Krznaric (1998) (постоянное факторное приближение).
  8. ^ Ченг, Голин и Цанг (1995).
  9. ^ Лингас (1987); Хит и Пеммараджу (1994).
  10. Ян, Сюй и Ю (1994).
  11. ^ Кейл (1994); Ченг, Голин и Цанг (1995); Ченг и Сюй (2001); Ху (2010).
  12. ^ Дикерсон и Монтегю (1996); Ченг, Като и Сугай (1996); Бейрути и Сноейинк (1998); Айххольцер, Ауренхаммер и Хайнц (1999).
  13. ^ Бозе, Деврой и Эванс (1996).
  14. ^ Цинь, Ван и Гун (1997); Кэпп и Джульстрем (1998).
  15. ^ Киода и др. (1997).
  16. ^ Джахани, Бигэм и Аскари (2010).
  17. ^ Хоффманн и Окамото (2004); Грантсон, Боргельт и Левкопулос (2005); Кнауэр и Спилнер (2006).
  18. ^ Ананьосту и Корней (1993); Мейер и Раппапорт (1992).
  19. ^ Эппштейн (1994).
  20. ^ Гудмундссон и Левкопулос (2007); Айххольцер и др. (2009).
  21. ^ Ленхарт и Лиотта (2002).

Ссылки

Внешние ссылки