怎样写出一个较好的高速排序程序

写出一个较好的高速排序程序

• 高速排序是经常使用的排序算法之中的一个，但要想写出一个又快又准的使用程序，就不是那么简单了

须要注意的事项

• 首先要写正确。通常使用递归实现。其递归相当于二叉树展开，因此假设要用迭代实现的话须要使用一个队列来保存兴许遍历信息。
• 高速排序须要找到一个pivot值，假设顺序选择pivot则易造成N^2的复杂度，假设使用随机数则效果最好，但开销又太大，採取三数中值法比較合适。三数中值法指的是选取第一个值，最后一个值，数组中间的值的中值。有文献表明能够提升5%的执行时间。
• 当数组长度较小时，如10个元素下面，最好使用插入排序或者选择排序完毕，以防止复杂度常数因子过大或多次函数调用带来的开销。而递归究竟层数组长度总是会变小的，因此这么做很有必要。
• 在合并前后两部分数组时，採用两边夹方法，在前后两部分各找到一个大于和小于的值再交换。相比通常情况下找到比pivot小的值就进行交换，能提高执行效率。

实现代码

• 代码例如以下。包含插入排序insert_sort，递归函数，三分中值函数三个辅助函数。
• 三分中值函数事实上採用的是插入排序。通过三次比較，确定中值。
• 插值算法使用暂时变量tmp避免了大量swap函数调用。

#include<iostream>
#include<iomanip>
#include<vector>
#include<cstdlib>
#include<ctime>
#include<algorithm>

using namespace std;

inline void swap(vector<int>& num, int p, int q){
    int t = num[p];
    num[p] = num[q];
    num[q] = t;
}

void insert_sort(vector<int>& num){
    int tmp, j;
    for (int i = 1; i < num.size(); i++){
        tmp = num[i];
        for (j = i - 1; j >= 0 && num[j] > tmp; j--)
            num[j + 1] = num[j];
        num[j + 1] =tmp;
    }
}

int quick_sort_sub(vector<int>& num, int p, int q){
    if (p >= q)
        return 0;

    // if 4 elements or less, use insert sort
    if (p + 10 > q){
        vector<int> tnum(num.begin() + p, num.begin() + q + 1);
        insert_sort(tnum);
        for (int i = 0; i < tnum.size(); i++)
            num[p + i] = tnum[i];
    }

    int idx = quick_three_partition(num, p, q);
    swap(num, idx, q);

    int pivot = num[q];
    int left = p, right = q - 1;
    while (1){
        while (num[left] < pivot)
            ++left;
        while (num[right] >= pivot)
            --right;
        if (left < right)
            swap(num, left, right);
        else
            break;
    }
    swap(num, left, q);
    quick_sort_sub(num, p, left - 1);
    quick_sort_sub(num, left + 1, q);

    return left;
}

void quick_sort(vector<int>& num){
    quick_sort_sub(num, 0, num.size() - 1);
}

int main(){
    const int n = 10;
    /*int num_array[n]= {2,1};
    vector<int> num(num_array, num_array + n);*/

    srand( time(NULL) );
    vector<int> num(n);
    for (auto& e : num)
        e = rand() % n;

    quick_sort(num);


    for (auto& e : num)
        cout << setw(4) << e << ‘ ‘;
    cout << endl;

    cout << "vector is sorted? : " << is_sorted(num.begin(), num.end()) << endl;

    return 0;

}