Shell Sort(希尔排序)

近日学习了Shell Sort，也就是希尔排序，也称递减增量排序算法。在1959年由DL．Shell提出于1959年提出，由此得名。

此版本算法是在插入排序(Insertion Sort)基础上，将数组分成了h份(gap).也就是在数组中每隔h个数取出一个数，为一个子数组。先在子数组上进行排序，然后不断减小h的大小，直到h == 1 时，也就是完全变成插入排序的时候，排序完成。

算法复杂度取决于h(步长/步进)的选择，在最差的时候，也就是h == 1的时候，希尔排序就变成了插入排序，复杂度为O(n^₂)。而且算法的性能不仅取决于h，还取决于h之间的数学性质，比如它们之间的公因式(Algorithm 4th).

但是，很多论文研究了各种不同的递增序列，都无法证明某个序列的性质是最好的。

以下引自Wiki：

步长的选择是希尔排序的重要部分。只要最终步长为1任何步长串行都可以工作。算法最开始以一定的步长进行排序。然后会继续以一定步长进行排序，最终算法以步长为1进行排序。当步长为1时，算法变为插入排序，这就保证了数据一定会被排序。

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

步长串行	最坏情况下复杂度

已知的最好步长串行是由Sedgewick提出的 (1, 5, 19, 41, 109,...)，该串行的项来自 9 * 4^i - 9 * 2^i + 1 和 4^i - 3 * 2^i + 1 这两个算式[1].这项研究也表明“比较在希尔排序中是最主要的操作，而不是交换。”用这样步长串行的希尔排序比插入排序和堆排序都要快，甚至在小数组中比快速排序还快，但是在涉及大量数据时希尔排序还是比快速排序慢。

另一个在大数组中表现优异的步长串行是(斐波那契数列除去0和1将剩余的数以黄金分区比的两倍的幂进行运算得到的数列)：（1, 9, 34, 182, 836, 4025, 19001, 90358, 428481, 2034035, 9651787, 45806244, 217378076, 1031612713, …）^[2]

此版本算法与复杂递增序列性能相近。

Source Code in Java：

public static void ShellSort(Comparable[] a) {
    int N = a.length;
    int h = 1;
    while(h < N / 3)
        h = 3 * h + 1;
    while(h >= 1) {
        for(int i = h; i < N; i++) {
            for(int j = i; j >= h && less(a[j], a[j-h]); j -= h) {
                exch(a, j, j-h);
            }
        }
        h /= 3;
    }
}

一般说来，h或(length - h)就是所分的组数：

所以也可以这么写：

Source Code in C++:

void ShellSort(int a[], int n)
{

    // 此版本下 h 从 n / 2 起步
    for (int h = n / 2; h > 0; h /= 2)
    {
        // 可以从前面开始数，有 h 组
        for (int i = 0 ;i < h; i++)
        {
            for (int j = i + h; j < n; j += h) 
            {
                // 如果a[j] < a[j-h]，则寻找a[j]位置，并将后面数据的位置都后移。
                // 这样做比exch效率稍高 in Algorithm 4th
                if (a[j] < a[j - h])
                {
                    int temp = a[j];
                    int k = j - h;
                    while (k >= 0 && a[k] > temp)
                    {
                        a[k + h] = a[k];
                        k -= h;
                    }
                    a[k + h] = temp;
                }
            }
        }

    }
}

图示实例如下：

  

将各自子序列用插入排序法排序。

完整版本图示：

　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　

                                                                            　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 Vane_Tse On the Road.   2014-06-19  11:20:07