希尔排序

对《大话数据结构》P389~P396—希尔排序，进行了自己的理解并完善了代码。

代码和解释如下（VS2012测试通过）：

#include <iostream>
using namespace std;
#define MAXSIZE 9//用于要排序数组个数最大值，可根据需要修改

//排序用的顺序表结构
typedef struct
{
    int r[MAXSIZE+1];//定义一个数组，用于存储要排序的数，r[0]作为临时变量
    int length;//用于记录顺序表的长度
}SqList;

//排序表的初始化
SqList *InitSqList(SqList *L)
{
    L=new SqList;
    L->length=MAXSIZE;//本例中长度是9
    cout<<"input list"<<endl;
    for(int i=1;i<=L->length;i++)
        cin>>L->r[i];
    return L;
}

//数组遍历输出
void PrintSqList(SqList *L)
{
    for(int i=1;i<=L->length;i++)
        cout<<L->r[i]<<" ";
    cout<<endl;
}

//希尔排序
void ShellSort(SqList *L)
{
    int i,j,k=0;
    int increment=L->length;//这是希尔排序的精华，将相隔某个增量的数组成子序列
    do
    {
        increment=increment/3+1;//增量选取是个难题，注意最后一个增量值必须等于1
        //每次循环结束后，相隔某个增量的数组成子序列必是顺序排列
        for(i=increment+1;i<=L->length;i++)
        {
            if (L->r[i]<L->r[i-increment])//相隔增量的某两个数比大小
            { 
                L->r[0]=L->r[i];//把r[i]挪到哨兵位置
                for(j=i-increment;j>0 && L->r[0]<L->r[j];j-=increment)//决定r[i]插入到哪个位置
                    L->r[j+increment]=L->r[j]; //把r[i]插入到合适的位置，不一定只减去一个增量，有可能是两个，三个...
                L->r[j+increment]=L->r[0]; //把哨兵位置的数插入到合适位置，不一定只与r[i]的相隔一个增量值
            }
        }
        cout<<++k<<"times ";
        PrintSqList(L);
    }
    while(increment>1);
}

int main()
{
    SqList *p=NULL;
    p=InitSqList(p);//初始化
    ShellSort(p);//排序
    cout<<"after sort"<<endl;
    PrintSqList(p);//输出
}

运行结果：

关于时间复杂度的分析，直接把书上的内容copy过来，不是很理解，欢迎讨论。