每日一“酷”之heapq

作用:heapq模块实现一个适用于Python列表的最小堆排序算法

　　堆（heap）是一个属性数据结构，其中子节点与父节点是一种有序关系。二叉堆是一种特殊的堆，二叉堆是完全二元树（二叉树）或者是近似完全二元树（二叉树）。可以使用以下如下方式组织的列表或数表示，即元素N的子元素位于2*N+1和2*N+2。这种布局允许原地重新组织堆，从而不必再增加或删除元素时分配大量内存。

      最大堆确保父节点大于或等于其两个子节点。最小堆要求父节点小雨或等于其子节点。Python的heqpq模块实现了一个最小堆。

1、示例数据

数据存储在：heapq_heapdata.py中

data = [19,9,4,10,11]

堆的输出存储在heapq_showtree.py 

import math
import cStringIO

def show_tree(tree,total_width=36,fill=' '):
    output = cStringIO.StringIO()
    last_row = -1
    for i,n in enumerate(tree):
        if i :
            row = int(math.floor(math.log(i+1,2)))
        else:
            row = 0
        if row != last_row:
            output.write('
')
        columns = 2 ** row
        col_width = int(math.floor((total_width * 1.0) / columns))
        output.write(str(n).center(col_width,fill))
        last_row = row
    print output.getvalue()
    print '-' * total_width
    print 
    return 

2、创建堆

创建堆的两种基本方式：heappush()和heapify() 

import heapq
from heapq_heapdata import data
from heapq_showtree import show_tree

heap = []
print 'random : ' , data
print 
for n in data:
    print 'add %3d:' % n
    heapq.heappush(heap,n)
    show_tree(heap)

运行结果：

使用heappush()时，从数据源增加新元素时会保持元素的堆顺序

如果数据已经存在内存中，使用heaoify()原地重新组合字列表中的元素会更高效

import heapq
from heapq_heapdata import data
from heapq_showtree import show_tree

heap = []
print 'random : ' , data
heapq.heapify(data)
print 'heapifed : '
show_tree(data)

运行结果：

如果按堆顺序一次一个元素构建列表，其结果与构建一个无序列表在调用heapify()是一样的。

3、访问堆内容

一旦堆已正确组织，就可以使用heappop()删除有最小值的元素

import heapq
from heapq_heapdata import data
from heapq_showtree import show_tree

print 'random    :', data
heapq.heapify(data)
print 'heapifed : '
show_tree(data)
print 
for i in xrange(2):
    smallest = heapq.heappop(data)
    print 'pop       %3d:' % smallest
    show_tree(data)

运行结果：

这个例子是由stdlib文档改写的，其中使用了heapify()和heappop()对一个数组列表排序。

如果希望在一个操作中删除现有元素并替换新值，可以使用heapreplace()。

import heapq
from heapq_heapdata import data
from heapq_showtree import show_tree

heapq.heapify(data)
print 'start:'
show_tree(data)
for n in [0,13]:
    smallest = heapq.heapreplace(data, n)
    print 'replace %2d with %3d' % (smallest,n)
    show_tree(data)

运行结果：

通过原地替换元素，这样可以维持一个固定大小的堆，如按优先级排序的作业队列。

4、堆的数据极值

Heapq还包括两个检查可迭代对象的函数，查找其中包含最大值或最小值的范围

import heapq
from heapq_heapdata import data
print 'all       :',data
print '3 largest :',heapq.nlargest(3,data)
print 'from sort :',list(reversed(sorted(data)[-3:]))
print '3 smallest:',heapq.nsmallest(3,data)
print 'from sort :',sorted(data)[:3]

运行结果：

只有当n值（n>1）相对小时使用nlargest()和nsmallest()才算高效，不过有些情况下这两个函数会很方便。