6.6. Kruskal’s Algorithm¶
Our next MCST algorithm is commonly referred to as Kruskal’s algorithm. Kruskal’s algorithm is also a simple, greedy algorithm. First partition the set of vertices into \(|\mathbf{V}|\) disjoint sets, each consisting of one vertex. Then process the edges in order of weight. An edge is added to the MCST, and two disjoint sets combined, if the edge connects two vertices in different disjoint sets. This process is repeated until only one disjoint set remains.
The edges can be processed in order of weight by using a min-heap. This is generally faster than sorting the edges first, because in practice we need only visit a small fraction of the edges before completing the MCST. This is an example of finding only a few smallest elements in a list.
The only tricky part to this algorithm is determining if two vertices
belong to the same equivalence class.
Fortunately, the ideal algorithm is available for the purpose —
the UNION/FIND.
Here is an implementation for Kruskal’s algorithm.
Class KruskalElem
is used to store the edges on the min-heap.
# Kruskal's MST algorithm
def Kruskal(G):
A = ParPtrTree(G.nodeCount()) # Equivalence array
E = [None] * G.edgeCount() # Minheap array
edgecnt = 0 # Count of edges
for i in range(G.nodeCount()): # Put edges in the array
for w in G.neighbors(i):
E[edgecnt] = (G.weight(i, w), i, w)
edgecnt += 1
H = MinHeap(E, edgecnt, edgecnt)
numMST = G.nodeCount() # Initially n disjoint classes
while numMST>1: # Combine equivalence classes
temp = H.removemin() # Next cheapest edge
if temp is None:
return # Must have disconnected vertices
_w, v, u = temp
if A.differ(v, u): # If in different classes
A.UNION(v, u) # Combine equiv classes
AddEdgetoMST(v, u) # Add this edge to MST
numMST -= 1 # One less MST
Kruskal’s algorithm is dominated by the time required to
process the edges.
The differ
and UNION
functions are nearly
constant in time if path compression and weighted union is used.
Thus, the total cost of the algorithm is
\(\Theta(|\mathbf{E}| \log |\mathbf{E}|)\) in the worst case,
when nearly all edges must be processed before all the edges of the
spanning tree are found and the algorithm can stop.
More often the edges of the spanning tree are the shorter ones,and
only about \(|\mathbf{V}|\) edges must be processed.
If so, the cost is often close to
\(\Theta(|\mathbf{V}| \log |\mathbf{E}|)\) in the average case.