排序算法之希尔排序

简介

之前介绍的几种排序算法，简单选择排序、直接插入排序，它们的时间复杂度都是$O(n^2)$。实际上在早期很长一段时间内，普遍认为排序算法的效率不会高于$O(n^2)$，直到希尔排序算法的出现才打破了这一认知。希尔排序是第一批冲破二次时间界($O(n^2)$)的算法之一。它通过比较相距一定距离的元素来排序，各趟比较所使用的距离随算法的进行而减小，直到比较相邻元素为止。因此，希尔排序也称缩减增量排序。

希尔排序通过使用一个增量序列$h_1$，$h_2$，…，$h_t$，将待排序元素序列分为若干个子序列。然后对这些子序列进行排序，使其$a[i]<=a[i+h_k]$，此时所以相隔$h_k$的元素都被排序，称为序列是$h_k$排序的。如下图。

下面我们继续通过代码进行逐步讲解，完整代码见https://github.com/kfcyh/sort/tree/master/shellsort。

代码讲解

void ShellSort(SqList L)
{
	int i, j;
	int increment = L->length;
	do
	{
		increment = increment / 3 + 1;	//增量序列
		for (i = increment + 1; i <= L->length; i++)
		{
			if (L->r[i] < L->r[i - increment])	//需将r[i]插入有序增量子表
			{
				L->r[0] = L->r[i];	//暂存元素
				for (j = i - increment; j > 0 && L->r[0] < L->r[j]; j -= increment)
					L->r[j + increment] = L->r[j];	//元素后移
				L->r[j + increment] = L->r[0];	//插入元素至正确位置
			}
		}
	} while (increment > 1);
}

由代码可知，希尔排序对增量子序列的排序采用插入排序的思想，下面我们通过对序列{9,1,5,8,3,7,4,6,2}进行排序来进行讲解。

1. 首先是一个$while$循环，中止条件是$increment=1$。然后初始化增量$increment=length/3+1$。内部的外循环$for$每次对增量为$increment$的子序列进行排序。$i$从$increment+1$开始，中止条件为$i>length$。第一次$while$循环，$increment=4$。

2. 第一次$for$循环，$i=5$，$r[i]=3<r[i-increment]=9$，交换元素。

3. 循环继续，$i=6$，$r[i]=7>r[i-increment]=1$，不交换元素。

4. 循环继续，$i=7$，$r[i]=4<r[i-increment]=5$，交换元素。

5. 同上继续循环至，$i=9$，$r[i]=2<r[i-increment]=9$，交换元素。注意此时还要继续比较$r[5]$与$r[1]$，$r[5]<r[1]$交换。

6. 第一轮$while$后，序列如下所示。

7. 由于此时$increment=4>1$，继续$do$循环。$increment=increment/3+1=2$，同上过程，$i$从3到9。当$i=3、4$时，不用交换，当$i=5$时，交换。

8. 当$i=6、7、8、9$时，不用交换。

9. 由于此时$increment=2>1$，继续$do$循环。$increment=increment/3+1=1$，此时为最后一轮$while$循环，同上过程。

10. 完成排序，结果如下图。

性能分析

在上述排序过程中，增量"$increment$"的选取十分关键，它选取什么最好目前还是一个数学难题，不过大量的研究表明，当增量序列为$dlta[k]=2^{t-k+1}$时可以获得不错效果。总结其时间复杂度为$O(n^{3/2})$，好于直接排序的$O(n^2)$，需要注意的是，增量序列的最后一个增量必须是1。而且希尔排序不是一个稳定的排序算法。