堆操作
实验目的
(一)建堆:将数组A[1..n]变成一个最大堆。
(二)堆排序:将一个堆中的元素按递减排序输出。
(三)用插入方法建堆:堆大小从1到n每次插入一个元素到堆中,直到n个元素入堆。
实验原理
(二叉)堆是一个数组,它可以被看成一个近似的完全二叉树。树上的每一个结点对应数组中的一个元素。除了最底层外,该树是完全充满的,而且是从左向右填充。表示堆的数组A包括两个属性:A.length(通常)给出数组元素的个数,A.heap-size表示有多少个堆元素存储在该数组中。也就是说,虽然A[1..A.length]可能都存有数据,但是只有A[1..A.heap-size]中存放的是堆的有效元素,这个,0 <=A.heap-size <= A.length 。树的根节点是A[1]。
二叉堆可以分为两种形式:最大堆和最小堆。在这两种堆中,结点的值都要满足堆的性质,但一些细节定义则有所差异。在最大堆中,最大堆性质是指除了根以外的所有结点i都要满足:A[PARENT(i)] >= A[i] ,也就是说,某个结点的值至多与其父结点一样大。因此,堆中的最大元素存放在根结点中;并且,在任一子树中,该子树所包含的所有结点的值都不大于该子树根节点的值。最小堆的组织方式正好相反:最小堆性质是指除了根以外的所有结点i都有A[PARENT(i)] <= A[i] ,最小堆中的最小元素存放在根结点中。
实验过程
(一)维护堆的性质
MAX-HEAPIFY是用于维护最大堆性质的重要过程。它的输入为一个数组A和一个下标i。在调用MAX-HEAPIFY的时候,对于违背了最大堆的性质的A[i]使其在最大堆中“逐级下降”,从而使得以下标i为根结点的子树重新遵循最大堆的性质。
MAX-HEAPIFY(A, i)
1 l = LEFT(i)
2 r = RIGHT(i)
3 if l <= A.heap-size andA[l] > A[i]
4 largest = l
5 else largest = i
6 if r <= A.heap-size andA[r] > A[largest]
7 largest = r
8 if largest != i
9 exchange A[i] with A[largest]
10 MAX-HEAPIFY(A, largest)
(二)建堆
我们可以用自底向上的方法利用过程MAX-HEAPIFY把一个大小为n = A.length的数组A[1..n]转换为最大堆。子数组A(|n/2|+1..n)中的元素都是树的叶节点。对每一个叶节点都调用MAX-HEAPIFY进行调整,从而满足最大堆的性质。
BUILD-MAX-HEAP
1 A.heap-size = A.length
2 for i = |A.length/2| downto1
3 MAX-HEAPIFY(A, i)
(三)堆排序
由于数组中的最大元素总是在根结点A[1]中,通过把它与A[n]进行互换,我们可以让该元素放到正确的位置。这个时候,我们从堆中去掉结点n(即通过减少A.heap-size的值来实现),剩余的节点中,原来根的孩子仍然是最大堆,而新的节点可能违背最大堆的性质,为了维护最大堆的性质,调用一次MAX-HEAPIFY(A, 1),从而在A[1..n-1]上构造一个新的最大堆。
HEAPSORT(A)
1 BUILD-MAX-HEAP(A)
2 for i = A.length downto 2
3 exchange A[1] with A[i]
4 A.heap-size = A.heap-size - 1
5 MAX-HEAPIFY(A, 1)
(四)插入方法建堆
MAX-HEAP-INSERT能够实现INSERT操作。它的输入是要被插入到最大堆A中的新元素的关键字。MAX-HEAP-INSERT首先通过增加一个关键字为-∞的叶节点来扩展最大堆,然后调用HEAP-INCREASE-KEY为新节点设置对应的关键字,同时保持最大堆的性质。
MAX-HEAP-INSERT(A, key)
1 A.heap-size = A.heap-size +1
2 A[A.heap-size] = -∞
3 HEAP-INCREASE-KEY(A,A.heap-size, key)
HEAP-INCREASE-KEY(A, i, key)
1 if key < A[i]
2 error“new key is smaller than current key”
3 A[i] = key
4 while i > 1 and A[PARENT(i)]< A[i]
5 exchange A[i] with A[PARENT(i)]
6 i = PARENT(i)
实验总结
(一)将整个实验过程完成后,深入理解了二叉堆的性质。堆的优点:在大多数计算机上,在实现PARENT(i),LEFT(i),RIGHT(i)过程中,可以通过较少的指令(左移,右移)进行计算。
(二)在建堆的过程中,有一个巧妙的地方:就是对于A[|A.length/2|..1]的叶节点进行最大堆的性质维护,而不是遍历整个堆。
(三)堆排序的实际应用场景:基于最大堆实现最大优先队列(记录将要执行的各个作业以及它们之间的相对优先级),从而应用到计算机系统的作业调度中。
附录(代码)
(一)建堆
#include <stdio.h>
#include <stdlib.h>
#define a_length 10
int a_heap_size;
void max_heapify(int a[], int i){
intl,r,largest;
int change;
l = 2 * i +1;
r = 2 * i +2;
if(l <=a_heap_size-1 && a[l] > a[i])
largest =l;
else largest= i;
if(r <=a_heap_size-1 && a[r] > a[largest])
largest =r;
if(largest !=i){
change = a[i];
a[i] = a[largest];
a[largest]= change;
max_heapify(a,largest);
}
}
void build_max_heap(int a[]){
int i;
a_heap_size =a_length;
for(i =a_length/2 - 1; i >= 0; i--)
max_heapify(a,i);
}
int main(){
inta[a_length] = {4,1,3,2,16,9,10,14,8,7};
int i;
printf("origin_heapa[]=");
for(i = 0; i< 10; i++)
printf("%d ",a[i]);
printf("\n");
build_max_heap(a);
printf(" max_heap a[]=");
for(i = 0; i< 10; i++)
printf("%d ",a[i]);
printf("\n");
}
(二)堆排序
#include <stdio.h>
#include <stdlib.h>
#define a_length 10
int a_heap_size;
void max_heapify(int a[], int i){
intl,r,largest;
int change;
l = 2 * i +1;
r = 2 * i +2;
if(l <=a_heap_size-1 && a[l] > a[i])
largest =l;
else largest= i;
if(r <=a_heap_size-1 && a[r] > a[largest])
largest =r;
if(largest !=i){
change = a[i];
a[i] = a[largest];
a[largest]= change;
max_heapify(a,largest);
}
}
void build_max_heap(int a[]){
int i,j;
a_heap_size =a_length;
for(i =a_length/2 - 1; i >= 0; i--)
max_heapify(a,i);
}
void heapsort(int a[]){
inti,exchange;
int j;
for(i =a_length-1; i >= 1; i--){
exchange = a[0];
a[0] = a[i];
a[i] = exchange;
a_heap_size= a_heap_size - 1;
max_heapify(a,0);
printf("\na_heap_size= %d i = %d \n",a_heap_size,i);
for(j = 0;j < 10; j++)
printf("%d ",a[j]);
printf("\n");
}
}
int main(){
int a[a_length] = {4,1,3,2,16,9,10,14,8,7};
int i;
printf("origin_heapa[]=");
for(i = 0; i< 10; i++)
printf("%d ",a[i]);
printf("\n");
build_max_heap(a);
printf(" max_heap a[]=");
for(i = 0; i< 10; i++)
printf("%d ",a[i]);
printf("\n");
heapsort(a);
printf("heap_sort a[]=");
for(i = 0; i< 10; i++)
printf("%d ",a[i]);
printf("\n");
}
(三)用插入法建堆
#include <stdio.h>
#include <stdlib.h>
#define a_length 10
int a_heap_size;
int insert_length;
int parent(int i){
int a;
a = (i-1)/2;
return a;
}
void heap_increase_key(int a[], int i, int key){
int exchange;
int a_end;
//if(key <a[i-1])
// printf("new key is smaller than currentkey!\n");
a_end = i - 1;
a[a_end] = key;
while(a_end> 0 && a[parent(a_end)] < a[a_end]){
exchange = a[a_end];
a[a_end] = a[parent(a_end)];
a[parent(a_end)]= exchange;
a_end = parent(a_end);
// printf("\n1111111111111\n");
}
}
void max_heap_insert(int a[], int key){
a_heap_size =a_heap_size + 1;
a[a_heap_size-1]= -1;
heap_increase_key(a,a_heap_size, key);
}
void build_max_heap(int a[]){
int i;
a_heap_size =1;
for(i = 1; i< insert_length; i++){
max_heap_insert(a,a[i]);
// printf("\n66666666666\n");
}
}
int main(){
inta[a_length];
int number,i;
printf("inputthe numbers:");
scanf("%d",&number);
//printf("\n");
//printf("number=%d\n",number);
insert_length= number;
printf("inputthe array:");
for(i = 0; i< number; i++)
scanf("%d",&a[i]);
printf("\n");
for(i = 0; i <number; i++)
printf("a[%d]=%d ",i,a[i]);
printf("\n");
build_max_heap(a);
printf("max_heapa[]= ");
for(i = 0; i< number; i++)
printf("%d ",a[i]);
printf("\n");
}