Triangulacja minimalnej wagi

Trzy możliwe triangulacje tego samego wielokąta. Centralna triangulacja ma mniejszą wagę (suma obwodów).

W geometrii obliczeniowej i informatyce problem triangulacji minimalnej wagi to problem znalezienia triangulacji o minimalnej całkowitej długości krawędzi . Oznacza to, że wielokąt wejściowy lub wypukła powłoka zbioru punktów wejściowych musi być podzielona na trójkąty, które stykają się krawędzią do krawędzi i wierzchołkiem do wierzchołka, w taki sposób, aby zminimalizować sumę obwodów trójkątów. Problem jest NP-trudny dla danych wejściowych zestawu punktów, ale można je przybliżyć z dowolnym pożądanym stopniem dokładności. W przypadku wejść wielokątnych można go rozwiązać dokładnie w czasie wielomianowym. Minimalna triangulacja wagowa była czasami nazywana triangulacją optymalną .

Historia

Problem triangulacji minimalnej wagi zbioru punktowego został postawiony przez Düppe i Gottschalk (1970) , którzy zasugerowali jego zastosowanie do konstrukcji trójkątnych nieregularnych modeli sieciowych obrysów terenu i zastosowali zachłanną heurystykę do jego przybliżenia. Shamos i Hoey (1975) przypuszczali, że triangulacja minimalnej wagi zawsze pokrywała się z triangulacją Delaunaya , ale zostało to szybko obalone przez Lloyda (1977) i rzeczywiście Kirkpatricka (1980) wykazało, że wagi dwóch triangulacji mogą różnić się o czynnik liniowy.

Problem triangulacji minimalnej wagi stał się znany, gdy Garey i Johnson (1979) umieścili go na liście otwartych problemów w swojej książce o NP-zupełności , a wielu późniejszych autorów opublikowało częściowe wyniki na jego temat. Wreszcie, Mulzer i Rote (2008) wykazali, że jest to NP-trudne, a Remy i Steger (2009) wykazali, że dokładne przybliżenia mogą być konstruowane wydajnie.

Złożoność

Wagę triangulacji zbioru punktów na płaszczyźnie euklidesowej definiuje się jako sumę długości jej krawędzi. Jego wariant decyzyjny polega na rozstrzygnięciu, czy istnieje triangulacja wag mniejszych od danej wagi; Mulzer i Rote (2008) udowodnili, że jest NP -trudny . Ich dowodem jest redukcja z PLANAR-1-IN-3-SAT, szczególny przypadek boolowskiego problemu spełnialności , w którym 3-CNF , którego wykres jest płaski jest akceptowane, gdy ma przypisanie prawdy, które spełnia dokładnie jeden literał w każdej klauzuli . Dowód wykorzystuje złożone gadżety i wymaga pomocy komputera w celu sprawdzenia prawidłowego działania tych gadżetów.

Nie wiadomo, czy problem decyzyjny triangulacji o minimalnej wadze jest NP-zupełny , ponieważ zależy to od znanego problemu otwartego, czy sumę rodników można obliczyć w czasie wielomianowym. Jednak Mulzer i Rote zauważają, że problem jest NP-zupełny, jeśli wagi krawędzi są zaokrąglone do wartości całkowitych.

być skonstruowana w czasie subwykładniczym za pomocą , który uwzględnia wszystkie możliwe proste separatory cykli punktów w triangulacji , rekurencyjnie znajduje optymalną triangulację po każdej stronie cyklu i wybiera separator cykli prowadzący do najmniejszej całkowitej masy. Całkowity czas dla tej metody wynosi .

Przybliżenie

