目录
1. 排序介绍
在以前的程序里,排序程序占用的时间大约是总程序时间的30%。可以看到排序是基础而又被广泛应用的算法。现在排序算法的时间比重比30%小,这并不意味着排序算法的重要性在降低,相反这是更有效的排序算法的功劳。本文带我们从基础到复杂,领略6种排序算法的思想,实现,时间空间复杂度和稳定性。
基于比较的排序(Comparison Sort)是算法中不可绕过的重要内容,核心地位不必赘述,如果从信息论的角度看,Comparison Sort就是找到一种恰当的策略,使得每次做完比较之后,能够获取的信息量最大。下面举一个栗子。
猜数字游戏,请用最少的次数猜出一个在1-100中的数字。
为了让获取的信息量最大,容易知道必须让一次比较确定下来的东西尽可能多,要达到合格目的最好的办法就是让比较后可能产生的两种结果的概率均等。在上面的栗子中,如果问这个数字是否大于50?这样的话一旦结果确定,就砍去了一半的可能性,获取的信息量最大。
类似的道理,对于12个球的称重问题:
有12个外观一模一样的球,其中有一个的质量和其余12个不一样,记为X,请用一个天平以最少的比较次数找出球X。
之所以不用6-6的比较方式,是因为这么做产生的结构是确定的可预测的,肯定是不平衡;相比之下如果用4-4的比较方法,产生的三种结果左倾,右倾,平衡的概率是完全相等的,这样就充分利用了这次试探的机会,获取了最大的信息量。 进一步的阅读请参照这里和这里。
接下来将要介绍若干种基本的排序算法,如下表所示。
1.1 Total Order & Partial Order
在介绍具体的算法之前,首先需要介绍一些基本的概念和Java的实现。 既然是排序,首先要明确顺序的确切含义。我们常说的顺序是一种二元关系,其中全序关系(Total Order)和偏序关系(Partial Order)是较为常见的两种。先盗一张wiki上的图。
如上所示就是一个按包含排序的偏序关系,我们发现虽然集合x,y和集合y,z都包含集合y,但是他们之间并不存在包含关系,也就是说,偏序其实并不要求所有的元素都可互相比较,其强调的只是一种继承关系,一般来说可以与有向无环图(DAG)一一对应起来。数学上形式化定义的非严格偏序如下所示:
- a ≤ a (自反性 reflexivity)
- if a ≤ b and b ≤ a then a = b (反对称性 antisymmetry)
- if a ≤ b and b ≤ c then a ≤ c (传递性 transitivity)
传递性这个比较有意思,我们小时候常玩的石头剪刀布,仔细想想,其实就是intransitivity的,对应于上图就是出现了环,因此不是DAG图。
而全序关系与偏序关系比较起来,少了自反省,多了完全性
- either a ≤ b or b ≤ a or both (完全性 totality)
这么一对比就看出了全序和偏序的区别了,值得注意的是,全序关系也是一种偏序关系。对于排列问题,由于要确定每个元素的大小关系,因此要求设计的Comparator满足全序关系。另外在本文介绍的排序如果不作特别说明都是按照从小到大的顺序。
对于一个排序算法而言,必须要有几个指标来考量算法的优劣:
- inplace? 直观上说,就是算法是否需要开辟额外的空间,这里的额外有严格的数学定义,即more than clogN space。由此可见,即使需要开辟额外的空间为linear(如mergesort),算法也不能称为inplace。
- stable? 对于相同的元素,排序后其相对应位置是否和排序前一致,如果一致则陈伟stable。
- worst, average, best cases。顾名思义,就是分别就三种情况讨论算法的时间复杂度。
1.2 Java Implementation
在Java中涉及元素的比较一般使用Comparable和Comparator接口,前者的定义非常简单,只需要实现compareTo函数即可,如下所示:
1
2
3
public interface Comparable<Key>{
public int compareTo(Key that);
}
在使用过程中直接调用即可
1
item1.compareTo(Item2);
与Comparable只传入一个参数不同,Comparator更像一个操作符,要求传入两个对象,然后比较这两个对象的大小,如下所示:
1
2
3
4
public interface Comparator<Key>{
public int compare(Key o1, Key o2);
public boolean equals(Object obj); //可以不实现
}
下面总结一下Comparable和Comparator的区别和联系:
- Comparable是排序接口,如果该类实现了Comparable接口,就意味着“该类支持排序”,相当于内部比较器;
- Comparator是比较器,若我们需要比较某个类的次序,可以建立一个“该类的比较器”进行排序,相当于外部比较器。本质是解耦了数据类型的定义和如何比较数据的定义。
如果一个类实现了上述的两个接口之一,就可以方便的利用Java框架的java.util.Arrays进行数组排序了。
- Array.sort(Comparable[] a),这个时候的a元素,需要实现Comparable接口;
- Array.sort(Object[] a,Comparator comparator),这个时候需要实现一个Comparator来定义如何比较a中的元素。
有了这两种比较方法,编程就变得非常方便,因为对于某一类型的数组,可以通过传入不同的Comparator来得到不同定义下的相对大小关系。值得注意的是,在调用上述两个函数过程中,对于Object对象,Java会用Merge Sort,对于Primitive原始类型,Java会使用Quick Sort。算法细节及原因将会在后文介绍。
2. 排序算法
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
public class Example{
public static void sort(Comparable[] a){
/* see 2.1, 2.2, 2.3, 2.4, 2.5, 2.6*/
}
public static boolean isSorted(Comparable[] a){
// test weather the array entry is in order
for (int i = 1; i < a.length; i++){
if(less(a[i],a[i-1])) return false;
}
return true;
}
private static boolean less(Comparable v, Comparable w){
return v.compareTo(w) < 0;
}
private static void exchange(Comparable[] a, int i, int j){
Comparable temp = a[i];
a[i] = a[j];
a[j] = temp;
}
private static void show(Comparable[] a){
// print the array, on a single line.
for(int i = 0; i < a.length; i++){
StdOut.print(a[i] + " ");
}
StdOut.println();
}
public static void main(String[] args) {
// Read strings from standard input, sort them, and print
String[] a = StdIn.readAllStrings();
sort(a);
assert isSorted(a);
show(a);
}
}
2.1 基本排序
2.1.1 选择排序
选择排序(Selection Sort): 每次迭代过程中选出最小的元素,并且与当前迭代的位置进行交换,循环往复直到数组末尾。
1
2
3
4
5
6
7
8
9
10
11
12
public class SelectionSort {
public static void sort(Comparable[] a){
int N = a.length;
for(int i = 0; i < N; i++){
int min = i;
for(int j = i+1; j < N; j++){
if(less(a[j],a[min])) min = j;
}
exchange(a,i,min);
}
}
}
就像每个人都有较为稳定的性格一样,每一种排序算法之所以不一样,是因为它们保持的循环不变式(Invariant)是不一样的。对于插入排序而言,假设当前迭代位置是↑,那么不难发现以下两个不变性质:
- ↑左侧的元素都是有序的;
- ↑右侧的元素皆不小于↑左侧的元素。
算法复杂度不难算出是O(N2),确切来说是:
- ~ N2/2 compares;
- ~ N exchanges。
虽然选择排序的比较次数(compares)不尽如意,即便是已经排好顺序的数据也一视同仁地比较N2/2次,但是交换次数(exchanges)是所有排序算法中最少的,没有比一步到位更快的了。另外,由于在交换过程中靠前的元素会被换到后面,因此选择排序是unstable的;由于仅需要常数的额外空间做交换中介,因此是原地排序(inplace)。
2.1.2 插入排序
插入排序(Insertion Sort): 类比理牌,每次我们抓了一张新牌,都会寻找合适的位置把牌插入,使得新得到的牌组依旧满足递增的性质,直到所有牌都插入完毕,这个过程就是插入排序.
1
2
3
4
5
6
7
8
9
10
public class InsertionSort {
public static void sort(Comparable[] a){
int N = a.length;
for(int i = 0; i < N; i++){
for(int j = i; j > 0 && less(a[j],a[j-1]); j--){
exchange(a,j,j-1);
}
}
}
}
插入排序的sort
average case: ~N2/4 compares and ~N2/4 exchanges; worst case: ~N2/2 compares and ~N2/2 exchanges;best case: N-1 compares and 0 exchanges。
由于排序的过程不难发现,插入排序满足以下两个Invariants:
- ↑左侧的都是有序的;
- ↑右侧都是位置的;
出入排序的复杂度分析比选择排序稍微复杂一点。
- Best Case,如果是Sorted,那么只需要经过N-1此的compares和0次的exchanges;
- Worst Case,如果正好是逆序的,那么易知需要经过N2/2次的compares和exchanges;
- Average Case,经过简单计算可知需要N2/4次的compares和exchanges。
插入排序看起来并没有什么特别之处,尤其在规模比较大的时候,quadratic的复杂度无法忍受。但是在工业界中,插入排序作为快速排序(Quick sort)的补充在规模较小时(一般小于8)得到较为广泛的应用,例如在STL的sort和stdlib的qsort。
此外,一个有意思的结论是,倘若数组已经处于部分有序(partially sorted)的话,那么插入排序能在线性时间内完成。
2.1.3 希尔排序
希尔排序(Shellsort): 在插入排序中,元素比较的步长h = 1,也就是和相邻的元素进行比较和交换。那么,是否可以从比较大的步长开始,每一次迭代都缩小步长,直到最后一次迭代步长h = 1,退化为插入排序。
对于每次步长为h的迭代,我们成为h-sorting,h逐渐减少,最后变成1-sorting。 希尔排序的一个核心问题是,如何选择h变化的序列使得算法复杂度最低?
高爷爷大笔一挥,给出了递推式h = 3h + 1,并且经过复杂证明可以知道在这个递推式下希尔排序的最坏复杂度是O(N1.5)。Java 代码如下。
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
public class ShellSort {
public static void sort(Comparable[] a){
int N = a.length;
int h = 1;
while (h < N/3) h = 3*h + 1; //1,4,13,39;
while(h >= 1){
for(int i = h; i < N; i++){
for(int j = i; j >= h && less(a[j],a[j-h]); j -= h){
exchange(a,j,j-h);
}
}
h = h/3;
}
}
}
有人可能会问,在h-sorting中为什么要用插入排序呢,用选择排序不行吗?下面给出不严谨的感性证明:
- 在h较大时,由于比较稀疏,因此h-sorting的子数组比较小,因此排序较快;
- 在h较小时,虽然子数组的规模增大了,但是由于之前的排序使得数组已经部分有序,这个时候插入排序能在最坏的情况下保证线性时间内解决问题,故依然比较快。
希尔排序另外一个有趣的性质是,如果序列先经过了h-sorting然后经过了g-sorting,那么序列依旧会保持h-sorting。推而广之,那么序列一定也已经是(h+g)-sorting了。
3.比较
in place: 没有用到额外的存储,如aux[N] array; stable: sort age, then sort name, the result is first sorted by age, within same age, sorted by name.
4. Sorting programming example: Collinear Points
Collinear Points: 在给定的点中,求出所有过至少四点的直线。
具体code见这里。