数据结构与算法——堆排序

优先队列可以用以O(NlogN)时间进行排序，基于该思想的算法叫做堆排序。

建立N个元素的二叉堆（比如是最小堆）需要花费O(N)时间，然后执行N次deleteMin操作，每次都是最小的元素离开堆。将每次离开的最小元素依次存放在另外一个数组，这样就得到了一个顺序序列。每次deleteMin操作花费O(logN)时间，总的花费时间为O(NlogN);因此堆排序的时间复杂度为O(N)+O(NlogN)，也就是O(NlogN)。

我们采取的方法是，将每次deleteMin操作删除的最小元素，堆缩小了1，我们将删除的元素放在堆的后面。

比如：堆中有6个元素，第一个deleteMin操作删除了a1，现在堆中剩了5个元素，因此把a1放在位置6上面。下一次deleteMin操作删除了a2,现在堆中剩了4个元素，因此把a1放在位置5上面。在最后一次deleteMin操作之后，该数组将以递减顺序包含这些元素。如果想将这些元素以递增的顺序排序，可以建立最大堆（父结点的值大于子结点的值）。

建立二叉堆的过程见博客：http://blog.csdn.net/linux_ever/article/details/50488816

源代码：

/*************************************************************************
	> File Name: heapsort.cpp
	> Author: 
	> Mail: 
	> Created Time: 2016年01月10日 星期日 20时34分07秒
 ************************************************************************/

#include <iostream>
#include <vector>
#include <algorithm>
#include <iterator>
using namespace std;


inline int leftChild(int i)
{
    return 2*i+1;//因为该堆是从数组下标0开始存放值的，和上次建立二叉堆的时候不一样
}

template<class T>
void percDown(vector<T> & a, int i, int n)
{
    int child;
    T tmp;

    for (tmp = a[i]; leftChild(i) < n; i = child){
        child = leftChild(i);
        if(child != n-1 && a[child] < a[child + 1])
            child++;

        if(tmp < a[child])
            a[i] = a[child];
        else
            break;
    }
    a[i] = tmp;
}

template<class T>
void heapsort(vector<T> & v)
{
    for (int i = v.size()/2; i >= 0; --i){//建立二叉堆
        percDown(v, i, v.size());
    }
    
    for (int j = v.size()-1; j > 0; j--){
        swap(v[0], v[j]);//删除最大值
        percDown(v, 0, j);
    }
}


int main()
{
    vector<int> v;
    v.push_back(97);
    v.push_back(53);
    v.push_back(28);
    v.push_back(33);
    v.push_back(23);
    v.push_back(77);
    v.push_back(89);
    v.push_back(-8);
    v.push_back(10);

    cout << "v: ";
    copy(v.begin(), v.end(), ostream_iterator<int> (cout, " "));
    cout << endl;

    heapsort(v);
    cout << "sorted v: ";
    copy(v.begin(), v.end(), ostream_iterator<int> (cout, " "));
    cout << endl;


    return 0;
}

下面是一个下滤的过程图