快速排序

韩冰

<>分而治之方法还可以用于实现另一种完全不同的排序方法，这种排序法称为快速排序（quick sort）。在这种方法中， n 个元素被分成三段（组）：左段l e f t，右段r i g h t和中段m i d d l e。中段仅包含一个元素。左段中各元素都小于等于中段元素，右段中各元素都大于等于中段元素。因此l e f t和r i g h t中的元素可以独立排序，并且不必对l e f t和r i g h t的排序结果进行合并。m i d d l e中的元素被称为支点( p i v o t )。图1 4 - 9中给出了快速排序的伪代码。

: ?9 G" }8 P' ^' C' g0 }
/ /使用快速排序方法对a[ 0 :n- 1 ]排序

& U8 `$ Q2 Y# U8 h. ^从a[ 0 :n- 1 ]中选择一个元素作为m i d d l e，该元素为支点

/ |0 W% }# \- M/ ]0 W2 @把余下的元素分割为两段left 和r i g h t，使得l e f t中的元素都小于等于支点，而right 中的元素都大于等于支点
0 c' k3 J- s) ?6 e' g- V! o
递归地使用快速排序方法对left 进行排序

递归地使用快速排序方法对right 进行排序
9 z5 a. u9 g/ v& p% w
+ M8 ?) M" b/ j4 [- @所得结果为l e f t + m i d d l e + r i g h t
3 K9 q( g5 N) T& ?: g
  C2 _4 X+ W, x. v图14-9 快速排序的伪代码
- L5 ]! e: l# [! k) R

考察元素序列[ 4 , 8 , 3 , 7 , 1 , 5 , 6 , 2 ]。假设选择元素6作为支点，则6位于m i d d l e；4，3，1，5，2位于l e f t；8，7位于r i g h t。当left 排好序后，所得结果为1，2，3，4，5；当r i g h t排好序后，所得结果为7，8。把right 中的元素放在支点元素之后， l e f t中的元素放在支点元素之前，即可得到最终的结果[ 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 ]。
0 j7 K1 D& I$ c  m# d1 a
! g' Q, T4 a1 H  D( z把元素序列划分为l e f t、m i d d l e和r i g h t可以就地进行（见程序1 4 - 6）。在程序1 4 - 6中，支点总是取位置1中的元素。也可以采用其他选择方式来提高排序性能，本章稍后部分将给出这样一种选择。

程序14-6 快速排序

template<CLASS T>
% s/ }, }0 T- G
, E" K. y/ O) T$ R- M3 Z! W" }void QuickSort(T*a, int n)

3 g+ G* r" p7 k! \4 Z{// 对a[0:n-1] 进行快速排序

{// 要求a[n] 必需有最大关键值

quickSort(a, 0, n-1);
8 ^. ]& n2 Z7 A6 C7 w  s) r
template<CLASS T>
9 b" S+ C7 @5 X# a
void quickSort(T a[], int l, int r)
, Y, b% K7 T$ h/ @$ e
' j7 d- N  W' w, s# X) j7 j{// 排序a [ l : r ]， a[r+1] 有大值
9 {9 g/ F) E4 B
1 h8 ^. D2 N0 W9 W6 a9 a0 Sif (l &gt;= r) return;
9 E- U8 L: _3 Q4 r4 |
% {. M' d& l' m. z: d5 j/ vint i = l, // 从左至右的游标

: y' N+ t; O% y! A- Gj = r + 1; // 从右到左的游标
0 z1 h- M& d8 d! F: C$ {
( Y0 Z4 C2 M$ z+ z. O+ \  A) WT pivot = a[l];

// 把左侧&gt;= pivot的元素与右侧&lt;= pivot 的元素进行交换
+ D3 |2 B3 Q: q& {0 V1 E
1 u0 y; d1 d+ P4 vwhile (true) {

4 d6 m7 Q0 y& _5 [2 V+ {do {// 在左侧寻找&gt;= pivot 的元素

. D) m/ V( {9 }2 M+ M9 ]i = i + 1;

} while (a &lt; pivot);

. ^! O. r7 M( \9 ^' n* wdo {// 在右侧寻找&lt;= pivot 的元素
$ }" E/ ~7 B" z' @7 J! `$ `' O
+ m9 C- R. q4 V9 q$ K7 Tj = j - 1;

} while (a[j] &gt; pivot);
( K0 V, g0 f, l
! ?  Y" k/ X* o+ E3 T4 u/ f! J3 uif (i &gt;= j) break; // 未发现交换对象
' J  z) ?5 _, Y- s% s
Swap(a, a[j]);

% ]. K9 b, T; X1 N/ Y/ ?- d& G, q; O}

1 a" J( q* g& W& N, q// 设置p i v o t

a[l] = a[j];
0 F5 s! i( E# B: Y! P
a[j] = pivot;
/ y; s3 G" U+ F; n/ K
+ j8 J. V7 V8 d9 I6 L( N8 oquickSort(a, l, j-1); // 对左段排序

quickSort(a, j+1, r); // 对右段排序
5 G( N# J9 s9 k/ R& M" T6 o
}
' n! F9 _( s5 u4 T: {# p
, `& M8 i  x/ r) e若把程序1 4 - 6中d o - w h i l e条件内的&lt;号和&gt;号分别修改为&lt; =和&gt; =，程序1 4 - 6仍然正确。实验结果表明使用程序1 4 - 6的快速排序代码可以得到比较好的平均性能。为了消除程序中的递归，必须引入堆栈。不过，消除最后一个递归调用不须使用堆栈。消除递归调用的工作留作练习（练习1 3）。程序1 4 - 6所需要的递归栈空间为O (n)。若使用堆栈来模拟递归，则可以把这个空间减少为O ( l o gn)。在模拟过程中，首先对left 和right 中较小者进行排序，把较大者的边界放入堆栈中。在最坏情况下l e f t总是为空，快速排序所需的计算时间为(n2 )。在最好情况下， l e f t和r i g h t中的元素数目大致相同，快速排序的复杂性为(nl o gn)。令人吃惊的是，快速排序的平均复杂性也是(nl o gn)。

