Eine Union-Find-Datenstruktur verwaltet die Partition einer Menge. Der abstrakte Datentyp wird durch die Menge der drei Operationen {Union, Find, Make-Set} gebildet. Union-Find-Strukturen dienen zur Verwaltung von Zerlegungen in disjunkte Mengen. Dabei bekommt jede Menge der Zerlegung ein kanonisches Element zugeordnet, welches als Name der Menge dient. Union vereinigt zwei solche Mengen, Find(x) bestimmt das kanonische Element derjenigen Menge, die x enthält, und Make-Set erzeugt eine Elementarmenge {x} mit dem kanonischen Element x.

Definition

Eine endliche Menge ${\displaystyle S}$ sei in die disjunkten Klassen ${\displaystyle X_{i}}$ partitioniert:

${\displaystyle S=X_{0}\uplus X_{1}\uplus X_{2}\uplus \ldots \uplus X_{k}}$ mit ${\displaystyle X_{i}\cap X_{j}=\varnothing \quad \forall i,j\in \lbrace 0,1,\ldots ,k\rbrace ,i\neq j}$.

Zu jeder Klasse ${\displaystyle X_{i}}$ wird ein Repräsentant ${\displaystyle r_{i}\in X_{i}}$ ausgewählt. Die zugehörige Union-Find-Struktur unterstützt die folgenden Operationen effizient:

• Init(${\displaystyle S}$): Initialisiert die Struktur und bildet für jedes ${\displaystyle x\in S}$ eine eigene Klasse mit ${\displaystyle x}$ als Repräsentant.
• Union(${\displaystyle r}$, ${\displaystyle s}$): Vereinigt die beiden Klassen, die zu den beiden Repräsentanten ${\displaystyle r}$ und ${\displaystyle s}$ gehören, und bestimmt ${\displaystyle r}$ zum neuen Repräsentanten der neuen Klasse.
• Find(${\displaystyle x}$): Bestimmt zu ${\displaystyle x\in S}$ den eindeutigen Repräsentanten, zu dessen Klasse ${\displaystyle x}$ gehört.

Implementierung

Eine triviale Implementierung speichert die Zugehörigkeiten zwischen den Elementen aus ${\displaystyle S}$ und den Repräsentanten ${\displaystyle r_{i}}$ in einem Array. Für kürzere Laufzeiten werden jedoch in der Praxis Mengen von Bäumen verwendet. Dabei werden die Repräsentanten in den Wurzeln der Bäume gespeichert, die anderen Elemente der jeweiligen Klasse in den Knoten darunter.

• Union(${\displaystyle x}$, ${\displaystyle y}$): Hängt die Wurzel des niedrigeren Baumes als neues Kind unter die Wurzel des höheren Baumes (gewichtete Vereinigung). Falls nun ${\displaystyle x}$ Kind von ${\displaystyle y}$ ist, werden ${\displaystyle x}$ und ${\displaystyle y}$ vertauscht.
• Find(${\displaystyle x}$): Wandert vom Knoten ${\displaystyle x}$ aus den Pfad innerhalb des Baumes nach oben bis zur Wurzel und gibt diese als Ergebnis zurück.

Heuristiken

Um die Operationen Find und Union zu beschleunigen, gibt es die Heuristiken union by size, union by rank (Rang- und Größenheuristik) und Pfadkompression, wobei sich union by size und union by rank gegenseitig ausschließen und union by rank oft in Kombination mit Pfadkompression verwendet wird.

function Union(x, y)
{
    // Knoten durch die Wurzeln des Baums ersetzen
    x = Find(x);
    y = Find(y);
    if (x = y) // Wenn x und y schon zur selben Menge gehören, wird die Funktion verlassen
    {
        return;
    }
    if (x.rank < y.rank) // Wenn der Rang von x kleiner als der Rang von y ist, wird y zur neuen Wurzel
    {
        x.parent = y;
    }
    else if (x.rank > y.rank) // Wenn der Rang von x größer als der Rang von y ist, wird x zur neuen Wurzel
    {
        y.parent = x;
    }
    else // Wenn die Ränge gleich sind, wird y zur neuen Wurzel und der Rang von y inkrementiert
    {
        x.parent = y;
        y.rank = y.rank + 1;
    }
}

function Find(x)
{
    if (x.parent ≠ x)
    {
        x.parent = Find(x.parent);
        return x.parent;
    }
    else
    {
        return x;
    }
}

Aufteilungsmethode (splitting)

Während des Durchlaufes lässt man jeden Knoten auf seinen bisherigen Großvater zeigen, falls vorhanden. Damit wird ein durchlaufender Pfad in zwei der halben Länge zerlegt.

