一、简介
堆数据结构的性质:堆:任何结点的值都小于或等于其孩子的值的完全二叉树为小根堆,任何结点的值都大于或等于其孩子的值的完全二叉树为大根堆。为了方便使用完全二叉树的性质,假定数组从下标1开始。这样:
leftChild = 2*i; rightChild = 2*i + 1; parent = i/2;
堆数据结构算法分析:
堆排序的最坏时间复杂度为O(nlogn)。堆序的平均性能较接近于最坏性能。由于建初始堆所需的比较次数较多,所以堆排序不适宜于记录数较少的文件。堆排序是就地排序,辅助空间为O(1)。堆排序是不稳定的。
堆数据结构实现如下所示:
#include"stdio.h" inline void swap(int &a,int &b) { int temp=a; a=b; b=temp; } void HeapAdjust(int array[],int i,int nLength)//自顶向下调整堆 { int nChild; int nTemp;//赋值为待调整的 节点 for(nTemp=array[i];2*i<nLength;i=nChild)//2*i<nLength说明还有左孩子 { nChild=2*i;//左孩子 /*一共两个子节点的话得到 较大的一个*/ if(nChild<nLength-1&&array[nChild+1]>array[nChild])//nChild<nLength-1 判断到头没有 ++nChild; /*如果较大子节点大于父节点 将子节点 调整到父节点*/ if(nTemp<array[nChild]) array[i]=array[nChild]; else break;//这个地方不加 会出错 第一个会输出第二个 array[nChild]=nTemp;//子节点 等于父节点 (不执行break) } } void HeapSort(int a[],int length) { /*初建堆 */ for(int i=length/2;i>0;--i)//从最后一个 非叶子节点调整 (这里的 i是下标) HeapAdjust(a,i,length); for(int i=length;i>1;--i) { swap(a[1],a[i]); /*第一个最大元素跟最后一个交换*/ HeapAdjust(a,1,i);//调整堆 (注意 length=i 由于堆是逐渐变小的) } } int main() { int a[10]={0,1,2,5,3,8,4,7,6}; HeapSort(a,8); for(int i=1;i<9;i++) printf("%d\n",a[i]); return 0; }
利用堆数据结构可实现优先级队列,优先队列是0个或多个元素的集合,每个元素都有一个优先权或值,对优先队列执行的操作有查找、插入一个新元素和删除。在最小优先队列(min priority queue)中,查找操作用来搜索优先权最小的元素,删除操作用来删除该元素;对于最大优先队列(max priority queue),查找操作用来搜索优先权最大的元素,删除操作用来删除该元素。优先权队列中的元素可以有相同的优先权,查找与删除操作可根据任意优先权进行。
优先级队列代码实现如下所示:
#include <iostream> using namespace std; template<class T> class priqueue { private: int n, maxsize; T *x; void swap(int &i, int &j)//根据坐标交换数组元素的值 { T t = i; i = j; j = t; } public: priqueue(int m)//初始化数组 { maxsize = m; x = new T[maxsize+1]; n = 0; } void insert(T t) { int i, p; x[++n] = t; //插入的元素放到最后 for (i = n; i > 1 && x[p=i/2] > x[i]; i = p) swap(x[p], x[i]); } T extractmin()//向下调整堆 { int i, c; T t = x[1]; x[1] = x[n--]; for (i = 1; (c=2*i) <= n; i = c) { if (c+1<=n && x[c+1]<x[c]) c++; if (x[i] <= x[c]) break; swap(x[c], x[i]); } return t; } void print(int n) { for (int i = 1; i < n; i++) //输出堆 cout << x[i] << " "; } }; template<class T> void pqsort(T v[], int n)//先初始化一个数组,然后插入建立一个堆 { priqueue<T> pq(n); int i; for (i = 0; i < n; i++) pq.insert(v[i]); cout<<"输出排序后的堆:"; pq.print(n); } int main() { const int n = 10; int i, v[n]; /*以下是通过向上调整堆 建立一个10个元素的堆*/ for (i = 0; i < n; i++) v[i] = n-i; pqsort(v, n); cout<<"\n执行插入和删除操作(输入0代表删除最小值,输入其他代表插入)"<<endl; priqueue<int> pq(100); int count=0; while (cin >> i) if (i == 0) { if(count) cout <<"删除的最小元素为:"<<pq.extractmin() << "\n"; else cout<<"请先插入元素"<<endl; } else { pq.insert(i); count++; } return 0; }
二、原理
(1)高效性。堆中所有结点和根结点之间相差的层数在logn之内。因此其执行效率很高,堆排序通过在同一数组中包含两种抽象结构(堆和元素序列)来避免使用额外空间。
(2)正确性,即堆数据结构的不变性质。
(3)抽象性。
(4)过程抽象。
(5)抽象数据类型。
三、习题
(1)习题1:实现基于堆的优先级队列,尽可能的提高运行速度。
为了提高向上调整堆的速度,在x[0]放置哨兵==当前插入的元素。省去了每次都判断i>1,向上调整堆结束:x[p] <= x[i]。代码如下:
void insert(T t) //向上调整堆 { int i, p; x[++n] = t; //插入的元素放到最后 x[0]=t; for (i = n; x[p=i/2] > x[i]; i = p) swap(x[p], x[i]); }
(2)习题4:如何使用优先级队列解决下列问题:
a、构件霍夫曼码;
b、计算大型浮点数集合的和;
c、在存有10亿个数的文件中找到最大的100W个数;
d、将多个较小的有序文件归并为一个较大的有序文件;
解答:
a、构造哈夫曼树时候,需要选取当前数组的两个最小值,删除两个最小值,并将计算之和插入原来数组,采用堆,初建堆,两次调用选取最小值的函数。计算之和之后,调用插入堆并调整堆;b、如果将较小浮点数和较大浮点数相加可能造成丢失精度。所以每次取最小的两个相加。然后将和插入数组集合。最后剩下一个就是所有浮点数的和;c、典型的topK;d、将所有小文件要插入的当前值组成一个堆。取堆最小值,插入排序数组。调整堆。然后插入该小文件下一个元素(无后继则不操作)。
(3)习题7:在一些计算机上,除以2求当前范围的中点是二分搜索程序中开销最大的部分。假设我们已经正确构建了待搜索的数组,说明如何使用乘以2的操作来替代除法。给出建立并搜索这样一个数组的算法。
修改后的二分搜索从i=1开始,每次迭代将i设置为2i或2i+1.元素x[1]包含中值,x[2]包含第一个四分位值,x[3]包含第三个四分位值,等等。一种能在O(n)时间内将n元有序数组调整为“堆搜索”顺序的算法。作为该方法的先驱,考虑把一个2^k-1元的有序数组a拷贝到一个“堆搜索”数组b中:a中奇数位的元素按顺序放到b的后半部分,模4余2位置的元素按顺序放到b中剩余部分的后半部分,等等。