【排序】希尔排序

希尔排序又称为“缩小增量排序”，是由希尔在1959年提出的。希尔排序是对插入排序的一种改进，在效率上较直接插入，冒泡，选择排序方法有较大改进。
基本思想：设定一个元素间隔增量gap，将参加排序的序列按这个间隔数gap从第1个元素开始依次分成若干个子序列。例如最开始gap=3，下面的序列可这样划分：
　　　　　　　原序列:{3 6 4 2 11 10}
　　　　　　　子序列1:3——2
　　　　　　　子序列2: 6——11
　　　　　　　子序列3: 4——10
将序列{3,6,4,2,11,10}划分为3个子序列。在子序列中采用其他排序方法(例如冒泡排序)。然后缩小增量gap进行划分，再分别对每个子序列进行排序。如此将“缩小增量gap–划分序列–将每个子序列排序”操作进行下去，直到间隔数增量gap=1为止。
由于排序时每一趟都是以不同的间隔对子序列进行排序，因此元素的移动在子序列中是跳跃的。间隔数gap越大，跳跃的跨度就越大。一般情况下，当增量gap缩小到1时，序列大都已经按值有序，不需要进行较多的元素移动就能达到排序的目的。
上述的元素序列{3,6,4,12,11,10,8,9}的希尔排序如下：
　　　　　　初始状态：{3,6,4,12,11,10,8,9}
　　　第1趟排序，gap=4：{3,6,4,9,11,10,8,12}
　　　第2趟排序，gap=2：{3,6,4,9,8,10,11,12}
　　　第3趟排序，gap=1：{3,4,6,8,9,10,11,12}
希尔排序的算法描述如下：

void sort(keytype k[],int n)
{
 	int i,j,flag,gap=n;
 	int tmp;
 	while(gap>1)
 	{
  		gap=gap/2;	//增量减小，每次减半 
  		do{		//子序列应用冒泡排序 
   			flag=0;
   			for(i=1;i<=n-gap;i++)
   			{
    				j=i+gap;
    				if(k[i]>k[j])
    				{
     					tmp=k[i];
     					k[i]=k[j];
     					k[j]=tmp;
     					flag=1;
    				}
   			}
  		}while(flag!=0);
 	}
}

上述算法描述希尔排序的过程。参数n为序列中元素的个数。将n赋值给变量gap，当gap>1时，循环地执行希尔排序的每一趟排序操作。在执行每一趟排序操作之前，gap的值都要减半，增量减小。在每趟的排序中使用的是冒泡排序，当然这里也可以使用其他的排序方法。
希尔排序算法的速度是一系列增量gap的函数，要对希尔排序的做法做出准确的分析并不容易，如何选取最合适的间隔数序列才能达到最优的排序效果是至今尚未解决的数学难题。
【实例】
编写一个C程序，实现数据序列{2,5,6,3,7,8,0,9,12,1}的希尔排序，要求从大到小，并输出排序后的数列元素。
【分析】
本序列中包含10个元素，因此间隔数gap的初始值为10/2=5。以后在进行每一趟排序之前，gap的值都要减半，直到gap=1为止。在这里要注意，本题要求将序列从大到小排列，因此要将每一趟的排序方法(上述算法中采用的是冒泡排序法)作适当的调整。

#include<stdio.h>
void sort(int k[],int n)
{
 	int i,j,flag,gap=n;
 	int tmp;
 	while(gap>1)
 	{
  		gap=gap/2;	//增量减小，每次减半 
  		do{		//子序列应用冒泡排序 
   			flag=0;
   			for(i=1;i<=n-gap;i++)
   			{
   	 			j=i+gap;
    				if(k[i]<k[j])
    				{
     					tmp=k[i];
     					k[i]=k[j];
     					k[j]=tmp;
    	 				flag=1;
    				}
   			}
  		}while(flag!=0);
 	}
}
main()
{
 	int i,a[11]={-111,2,5,6,3,7,8,0,9,12,1};
 	//初始化序列,a[0]可任意置数
 	printf("原序列为：\n");//显示原序列之中的元素
 	for(i=1;i<=10;i++)
  		printf("%d ",a[i]);
 	sort(a,10);//执行希尔排序 
 	printf("\n排序后的序列为：\n");
 	for(i=1;i<=10;i++)
  		printf("%d ",a[i]); //输出排序后结果
 	return 0;  
}

运行结果为：