Laufzeiten

Union-Find-Datenstrukturen ermöglichen die Ausführung der obigen Operationen mit den folgenden Laufzeiten (${\displaystyle n=\left|S\right|}$):

Triviale Implementierung mit einem Array A (A[i] = x: Knoten i liegt in der Menge x), worst-case: Find ${\displaystyle O(1)}$, Union: ${\displaystyle O(n)}$

Implementierung mit Bäumen:

• ohne Heuristiken: Find ${\displaystyle O(n)}$, Union ${\displaystyle O(n)}$
• mit Union-By-Size, worst-case: Find ${\displaystyle O(\log(n))}$, Union ${\displaystyle O(1)}$
• mit Union-By-Size, Pfadkompression, worst-case: Find ${\displaystyle O(\log(n))}$, Union ${\displaystyle O(1)}$
• mit Union-By-Size, Pfadkompression, Folge von m Find- und n-1 Union-Operationen: ${\displaystyle O\left(n+m\cdot \alpha (n)\right)}$

Dabei ist ${\displaystyle \alpha (n)}$ die Umkehrfunktion der Ackermannfunktion ${\displaystyle A(n,n)}$. Sie wächst sehr langsam und ist für alle „praktischen“ Eingaben (${\displaystyle n<10^{80}}$) kleiner als 5. Im Jahr 1977 hat Robert Tarjan gezeigt, dass für eine Sequenz m Find- und n-1 Union-Operationen jede Datenstruktur im Pointer-Maschinen Modell eine Laufzeit von ${\displaystyle \Omega \left(n+m\cdot \alpha (n)\right)}$ benötigt. Fredman und Saks haben 1989 gezeigt, dass diese untere Schranke auch im allgemeineren Cell-Probe Modell gilt. Daraus folgt, dass die Datenstruktur mit Union-By-Size und Pfadkompression eine asymptotisch optimale Laufzeit besitzt, und es keine Union-Find-Datenstruktur geben kann, die sowohl Find als auch Union amortisiert in O(1) ausführt.

Anwendungen

Union-Find-Strukturen eignen sich gut, um Zusammenhangskomponenten von Graphen zu verwalten.

Union-Find-Strukturen können auch genutzt werden, um effizient herauszufinden, ob ein Graph Zyklen enthält.

Das folgende Computerprogramm in der Programmiersprache C++ zeigt die Implementierung eines ungerichteten Graphen mit einem Array von Kanten. Der ungerichtete Graph wird als Datentyp UndirectedGraph deklariert. Bei der Ausführung des Programms wird die Methode main verwendet, das auf der Konsole ausgibt, ob der gegebene Graph Zyklen enthält. Die Funktion unionSubsets verwendet union by rank, um zwei Teilmengen von Kanten des Graphen zu vereinigen.

#include <iostream>
using namespace std;

// Deklariert den Datentyp für die Kanten des Graphen
struct Edge
{
    int startIndex, endIndex;
};

// Deklariert den Datentyp für den ungerichteten Graphen
struct UndirectedGraph
{
    int numberOfVertices, numberOfEdges;
    Edge* edges; // Pointer auf das Array für die Kanten
};

// Deklariert den Datentyp für Teilmengen der Kantenmenge des ungerichteten Graphen
struct subset
{
    int parent;
    int rank;
};

// Diese Funktion erzeugt einen ungerichteten Graphen
struct UndirectedGraph* createGraph(int numberOfVertices, int numberOfEdges)
{
    struct UndirectedGraph* undirectedGraph = new UndirectedGraph();
    undirectedGraph->numberOfVertices = numberOfVertices;
    undirectedGraph->numberOfEdges = numberOfEdges;
    undirectedGraph->edges = new Edge[numberOfEdges];
    return undirectedGraph;
}

// Diese rekursive Funktion gibt den Index der Wurzel der Teilmenge mit dem Index i zurück
int find(subset subsets[], int i)
{
    // Setzt Index der Wurzel auf den Index der Wurzel der Teilmenge mit dem Index i
    if (subsets[i].parent != i)
    {
        subsets[i].parent = find(subsets, subsets[i].parent); // Rekursiver Aufruf der Funktion
    }
    return subsets[i].parent;
}

