希尔排序

对于增量gap的理解

a[0]	a[1]	a[2]	a[3]	a[4]	a[5]	a[6]	a[7]	a[8]	a[9]
9	8	7	6	5	4	3	2	1	0

对于上述一个数组来说：
gap=5(初始值=length/2，其后为gap/=2)
gap=2
gap=1

例如gap=5,有：
①gap=组数：将该数组分为5组:(a[0]\a[5]) (a[1]\a[6]) (a[2]\a[7]) (a[3]\a[8]) (a[4]\a[9])
②gap=某组内相邻两数之间的间隔
③[0,gap)即0、1、2、3、4为每组第一个元素的下标。

实现代码（c++）

//希尔排序
void Shell_sort(int a[],int length)
{
	int gap=length/2;
	while(gap>0)
	{
		for(int k=0;k<gap;k++)//每个组
		{
			for(int i=k+gap;i<length;i+=gap)//某个组内，进行插入排序。k+gap是第二个元素。
			{
		    	int current=a[i];
		    	int position=i;//这里这个初始值很重要

		     	for(int j=i-gap;j>=0;j-=gap)//用while的话看起来会更清楚
				{
		 	    	if(a[j]>current)
					{
			     		position=j;
			     		a[j+gap]=a[j];
					}
				}
			a[position]=current;
			}

		}
		gap/=2;
	}
}

第二次再实现该排序算法的时候出现了错误。找了很久了原因，最后发现是每组第一个元素的角标会随着循环的进行而变化。导致在插入时，只与前面的第一个元素比较了一次，而不是一直比较到第一个元素。
修改后的代码如下（错误在注释中）

void Shell_sort2(int a[],int length)
{
	int gap=length/2;
	while(gap>0)
	{
		for(int k=0;k<gap;k++)//每个组
		{
			for(int i=k+gap;i<length;i+=gap)//某个组内，进行插入排序。k+gap是第二个元素。
			{
		    	int current=a[i];
		    	int position=i-gap;//这里这个初始值很重要

				while(position>=k&&a[position]>current)//这里是position>=k(k表示每组第一个元素的角标).
				{                                      //一开始用的是position>=i-gap，但是i值在循环过程中会改变，所以i-gap
					                                   //就不一定是第一个元素的角标。
					a[position+gap]=a[position];
					position-=gap;
				}
			a[position+gap]=current;
			}

		}
		gap/=2;
	}
}

第三次再实现该排序算法时，发现前两次编写的代码都执行了四次循环，而实际上值需要用三次循环就好了。

void Shell_sort(int *a,int n)
{
	int gap=n/2;
	while(gap)
	{
		for(int i=gap;i<n;i++)//每组的第二个元素。第一组a[gap]、第二组a[gap+1]
		{
			int current=a[i];
			int pre_index=i-gap;
			while(pre_index>=0&&a[pre_index]>current)
			{
				a[pre_index+gap]=a[pre_index];
				pre_index-=gap;
			}
			a[pre_index+gap]=current;
		}
		gap/=2;
	}
}

总结

①稳定性：不稳定。
②平均时间复杂度O(n logn)。

声明

参考源