Kilku autorów udowodniło wyniki odnoszące triangulację minimalnej wagi do innych triangulacji pod względem współczynnika aproksymacji , najgorszego przypadku stosunku całkowitej długości krawędzi alternatywnej triangulacji do całkowitej długości triangulacji minimalnej wagi. wiadomo, że triangulacja Delaunaya ma współczynnik aproksymacji a triangulacja zachłanna (triangulacja utworzona przez dodanie krawędzi w kolejności od najkrótszej do najdłuższej każdy krok, w tym krawędź, gdy nie przecina wcześniejszej krawędzi) ma współczynnik aproksymacji równy . Niemniej jednak, w przypadku losowo rozłożonych zbiorów punktów, zarówno triangulacja Delaunaya, jak i zachłanna mieszczą się w logarytmicznym współczynniku minimalnej wagi.

Wynik twardości Mulzera i Rote implikuje również NP-twardość znalezienia przybliżonego rozwiązania ze względnym błędem przybliżenia co najwyżej O(1/n 2 ). Zatem w pełni wielomianowy schemat aproksymacji dla triangulacji minimalnej wagi jest mało prawdopodobny. Możliwy jest jednak quasi-wielomianowy schemat aproksymacji: dla dowolnej stałej ε 0 rozwiązanie o współczynniku aproksymacji 1 + ε można znaleźć w quasi-wielomianowym czasie exp(O((log n ) 9 ).

Heurystyka

Ze względu na trudność w znalezieniu dokładnych rozwiązań triangulacji minimalnej wagi, wielu autorów badało heurystyki, które w niektórych przypadkach mogą znaleźć rozwiązanie, chociaż nie można udowodnić, że działają we wszystkich przypadkach. W szczególności wiele z tych badań koncentrowało się na problemie znajdowania zestawów krawędzi, które z pewnością należą do triangulacji o minimalnej masie. Jeśli znaleziony w ten sposób podwykres triangulacji o minimalnej wadze okaże się spójny, to pozostałą nietriangulowaną przestrzeń można postrzegać jako tworzącą prosty wielokąt, a całą triangulację można znaleźć za pomocą programowania dynamicznego znaleźć optymalną triangulację tej wielokątnej przestrzeni. To samo podejście do programowania dynamicznego można również rozszerzyć na przypadek, w którym podgraf ma ograniczoną liczbę połączonych elementów, i to samo podejście polegające na znalezieniu połączonego wykresu, a następnie zastosowaniu programowania dynamicznego w celu wypełnienia wielokątnych luk otaczających krawędzie grafu jako część heurystyki do znajdowania triangulacji o małej wadze, ale nie o minimalnej wadze.

Wykres utworzony przez połączenie dwóch punktów, gdy są one najbliższymi sąsiadami, jest z konieczności podgrafem triangulacji minimalnej wagi. Jednak ten graf wzajemnego najbliższego sąsiada jest pasujący , a zatem nigdy nie jest połączony. Pokrewny kierunek badań znajduje duże podgrafy triangulacji o minimalnej masie przy użyciu β -szkieletów opartych na okręgu , wykresy geometryczne utworzone przez włączenie krawędzi między dwoma punktami u i v , gdy nie istnieje trzeci punkt w tworzący kąt uwv większy niż jakiś parametr θ. Wykazano, że dla wystarczająco małych wartości θ tak utworzony wykres jest podgrafem triangulacji o minimalnej masie. Jednak wybór θ potrzebny, aby zapewnić, że jest mniejszy niż kąt {{{1}}} , dla którego β -szkielet jest zawsze spójny.

Bardziej wyrafinowaną technikę zwaną szkieletem LMT zaproponowali Dickerson i Montague (1996) . Jest tworzony w procesie iteracyjnym, w którym utrzymywane są dwa zestawy krawędzi, zestaw krawędzi, o których wiadomo, że należą do triangulacji o minimalnej wadze, oraz zestaw krawędzi, które są kandydatami do przynależności do niej. Początkowo zestaw znanych krawędzi jest inicjalizowany do otoczki wypukłej wejścia, a wszystkie pozostałe pary wierzchołków tworzą krawędzie kandydujące. Następnie w każdej iteracji procesu konstrukcyjnego krawędzie kandydujące są usuwane, gdy nie ma pary trójkątów utworzonych przez pozostałe krawędzie tworzące czworobok, dla którego kandydująca krawędź jest najkrótszą przekątną, a kandydujące krawędzie są przesuwane do zbioru znanych krawędzi gdy nie ma innej kandydującej krawędzi, która je przecina. Szkielet LMT jest zdefiniowany jako zestaw znanych krawędzi powstałych po tym, jak ten proces przestanie wprowadzać więcej zmian. Gwarantuje, że będzie to podwykres triangulacji o minimalnej wadze, może być konstruowany wydajnie, aw eksperymentach na zbiorach do 200 punktów często był połączony. Jednak wykazano, że średnio dla dużych zbiorów punktowych ma on liniową liczbę połączonych elementów.

Inne heurystyki, które zostały zastosowane do problemu triangulacji minimalnej wagi, obejmują algorytmy genetyczne rozgałęzione i związane oraz algorytmy optymalizacji kolonii mrówek .

Wariacje

Triangulację wielokątów o minimalnej wadze można skonstruować w czasie sześciennym przy użyciu podejścia do programowania dynamicznego , opisanego niezależnie przez Gilberta (1979) i Klincska (1980) . W tej metodzie wierzchołki są numerowane kolejno wokół granicy wielokąta, a dla każdej przekątnej od wierzchołka i do wierzchołka j , która leży w wielokącie, optymalna triangulacja jest obliczana poprzez uwzględnienie wszystkich możliwych trójkątów ijk w wielokącie, dodanie wag optymalnych triangulacji poniżej przekątnych ik i jk oraz wybierając wierzchołek k , który prowadzi do minimalnej masy całkowitej. Oznacza to, że jeśli MWT( i , j ) oznacza wagę triangulacji minimalnej wagi wielokąta poniżej krawędzi ij , to ogólny algorytm wykonuje następujące kroki:

  • Dla każdej możliwej wartości i od n − 1 do 1 wykonaj:
    • Dla każdej możliwej wartości j , od i + 1 do n , wykonaj:
      • Jeśli ij jest krawędzią wielokąta, ustaw MWT( i , j ) = długość ( ij )
      • W przeciwnym razie, jeśli ij nie jest wewnętrzną przekątną wielokąta, ustaw MWT( i , j ) = +∞
      • Inaczej ustaw MWT( i , j ) = długość ( ij ) + min i < k < j MWT( i , k ) + MWT( k,j )

Po zakończeniu tej iteracji MWT(1, n ) zapisze całkowitą wagę triangulacji minimalnej wagi. Samą triangulację można uzyskać rekurencyjnie przeszukując tę ​​tablicę, zaczynając od MWT(1, n ), na każdym kroku wybierając trójkąt ijk prowadzący do minimalnej wartości dla MWT( i , j ) i rekurencyjnie przeszukując MWT( i , k ) i MWT( j , k ).

Podobne metody programowania dynamicznego można również dostosować do danych wejściowych zbioru punktów, w których wszystkie oprócz stałej liczby punktów leżą na wypukłym korpusie danych wejściowych , oraz do zbiorów punktów, które leżą na stałej liczbie zagnieżdżonych wielokątów wypukłych lub na stałej liczbie linii żadne dwa z nich nie krzyżują się w wypukłej otoczce.

Możliwe jest również sformułowanie wersji problemów ze zbiorem punktów lub triangulacją wielokątów, w których można dodawać punkty Steinera , dodatkowe wierzchołki, w celu zmniejszenia całkowitej długości krawędzi wynikowych triangulacji. W niektórych przypadkach wynikowa triangulacja może być krótsza niż triangulacja minimalnej wagi nawet o współczynnik liniowy. Możliwe jest przybliżenie minimalnej wagi triangulacji Steinera punktu ustawionego w ramach stałego współczynnika optymalnego, ale stały współczynnik przybliżenia jest duży.

Powiązane problemy, które również zostały zbadane, obejmują konstrukcję pseudotriangulacji o minimalnej wadze i charakterystykę wykresów triangulacji o minimalnej wadze.

Notatki

Linki zewnętrzne