定理2-1 快速排序的平均复杂性为(nl o gn)。
$ p3 d7 E, }4 I
; I+ ~% G5 S# `! G证明用t (n) 代表对含有n 个元素的数组进行排序的平均时间。当n≤1时，t (n)≤d，d为某一常数。当n &lt;1时，用s 表示左段所含元素的个数。由于在中段中有一个支点元素，因此右段中元素的个数为n-s- 1。所以左段和右段的平均排序时间分别为t (s), t (n-s- 1 )。分割数组中元素所需要的时间用cn 表示，其中c 是一个常数。因为s 有同等机会取0 ~n- 1中的任何一个值.

0 _. O. {, x' X* c, `9 T6 Z; a6 I如对（2 - 8）式中的n 使用归纳法，可得到t (n)≤kn l o ge n，其中n&gt; 1且k=2(c+d)，e～2 . 7 1 8为自然对数的基底。在归纳开始时首先验证n= 2时公式的正确性。根据公式（ 1 4 - 8），可以得到t( 2 )≤2c+ 2d≤k nl o ge 2。在归纳假设部分，假定t(n)≤kn l o ge n（当2≤n＜m 时，m 是任意一个比2大的整数＝.

) c, f* U& r  y6 F图1 4 - 1 0对本书中所讨论的算法在平均条件下和最坏条件下的复杂性进行了比较。


方法最坏复杂性平均复杂性
. d) V4 J+ T4 v- b* {% ]' L
冒泡排序n2 n2
6 a  Q3 e, {8 D  t8 B
计数排序n2 n2

插入排序n2 n2

+ C& s0 @7 C" w* r选择排序n2 n2
4 w0 J/ H4 b5 G2 M8 m/ @
1 ]' ~' t, F/ H/ F0 `堆排序nl o gn nl o gn

6 U4 _8 ^, k  K" p归并排序nl o gn nl o gn
2 `& U! S  Y: s2 m/ J7 I
快速排序n2 nl o gn
! x& N8 c. _4 I; D$ _
图14-10 各种排序算法的比较
3 S) o3 }' l, P; R! y% M+ S9 e

中值快速排序（ median-of-three quick sort）是程序1 4 - 6的一种变化，这种算法有更好的平均性能。注意到在程序1 4 - 6中总是选择a [ 1 ]做为支点，而在这种快速排序算法中，可以不必使用a [ 1 ]做为支点，而是取{a[1],a[(1+r)/2],a[r]} 中大小居中的那个元素作为支点。例如，假如有三个元素，大小分别为5，9，7，那么取7为支点。为了实现中值快速排序算法，一种最简单的方式就是首先选出中值元素并与a[1] 进行交换，然后利用程序1 4 - 6完成排序。如果a [ r ]是被选出的中值元素，那么将a[1] 与a[r] 进行交换，然后将a [ 1 ]（即原来的a [ r ]）赋值给程序1 4 - 6中的变量p i v o t，之后继续执行程序1 4 - 6中的其余代码。

图2 - 11中分别给出了根据实验所得到的归并排序、堆排序、插入排序、快速排序的平均时间。对于每一个不同的n, 都随机产生了至少1 0 0组整数。随机整数的产生是通过反复调用s t d l i b . h库中的r a n d o m函数来实现的。如果对一组整数进行排序的时间少于1 0个时钟滴答，则继续对其他组整数进行排序，直到所用的时间不低于1 0个时钟滴答。在图2 - 11中的数据包含产生随机整数的时间。对于每一个n，在各种排序法中用于产生随机整数及其他开销的时间是相同的。因此，图2 - 11中的数据对于比较各种排序算法是很有用的。

对于足够大的n，快速排序算法要比其他算法效率更高。从图中可以看到快速排序曲线与插入排序曲线的交点横坐标比2 0略小，可通过实验来确定这个交点横坐标的精确值。可以分别用n = 1 5 , 1 6 , 1 7 , 1 8 , 1 9进行实验，以寻找精确的交点。令精确的交点横坐标为nBr e a k。当n≤nBreak 时，插入排序的平均性能最佳。当n＞nBreak 时，快速排序性能最佳。当n＞nBreak 时，把插入排序与快速排序组合为一个排序函数，可以提高快速排序的性能，实现方法是把程序1 4 - 6中的以下语句：

7 B, p5 a  ?; o, yif(l &gt;= r)r e t u r n ;

替换为

if (r-1<NBREAK) {InsertionSort(a,l,r); return;}
0 e' p6 H  H5 u
这里I n s e r t i o n S o r t ( a , l , r )用来对a [ 1 : r ]进行插入排序。测量修改后的快速排序算法的性能留作练习（练习2 0）。用更小的值替换n B r e a k有可能使性能进一步提高（见练习2 0）。

. n7 N2 [% M+ S1 ?大多数实验表明，当n＞c时（c为某一常数），在最坏情况下归并排序的性能也是最佳的。而当n≤c时，在最坏情况下插入排序的性能最佳。通过将插入排序与归并排序混合使用，可以提高归并排序的性能（练习2 1）。</P>