Algorithmus von Christofides
# Algorithmus von Christofides

Der Algorithmus von Christofides oder der Algorithmus von Christofides und Serdyukov ist ein Algorithmus, der zur Approximation des metrischen Problem des Handlungsreisenden dient. Er wurde 1976 unabhängig von Nicos Christofides und Anatoliy I. Serdyukov entdeckt und war lange Zeit die beste Approximation des Problems für euklidische Graphen. 1996 stellten Arora und Mitchell für diese jedoch einen besseren Approximationsalgorithmus vor.

Formal geht man ähnlich wie bei der Minimum-Spanning-Tree-Heuristik vor:

* Erzeuge einen minimalen aufspannenden Baum $T$ für den zugrunde liegenden Graphen $G=\left(V,E\right)$ mit Kantengewichten.
* Suche ein (bezüglich Kantengewicht) minimales perfektes Matching $M$ im Graphen zwischen den Knoten, die ungeraden Grad in dem gerade erzeugten Baum $T$ besitzen.
* Füge diese Kanten zu $T$ hinzu. Dabei können Multikanten auftreten. Der entstehende Graph $T\cup M$ ist dann eulersch.
* Konstruiere eine Eulertour in $T\cup M$.
* Konstruiere einen Hamiltonkreis in $T\cup M$. Wähle dazu einen beliebigen Startknoten und gehe die Eulertour ab. Ersetze dabei die bereits besuchten Knoten durch direkte Verbindungen (bzw. Abkürzungen) zum nächsten noch nicht besuchten Knoten.

## Gütegarantie

Es lässt sich zeigen, dass die Christofides-Heuristik eine $1,5$-Approximation ist. Das heißt, die so entstandene Rundreise ist maximal um die Hälfte länger als die optimale Tour. Der Beweis beruht dabei auf einer wiederholten Anwendung der Dreiecksungleichung.

* Die Summe der Kantengewichte im Minimum-Spanning-Tree (MST) ist sowieso kleiner gleich der optimalen Lösung, da jede Lösung des Traveling Salesman Problem (TSP) einen Spannbaum enthält.
* Bezüglich des Matchings gilt folgendes: Sei Fehler beim Parsen (SVG (MathML kann über ein Browser-Plugin aktiviert werden): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „http://localhost:6011/de.wikipedia.org/v1/“:): i_{1},\ldots ,i_{n} die Folge der Knoten vormals ungeraden Grades in der optimalen Lösung; dazwischen liegen irgendwelche anderen Knoten: $a-i_{1}-b-i_{2}-c-\cdots -i_{n}$. Betrachte die beiden Matchings $\left\{\left\{i_{1},i_{2}\right\},\left\{i_{3},i_{4}\right\},\ldots \left\{i_{n-1},i_{n}\right\}\right\}$ sowie $\left\{\left\{i_{2},i_{3}\right\},\left\{i_{4},i_{5}\right\},\ldots \left\{i_{n},i_{1}\right\}\right\}$. Dann gilt aufgrund der Dreiecksungleichung, dass $c(i_{1},i_{2})\leq c(i_{1},b)+c(b,i_{2}),c(i_{2},i_{3})\leq c(i_{2},c)+c(c,i_{3}),\ldots$ Also sind die Gesamtkosten der optimalen Lösung größer gleich derer zweier beliebiger Matchings, insbesondere also zwei Mal des minimalen Matchings. Dann ist ein minimales Matching auch nur maximal halb so groß wie die optimale Lösung. So lässt sich die Summe der Kantengewichte entlang der Eulertour in $T\cup M$ (d. h. die Summe der Gewichte aller Kanten in $T\cup M$) nach oben hin abschätzen.
* Schließlich lässt sich die Summe der Kanten in dem aus der Eulertour erzeugten Hamiltonkreis durch erneutes Anwenden der Dreiecksungleichung nach oben hin durch die Summe der Kanten in der Eulertour abschätzen (denn die Direktkanten können nicht länger sein als die Verbindung über einen schon früher besuchten Knoten), also transitiv durch das 1,5-Fache der optimalen Lösung.

## Beispiel

* Ausgangslage: metrischer Graph $G=\left(V,E\right)$ mit Kantengewichten
* Minimalen Spannbaum $T$ berechnen.
* Die Menge der Knoten mit ungeradem Grad im Spannbaum bestimmen ($V'$).
* $G$ auf die Knoten aus Fehler beim Parsen (SVG (MathML kann über ein Browser-Plugin aktiviert werden): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „http://localhost:6011/de.wikipedia.org/v1/“:): V'
* reduzieren ($G|_{V'}$).
* Matching $M$ mit minimalem Gewicht auf $G|_{V'}$ bestimmen.
* Matching und Spannbaum vereinigen ($T\cup M$). Dieser Schritt sorgt dafür, dass Knoten mit vormals ungeradem Grad nun einen geraden Grad aufweisen. Dies ist eine notwendige Bedingung für die Berechnung der Euler-Tour im nächsten Schritt.
* Euler-Tour auf $T\cup M$ berechnen (A-B-C-A-D-E-A).
* Wiederholt vorkommende Knoten entfernen und durch Direktverbindung ersetzen (A-B-C-D-E-A). In metrischen Graphen führt dies nicht zu einer längeren Strecke.
* Diese Tour ist die Ausgabe des Algorithmus.
