1、算法出现的背景

之前讲的,当我们排序的是一个近乎有序的序列时,快速排序会退化到一个O(n^2)级别的排序算法,

而对此的改进就是引入了随机化快速排序算法;但是当我们排序的是一个数值重复率非常高的序列时,

此时随机化快速排序算法就不再起作用了,而将会再次退化为一个O(n^2)级别的排序算法,那为什么

会出现这种情况呢?且听下面的分析:

 

如上图所示就是之前分析的快速排序算法的partition的操作原理,我们通过判断此时i索引指向的数组

元素e>v还是<v,将他放在橙色或者是紫色两个不同的位置,然后将整个数组分成两个部分递归下去;

但是这里其实我们是没有考虑=v的情况,其实隐含的意思就是下面的两种情况:

 

         

其实从这里就可以看出来了,不管是>=v还是<=v,当我们的序列中存在大量重复的元素时,

排序完成之后就会将整个数组序列分成两个极度不平衡的部分,所以又退化到了O(n^2)级别

的时间复杂度,这是因为对于每一个"基准"元素来说,重复的元素太多了,如果我们选的"基准"

元素稍微有一点的不平衡,那么就会导致两部分的差距非常大;即时我们的"基准"元素选在了

一个平衡的位置,但是由于等于"基准"元素的元素也非常多,也会使得序列被分成两个及其不平

衡的部分,那么在这种情况下快速排序就又会退化成O(n^2)级别的排序算法。如何解决呢?

这就要用到今天讲的双路快速排序算法的原理了。

 

 

 

2、双路快速排序算法的原理

之前说的快速排序算法是将>v和<v两个部分元素都放在索引值i所指向的位置的左边部分,而我们

的双路快速排序算法则不同,他使用两个索引值(i、j)用来遍历我们的序列,将<v的元素放在索

引i所指向位置的左边,而将>v的元素放在索引j所指向位置的右边,这也正是双路排序算法的

partition原理:

 

基本思想: 

 首先从左边的i索引往右边遍历,如果i指向的元素<v,

                                                                那直接将i++移动到下一个位置,直道i指向的元素>=v则停止

 

 然后使用j索引从右边开始往左边遍历,如果j指向的元素>v,

                                                                那直接将j--移动到下一个位置,直道j指向的元素<=v则停止

 

 此时i之前的元素都已经归并为<v的部分了,而j之后的元素也

                                                               都已经归并为>v的部分了,此时只需要将arr[i]和arr[j]交换位置即可

 

 这样就可以避免出现=v的元素全部集中在某一个部分,这正

                                                               是双路排序算法的一个核心

 

  将i++,j--开始遍历后后面的元素

 

3、双路快速排序算法的实现(基于C++)

 1 /******************************** 双路快速排序算法实现 ***************************************/
 2 template<typename T>
 3 T __partition2 (T arr[], int left, int right)
 4 {
 5     std::swap(arr[left], arr[std::rand() % (right - left + 1) + left]);   // 随机产生"基准"元素所在位置,并与第一个元素交换位置
 6     T v = arr[left];                       // 将第一个元素作为"基准"元素
 7 
 8     // 使用i索引从左到右遍历,使用j索引从右到左遍历
 9     int i = left + 1;       // 索引值i初始化为第二个元素位置
10     int j = right;          // 索引值j初始化为最后一个元素位置
11     while (true) {
12         while ((i < right) && (arr[i] < v)) i++;    // 使用索引i从左往右遍历直到 arr[i] < v
13         while ((j > left + 1) && (arr[j] > v)) j--; // 使用索引j从右往左遍历直到 arr[j] > v
14         if (i >= j) break;                          // 退出循环的条件
15         std::swap(arr[i], arr[j]);                  // 将 arr[i] 与 arr[j] 交换位置
16         i++;                                        // i++    j--
17         j--;
18     }
19 
20     std::swap(arr[left], arr[j]);             // 最后将"基准"元素v放置到合适的位置
21 
22     return j;
23 }
24 
25 template<typename T>
26 void __quickSort2 (T arr[], int left, int right)
27 {
28     if (right - left <= 40) {           // 当递归到序列数据量较小时使用插入排序算法
29         __insertSortMG<T>(arr, left, right);
30         return;
31     }
32 
33     int p = __partition2<T>(arr, left, right);  // 对arr[left...right]区间元素进行partition操作,找到"基准"元素
34     __quickSort2<T>(arr, left, p - 1);          // 对基准元素之前的序列递归调用__quickSort函数
35     __quickSort2<T>(arr, p + 1, right);         // 对基准元素之后的序列递归调用__quickSort函数
36 }
37 
38 template<typename T> 
39 void quickSort2 (T arr[], int count)
40 {
41     std::srand(std::time(NULL));         // 种下随机种子
42     __quickSort2<T>(arr, 0, count - 1);  // 调用__quickSort函数进行快速排序
43 }
44 /*********************************************************************************************/

 

4、性能测试

一般序列:

 

近乎有序的序列:

 

重复率非常高的序列: