FibonacciHeaps 2x2
FibonacciHeaps 2x2
F IBONACCI H EAPS operation linked list binary heap binomial heap Fibonacci heap †
† amortized
Theorem. [Fredman-Tarjan 1986] Starting from an empty Fibonacci heap, Theorem. [Fredman-Tarjan 1986] Starting from an empty Fibonacci heap,
any sequence of m INSERT, EXTRACT-MIN, and DECREASE-KEY operations any sequence of m INSERT, EXTRACT-MIN, and DECREASE-KEY operations
involving n INSERT operations takes O(m + n log n) time. involving n INSERT operations takes O(m + n log n) time.
History.
Fibonacci Heaps and Their Uses in Improved Network
Optimization Algorithms
・Ingenious data structure and application of amortized analysis.
MICHAEL L. FREDMAN
・Original motivation: improve Dijkstra's shortest path algorithm
University of California, San Diego, L.a Jolla, California
from O(m log n) to O(m + n log n).
AND
ROBERT ENDRE TARJAN ・Also improved best-known bounds for all-pairs shortest paths,
AT&T Bell Laboratories, Murray HilI, New Jersey
assignment problem, minimum spanning trees.
Abstract. In this paper we develop a new data structure for implementing heaps (priority queues). Our
structure, Fibonacci heaps (abbreviated F-heaps), extends the binomial queues proposed by Vuillemin
and studied further by Brown. F-heaps support arbitrary deletion from an n-item heap in qlogn)
amortized time and all other standard heap operations in o( 1) amortized time. Using F-heaps we are
able to obtain improved running times for several network optimization algorithms. In particular, we
obtain the following worst-case bounds, where n is the number of vertices and m the number of edges
in the problem graph:
( 1) O(n log n + m) for the single-source shortest path problem with nonnegative edge lengths, improved
from O(mlogfmh+2)n);
(2) O(n*log n + nm) for the all-pairs shortest path problem, improved from O(nm lo&,,,+2,n);
(3) O(n*logn + nm) for the assignment problem (weighted bipartite matching), improved from
O(nmlog0dn+2)n);
(4) O(mj3(m, n)) for the minimum spanning tree problem, improved from O(mloglo&,,.+2,n), where
j3(m, n) = min {i 1log% 5 m/n). Note that B(m, n) 5 log*n if m 2 n.
Of these results, the improved bound for minimum spanning trees is the most striking, although all the
results give asymptotic improvements for graphs of appropriate densities.
Categories and Subject Descriptors: E.l [Data]: Data Structures--trees; graphs; F.2.2 [Analysis of
Algorithms and Problem Complexity]: Nonnumerical Algorithms and Problems-computations on
discrete structures; sorting and searching; G.2.2 [Discrete Mathematics]: Graph Theory-graph algo- 3 4
rithms; network problems; trees
General Terms: Algorithms, Theory
Additional Key Words and Phrases: Heap, matching, minimum spanning tree, priority queue, shortest
Fibonacci heaps
Basic idea.
F IBONACCI H EAPS ・Similar to binomial heaps, but less rigid structure.
・Binomial heap: eagerly consolidate trees after each INSERT;
‣ structure implement DECREASE-KEY by repeatedly exchanging node with its parent.
‣ insert
‣ extract the minimum
‣ decrease key
‣ bounding the rank
‣ meld and delete
SECTION 19.1
・Fibonacci heap: lazily defer consolidation until next EXTRACT-MIN;
implement DECREASE-KEY by cutting off node and splicing into root list.
Remark. Height of Fibonacci heap is Θ(n) in worst case, but it doesn't use
sink or swim operations.
heap-ordered tree
min
17 24 23 7 3 root 17 24 23 7 3
30 26 46 30 26 46
18 52 41 18 52 41
min min
17 24 23 7 3 17 24 23 7 3 rank = 3
30 26 46 30 26 46
18 52 41 18 52 41
children are in a
heap H heap H
35 35 circular doubly-linked list
39 44 39 44
9 10
・Add or remove a node from the root list. rank(H) max rank of any node in heap H
・Remove a subtree and merge into root list. trees(H) number of trees in heap H
・Link the root of a one tree to root of another tree. marks(H) number of marked nodes in heap H
30 26 46 30 26 46
18 52 41 18 52 41
heap H heap H
35 35
39 44 39 44
11 12
Fibonacci heap: potential function
Potential function.
F IBONACCI H EAPS
30 26 46
18 52 41
heap H
35
39 44
13
insert 21 insert 21
21
min min
17 24 23 7 3 17 24 23 7 21 3
30 26 46 30 26 46
18 52 41 18 52 41
heap H heap H
35 35
39 44 39 44
15 16
Fibonacci heap: insert analysis
‣ preliminaries
Amortized cost. ĉi = ci + ∆Φ = O(1).
‣ insert
‣ extract the minimum
Φ(H) = trees(H) + 2 ⋅ marks(H) ‣ decrease key
‣ bounding the rank
min
‣ meld and delete
17 24 23 7 21 3
SECTION 19.2
30 26 46
18 52 41
heap H
35
39 44
17
Useful primitive. Combine two trees T1 and T2 of rank k. ・Delete min; meld its children into root list; update min.
・Make larger root be a child of smaller root. ・Consolidate trees so that no two roots have same rank.
・Resulting tree T ' has rank k + 1.
15 3 3
min
15 18 52 41 7 24 23 17 3
56 24 33 18 52 41
77 39 44 56 24 33 39 44
30 26 46
18 52 41
77
still heap-ordered 35 39 44
19 20
Fibonacci heap: extract the minimum Fibonacci heap: extract the minimum
・Delete min; meld its children into root list; update min. ・Delete min; meld its children into root list; update min.
・Consolidate trees so that no two roots have same rank. ・Consolidate trees so that no two roots have same rank.
current
min 7 24 23 17 18 52 41 min 7 24 23 17 18 52 41
30 26 46 39 44 30 26 46 39 44
35 35
21 22
Fibonacci heap: extract the minimum Fibonacci heap: extract the minimum
・Delete min; meld its children into root list; update min. ・Delete min; meld its children into root list; update min.
・Consolidate trees so that no two roots have same rank. ・Consolidate trees so that no two roots have same rank.
rank rank
0 1 2 3 0 1 2 3
current current
min 7 24 23 17 18 52 41 min 7 24 23 17 18 52 41
30 26 46 39 44 30 26 46 39 44
35 35
23 24
Fibonacci heap: extract the minimum Fibonacci heap: extract the minimum
・Delete min; meld its children into root list; update min. ・Delete min; meld its children into root list; update min.
・Consolidate trees so that no two roots have same rank. ・Consolidate trees so that no two roots have same rank.
rank rank
0 1 2 3
link 23 to 17 0 1 2 3
current current
min 7 24 23 17 18 52 41 min 7 24 23 17 18 52 41
30 26 46 39 44 30 26 46 39 44
35 35
25 26
Fibonacci heap: extract the minimum Fibonacci heap: extract the minimum
・Delete min; meld its children into root list; update min. ・Delete min; meld its children into root list; update min.
・Consolidate trees so that no two roots have same rank. ・Consolidate trees so that no two roots have same rank.
rank rank
link 17 to 7 0 1 2 3
link 24 to 7 0 1 2 3
current current
min
min 7 24 17 18 52 41 24 7 18 52 41
30 26 46 23 39 44 26 46 17 30 39 44
35 35 23
27 28
Fibonacci heap: extract the minimum Fibonacci heap: extract the minimum
・Delete min; meld its children into root list; update min. ・Delete min; meld its children into root list; update min.
・Consolidate trees so that no two roots have same rank. ・Consolidate trees so that no two roots have same rank.
rank rank
0 1 2 3 0 1 2 3
current current
min min
7 18 52 41 7 18 52 41
24 17 30 39 44 24 17 30 39 44
26 46 23 26 46 23
35 35
29 30
Fibonacci heap: extract the minimum Fibonacci heap: extract the minimum
・Delete min; meld its children into root list; update min. ・Delete min; meld its children into root list; update min.
・Consolidate trees so that no two roots have same rank. ・Consolidate trees so that no two roots have same rank.
rank rank
0 1 2 3
link 41 to 18 0 1 2 3
current current
min min
7 18 52 41 7 18 52 41
24 17 30 39 44 24 17 30 39 44
26 46 23 26 46 23
35 35
31 32
Fibonacci heap: extract the minimum Fibonacci heap: extract the minimum
・Delete min; meld its children into root list; update min. ・Delete min; meld its children into root list; update min.
・Consolidate trees so that no two roots have same rank. ・Consolidate trees so that no two roots have same rank.
rank rank
0 1 2 3 0 1 2 3
current current
min min
7 52 18 7 52 18
24 17 30 41 39 24 17 30 41 39
26 46 23 44 26 46 23 44
35 35
33 34
Fibonacci heap: extract the minimum Fibonacci heap: extract the minimum analysis
・Delete min; meld its children into root list; update min. Actual cost. ci = O(rank(H)) + O(trees(H)).
・Consolidate trees so that no two roots have same rank. ・O(rank(H)) to meld min's children into root list. ≤ rank(H) children
min
7 52 18 Amortized cost. O(log n).
・ĉi = ci + ∆Φ = O(rank(H)) + O(rank(H')).
・The rank of a Fibonacci heap with n elements is O(log n).
24 17 30 41 39
Fibonacci lemma
(stay tuned)
35
35 36
Fibonacci heap vs. binomial heaps
Observation. If only INSERT and EXTRACT-MIN operations, then all trees are
binomial trees. F IBONACCI H EAPS
we link only trees of equal rank
‣ preliminaries
B0 B1 B2 B3 ‣ insert
‣ extract the minimum
‣ decrease key
‣ bounding the rank
‣ meld and delete
Binomial heap property. This implies rank(H) ≤ log2 n.
SECTION 19.3
37
Intuition for deceasing the key of node x. Intuition for deceasing the key of node x.
・If heap-order is not violated, decrease the key of x. ・If heap-order is not violated, decrease the key of x.
・Otherwise, cut tree rooted at x and meld into root list. ・Otherwise, cut tree rooted at x and meld into root list.
6 6 7
8 29 10 44 8 29 10 44 45 32
30 23 22 48 31 17 23 22 48 31 17 55
45 32 24 50 24 50
55
39 40
Fibonacci heap: decrease key Fibonacci heap: decrease key
Intuition for deceasing the key of node x. Intuition for deceasing the key of node x.
・ If heap-order is not violated, decrease the key of x. ・If heap-order is not violated, decrease the key of x.
・Otherwise, cut tree rooted at x and meld into root list. ・Otherwise, cut tree rooted at x and meld into root list.
・Problem: number of nodes not exponential in rank.
decrease-key of 22 to 4
decrease-key of 48 to 3
decrease-key of 31 to 2 6 6 3 2 1
7 5 7 5 4
decrease-key of 17 to 1
8 29 10 44 45 32 24 8 29 10 44 45 32 24 50
rank = 4, nodes = 5
22 48 31 17 55 55
50
41 42
Intuition for deceasing the key of node x. Case 1. [heap order not violated]
・If heap-order is not violated, decrease the key of x. ・Decrease key of x.
・Otherwise, cut tree rooted at x and meld into root list. ・Change heap min pointer (if necessary).
・Solution: as soon as a node has its second child cut,
cut it off also and meld into root list (and unmark it).
decrease-key of x from 46 to 29
min min
7 18 38 7 18 38
marked node:
24 17 23 21 39 41 24 17 23 21 39 41
one child already cut
26 46 30 52 26 46 30 52
35 88 72 35 88 72
43 44
Fibonacci heap: decrease key Fibonacci heap: decrease key
Case 1. [heap order not violated] Case 2a. [heap order violated]
・ Decrease key of x. ・Decrease key of x.
・Change heap min pointer (if necessary). ・Cut tree rooted at x, meld into root list, and unmark.
・If parent p of x is unmarked (hasn't yet lost a child), mark it;
Otherwise, cut p, meld into root list, and unmark
(and do so recursively for all ancestors that lose a second child).
7 18 38 7 18 38
24 17 23 21 39 41 24 17 23 21 39 41
p
26 29 30 52 26 29 30 52
x x
35 88 72 35 88 72
45 46
Case 2a. [heap order violated] Case 2a. [heap order violated]
・Decrease key of x. ・Decrease key of x.
・Cut tree rooted at x, meld into root list, and unmark. ・Cut tree rooted at x, meld into root list, and unmark.
・If parent p of x is unmarked (hasn't yet lost a child), mark it; ・If parent p of x is unmarked (hasn't yet lost a child), mark it;
Otherwise, cut p, meld into root list, and unmark Otherwise, cut p, meld into root list, and unmark
(and do so recursively for all ancestors that lose a second child). (and do so recursively for all ancestors that lose a second child).
7 18 38 7 18 38
24 17 23 21 39 41 24 17 23 21 39 41
p p
26 15 30 52 26 15 30 52
x x
35 88 72 35 88 72
47 48
Fibonacci heap: decrease key Fibonacci heap: decrease key
Case 2a. [heap order violated] Case 2a. [heap order violated]
・ Decrease key of x. ・Decrease key of x.
・Cut tree rooted at x, meld into root list, and unmark. ・Cut tree rooted at x, meld into root list, and unmark.
・If parent p of x is unmarked (hasn't yet lost a child), mark it; ・If parent p of x is unmarked (hasn't yet lost a child), mark it;
Otherwise, cut p, meld into root list, and unmark Otherwise, cut p, meld into root list, and unmark
(and do so recursively for all ancestors that lose a second child). (and do so recursively for all ancestors that lose a second child).
x 15 7 18 38 x 15 7 18 38
72 24 17 23 21 39 41 72 24 17 23 21 39 41
p p
26 30 52 26 30 52
35 88 35 88
49 50
Case 2b. [heap order violated] Case 2b. [heap order violated]
・Decrease key of x. ・Decrease key of x.
・Cut tree rooted at x, meld into root list, and unmark. ・Cut tree rooted at x, meld into root list, and unmark.
・If parent p of x is unmarked (hasn't yet lost a child), mark it; ・If parent p of x is unmarked (hasn't yet lost a child), mark it;
Otherwise, cut p, meld into root list, and unmark Otherwise, cut p, meld into root list, and unmark
(and do so recursively for all ancestors that lose a second child). (and do so recursively for all ancestors that lose a second child).
15 7 18 38 15 7 18 38
72 24
24 17 23 21 39 41 72 24
24 17 23 21 39 41
p 26 p 26
30 52 30 52
x 35 88 x 5 88
51 52
Fibonacci heap: decrease key Fibonacci heap: decrease key
Case 2b. [heap order violated] Case 2b. [heap order violated]
・ Decrease key of x. ・Decrease key of x.
・Cut tree rooted at x, meld into root list, and unmark. ・Cut tree rooted at x, meld into root list, and unmark.
・If parent p of x is unmarked (hasn't yet lost a child), mark it; ・If parent p of x is unmarked (hasn't yet lost a child), mark it;
Otherwise, cut p, meld into root list, and unmark Otherwise, cut p, meld into root list, and unmark
(and do so recursively for all ancestors that lose a second child). (and do so recursively for all ancestors that lose a second child).
x min x min
15 5 7 18 38 15 5 7 18 38
72 24
24 17 23 21 39 41 72 24
24 17 23 21 39 41
second child cut
p 26 p 26
30 52 30 52
88 88
53 54
Case 2b. [heap order violated] Case 2b. [heap order violated]
・Decrease key of x. ・Decrease key of x.
・Cut tree rooted at x, meld into root list, and unmark. ・Cut tree rooted at x, meld into root list, and unmark.
・If parent p of x is unmarked (hasn't yet lost a child), mark it; ・If parent p of x is unmarked (hasn't yet lost a child), mark it;
Otherwise, cut p, meld into root list, and unmark Otherwise, cut p, meld into root list, and unmark
(and do so recursively for all ancestors that lose a second child). (and do so recursively for all ancestors that lose a second child).
x min p x min p
15 5 26 7 18 38 15 5 26 7 18 38
72 88 24
24 17 23 21 39 41 72 88 p' 24
24 17 23 21 39 41
55 56
Fibonacci heap: decrease key Fibonacci heap: decrease key analysis
Case 2b. [heap order violated] Actual cost. ci = O(c), where c is the number of cuts.
・ Decrease key of x. ・O(1) time for changing the key.
・Cut tree rooted at x, meld into root list, and unmark. ・O(1) time for each of c cuts, plus melding into root list.
・If parent p of x is unmarked (hasn't yet lost a child), mark it;
Otherwise, cut p, meld into root list, and unmark Change in potential. ∆Φ = O(1) – c.
(and do so recursively for all ancestors that lose a second child). ・trees(H') = trees(H) + c.
decrease-key of x from 35 to 5
・marks(H') ≤ marks(H) – c + 2.
each cut (except first) unmarks a node
last cut may or may not mark a node
・ΔΦ ≤ c + 2 ⋅ (-c + 2) = 4 – c.
x min p p' p''
15 5 26 24 7 18 38
Amortized cost. ĉi = ci + ∆Φ = O(1).
but don't
72 88 mark parent 17 23 21 39 41
if it's a root
30 52
Φ(H) = trees(H) + 2 ⋅ marks(H)
57 58
Analysis summary
Insert. O(1).
F IBONACCI H EAPS Delete-min. O(rank(H)) amortized.
Decrease-key. O(1) amortized.
‣ preliminaries
‣ insert
‣ extract the minimum Fibonacci lemma. Let H be a Fibonacci heap with n elements.
‣ decrease key Then, rank(H) = O(log n).
SECTION 19.4
60
Bounding the rank Bounding the rank
Lemma 1. Fix a point in time. Let x be a node of rank k, and let y1, …, yk Lemma 1. Fix a point in time. Let x be a node of rank k, and let y1, …, yk
denote its current children in the order in which they were linked to x. denote its current children in the order in which they were linked to x.
Then: Then:
x x
0 B7 i = 1 0 B7 i = 1
rank(yi ) rank(yi )
i 2 B7 i 2 i 2 B7 i 2
y1 y2 … yk y1 y2 … yk
61 62
Lemma 1. Fix a point in time. Let x be a node of rank k, and let y1, …, yk Lemma 2. Let sk be minimum number of elements in any Fibonacci heap of
denote its current children in the order in which they were linked to x. rank k. Then sk ≥ Fk+2, where Fk is the kth Fibonacci number.
Then:
x
0 B7 i = 1 Pf. [by strong induction on k]
rank(yi )
i 2 B7 i 2 ・Base cases: s0 = 1 and s1 = 2.
y1 y2 … yk ・Inductive hypothesis: assume si ≥ Fi+2 for i = 0, …, k – 1.
・As in Lemma 1, let let y1, …, yk denote its current children in the order in
which they were linked to x.
Def. Let Tk be smallest possible tree of rank k satisfying property.
sk ≥ 1 + 1 + (s0 + s1 + … + sk–2) (Lemma 1)
T4 T5 T6
≥ (1 + F1) + F2 + F3 + … + Fk (inductive hypothesis)
F6 = 8 F7 = 13 F8 = F6 + F7 = 8 + 13 = 21
63 64
Bounding the rank Fibonacci fact 1
Fibonacci lemma. Let H be a Fibonacci heap with n elements. Def. The Fibonacci sequence is: 0, 1, 1, 2, 3, 5, 8, 13, 21, …
Then, rank(H) ≤ logφ n, where φ is the golden ratio = (1 + √5) / 2 ≈ 1.618.
0 B7 k = 0
Fk = 1 B7 k = 1
Pf.
Fk + Fk B7 k 2
・Let H is a Fibonacci heap with n elements and rank k.
1 2
65 66
Def. The Fibonacci sequence is: 0, 1, 1, 2, 3, 5, 8, 13, 21, … Fibonacci numbers arise both in nature and algorithms.
0 B7 k = 0
Fk = 1 B7 k = 1
Fk 1 + Fk 2 B7 k 2
≥ φk – 1 + φk – 2 (inductive hypothesis)
= φk – 2 (1 + φ) (algebra) pinecone
= φk – 2 φ2 (φ2 = φ + 1)
cauliflower
= φk. ▪ (algebra)
67 68
Fibonacci heap: meld
23 24 17 7 3 21
SECTION 19.2, 19.3
30 26 46
18 52 41
heap H1 heap H2
35
39 44
70
Meld. Combine two Fibonacci heaps (destroying old heaps). Actual cost. ci = O(1).
Change in potential. ∆Φ = 0.
Recall. Root lists are circular, doubly-linked lists. Amortized cost. ĉi = ci + ∆Φ = O(1).
min min
23 24 17 7 3 21 23 24 17 7 3 21
30 26 46 30 26 46
18 52 41 18 52 41
35 heap H 35 heap H
39 44 39 44
71 72
Fibonacci heap: delete Priority queues performance cost summary
† amortized
73 74
Heaps of heaps
・b-heaps.
P RIORITY Q UEUES ・Fat heaps.
・2-3 heaps.
‣ binary heaps ・Leaf heaps.
‣ d-ary heaps ・Thin heaps.
・Skew heaps.
‣ binomial heaps
・Splay heaps.
‣ Fibonacci heaps ・Weak heaps.
‣ advanced topics ・Leftist heaps.
・Quake heaps.
・Pairing heaps.
・Violation heaps.
・Run-relaxed heaps.
・Rank-pairing heaps.
・Skew-pairing heaps.
・Rank-relaxed heaps.
・Lazy Fibonacci heaps. 76
Brodal queues Strict Fibonacci heaps
Q. Can we achieve same running time as for Fibonacci heap but with Q. Can we achieve same running time as for Fibonacci heap but with
worst-case bounds per operation (instead of amortized)? worst-case bounds per operation (instead of amortized) in pointer model?
Theory. Same amortized running times as Fibonacci heaps for all operations
except DECREASE-KEY.
・O(log n) amortized. [Fredman et al. 1986]
・Ω(log log n) lower bound on amortized cost. [Fredman 1999]
・ 2 O(log log n) amortized. [Pettie 2005]
79 80
Pairing heaps Priority queues performance cost summary
RESEARCHCONlRlWlIONS
INSERT O(1) O(log n) O(log n) O(1) O(1) O(1)
Algorithms and
Data Structures Pairing Heaps:
G. Scott Graham
Editor
Experiments and Analysis EXTRACT-MIN O(n) O(log n) O(log n) O(log n) O(log n) O(log n)
JOHN T. STASKO and JEFFREY SCOTT VlllER MELD O(1) O(n) O(log n) O(1) O(1) O(1)
ABSTRACT: The pairing heap has recently been and practical importance from their use in solving a
introduced as a new data structure for priority queues.
Pairing heaps are extremely simple to implement and
seem to be very efficient in practice, but they are difficult
wide range of combinatorial problems, including job
scheduling, minimal spanning tree, shortest path,
and graph traversal.
FIND-MIN O(n) O(1) O(log n) O(1) O(1) O(1)
to analyze theoretically, and open problems remain. It Priority queues support the operations insert,
has been conjectured that they achieve the same find-min, and delete-min; additional operations often
amortized time bounds as Fibonacci heaps, namely, include decrease-key and delete. The insert(t, v) opera-
O(log n) time for delete and delete-min and O(1) for
all other operations, where n is the size of the priority
tion adds item t with key value v to the priority
queue. The find-min operation returns the item
† amortized
queue at the time of the operation. We provide empirical with minimum key value. The delete-min operation
evidence that supports this conjecture. The most returns the item with minimum key value and
promising algorithm in our simulations is a new variant removes it from the priority queue. The decrease-
of the twopass method, called auxiliary twopass. We key(t, d) operation reduces item t’s key value by d.
prove that, assuming no decrease-key operations are The delete(t) operation removes item t from the
performed, it achieves the same amortized time bounds as priority queue. The decrease-key and delete opera-
Fibonacci heaps. tions require that a pointer to the location in the
priority queue of item t be supplied explicitly, since
1. INTRODUCTION priority queues do not support searching for arbi-
A priority queue is an abstract data type for main- trary items by value. Some priority queues also sup-
taining and manipulating a set of items based on port the merge operation, which combines two item-
priority [I]. Prio’rity queues derive great theoretical disjoint priority queues. 81 82
Support was provided in part by NSF research grant DCR-84-03613, an NSF
We will concentrate on the insert, delete-min, and
Presidential Young Investigator Award, an IBM Faculty Development Award, decrease-key operations because they are the opera-
and a Guggenheim Fellowship.
tions that primarily distinguish priority queues from
Part of this research was performed at Mathematical Sciences Research Insti-
tute. Berkeley, Calif., and the Institut National de Recherche en Informatique other set manipulation algorithms and because they
et en Automatique, Rocquencourt. France. are the critical operations as far as the time bounds
0 1987 ACM OOOl-0782/87/0300-0234 75a: are concerned.
Priority queues with integer priorities Priority queues with integer priorities
Assumption. Keys are integers between 0 and C. Assumption. Keys are integers between 0 and C.
Theorem. [Thorup 2004] There exists a priority queue that supports INSERT, Theorem. [Thorup 2004] There exists a priority queue that supports INSERT,
FIND-MIN, and DECREASE-KEY in constant time and EXTRACT-MIN and DELETE-KEY FIND-MIN, and DECREASE-KEY in constant time and EXTRACT-MIN and DELETE-KEY
in either O(log log n) or O(log log C) time. in either O(log log n) or O(log log C) time.
ARTICLE IN PRESS Corollary 1. Can implement Dijkstra's algorithm in either O(m log log n) or
O(m log log C) time.
Journal of Computer and System Sciences 69 (2004) 330–353
http://www.elsevier.com/locate/jcss
Abstract
We consider Fibonacci heap style integer priority queues supporting find-min, insert, and decrease key
operations in constant time. We present a deterministic linear space solution that with n integer keys
supports delete in Oðlog log nÞ time. If the integers are in the range ½0; NÞ; we can also support delete in
Oðlog log NÞ time.
Even for the special case of monotone priority queues, where the minimum has to be non-decreasing, the
best previous bounds on delete were Oððlog nÞ1=ð3$eÞ Þ and Oððlog NÞ1=ð4$eÞ Þ: These previous bounds used
both randomization and amortization. Our new bounds are deterministic, worst-case, with no restriction to
monotonicity, and exponentially faster. 83 84
As a classical application, for a directed graph with n nodes and m edges with non-negative integer
weights, we get single source shortest paths in Oðm þ n log log nÞ time, or Oðm þ n log log CÞ if C is the
maximal edge weight. The latter solves an open problem of Ahuja, Mehlhorn, Orlin, and Tarjan from 1990.
Soft heaps Soft heaps
Goal. Break information-theoretic lower bound by allowing priority queue to Goal. Break information-theoretic lower bound by allowing priority queue to
corrupt 10% of the keys (by increasing them). corrupt 10% of the keys (by increasing them).
Representation.
・Set of binomial trees (with some subtrees missing).
・Each node may store several elements.
・Each node stores a value that is an upper bound on the original keys.
・Binomial trees are heap-ordered with respect to these values.
elements inserted 0 11 22 44 55
88 11 22 99 44 33 77 66 0 55
57 66 77 88 99
soft heap
corrupted
85 86
Goal. Break information-theoretic lower bound by allowing priority queue to Goal. Break information-theoretic lower bound by allowing priority queue to
corrupt 10% of the keys (by increasing them). corrupt 10% of the keys (by increasing them).
Theorem. [Chazelle 2000] Starting from an empty soft heap, any sequence Q. Brilliant. But how could it possibly be useful?
of n INSERT, MIN, EXTRACT-MIN, MELD, and DELETE operations takes O(n) time Ex. Linear-time deterministic selection. To find kth smallest element:
and at most 10% of its elements are corrupted at any given time. ・Insert the n elements into soft heap.
・Extract the minimum element n / 2 times.
The Soft Heap: An Approximate Priority Queue with
・The largest element deleted ≥ 4n / 10 elements and ≤ 6n / 10 elements.
Optimal Error Rate ・Can remove ≥ 5n / 10 of elements and recur.
BERNARD CHAZELLE
・T(n) ≤ T(3n / 5) + O(n) ⇒ T(n) = O(n). ▪
Princeton University, Princeton, New Jersey, and NEC Research Institute
Abstract. A simple variant of a priority queue, called a soft heap, is introduced. The data structure
supports the usual operations: insert, delete, meld, and findmin. Its novelty is to beat the logarithmic
bound on the complexity of a heap in a comparison-based model. To break this information-theoretic
barrier, the entropy of the data structure is reduced by artificially raising the values of certain keys.
Given any mixed sequence of n operations, a soft heap with error rate ! (for any 0 ! ! " 1/2) ensures
that, at any time, at most !n of its items have their keys raised. The amortized complexity of each
operation is constant, except for insert, which takes O(log 1/!) time. The soft heap is optimal for any
value of ! in a comparison-based model. The data structure is purely pointer-based. No arrays are
used and no numeric assumptions are made on the keys. The main idea behind the soft heap is to
move items across the data structure not individually, as is customary, but in groups, in a
data-structuring equivalent of “car pooling.” Keys must be raised as a result, in order to preserve the
heap ordering of the data structure. The soft heap can be used to compute exact or approximate
medians and percentiles optimally. It is also useful for approximate sorting and for computing
minimum spanning trees of general graphs.
Categories and Subject Descriptors: E.1 [Data Structures]: Nonnumerical Algorithms and Problems 87 88
General Terms: Theory
Additional Key Words and Phrases: Amortization, heap, priority queue, soft heap
Soft heaps
Theorem. [Chazelle 2000] There exists an O(m α(m, n)) time deterministic
algorithm to compute an MST in a graph with n nodes and m edges.
BERNARD CHAZELLE
Princeton University, Princeton, New Jersey, and NEC Research Institute
Abstract. A deterministic algorithm for computing a minimum spanning tree of a connected graph is
presented. Its running time is O(m ! (m, n)), where ! is the classical functional inverse of
Ackermann’s function and n (respectively, m) is the number of vertices (respectively, edges). The
algorithm is comparison-based: it uses pointers, not arrays, and it makes no numeric assumptions on
the edge costs.
Categories and Subject Descriptors: F.2.2 [Analysis of Algorithms and Problem Complexity]:
Nonnumerical Algorithms and Problems.
General Terms: Theory
Additional Key Words and Phrases: Graphs, matroids, minimum spanning trees 89
1. Introduction
The history of the minimum spanning tree (MST) problem is long and rich, going
as far back as Borůvka’s work in 1926 [Borůvka 1926; Graham and Hell 1985;
Nešetřil 1997]. In fact, MST is perhaps the oldest open problem in computer
science. According to Nešetřil [1997], “this is a cornerstone problem of combina-
torial optimization and in a sense its cradle.” Textbook algorithms run in O(m
log n) time, where n and m denote, respectively, the number of vertices and
edges in the graph. Improvements to O(m log log n) were given independently
by Yao [1975] and by Cheriton and Tarjan [1976]. In the mid-eighties, Fredman
and Tarjan [1987] lowered the complexity to O(m " (m, n)), where " (m, n) is
the number of log-iterations necessary to map n to a number less than m/n. In
the worst case, m ! O(n) and the running time is O(m log* m). Soon after, the
A preliminary version of this paper appeared as CHAZELLE, B. 1997. A faster deterministic algorithm
for minimum spanning trees. In Proceedings of the 38th Annual IEEE Symposium on Foundations of
Computer Science. IEEE Computer Society Press, Los Alamitos, Calif., pp. 22–31.
This work was supported in part by the National Science Foundation (NSF) Grants CCR 93-01254
and CCR 96-23768, ARO Grant DAAH04-96-1-0181, and NEC Research Institute.
Author’s address: Department of Computer Science, Princeton University, 35 Olden Street, Prince-
ton, NJ 083-44-2087, e-mail: [email protected] and NEC Research Institute, e-mail:
[email protected].
Permission to make digital / hard copy of part or all of this work for personal or classroom use is
granted without fee provided that the copies are not made or distributed for profit or commercial
advantage, the copyright notice, the title of the publication, and its date appear, and notice is given
that copying is by permission of the Association for Computing Machinery (ACM), Inc. To copy
otherwise, to republish, to post on servers, or to redistribute to lists, requires prior specific permission
and / or a fee.
© 2000 ACM 0004-5411/00/1100-1028 $05.00
Journal of the ACM, Vol. 47, No. 6, November 2000, pp. 1028 –1047.