[数据结构与算法] - 排序 - 希尔排序

希尔排序

思路分析：

希尔排序是把记录按下标的一定增量分组，对每组使用直接插入排序算法排序；随着增量逐渐减少，每组包含的关键词越来越多，当增量减至1时，整个文件恰被分成一组，算法便终止。

希尔排序在数组中采用跳跃式分组的策略，通过某个增量将数组元素划分为若干组，然后分组进行插入排序，随后逐步缩小增量，继续按组进行插入排序操作，直至增量为1。希尔排序通过这种策略使得整个数组在初始阶段达到从宏观上看基本有序，小的基本在前，大的基本在后。然后缩小增量，到增量为1时，其实多数情况下只需微调即可，不会涉及过多的数据移动。

来看下希尔排序的基本步骤，在此选择增量gap=length/2，缩小增量继续以gap = gap/2的方式，这种增量选择可以用一个序列来表示，{n/2,(n/2)/2…1}，称为增量序列。希尔排序的增量序列的选择与证明是个数学难题，选择的这个增量序列是比较常用的，也是希尔建议的增量，称为希尔增量，但其实这个增量序列不是最优的。此处做示例使用希尔增量。

步骤：

1. 分组预排序，整体接近有序
2. 直接插入排序
   
   
   

处理经过gap分组后数据的处理方式：

方式一：先处理紫色组，再处理红色组，最后处理蓝色组

下面的代码的一部分有直接插入排序的影子。

void ShellSort(int* a, int n)
{
	//按gap分组数据进行预排序
	int gap = 3;
	//一组一组排
	for (int j = 0; j < gap; ++j)
	{
		//同组的排序
		//下面的代码类似与直接插入排序，只不过1变为了gap
		for (int i = j; i < n - gap; i+=gap)
		{
			int end = i;
			int x = a[end + gap];
			while (end >= 0)
			{
				if (a[end] > x)
				{
					a[end + gap] = a[end];
					end = end - gap;
				}
				else
				{
					break;
				}
			}
			a[end + gap] = x;
		}
	}
	
}

方式二：把gap分成的组一起处理

void ShellSort(int* a, int n)
{
	//按gap分组数据进行预排序
	int gap = 3;
	//把红色、蓝色、紫色一起炖
	for (int i = 0; i < n - gap; ++i)
	{
		int end = i;
		int x = a[end + gap];
		while (end >= 0)
		{
			if (a[end] > x)
			{
				a[end + gap] = a[end];
				end = end - gap;
			}
			else
			{
				break;
			}
		}
		a[end + gap] = x;
	}
}

总代码

1. 这里我们处理分组后的数据使用的是一锅炖的方式。
2. 然后我们要手动控制gap变量，使希尔排序的时间复杂度低些
3. 为什么while循环里gap要大于1？
   因为在gap等于1之前，代码 gap /= 2;可能已经把gap置为1，那么如果循环条件变为大于等于1，就会在重复循环一遍。
4. 为什么最后排序的时候gap要为1？
   因为我们之前是进行分组预排序，最后是要进行直接插入排序的，然后我们得知直接插入排序中，数组越是有序，它的时间复杂度越接近O(N)，这也是之前预排序的意义。

//希尔排序
//时间复杂度最好是O(N)
//最坏：F(N,gap) = (1 + 2 + ... + N/gap)*gap
//gap越大，预排越快，预排后越不接近有序
//gap越小，预排越慢,预排后越接近有序
void ShellSort(int* a, int n)
{
	//按gap分组数据进行预排序
	//动手控制gap
	//多次预排序(gap > 1) + 直接插入(gap == 1)
	int gap = n;
	while (gap > 1)
	{
		gap /= 2; //一种控制方式
		//gap = gap / 3; //另一种控制方式，但是最后一次gap可能达不到1
		//gap = gap / 3 + 1; //加1是为了最后gap可以达到1
		//多组一起炖
		for (int i = 0; i < n - gap; ++i)
		{
			int end = i;
			int x = a[end + gap];
			while (end >= 0)
			{
				if (a[end] > x)
				{
					a[end + gap] = a[end];
					end = end - gap;
				}
				else
				{
					break;
				}
			}
			a[end + gap] = x;
		}
	}
}