希尔排序

希尔排序

希尔排序，也称为递减增量排序算法，是插入排序的一种更高效的改进版本。希尔排序是非稳定排序算法。

希尔排序是基于插入排序的以下两点性质而提出改进方法的：

• 插入排序对几乎已经排好序的数据操作时，效率高，即可以达到线性排序的效率
• 但插入排序一般来说是低效的，因为插入排序每次只能将数据移动一位(in-place)

算法实现

原始的算法实现在最坏的情况下需要进行O(n^2)的比较和交换。V. Pratt的书对算法进行了少量修改，可以使得性能提升至O(n log2 n)。这比最好的比较算法的O(n log n)要差一些。

希尔排序通过将比较的全部元素分为几个区域来提升插入排序的性能。这样可以让一个元素可以一次性地朝最终位置前进一大步。然后算法再取越来越小的步长进行排序，算法的最后一步就是普通的插入排序，但是到了这步，需排序的数据几乎是已排好的了（此时插入排序较快）。

假设有一个很小的数据在一个已按升序排好序的数组的末端。如果用复杂度为O(n^2)的排序（冒泡排序或插入排序），可能会进行n次的比较和交换才能将该数据移至正确位置。而希尔排序会用较大的步长移动数据，所以小数据只需进行少数比较和交换即可到正确位置。

一个更好理解的希尔排序实现：将数组列在一个表中并对列排序（用插入排序）。重复这过程，不过每次用更长的列来进行。最后整个表就只有一列了。将数组转换至表是为了更好地理解这算法，算法本身仅仅对原数组进行排序（通过增加索引的步长，例如是用i += step_size而不是i++）。

例如，假设有这样的一组数 [ 13 14 94 33 82 25 59 94 65 23 45 27 73 25 39 10 ]，如果我们以步长为5开始进行排序，我们可以通过将这列表放在有5列的表中来更好地描述算法，这样他们就应该看起来是这样： 

　　　　　　　　13  14  94  33  82

　　　　　　　　25  59  94  65  23

　　　　　　　　45  27  73  25  39

　　　　　　　　10

然后，我们对列进行排序

　　　　　　　　10  14  73  25  23

　　　　　　　　13  27  94  33  39

　　　　　　　　25  59  94  65  82

　　　　　　　　45

将上述的四行数，依序连接在一起我们得到[10 14 73 25 23 13 27 94 33 39 25 59 94 65 82 45 ]。这时10已经移至正确位置了，然后再以3为步长进行排序： 

　　　　　　　　10  14  73

　　　　　　　　25  23  13

　　　　　　　　27  94  33

　　　　　　　　39  25  59

　　　　　　　　94  65  82

　　　　　　　　45

排序之后变为：

　　　　　　　　10  14  13

　　　　　　　　25  23  33

　　　　　　　　27  25  59

　　　　　　　　39  65  73

　　　　　　　　45  94  82

　　　　　　　　94

最后以1步长进行排序（此时就是简单的插入排序了）。

步长序列

步长的选择是希尔排序的重要部分。只要最终步长为1任何步长序列都可以工作。

Donald Shell最初建议的步长选择为N/2并且对步长取半直到步长为1，虽然这样取可以比O(n^2)类的算法（插入排序）更好，但这样仍然有减少平均时间和最差时间的余地。可能希尔排序最重要的地方在于当用较小步长排序后，以前用的较大步长仍然是有序的。比如，如果一个数列以步长5进行了排序然后再以步长3进行排序，那么该数列不仅是以步长3有序，而且是以步长5有序。如果不是这样，那么算法在迭代过程中会打乱以前的顺序，那就不会以如此短的时间完成排序了。 

步长序列	最坏情况下的复杂度
n/2^i	O(n^2)
2^k - 1	O(n^3/2)
2^i*3^j	O(nlog2n)

 

 

 

 

已知的最好步长序列是由Sedgewick提出的(1, 5, 19, 41, 109,...)，该序列的项来自和这两个算式。

C语言示例

 1 void shell_sort(int arr[], int len) {

 2     int gap, i, j;

 3     int temp;

 4     for (gap = len >> 1; gap > 0; gap >>= 1) {

 5         for (i = gap; i < len; i++) {

 6             temp = arr[i];

 7             for (j = i - gap; j >= 0 && arr[j] > temp; j -= gap) {

 8                 arr[j + gap] = arr[j];

 9             }

10             arr[j + gap] = temp;

11         }

12     }

13 }