八大排序算法之二希尔排序（Shell Sort）

希尔排序是1959 年由D.L.Shell 提出来的，相对直接排序有较大的改进。希尔排序又叫缩小增量排序

基本思想：

先将整个待排序的记录序列分割成为若干子序列分别进行直接插入排序，待整个序列中的记录“基本有序”时，再对全体记录进行依次直接插入排序。

操作方法：

1. 选择一个增量序列t₁，t₂，…，t_k，其中t_i>t_j，其中i<j,t_k=1；
2. 按增量序列个数k，对序列进行k 趟排序；
3. 每趟排序，根据对应的增量t_i，将待排序列分割成若干长度为m 的子序列，分别对各子表进行直接插入排序。仅增量因子为1 时，整个序列作为一个表来处理，表长度即为整个序列的长度。

算法实现：

我们简单处理增量序列：增量序列d = {n/2 ,n/4, n/8 .....1} n为要排序数的个数

即：先将要排序的一组记录按某个增量d（n/2,n为要排序数的个数）分成若干组子序列，每组中记录的下标相差d.对每组中全部元素进行直接插入排序，然后再用一个较小的增量（d/2）对它进行分组，在每组中再进行直接插入排序。继续不断缩小增量直至为1，最后使用直接插入排序完成排序。

 同样使用杭电acm 1040 来测试： AC啦

#include<iostream>
using namespace std;
const int N =1005;
void sort(int a[],int num);
void print(int a[],int num);
void shellSort(int a,int num);
int main()
{
    int s[N];
    int T;
    cin>>T;
    while(T--)
    {
        int n;
        cin>>n;
        for(int j=0;j<n;j++)
           cin>>s[j];
        sort(s,n);
        print(s,n);       
        cout<<endl;        
    }
    return 0;
}
void sort(int a[],int num)//插入排序代码 
{
    for(int i=1;i<num;i++)
      if(a[i]<a[i-1])
      {
          int x=a[i];
          int j=i;
          while(j>0&&x<a[j-1])
          {
               a[j]=a[j-1];//后移 
             j--;    
        }
        a[j]=x;
      }
}
void shellSort(int a[],int n)//希尔排序的核心代码 
{
    int d=n/2;
    while(d>=1)
    {
        sort(a,d);//同样使用了插入排序 
        d=d/2;
    }
}
void print(int a[],int n)//输出数组元素 
{
    cout<<a[0];
    for(int k=1;k<n;k++)
    {
        cout<<" "<<a[k];
    }
}

希尔排序时效分析很难，关键码的比较次数与记录移动次数依赖于增量因子序列d的选取，特定情况下可以准确估算出关键码的比较次数和记录的移动次数。目前还没有人给出选取最好的增量因子序列的方法。增量因子序列可以有各种取法，有取奇数的，也有取质数的，但需要注意：增量因子中除1 外没有公因子，且最后一个增量因子必须为1。希尔排序方法是一个不稳定的排序方法。

void shellsort2(int a[], int n)
{
    int j, gap;
    
    for (gap = n / 2; gap > 0; gap /= 2)
        for (j = gap; j < n; j++)//从数组第gap个元素开始
            if (a[j] < a[j - gap])//每个元素与自己组内的数据进行直接插入排序
            {
                int temp = a[j];
                int k = j - gap;
                while (k >= 0 && a[k] > temp)
                {
                    a[k + gap] = a[k];
                    k -= gap;
                }
                a[k + gap] = temp;
            }
}