Towards a Final Analysis of Pairing Heaps

Total Page:16

File Type:pdf, Size:1020Kb

Towards a Final Analysis of Pairing Heaps Towards a Final Analysis of Pairing Heaps Seth Pettie£ Max Planck Institut f¨ur Informatik Abstract perimentally that the cost of decreasekey is virtually, if not asymptotically, constant, and Moret & Shapiro [10] con- Fredman, Sedgewick, Sleator, and Tarjan proposed the cluded that in Prim’s minimum spanning tree algorithm pairing heap as a self-adjusting, streamlined version of the pairing heaps are faster than the alternatives. These in- Fibonacci heap. It provably supports all priority queue cluded binary heaps, splay heaps, Fibonacci heaps, and oth- operations in logarithmic time and is known to be extremely ers. The pairing heap is now included in implementations efficient in practice. However, despite its simplicity and of the GNU C++ library and the LEDA library [9]. empirical superiority, the pairing heap is one of the few Fredman [4] proved the remarkable result that on a spe- popular data structures whose basic complexity remains cific distribution of operation sequences, no (generalized) open. In this paper we prove that pairing heaps support pairing heap can perform optimally. Specifically, there Ò Ò ÐÓ Ò the deletemin operation in optimal logarithmic time and are sequences of Ò insertions, deletions, and de- Ò ÐÓ Ò ÐÓ ÐÓ Òµ all other operations (insert, meld, and decreasekey) in time ª´ Ô creasekeys that take time to execute. ¾ ÐÓ ÐÓ Ò Ç ´¾ µ. This result gives the first sub-logarithmic He concluded that the complexity of the decreasekey op- Òµ Ç ´ÐÓ Òµ time bound for decreasekey and comes close to the lower eration lies somewhere between ª´ÐÓ ÐÓ and . Òµ bound of ª´ÐÓ ÐÓ established by Fredman. The difference between optimality and suboptimality is not as stark as one might think. Because the user of Pairing heaps have a well known but poorly understood re- a priority queue has already resigned himself to paying lationship to splay trees and, to date, the transfer of ideas Ò ÐÓ Òµ Ò ¢´ time for deletemins, a non-constant cost for has flowed in one direction: from splaying to pairing. One decreasekey will only be detectable for sufficiently dense contribution of this paper is a new analysis that reasons ex- input sequences. (Density is the ratio of decreasekeys plicitly with information-theoretic measures. Whether these to deletemins.) For example, if Fredman’s lower bound ideas could contribute to the analysis of splay trees is an turns out to be tight, this would imply that in the stan- open question. dard graph applications (shortest paths, weighted match- ings, minimum spanning trees, etc.) pairing heaps are no worse than Fibonacci heaps when the number of edges is ´Ò ÐÓ Ò ÐÓ ÐÓ Òµ 1 Introduction Ç . Fredman’s lower bound only estab- lishes that the pairing heap is suboptimal in the density ´ÐÓ Ò ÐÓ ÐÓ Òµ Ó´ÐÓ Ò ÐÓ ÐÓ Òµ Twenty years ago Fredman et al. [5] proposed the pair- range to ; for other ing heap, a simple priority queue data structure that they ranges its complexity remains open. conjectured to be as efficient as Fibonacci heaps [6]. In In this paper we provide the first proof that pairing heaps their original analysis Fredman et al. proved that all priority support decreasekeys in sub-logarithmic time. queue operations (insert, deletemin, meld, and decreasekey) ´ÐÓ Òµ Ò were supported in Ç amortized time, where is the Theorem 1 The pairing heap executes a mixed sequence size of the heap. Although logarithmic time is optimal for Ò of Ñ priority queue operations, including deletemins, in deletemin, theoretically efficient priority queues like the Fi- Ô ¾ ÐÓ ÐÓ Ò Ç ´Ò ÐÓ Ò · Ѿ µ time. ´½µ bonacci heap enjoy Ç amortized time for insert, meld, and decreasekey. Ô ¾ ÐÓ ÐÓ Ò In the intervening years pairing heaps have become the The ¾ function lies in the asymptotic spectrum ¯ ´ ´ÐÓ ÐÓ Òµµ Ó´ÐÓ Òµ priority queue of choice in applications requiring the de- between poly and . It is an ugly creasekey operation. Stasko & Vitter [8] demonstrated ex- bound and undoubtedly not tight. Nonetheless, we can now £ claim with certainty that the pairing heap is asymptotically Supported by an Alexander von Humboldt Postdoctoral Fellowship. Ô ¾ ÐÓ ÐÓ Ò Ç ´ÐÓ Ò¾ Email: [email protected]. optimal for any density µ. The Pairing Heap. The structure consists of a single rectly. rooted tree where the children of a node are assigned some left-to-right ordering. Each node is identified with a key and Overview. In Section 2 we discuss the original analysis of the key of a parent is no larger than the key of any child. (It pairing heaps in order to provide a reference point for our is heap ordered.) During the execution of an operation there approach. In Section 3 we introduce some key ideas and, as Ô Ü Ü ¾ may be multiple rooted trees. If ½ are two roots then ´ ÐÓ Òµ a warmup, sketch a proof that decreasekeys take Ç Ü Ü Ü Ü ¾ ½ ¾ linking ½ and either makes the leftmost child of or time. The full proof of Theorem 1 appears in Section 4. vice versa, depending on which node’s key is smaller. The In Section 5 we discuss some open problems and make a operations are implemented as follows. An insert creates conjecture on the actual complexity of pairing heaps. a new node with a given key and links it with the existing root. A meld of two given heaps simply links their roots. A decreasekey operation is given both a node and the heap to 2 Pairing Heaps, Entropy, and Potentials which it belongs. It detaches the subtree rooted at the given node, reduces the node’s key appropriately, then links the The original analysis of pairing heaps [5] was a stroke of given node with the existing heap root. elegance. However, it was achieved by drawing a slightly By the heap-order property the minimum node is neces- unnatural correspondence between pairing heaps and splay sarily at the root. A deletemin proceeds by detaching ev- trees. In particular, Fredman et al. [5] map a general rooted ery child of the root then linking children in a prescribed tree to a binary tree: a node’s left child in the binary rep- fashion until one tree remains. Several linking strategies resentation corresponds to its leftmost child in the general have been proposed [5, 8, 3]. In this paper we consider the tree and its right child corresponds to its right sibling in the Ü Ü Ü ¾ standard 2-pass version given in [5]. Let ½ be general tree. They observe that in the binary representation the children of the root in left-to-right order. In the cou- a deletemin operation looks very similar to a splay up the ´Ü Ü µ ´Ü Ü µ ´Ü Ü µ ¾ ¿ pling pass we link the pairs ½ . right shoulder of the tree. By borrowing the same poten- Ý Ý Ý ¾ Let ½ be the winners of these linkings, where tial function and analysis from [11] they derive amortized Ý Ü Ý if is odd. In the accumulation pass we link logarithmic bounds on deletemins and all other operations. Ý Ý ½ ¾ with , the winner of this comparison with ,andso Ý Ý Ý ·½ on; in general, is linked with the winner of . Previous Work. Fredmanetal.[5]provedthatthe ´ÐÓ Òµ 2-pass version of deletemin takes Ç amortized time. They also considered arbitrary coupling, multi- pass,andbackward 2-pass.1 It was shown that arbi- Ô Òµ trary coupling takes ¢´ time and that multipass takes Ô ´ÐÓ Ò ÐÓ ÐÓ Ò ÐÓ ÐÓ ÐÓ Òµ Ó´ Òµ Ç time. No time bound is known for backward 2-pass. Stasko and Vitter [8] Figure 1. The general trees are marked with and Elmasry [2] proposed several variants of these linking circular nodes, the binary representation with strategies. In [8] it was observed that with a small inser- square nodes. On the left is a tree with en- tion buffer the amortized cost of insertion can be reduced tropy zero and on the right is a tree with en- ´½µ ÐÓ´ µ to Ç . Iacono [7] proved that even the standard 2-pass tropy . ´½µ version supports insert and meld in Ç time. Bringing the ´½µ insertion cost to Ç is significant in applications where the heap is likely to end in a non-empty state, for instance, The analysis of [5] uses the potential function ¨ È ×´Üµ ״ܵ when computing single-source, single-destination shortest ÐÓ where is the size of the subtree rooted Ü paths with Dijkstra’s algorithm. Iacono [7] and Elmasry at Ü in the binary representation. One property of this func- [2] studied the distribution-sensitive properties of pairing tion is that general trees with completely different charac- heaps. Before the pairing heap was introduced Sleator and teristics can get assigned the same potential. In Figure 1 Tarjan [12] invented the skew-heap, a self-adjusting data we have, among rooted 5-node trees, the two with mini- structure that supported insert, deletemin, and meld in log- mum and maximum entropy, respectively. Here entropy is arithmic time; however, it did not support decreasekeys di- the logarithm of the number of linear extentions consistent with the heap-order. Since both trees get assigned a poten- 1Arbitrary coupling partitions the roots into couplets arbitrarily and ac- tial of ÐÓ´ µ one wonders if there is an interpretation of cumulates the couplet winners in no specific order. Multipass repeatedly executes the coupling pass of 2-pass and backward 2-pass accumulates the this potential that makes sense in the context of a priority Ý Ý ½ roots from left-to-right rather than right-to-left.
Recommended publications
  • Lecture 04 Linear Structures Sort
    Algorithmics (6EAP) MTAT.03.238 Linear structures, sorting, searching, etc Jaak Vilo 2018 Fall Jaak Vilo 1 Big-Oh notation classes Class Informal Intuition Analogy f(n) ∈ ο ( g(n) ) f is dominated by g Strictly below < f(n) ∈ O( g(n) ) Bounded from above Upper bound ≤ f(n) ∈ Θ( g(n) ) Bounded from “equal to” = above and below f(n) ∈ Ω( g(n) ) Bounded from below Lower bound ≥ f(n) ∈ ω( g(n) ) f dominates g Strictly above > Conclusions • Algorithm complexity deals with the behavior in the long-term – worst case -- typical – average case -- quite hard – best case -- bogus, cheating • In practice, long-term sometimes not necessary – E.g. for sorting 20 elements, you dont need fancy algorithms… Linear, sequential, ordered, list … Memory, disk, tape etc – is an ordered sequentially addressed media. Physical ordered list ~ array • Memory /address/ – Garbage collection • Files (character/byte list/lines in text file,…) • Disk – Disk fragmentation Linear data structures: Arrays • Array • Hashed array tree • Bidirectional map • Heightmap • Bit array • Lookup table • Bit field • Matrix • Bitboard • Parallel array • Bitmap • Sorted array • Circular buffer • Sparse array • Control table • Sparse matrix • Image • Iliffe vector • Dynamic array • Variable-length array • Gap buffer Linear data structures: Lists • Doubly linked list • Array list • Xor linked list • Linked list • Zipper • Self-organizing list • Doubly connected edge • Skip list list • Unrolled linked list • Difference list • VList Lists: Array 0 1 size MAX_SIZE-1 3 6 7 5 2 L = int[MAX_SIZE]
    [Show full text]
  • Advanced Data Structures
    Advanced Data Structures PETER BRASS City College of New York CAMBRIDGE UNIVERSITY PRESS Cambridge, New York, Melbourne, Madrid, Cape Town, Singapore, São Paulo Cambridge University Press The Edinburgh Building, Cambridge CB2 8RU, UK Published in the United States of America by Cambridge University Press, New York www.cambridge.org Information on this title: www.cambridge.org/9780521880374 © Peter Brass 2008 This publication is in copyright. Subject to statutory exception and to the provision of relevant collective licensing agreements, no reproduction of any part may take place without the written permission of Cambridge University Press. First published in print format 2008 ISBN-13 978-0-511-43685-7 eBook (EBL) ISBN-13 978-0-521-88037-4 hardback Cambridge University Press has no responsibility for the persistence or accuracy of urls for external or third-party internet websites referred to in this publication, and does not guarantee that any content on such websites is, or will remain, accurate or appropriate. Contents Preface page xi 1 Elementary Structures 1 1.1 Stack 1 1.2 Queue 8 1.3 Double-Ended Queue 16 1.4 Dynamical Allocation of Nodes 16 1.5 Shadow Copies of Array-Based Structures 18 2 Search Trees 23 2.1 Two Models of Search Trees 23 2.2 General Properties and Transformations 26 2.3 Height of a Search Tree 29 2.4 Basic Find, Insert, and Delete 31 2.5ReturningfromLeaftoRoot35 2.6 Dealing with Nonunique Keys 37 2.7 Queries for the Keys in an Interval 38 2.8 Building Optimal Search Trees 40 2.9 Converting Trees into Lists 47 2.10
    [Show full text]
  • COS 423 Lecture 6 Implicit Heaps Pairing Heaps
    COS 423 Lecture 6 Implicit Heaps Pairing Heaps ©Robert E. Tarjan 2011 Heap (priority queue): contains a set of items x, each with a key k(x) from a totally ordered universe, and associated information. We assume no ties in keys. Basic Operations : make-heap : Return a new, empty heap. insert (x, H): Insert x and its info into heap H. delete-min (H): Delete the item of min key from H. Additional Operations : find-min (H): Return the item of minimum key in H. meld (H1, H2): Combine item-disjoint heaps H1 and H2 into one heap, and return it. decrease -key (x, k, H): Replace the key of item x in heap H by k, which is smaller than the current key of x. delete (x, H): Delete item x from heap H. Assumption : Heaps are item-disjoint. A heap is like a dictionary but no access by key; can only retrieve the item of min key: decrease-ke y( x, k, H) and delete (x, H) are given a pointer to the location of x in heap H Applications : Priority-based scheduling and allocation Discrete event simulation Network optimization: Shortest paths, Minimum spanning trees Lower bound from sorting Can sort n numbers by doing n inserts followed by n delete-min’s . Since sorting by binary comparisons takes Ω(nlg n) comparisons, the amortized time for either insert or delete-min must be Ω(lg n). One can modify any heap implementation to reduce the amortized time for insert to O(1) → delete-min takes Ω(lg n) amortized time.
    [Show full text]
  • Pairing Heaps: the Forward Variant
    Pairing heaps: the forward variant Dani Dorfman Blavatnik School of Computer Science, Tel Aviv University, Israel [email protected] Haim Kaplan1 Blavatnik School of Computer Science, Tel Aviv University, Israel [email protected] László Kozma2 Eindhoven University of Technology, The Netherlands [email protected] Uri Zwick3 Blavatnik School of Computer Science, Tel Aviv University, Israel [email protected] Abstract The pairing heap is a classical heap data structure introduced in 1986 by Fredman, Sedgewick, Sleator, and Tarjan. It is remarkable both for its simplicity and for its excellent performance in practice. The “magic” of pairing heaps lies in the restructuring that happens after the deletion of the smallest item. The resulting collection of trees is consolidated in two rounds: a left-to-right pairing round, followed by a right-to-left accumulation round. Fredman et al. showed, via an elegant correspondence to splay trees, that in a pairing heap of size n all heap operations take O(log n) amortized time. They also proposed an arguably more natural variant, where both pairing and accumulation are performed in a combined left-to-right round (called the forward variant of pairing heaps). The analogy to splaying breaks down in this case, and the analysis of the forward variant was left open. In this paper we show that inserting an item and√ deleting the minimum in a forward-variant pairing heap both take amortized time O(log n · 4 log n). This is the first improvement over the √ O( n) bound showed by Fredman et al. three decades ago.
    [Show full text]
  • Advanced Data Structures and Implementation
    www.getmyuni.com Advanced Data Structures and Implementation • Top-Down Splay Trees • Red-Black Trees • Top-Down Red Black Trees • Top-Down Deletion • Deterministic Skip Lists • AA-Trees • Treaps • k-d Trees • Pairing Heaps www.getmyuni.com Top-Down Splay Tree • Direct strategy requires traversal from the root down the tree, and then bottom-up traversal to implement the splaying tree. • Can implement by storing parent links, or by storing the access path on a stack. • Both methods require large amount of overhead and must handle many special cases. • Initial rotations on the initial access path uses only O(1) extra space, but retains the O(log N) amortized time bound. www.getmyuni.com www.getmyuni.com Case 1: Zig X L R L R Y X Y XR YL Yr XR YL Yr If Y should become root, then X and its right sub tree are made left children of the smallest value in R, and Y is made root of “center” tree. Y does not have to be a leaf for the Zig case to apply. www.getmyuni.com www.getmyuni.com Case 2: Zig-Zig L X R L R Z Y XR Y X Z ZL Zr YR YR XR ZL Zr The value to be splayed is in the tree rooted at Z. Rotate Y about X and attach as left child of smallest value in R www.getmyuni.com www.getmyuni.com Case 3: Zig-Zag (Simplified) L X R L R Y Y XR X YL YL Z Z XR ZL Zr ZL Zr The value to be splayed is in the tree rooted at Z.
    [Show full text]
  • Quake Heaps: a Simple Alternative to Fibonacci Heaps
    Quake Heaps: A Simple Alternative to Fibonacci Heaps Timothy M. Chan Cheriton School of Computer Science, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada, [email protected] Abstract. This note describes a data structure that has the same the- oretical performance as Fibonacci heaps, supporting decrease-key op- erations in O(1) amortized time and delete-min operations in O(log n) amortized time. The data structure is simple to explain and analyze, and may be of pedagogical value. 1 Introduction In their seminal paper [5], Fredman and Tarjan investigated the problem of maintaining a set S of n elements under the operations – insert(x): insert an element x to S; – delete-min(): remove the minimum element x from S, returning x; – decrease-key(x, k): change the value of an element x toasmallervaluek. They presented the first data structure, called Fibonacci heaps, that can support insert() and decrease-key() in O(1) amortized time, and delete-min() in O(log n) amortized time. Since Fredman and Tarjan’s paper, a number of alternatives have been pro- posed in the literature, including Driscoll et al.’s relaxed heaps and run-relaxed heaps [1], Peterson’s Vheaps [9], which is based on AVL trees (and is an instance of Høyer’s family of ranked priority queues [7]), Takaoka’s 2-3 heaps [11], Kaplan and Tarjan’s thin heaps and fat heaps [8], Elmasry’s violation heaps [2], and most recently, Haeupler, Sen, and Tarjan’s rank-pairing heaps [6]. The classical pair- ing heaps [4, 10, 3] are another popular variant that performs well in practice, although they do not guarantee O(1) decrease-key cost.
    [Show full text]
  • Heaps Simplified
    Heaps Simplified Bernhard Haeupler2, Siddhartha Sen1;4, and Robert E. Tarjan1;3;4 1 Princeton University, Princeton NJ 08544, fsssix, [email protected] 2 CSAIL, Massachusetts Institute of Technology, [email protected] 3 HP Laboratories, Palo Alto CA 94304 Abstract. The heap is a basic data structure used in a wide variety of applications, including shortest path and minimum spanning tree algorithms. In this paper we explore the design space of comparison-based, amortized-efficient heap implementations. From a consideration of dynamic single-elimination tournaments, we obtain the binomial queue, a classical heap implementation, in a simple and natural way. We give four equivalent ways of representing heaps arising from tour- naments, and we obtain two new variants of binomial queues, a one-tree version and a one-pass version. We extend the one-pass version to support key decrease operations, obtaining the rank- pairing heap, or rp-heap. Rank-pairing heaps combine the performance guarantees of Fibonacci heaps with simplicity approaching that of pairing heaps. Like pairing heaps, rank-pairing heaps consist of trees of arbitrary structure, but these trees are combined by rank, not by list position, and rank changes, but not structural changes, cascade during key decrease operations. 1 Introduction A meldable heap (henceforth just a heap) is a data structure consisting of a set of items, each with a real-valued key, that supports the following operations: – make-heap: return a new, empty heap. – insert(x; H): insert item x, with predefined key, into heap H. – find-min(H): return an item in heap H of minimum key.
    [Show full text]
  • Lecture 6 1 Overview (BST) 2 Splay Trees
    ICS 621: Analysis of Algorithms Fall 2019 Lecture 6 Prof. Nodari Sitchinava Scribe: Evan Hataishi, Muzamil Yahia (from Ben Karsin) 1 Overview (BST) In the last lecture we reviewed on Binary Search Trees (BST), and used Dynamic Programming (DP) to design an optimal BST for a given set of keys X = fx1; x2; : : : ; xng with frequencies Pn F = ff1; f2; : : : ; fng i.e. f(xi) = fi. By optimal we mean the minimal cost of m = i=1 fi queries Pn 2 given by Cm = i=1 fi · ci where ci is the cost for searching xi. Knuth designed O(n ) algorithm to construct an optimal BST on X given F [3]. In this lecture, we look at self-balancing data structures that optimize themselves without knowing in advance what those frequencies are and also allow update operations on the BST, such as insertions and deletions. Pn Note that if we define p(x) = f(x)=m, then Cm = i=1 fi · ci is bounded from above by m(Hn + 2) Hn (a proof is given by Knuth [4, Section 6.2.2]), and bounded from below by m log 3 (proved by Pn 1 Mehlhorn [5]) where Hn is the entropy of n probabilities given by Hn = p(xi) log and i=1 p(xi) therefore Cm 2 Θ(mHn). Since entroby is bounded by Hn ≤ log n (see Lemma 2.10.1 from [1]), we have Cm = O(m log n). 2 Splay Trees Splay trees, detailed in [6] are self-adjusting BSTs that: • implement dictionary APIs, • are simply structured with limited overhead, • are c-competitive 1 with best offline BST, i.e., when all queries are known in advance, • are conjectured to be c-competitive with the best online BST, i.e., when future queries are not know in advance 2.1 Splay Tree Operations Simply, Splay trees are BSTs that move search targets to the root of the tree.
    [Show full text]
  • Improved Bounds for Multipass Pairing Heaps and Path-Balanced Binary Search Trees
    Improved Bounds for Multipass Pairing Heaps and Path-Balanced Binary Search Trees Dani Dorfman Blavatnik School of Computer Science, Tel Aviv University, Israel [email protected] Haim Kaplan1 Blavatnik School of Computer Science, Tel Aviv University, Israel [email protected] László Kozma2 Eindhoven University of Technology, The Netherlands [email protected] Seth Pettie3 University of Michigan [email protected] Uri Zwick4 Blavatnik School of Computer Science, Tel Aviv University, Israel [email protected] Abstract We revisit multipass pairing heaps and path-balanced binary search trees (BSTs), two classical algorithms for data structure maintenance. The pairing heap is a simple and efficient “self- adjusting” heap, introduced in 1986 by Fredman, Sedgewick, Sleator, and Tarjan. In the multi- pass variant (one of the original pairing heap variants described by Fredman et al.) the minimum item is extracted via repeated pairing rounds in which neighboring siblings are linked. Path-balanced BSTs, proposed by Sleator (cf. Subramanian, 1996), are a natural alternative to Splay trees (Sleator and Tarjan, 1983). In a path-balanced BST, whenever an item is accessed, the search path leading to that item is re-arranged into a balanced tree. Despite their simplicity, both algorithms turned out to be difficult to analyse. Fredman et al. showed that operations in multipass pairing heaps take amortized O(log n·log log n/ log log log n) time. For searching in path-balanced BSTs, Balasubramanian and Raman showed in 1995 the same amortized time bound of O(log n · log log n/ log log log n), using a different argument.
    [Show full text]
  • Masterarbeit Hartmann
    Master's thesis Efficiency of Self-Adjusting Heap Variants LM Hartmann Freie Universit¨atBerlin October 2020 E-Mail address: [email protected] Student Number: 5039594 Target Degree: Master of Science Subject: Computer Science Advisor: Prof. Dr. L´aszl´oKozma 2nd Advisor: Prof. Dr. Wolfgang Mulzer This copy is an edited version of the thesis that was submitted originally. The following changes were made: Minor spelling corrections in the \Acknowledgements" section; reformatting from US letter format to A4 format. Statement of authorship I declare that this document and the accompanying code has been composed by myself, and de- scribes my own work, unless otherwise acknowledged in the text. It has not been accepted in any previous application for a degree. All verbatim extracts have been distinguished by quotation marks, and all sources of information have been specifically acknowledged. Maria Hartmann Date i Acknowledgements This has been a challenging time to write my thesis, and I would like to thank all those who have supported me along the way. First and foremost, I wish to thank Prof. Dr. L´aszl´oKozma, who has been a wonderfully encouraging and supportive thesis supervisor. He has also suggested that I flaunt my Latin skills more, so Gratias tibi propter opportunitatem auxiliumque volo. Furthermore, I would like to thank the friends and family members who have listened to me during moments of frustration, and whose steadfast encouragement has been immensely helpful. I am also grateful to everyone who has read previous drafts of my thesis and offered suggestions. ii Contents Abstract vi 1 Introduction and Outline 1 2 Description of Heaps 3 2.1 Pairing Heap Variants .
    [Show full text]
  • Pairing Heap
    Pairing Heap Hauke Brinkop and Tobias Nipkow February 23, 2021 Abstract This library defines three different versions of pairing heaps: afunc- tional version of the original design based on binary trees [1], the ver- sion by Okasaki [2] and a modified version of the latter that is free of structural invariants. The amortized complexities of these implementations are analyzed in the AFP article Amortized Complexity. Contents 1 Pairing Heap in Binary Tree Representation 1 1.1 Definitions ............................. 2 1.2 Correctness Proofs ........................ 2 1.2.1 Invariants ......................... 2 1.2.2 Functional Correctness .................. 3 2 Pairing Heap According to Okasaki 5 2.1 Definitions ............................. 5 2.2 Correctness Proofs ........................ 6 2.2.1 Invariants ......................... 6 2.2.2 Functional Correctness .................. 6 3 Pairing Heap According to Oksaki (Modified) 8 3.1 Definitions ............................. 8 3.2 Correctness Proofs ........................ 9 3.2.1 Invariants ......................... 9 3.2.2 Functional Correctness .................. 10 1 Pairing Heap in Binary Tree Representation theory Pairing-Heap-Tree imports HOL−Library:Tree-Multiset HOL−Data-Structures:Priority-Queue-Specs begin 1 1.1 Definitions Pairing heaps [1] in their original representation as binary trees. fun get-min :: 0a :: linorder tree ) 0a where get-min (Node - x -) = x fun link :: ( 0a::linorder) tree ) 0a tree where link (Node hsx x (Node hsy y hs)) = (if x < y then Node (Node
    [Show full text]
  • COS 423 Lecture 7 Rank-Pairing Heaps
    COS 423 Lecture 7 Rank -Pairing Heaps ©Robert E. Tarjan 2011 Inefficiency in pairing heaps: Links during inserts , melds , decrease-keys can combine two trees of very different sizes, increasing the potential by Θ(lg n). To avoid such bad links, we use ranks : Each node has a non-negative integer rank Missing nodes have rank –1 Rank of root = 1 + rank of left child rank of tree = rank of root Store ranks with nodes. Link only roots of equal rank. Rank of winning root increases by one: r r r + 1 10 + 8 8 r 10 A B A B A heap is a set of half trees, not just one half tree (can’t link half trees of different ranks) Representation of a heap: a circular singly-linked list of roots of half trees, with the root of minimum key (the min -root ) first Circular linking → catenation takes O(1) time Node ranks depend on how operations are done find-min : return item in min-root make-heap : return a new, empty list insert : create new one-node half tree of rank 0, insert in list, update min-root 5 + 4 9 6 8 4 5 9 6 8 meld : catenate lists, update min-root 4 9 6 + 2 8 5 7 2 9 6 4 8 5 7 delete-min : Delete min-root. Cut edges along right path down from new root to give new half-trees. Link roots of equal rank until no two roots have equal rank. Form list of remaining half trees. To do links, use array of buckets, one per rank.
    [Show full text]