希尔排序(shell sort)

希尔排序,插入排序的一种变体。算法复杂度:

最好情况:O(n);平均:O(n^1.5);最坏:O(n^1.5或n^2)。

算法简介
回顾插入排序,元素朝着其正确位置一次一个位置的交换swap过去的(插入排序是从右向左比较哦),一个元素离其正确位置越远,其要移动的次数就越多,开销自然也就越大了。一个数组越接近有序,交换操作就越少,算法也就执行得越快。理想情况,当一个数组已经排序好了,那么插入算法能达到其最最理想的情况:O(n)。
希尔排序的基本思想就是:一个数据越接近有序,插入排序的效率就越高。我们定义一个step,将原数组分成step 个子数组(第一个子数组为:A[0],A[step], A[step * 2]...;第二个子数组为:A[0],A[0 + step],A[10+ step*2]...;…;第step个子数组为:A[step -1], A[step -1 + step],A[step-1+step*2]...)。将这些子数组进行插入排序,那么整个数组就越接近有序了。然后step= step / 2(减少步长,减少方式多样,这里/2),再对数据进行划分排序。。。直到step = 1进行最后一次插入排序。子数组的合并也有很多种方式:1,可以A[0]= A[0],A[1] = A[step],即先把一个子数组中的数据放完了,再来放第二个子数组中的数据;2,可以A[0] = A[0], A[1] = A[1],即在第一个子数组中取一个数放在数组中,再在第二个子数组中取一个数放在数组中个,再在第三个子数组中取一个数放在数组中。。。
关于step的选择(下面的内容摘自维基百科):

步长的选择是希尔排序的重要部分。只要最终步长为1任何步长序列都可以工作。算法最开始以一定的步长进行排序。然后会继续以一定步长进行排序,最终算法以步长为1进行排序。当步长为1时,算法变为插入排序,这就保证了数据一定会被排序。

Donald Shell 最初建议步长选择为\frac{n}{2}并且对步长取半直到步长达到 1。虽然这样取可以比\mathcal{O}(n^2)类的算法(插入排序)更好,但这样仍然有减少平均时间和最差时间的余地。 可能希尔排序最重要的地方在于当用较小步长排序后,以前用的较大步长仍然是有序的。比如,如果一个数列以步长5进行了排序然后再以步长3进行排序,那么该数列不仅是以步长3有序,而且是以步长5有序。如果不是这样,那么算法在迭代过程中会打乱以前的顺序,那就不会以如此短的时间完成排序了。

步长序列 最坏情况下复杂度
{n/2^i} \mathcal{O}(n^2)
2^k - 1 \mathcal{O}(n^{3/2})
2^i 3^i \mathcal{O}( n\log^2 n )

已知的最好步长序列由Marcin Ciura设计(1,4,10,23,57,132,301,701,1750,…) 这项研究也表明“比较在希尔排序中是最主要的操作,而不是交换。”用这样步长序列的希尔排序比插入排序堆排序都要快,甚至在小数组中比快速排序还快,但是在涉及大量数据时希尔排序还是比快速排序慢。

另一个在大数组中表现优异的步长序列是(斐波那契数列除去0和1将剩余的数以黄金分割比的两倍的进行运算得到的数列):(1, 9, 34, 182, 836, 4025, 19001, 90358, 428481, 2034035, 9651787, 45806244, 217378076, 1031612713, …)


代码
package sorting;

public class Data {
	
	/**
	 * generate an unsorted array of length n
	 * @param n length
	 * @param max the maximum integer element of this array
	 * @return an unsorted array consists of elements range from 1 to max
	 */
	public static int[] getArray(int n, int max){
		int[] result = new int[n];
		for(int i =0; i < n; i++){
			result[i] = (int)(Math.random() * max + 1);
		}
		return result;
	}
	/**
	 * print the array
	 * @param arg array
	 */
	public static void printArray(int[] arg){
		StringBuffer temp = new StringBuffer();
		for(int i = 0; i < arg.length; i++){
			temp.append(arg[i] + "\t");
		}
		System.out.println(temp);
	}
}
package sorting;

public class ShellSort {
	private int step = 4;

	public ShellSort(int step) {
		this.step = step;
	}

	public ShellSort() {

	}

	/**
	 * @param args
	 */
	public static void main(String[] args) {
		// TODO Auto-generated method stub
		ShellSort ss = new ShellSort();
		int[] data = Data.getArray(20, 100);
		System.out.print("Source data:\t");
		Data.printArray(data);
		ss.shellSort(data);
		System.out.print("Sorted data:\t");
		Data.printArray(data);

	}

	public void shellSort(int[] arg) {
		int curStep = step;
		int i = 0, j = 0, k = 0;
		int key = 0;
		for (; curStep > 0; curStep = curStep / 2) {
			for (i = 0; i < curStep; i++) {
				for (k = i + curStep; k < arg.length; k += curStep) {
					key = arg[k];
					for (j = k - curStep; j >= i; j -= curStep) {
						if (arg[j] > key) {
							arg[j + curStep] = arg[j];
							arg[j] = key;
						} else {
							break;
						}
					}
				}
			}
		}
	}

}


  • 2
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值