堆排序
堆排序是利用堆的性质进行的一种选择排序。下面先讨论一下堆。
1.堆
堆实际上是一棵完全二叉树,其任何一非叶节点满足性质:
Key[i]<=key[2i+1]&&Key[i]<=key[2i+2]或者Key[i]>=Key[2i+1]&&key>=key[2i+2]
即任何一非叶节点的关键字不大于或者不小于其左右孩子节点的关键字。
堆分为大顶堆和小顶堆,满足Key[i]>=Key[2i+1]&&key>=key[2i+2]称为大顶堆,满足 Key[i]<=key[2i+1]&&Key[i]<=key[2i+2]称为小顶堆。由上述性质可知大顶堆的堆顶的关键字肯定是所有关键字中最大的,小顶堆的堆顶的关键字是所有关键字中最小的。
2.堆排序的思想
利用大顶堆(小顶堆)堆顶记录的是最大关键字(最小关键字)这一特性,使得每次从无序中选择最大记录(最小记录)变得简单。
其基本思想为(大顶堆):
1)将初始待排序关键字序列(R1,R2....Rn)构建成大顶堆,此堆为初始的无序区;
2)将堆顶元素R[1]与最后一个元素R[n]交换,此时得到新的无序区(R1,R2,......Rn-1)和新的有序区(Rn),且满足R[1,2...n-1]<=R[n];
3)由于交换后新的堆顶R[1]可能违反堆的性质,因此需要对当前无序区(R1,R2,......Rn-1)调整为新堆,然后再次将R[1]与无序区最后一个元素交换,得到新的无序区(R1,R2....Rn-2)和新的有序区(Rn-1,Rn)。不断重复此过程直到有序区的元素个数为n-1,则整个排序过程完成。
操作过程如下:
1)初始化堆:将R[1..n]构造为堆;
2)将当前无序区的堆顶元素R[1]同该区间的最后一个记录交换,然后将新的无序区调整为新的堆。
因此对于堆排序,最重要的两个操作就是构造初始堆和调整堆,其实构造初始堆事实上也是调整堆的过程,只不过构造初始堆是对所有的非叶节点都进行调整。
下面举例说明:
给定一个整形数组a[]={16,7,3,20,17,8},对其进行堆排序。
首先根据该数组元素构建一个完全二叉树,得到
20和16交换后导致16不满足堆的性质,因此需重新调整
这样就得到了初始堆。
此时3位于堆顶不满堆的性质,则需调整继续调整
public class SortArray {
/*堆排序(大顶堆) 2011.9.14*/
public static void swap(int[] arr, int a, int b) {
int tmp = arr[a];
arr[a] = arr[b];
arr[b] = tmp;
}
public static void heapAdjust(int[] arr, int i, int size) {
int lchild = 2 * i + 1;//i的左孩子节点序号
int rchild = 2 * i + 2;//i的右孩子节点序号
int max = i;//临时变量
if (i < size / 2) { //如果i是叶节点就不用进行调整
if (lchild < size && arr[lchild] > arr[max])
max = lchild;
if (rchild < size && arr[rchild] > arr[max])
max = rchild;
if (max != i) {
swap(arr, max, i);
heapAdjust(arr, max, size);
}
}
}
public static void buildHeap(int[] arr) {
int size = arr.length;
for (int i = size / 2 - 1; i >= 0; i--) {
heapAdjust(arr, i, size);
}
}
public static void heapSort(int[] arr) //堆排序
{
int i;
int size = arr.length;
buildHeap(arr);
for (i = size - 1; i >= 0; i--) {
swap(arr, 0, i); //交换堆顶和最后一个元素,即每次将剩余元素中的最大者放到最后面
heapAdjust(arr, 0, i); //重新调整堆顶节点成为大顶堆
}
}
public static void main(String[] args) {
int[] arr = {4, 6, 3, 9, 1};
heapSort(arr);
for (int i = 0; i < arr.length; i++) {
System.out.print(arr[i] + " ");
}
}
}
c++版本:
/*堆排序(大顶堆) 2011.9.14*/
#include <iostream>
#include<algorithm>
using namespace std;
//----------------duipaixu--------------------------------------
void buildheap(int a[], int k, int size)
{
int l = 2 * k + 1;
int r = 2 * k + 2;
if (k<size/2)
{
int max = k;
if (l<size&&a[l]>a[k])
max = l;
if (r<size&&a[r]>a[max])
max = r;
if (max != k)
{
int tmp = a[k];
a[k] = a[max];
a[max] = tmp;
buildheap(a, max, size);
}
}
}
void buildall(int a[], int size)
{
for (int i = (size) / 2 - 1; i >= 0; i--)
{
buildheap(a, i, size);
}
}
void duipx(int a[], int size)
{
buildall(a, size);
for (int i = size-1; i >=0; i--)
{
int tmp = a[i];
a[i] = a[0];
a[0] = tmp;
buildheap(a, 0, i);
}
}
int main()
{
int a[] = { 5, 4, 7, 1, 9, 3,4 };
//kp(a, 0, 6);
duipx(a, 7);
for (int i = 0; i < 7; i++)
{
cout << a[i] << " ";
}
system("pause");
return 0;
}