// Diese Methode bildet die Vereinigungsmenge der zwei Teilmengen mit den Indexen index1 und index2
void unionSubsets(subset subsets[], int index1, int index2)
{
    int newIndex1 = find(subsets, index1); // Index der Teilmenge mit dem Index index1
    int newIndex2 = find(subsets, index2); // Index der Teilmenge mit dem Index index2
     // Hängt den Teilbaum mit dem niedrigeren Rang unter die Wurzel des Baums mit dem höheren Rang
    if (subsets[newIndex1].rank < subsets[newIndex2].rank)
    {
        subsets[newIndex1].parent = newIndex2;
    }
    else if (subsets[newIndex1].rank > subsets[newIndex2].rank)
    {
        subsets[newIndex2].parent = newIndex1;
    }
    else // Wenn die Teilbäume denselben Rang haben, wird der Rang des einen Baums erhöht und der andere Baum unter die Wurzel des anderen Baums gehängt
    {
        subsets[newIndex2].parent = newIndex1;
        subsets[newIndex1].rank++;
    }
}

// Diese Methode prüft, ob der Graph Zyklen enthält
bool containsCycle(struct UndirectedGraph* graph)
{
    int numberOfVertices = graph->numberOfVertices;
    int numberOfEdges = graph->numberOfEdges;
    Edge* edges = graph->edges; // Pointer auf das Array für die Kanten
    subset* subsets = new subset[numberOfVertices]; // Deklariert ein Array für die Teilmengen der Kantenmenge
    for (int i = 0; i < numberOfVertices; i++) // for-Schleife, die Teilmengen mit einzelnen Kanten erzeugt
    {
        subsets[i].parent = i;
        subsets[i].rank = 0;
    }
    for (int i = 0; i < numberOfEdges; i++) // foreach-Schleife, die alle Kanten durchläuft
    {
        Edge nextEdge = edges[i++]; // Weist die verbleibende Kante mit dem kleinsten Kantengewicht zu und erhöht den Kantenindex für die nächste Iteration um 1
        int index1 = find(subsets, nextEdge.startIndex); // Index der Wurzel der Teilmenge mit dem Index nextEdge.startIndex
        int index2 = find(subsets, nextEdge.endIndex); // Index der Wurzel der Teilmenge mit dem Index nextEdge.endIndex
        if (index1 == index2) // Wenn die Indexe der Wurzeln, also auch die Mengen übereinstimmen, enthält der Graph einen Zyklus
        {
            return true;
        }
        unionSubsets(subsets, index1, index2);
    }
    return false;
}

// Hauptfunktion die das Programm ausführt
int main()
{
    int numberOfVertices = 3;
    int numberOfEdges = 3;

    // Erzeugt den ungerichteten Graphen mit den gegebenen Kanten
    struct UndirectedGraph* undirectedGraph = createGraph(numberOfVertices, numberOfEdges);
    // Initialisiert das Array mit 3 Kanten
    Edge* edges = undirectedGraph->edges;
    edges[0].startIndex = 0;
    edges[0].endIndex = 1;
    edges[1].startIndex = 1;
    edges[1].endIndex = 2;
    edges[2].startIndex = 0;
    edges[2].endIndex = 2;
    if (containsCycle(undirectedGraph)) // Aufruf der Methode
    {
        cout << "Der Graph enthält Zyklen." << endl; // Ausgabe auf der Konsole
    }
    else
    {
        cout << "Der Graph enthält keine Zyklen." << endl; // Ausgabe auf der Konsole
    }
}

Der Algorithmus von Kruskal erfordert beispielsweise eine solche Datenstruktur, um effizient implementiert werden zu können.

Weblinks

• Demonstration von Union-Find

Einzelnachweise

1. ↑ Disjoint Set | Union-Find Algorithm – Union by rank and path compression. In: algorithms.tutorialhorizon.com (englisch)
2. ↑ Disjoint-set data structure (Memento des Originals vom 2. Dezember 2016 im Internet Archive)  Info: Der Archivlink wurde automatisch eingesetzt und noch nicht geprüft. Bitte prüfe Original- und Archivlink gemäß Anleitung und entferne dann diesen Hinweis.@1@2Vorlage:Webachiv/IABot/www.mathblog.dk. In: mathblog.dk
3. ↑ Robert E. Tarjan: A class of algorithms which require nonlinear time to maintain disjoint sets. In: Journal of Computer and System Sciences. 18. Jahrgang, Nr. 2, 1979, S. 110–127, doi:10.1016/0022-0000(79)90042-4 (englisch). 
4. ↑ M. Fredman, M. Saks, Saks: The cell probe complexity of dynamic data structures. Proceedings of the Twenty-First Annual ACM Symposium on Theory of Computing. Mai 1989, S. 345–354, doi:10.1145/73007.73040 (englisch). 
5. ↑ Union-Find Algorithm | Set 2 (Union By Rank and Path Compression). In: GeeksforGeeks (englisch)