9.2希尔排序

希尔排序

void Shell_sort( ElementType A[], int N )
{
    for( D=N/2; D>0; D/=2 ){// 希尔增量序列
        for( P=D; P<N; P++ ) {//插入序列
            Tmp = A[P];
            for( i=P; i>=D && A[i-D]>Tmp; i-=D )
                A[i] = A[i-D];
            A[i] = Tmp;
        }
    }
}

• 最坏情况：$T=\Theta (N^2)$

希尔增量序列

• 原始希尔排列 $D_M=[N/2],D_k=[D_{k+1}/2]$
• Hibbard增量序列
  • $D_k=2^k-1$ 相邻元素互质
  • 最坏情况：$T=\Theta (N^{3/2})$
  • 猜想：$T_{avg}=O(N^{5/4})$
• Sedgewick增量序列
  • {1,5,19,41,109,…}
  • 猜想：$T_{avg}=O(N^{7/6}),T_{worst}=O(N^{4/3})$

void ShellSort( ElementType A[], int N )
{ /* 希尔排序 - 用Sedgewick增量序列 */
     int Si, D, P, i;
     ElementType Tmp;
     /* 这里只列出一小部分增量 */
     int Sedgewick[] = {929, 505, 209, 109, 41, 19, 5, 1, 0};
      
     for ( Si=0; Sedgewick[Si]>=N; Si++ ) 
         ; /* 初始的增量Sedgewick[Si]不能超过待排序列长度 */
 
     for ( D=Sedgewick[Si]; D>0; D=Sedgewick[++Si] )
         for ( P=D; P<N; P++ ) { /* 插入排序*/
             Tmp = A[P];
             for ( i=P; i>=D && A[i-D]>Tmp; i-=D )
                 A[i] = A[i-D];
             A[i] = Tmp;
         }
}