二项堆python实现——eager binomial heap

关于eager binomial heaps的内容详见算法设计搬运(3)——Heaps,更多细节在后面的Lazy Binomial Heaps也有。比较懒用的双向链表写的。双向链表可以实现 O(\log n)的 decreaseKey 和 deleteMin,可以参考石溪的ppt,尽管它是用array来存heap指针。有的ppt就错得离谱,单向链表decreaseKey就不只是在树的一个路径上比较了,没法O(1)找它的left。

资源已上传,新手上路,就只当为熟悉python指针练手...

# author: zyx
# eager binomial heaps
# merge,insert,find_min,extractMin

class Binomial_Heap():
    def __init__(self):
        self.head = None
        self.min = None

    def merge(self,H2):
        if not self.min:
            self.min = H2.min
        elif not H2.min:
            return
        else:
            if self.min.data > H2.min.data:
                self.min = H2.min
        h1 = self.head
        h2 = H2.head
        unique_order_heap = Binomial_Heap()
        c = unique_order_heap.head
        while h1 and h2:
            if h1.order == h2.order:
                if h1.data <= h2.data:
                    h1.order+=1
                    h2.pre = None
                    d = h2
                    h2 = h2.sibling
                    d.sibling = h1.Lchild
                    d.parent = h1
                    h1.Lchild = d
                    if unique_order_heap.head == None:
                        unique_order_heap.head = h1
                        c = h1
                        h1.pre = None
                        h1 = h1.sibling
                    else:
                        h1.pre = c
                        c.sibling = h1
                        h1 = h1.sibling
                        c = c.sibling
                else:
                    h2.order += 1
                    h1.pre = None
                    d = h1
                    h1 = h1.sibling
                    d.sibling = h2.Lchild
                    d.parent = h2
                    h2.Lchild = d
                    if unique_order_heap.head == None:
                        unique_order_heap.head = h2
                        c = h2
                        h2.pre = None
                        h2 = h2.sibling
                    else:
                        h2.pre = c
                        c.sibling = h2
                        h2 = h2.sibling
                        c = c.sibling
            elif h1.order < h2.order:
                if unique_order_heap.head == None:
                    unique_order_heap.head = h1
                    c = unique_order_heap.head
                    h1.pre = None
                    h1 = h1.sibling
                else:
                    if h1.order == c.order:
                        if c.data <= h1.data:
                            c.order += 1
                            d = h1
                            h1.pre = None
                            h1 = h1.sibling
                            d.sibling = c.Lchild
                            c.Lchild = d
                            d.parent = c
                        else:
                            if c.pre == None: # single B
                                unique_order_heap.head = h1
                                h1.pre = None
                            else:
                                c.pre.sibling = h1
                                h1.pre = c.pre
                            h1.order += 1
                            c.sibling = h1.Lchild
                            h1.Lchild = c
                            c.parent = h1
                            c = h1
                            h1 = h1.sibling
                    else:
                        c.sibling = h1
                        h1.pre = c
                        h1 = h1.sibling
                        c = c.sibling
            else:
                if unique_order_heap.head == None:
                    unique_order_heap.head = h2
                    c = unique_order_heap.head
                    h2.pre = None
                    h2 = h2.sibling
                else:
                    if h2.order == c.order:
                        if c.data <= h2.data:
                            d = h2
                            c.order += 1
                            h2 = h2.sibling
                            d.sibling = c.Lchild
                            c.Lchild = d
                            d.parent = c
                        else:
                            if c.pre == None: # single B
                                unique_order_heap.head = h2
                                h2.pre = None
                            else:
                                c.pre.sibling = h2
                                h2.pre = c.pre
                            h2.order += 1
                            c.sibling = h2.Lchild
                            h2.Lchild = c
                            c.parent = h2
                            c = h2
                            h2 = h2.sibling
                    else:
                        c.sibling = h2
                        h2.pre = c
                        h2 = h2.sibling
                        c = c.sibling
        if c:
            c.sibling = None
        while h1:
            if c:
                if h1.order == c.order:
                    if c.data <= h1.data:
                        c.order += 1
                        d = h1
                        h1.pre = None
                        h1 = h1.sibling
                        d.sibling = c.Lchild
                        c.Lchild = d
                        d.parent = c
                    else:
                        if c.pre == None:  # single B
                            unique_order_heap.head = h1
                            h1.pre = None
                        else:
                            c.pre.sibling = h1
                            h1.pre = c.pre
                        h1.order += 1
                        c.sibling = h1.Lchild
                        h1.Lchild = c
                        c.parent = h1
                        c = h1
                        h1 = h1.sibling
                else:
                    c.sibling = h1
                    h1.pre = c
                    break
            else:
                unique_order_heap.head = h1
                h1.pre = None
                break
        while h2:
            if c:
                if h2.order == c.order:
                    if c.data <= h2.data:
                        d = h2
                        c.order += 1
                        h2 = h2.sibling
                        d.sibling = c.Lchild
                        c.Lchild = d
                        d.parent = c
                    else:
                        if c.pre == None:  # single B
                            unique_order_heap.head = h2
                            h2.pre = None
                        else:
                            c.pre.sibling = h2
                            h2.pre = c.pre
                        h2.order += 1
                        c.sibling = h2.Lchild
                        h2.Lchild = c
                        c.parent = h2
                        c = h2
                        h2 = h2.sibling
                else:
                    c.sibling = h2
                    h2.pre = c
                    break
            else:
                unique_order_heap.head = h2
                h2.pre = None
                break
        self.head = unique_order_heap.head

    def find_Min(self):
        return self.min

    def insert(self,item):
        node = Binomial_Trees(item)
        H2 = Binomial_Heap()
        H2.head = node
        H2.min = node
        self.merge(H2)

    def extractMin(self):
        m = self.min
        if m.pre:
            m.pre.sibling = m.sibling
        else:
            self.head = m.sibling
        if m.sibling:
            m.sibling.pre = m.pre
        new_heap = Binomial_Heap()
        to_transform = m.Lchild
        minimum = to_transform
        while to_transform:
            if minimum.data > to_transform.data:
                minimum = to_transform
            to_transform.parent = None
            if not new_heap.head:
                new_heap.head = to_transform
                to_transform = to_transform.sibling
                new_heap.head.sibling = None
            else:
                d = to_transform
                to_transform = to_transform.sibling
                d.sibling = new_heap.head
                new_heap.head.pre = d
                new_heap.head = d
        new_heap.min = minimum
        cur = self.head
        minimum = self.head
        while cur:
            if minimum.data < cur.data:
                minimum = cur
            cur = cur.sibling
        self.min = minimum
        self.merge(new_heap)
        return m.data

class Binomial_Trees():
    def __init__(self,item):
        self.order = 0
        self.parent = None
        self.sibling = None
        self.Lchild = None
        self.data = item
        self.pre = None

if __name__ == '__main__':
    a = Binomial_Heap()
    a.insert(1)
    a.insert(2)
    a.insert(3)
    a.insert(4)
    a.insert(-1)
    a.extractMin()
    a.extractMin()
    b = Binomial_Heap()
    b.insert(2)
    b.insert(5)
    b.insert(0)
    b.insert(3)
    a.merge(b)
    a.insert(-1)
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 217,657评论 6 505
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,889评论 3 394
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 164,057评论 0 354
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,509评论 1 293
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,562评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,443评论 1 302
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,251评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,129评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,561评论 1 314
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,779评论 3 335
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,902评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,621评论 5 345
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,220评论 3 328
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,838评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,971评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,025评论 2 370
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,843评论 2 354