排序 —— 希尔排序（Shell sort）

希尔排序（Shell sort）的名称源于它的发明者 Donald Shell，该算法是冲破二次时间屏障（冒泡和插入排序，基于相邻元素的交换）的第一批算法。希尔排序改进了冒泡和插入排序的相邻元素才进行交换，而是比较相距一段距离的元素来工作，各趟比较所用的距离随着算法的进行而减少，直到只比较相邻元素的最后一趟排序为止。正是因为这样的工作机制，希尔排序有时也称为缩小增量排序（diminishing increment sort）。

• C 语言版

  typedef int ElementType;
  
  void ShellSort (ElementType A[], int N) {
  	
  	int Increment, i, j;
  	ElementType Tmp;
  	
  	for (Increment = N/2; Increment >= 1; Increment /= 2) {
  		
  		for (i = Increment; i < N; ++i) {
  			Tmp = A[i];
  			for (j = i; j >= Increment; j -= Increment) {
  				if (Tmp < A[j-Increment])
  					A[j] = A[j-Increment];
  				else break;
  						// 尤其注意这里的 break；
  			}
  			A[j] = Tmp;
  		}
  	}
  }
  

• python

  def shell_sort(seq):
      n = len(seq)
      gap = n // 2
      while gap >= 1:
          for i in range(gap, n):
              temp = seq[i]
              j = i
              while j >= gap and seq[j-gap] > temp:
                  seq[j] = seq[j-gap]
                  # print(seq)
                  j -= gap
              seq[j] = temp
          gap //= 2
  

• 实现思路：

  • 三层循环：
    • 最外层是遍历逐渐缩小的增量，直至为 1；
    • 内层遍历待比较的元素，起点；
    • 最内层在增量上进行移动，先对起始元素进行保存；
      • 最内层循环是十分严苛的，如果向左逐gap的遍历时，j-gap <= temp 时立刻退出循环的；
      • 只有 j-gap > temp 才会向左继续走下去，并进行交换，交换是把左边的较大的值覆盖右端的值，所以在最内层循环之前要对最右侧的值，也就是最内层遍历的起点进行保存；

1. 一种先验的增量

已知的最好步长序列由Marcin Ciura设计（1，4，10，23，57，132，301，701，1750，…）这项研究也表明“比较在希尔排序中是最主要的操作，而不是交换。”用这样步长序列的希尔排序比插入排序和堆排序都要快，甚至在小数组中比快速排序还快，但是在涉及大量数据时希尔排序还是比快速排序慢。

def shell_sort_fixed_gaps(seq):
    gaps = [701, 301, 132, 57, 23, 10, 4, 1]
    n = len(seq)
    for gap in gaps:
        for i in range(gap, n):
            temp = seq[i]
            j = i
            while j >= i and seq[j-gap] > temp:
                seq[j] = seq[j-gap]
                j -= gap
            seq[j] = temp