斐波那契堆

版权声明:本文为博主原创文章,未经允许不得转载。 https://blog.csdn.net/Koala_Tree/article/details/77971184

斐波那契堆


作者大树先生
博客http://blog.csdn.net/koala_tree
GitHubhttps://github.com/koalatree
2017 年 09 月 13 日


自《算法导论》.

斐波那契堆有两种用途:第一种,支持一系列操作,这些操作构成了所谓的“可合并堆”。第二种,斐波那契堆的一些操作可以在常数摊还时间内完成。

可合并堆的两种实现方式下各操作的运行时间。在操作时堆中的项数用n表示。

这里写图片描述

一、斐波那契堆结构

一个斐波那契堆是一系列具有最小堆序的有根树的集合。也就是说,每棵树均遵循最小堆性质:每个结点的关键字大于或等于它的父结点的关键字。

这里写图片描述

结点属性

  • x.p:指向其父结点的指针;
  • x.child:指向其某一个孩子结点的指针;x结点所有的孩子链接成一个环形双向链表,称为孩子链表;
  • y.left、y.right:每个孩子y包含,分别指向y的左兄弟和右兄弟,若只有一个孩子,则y.left=y.right=y;
  • x.degree:结点x的孩子链表中的孩子树木;
  • x.mark:结点x自从上一次成为另一个结点的孩子后,是否失去过孩子;
  • H.min:用来访问给定的斐波那契堆,指向具有最小关键字的树的根结点,称为最小结点;
  • H.n:表示H中当前结点数目。

势函数

使用势方法来分析斐波那契堆操作的性能。

  • t(H):表示H中跟链表中树的数目;
  • m(H):表示H中已标记的结点数目;
  • 势函数Φ(H)
    Φ(H)=t(H)+2m(H)

最大度数

在一个n个结点的斐波那契堆中任何结点的最大度数都有上界D(n)lgn.

二、可合并堆操作

1. 创建一个新的斐波那契堆

过程分配并返回一个斐波那契堆对象H,其中H.n=0和H.min = NIL,H中不存在树。

由于t(H)=0,m(H)=0,空斐波那契堆的势为Φ(H)=0,因此创建新堆的摊还代价等于他的实际代价O(1)

2. 插入一个结点

Fib_heap_insert(H, x)
    x.degree = 0
    x.p = NIL
    x.child = NIL
    x.mark = FALSE
    if H.min == NIL
        create a root list for H containing just x
        H.min = x
    else insert x into H's root list
        if x.key < H.min.key
            H.min = x
    H.n = H.n + 1

这里写图片描述

t(H)=t(H)+1,m(H)=m(H),则增加的势能为1,实际代价为O(1),摊还代价为O(1)+1=O(1).

3. 寻找最小结点

通过指针H.min得到。可以在O(1)的实际代价内找到最小结点。

4. 两个斐波那契堆的合并

合并斐波那契堆H1,H2,销毁这两个堆。

Fib_heap_union(H1, H2)
    H = Make_fib_heap()
    H.min = H1.min
    concatenate the root list of H2 with the root list of H
    if (H1.min == NIL) or (H2.min /= NIL and H2.min.key < H1.min.key)
        H.min = H2.min
    H.n = H1.n + H2.n
    return H

势函数变化为0,所以摊还代价等于实际代价O(1)

5. 抽取最小结点

Fib_heap_extract_min(H)
    z = H.min
    if z /= NIL
        for each child x of z
            add x to the root list of H
            x.p = NIL
        remove z from the root list of H
        if z == z.right
            H.min = NIL
        else H.min = z.right
            Consolidate(H)
        H.n = H.n - 1
    return z

以上代码用到的合并(Consolidating)H根链表的操作,通过调用Consolidating(H)来减少斐波那契堆中树的数目。

过程重复执行以下步骤:

  • 在根链表中找到两个具有相同度数的根x和y,其中,x.key <= y.key;
  • 把y链接到x:从根链表中移除y,调用Fib_heap_link过程,使y称为x的孩子,该过程将x.degree属性增加1,并清除y上的标记。

使用辅助数组A[0..D(H.n)]记录根结点对应的度数的轨迹。如A[i] = y,那么当前的y是一个具有y.degree = i的根。

Consolidate(H)
    let A[0..D(H.n)] be a new array
    for i = 0 to D(H.n)
        A[i] = NIL
    for each node w in the root list of H
        x = w
        d = x.degree
        while A[d] /= NIL
            y = A[d]  #another node with the same degree as x
            if x.key > y.key
                exchange x with y
            Fib_heap_link(H, y, x)
            A[d] = NIL
            d = d + 1
        A[d] = x
    H.min = NIL
    for i = 0 to D(H.n)
        if A[i] /= NIL
            if H.min == NIL
                create a root list for H containing just A[i]
                H.min = A[i]
            else insert A[i] into H's root list
                if A[i].key < H.min.key
                    H.min = A[i]
Fib_heap_link(H, y, x)
    remove y from the root list of H
    make y a child of x, incrementing x.degree
    y.mark = FALSE

这里写图片描述

这里写图片描述

抽取最小结点的摊还代价为O(D(n))=O(lgn).

三、关键字减值和删除一个结点

1. 关键字减值

摊还时间:O(1).

Fib_heap_decrease_key(H, x, k)
    error "new key is greater than current key"
    x.key = k
    y = x.p
    if y /= NIL and x.key < y.key
        Cut(H, x, y)
        Cascading_cut(H, y)
    if x.key < H.min.key
        H.min = x
Cut(H, x, y)
    remove x from the child list of y, decrementing y.degree
    add x to the root list of H
    x.p = NIL
    x.mark = FALSE
Cascading_cut(H, y)
    z = z.p
    if z /= NIL
        if y.mark == FALSE
            y.mark = TRUE
        else Cut(H, y, z)
            Cascading_cut(H, z)

斐波那契堆中规定,某个结点x一旦失掉第二个孩子,就切断x与其父结点的链接,是它称为一个新的跟。

所以,如果切掉的结点是其父结点的第二个孩子,则需要进行一次级联切断(cascading cut)。

过程示意图:

这里写图片描述

2. 删除一个结点

摊还时间:O(D(n))=O(lgn).

Fib_heap_delete(H, x)
    Fib_heap_decrease_key(H, x, -inf)
    Fib_heap_extract_min(H)

O(1)+O(D(n))=O(D(n))=O(lgn).

展开阅读全文

没有更多推荐了,返回首页