【算法】三、堆排序

堆排序是一种选择排序，其时间复杂度为O(nlogn)。

堆的定义

　　n个元素的序列{k₁，k₂，…,k_n}当且仅当知足下列关系之一时，称之为堆。

　　情形1：k_i <= k_2i 且k_i <= k_2i+1 （最小化堆或小顶堆）

　　情形2：k_i >= k_2i 且k_i >= k_2i+1 （最大化堆或大顶堆）

　　其中i=1,2,…,n/2向下取整;

     

                 

 

　　若将和此序列对应的一维数组（即以一维数组做此序列的存储结构）当作是一个彻底二叉树，则堆的含义代表，彻底二叉树中全部非终端结点的值均不大于（或不小于）其左、右孩子结点的值。

　　由此，若序列{k₁，k₂，…,k_n}是堆，则堆顶元素（或彻底二叉树的根）必为序列中n个元素的最小值（或最大值）。

　　例如，下列两个序列为堆，对应的彻底二叉树如图：

　　

 

　　若在输出堆顶的最小值以后，使得剩余n-1个元素的序列重又建成一个堆，则获得n个元素的次小值。如此反复执行，便能获得一个有序序列，这个过程称之为堆排序。

　　堆排序（Heap Sort）只须要一个记录元素大小的辅助空间（供交换用），每一个待排序的记录仅占有一个存储空间。

 

堆的存储

　　通常用数组来表示堆，若根结点存在序号0处， i结点的父结点下标就为(i-1)/2。i结点的左右子结点下标分别为2*i+1和2*i+2。

　　（注：若是根结点是从1开始，则左右孩子结点分别是2i和2i+1。）

　　如第0个结点左右子结点下标分别为1和2。

　　如最大化堆以下：

 　　

 

　　左图为其存储结构，右图为其逻辑结构。

堆排序的实现

　　实现堆排序须要解决两个问题：

　　　　1.如何由一个无序序列建成一个堆？

　　　　2.如何在输出堆顶元素以后，调整剩余元素成为一个新的堆？

 

　　先考虑第二个问题，通常在输出堆顶元素以后，视为将这个元素排除，而后用表中最后一个元素填补它的位置，自上向下进行调整：首先将堆顶元素和它的左右子树的根结点进行比较，把最小的元素交换到堆顶；而后顺着被破坏的路径一路调整下去，直至叶子结点，就获得新的堆。

　　咱们称这个自堆顶至叶子的调整过程为“筛选”。

　　从无序序列创建堆的过程就是一个反复“筛选”的过程。

构造初始堆

　　初始化堆的时候是对全部的非叶子结点进行筛选。

　　最后一个非终端元素的下标是[n/2]向下取整，因此筛选只须要从第[n/2]向下取整个元素开始，从后往前进行调整。

　　好比，给定一个数组，首先根据该数组元素构造一个彻底二叉树。

　　而后从最后一个非叶子结点开始，每次都是从父结点、左孩子、右孩子中进行比较交换，交换可能会引发孩子结点不知足堆的性质，因此每次交换以后须要从新对被交换的孩子结点进行调整。

进行堆排序

　　有了初始堆以后就能够进行排序了。

　　堆排序是一种选择排序。创建的初始堆为初始的无序区。

　　排序开始，首先输出堆顶元素（由于它是最值），将堆顶元素和最后一个元素交换，这样，第n个位置（即最后一个位置）做为有序区，前n-1个位置还是无序区，对无序区进行调整，获得堆以后，再交换堆顶和最后一个元素，这样有序区长度变为2。。。

　　不断进行此操做，将剩下的元素从新调整为堆，而后输出堆顶元素到有序区。每次交换都致使无序区-1，有序区+1。不断重复此过程直到有序区长度增加为n-1，排序完成。

堆排序实例

 　　首先，创建初始的堆结构如图：

　　

　　而后，交换堆顶的元素和最后一个元素，此时最后一个位置做为有序区（有序区显示为黄色），而后进行其余无序区的堆调整，从新获得大顶堆后，交换堆顶和倒数第二个元素的位置……

　　

　　重复此过程：

　　

 

　　最后，有序区扩展完成即排序完成：

　　

 

　　由排序过程可见，若想获得升序，则创建大顶堆，若想获得降序，则创建小顶堆。

代码

　　假设排列的元素为整型，且元素的关键字为其自己。

　　由于要进行升序排列，因此用大顶堆。

　　根结点从0开始，因此i结点的左右孩子结点的下标为2i+1和2i+2。

 

void HeapAdjust(int H[], int start, int end) {

        int temp = H[start];
        for (int i = 2 * start + 1; i <= end; i = 2 * i + 1) {
            //由于假设根结点的序号为0而不是1，因此i结点左孩子和右孩子分别为2i+1和2i+2
            if (i < end && H[i] < H[i + 1])//左右孩子的比较
            {
                i++;//i为较大的记录的下标
            }

            if (temp > H[i])//左右孩子中获胜者与父亲的比较
            {
                break;
            }

            //将孩子结点上位，则以孩子结点的位置进行下一轮的筛选
            H[start] = H[i];
            start = i;

        }
        H[start] = temp; //插入最开始不和谐的元素
    }

    void HeapSort(int A[], int n) {
        //先创建大顶堆
        for (int i = n / 2 - 1; i >= 0; i--) {
            HeapAdjust(A, i, n - 1);
        }
        //进行排序
        for (int i = n; i > 0; i--) {
            //最后一个元素和第一元素进行交换
            int temp = A[i];
            A[i] = A[0];
            A[0] = temp;

            //而后将剩下的无序元素继续调整为大顶堆
            HeapAdjust(A, 0, i - 1);
        }
    }

 

堆排序分析

　　堆排序方法对记录数较少的文件并不值得提倡，但对n较大的文件仍是颇有效的。由于其运行时间主要耗费在建初始堆和调整建新堆时进行的反复“筛选”上。

　　堆排序在最坏的状况下，其时间复杂度也为O(nlogn)。相对于快速排序来讲，这是堆排序的最大优势。此外，堆排序仅需一个记录大小的供交换用的辅助存储空间。

参考资料：

　　严蔚敏《数据结构》

　　http://www.cnblogs.com/dolphin0520/archive/2011/10/06/2199741.html

　　http://blog.csdn.net/morewindows/article/details/6709644

出处：http://www.cnblogs.com/mengdd/archive/2012/11/30/2796845.html  （对原文代码有改动。本文代码为java代码。通过了博主测试，运行没有问题。）