希尔排序,插入排序的一种变体。算法复杂度:
最好情况:O(n);平均:O(n^1.5);最坏:O(n^1.5或n^2)。
算法简介
步长的选择是希尔排序的重要部分。只要最终步长为1任何步长序列都可以工作。算法最开始以一定的步长进行排序。然后会继续以一定步长进行排序,最终算法以步长为1进行排序。当步长为1时,算法变为插入排序,这就保证了数据一定会被排序。
Donald Shell 最初建议步长选择为并且对步长取半直到步长达到 1。虽然这样取可以比
类的算法(插入排序)更好,但这样仍然有减少平均时间和最差时间的余地。 可能希尔排序最重要的地方在于当用较小步长排序后,以前用的较大步长仍然是有序的。比如,如果一个数列以步长5进行了排序然后再以步长3进行排序,那么该数列不仅是以步长3有序,而且是以步长5有序。如果不是这样,那么算法在迭代过程中会打乱以前的顺序,那就不会以如此短的时间完成排序了。
步长序列 | 最坏情况下复杂度 |
---|---|
![]() | ![]() ![]() |
![]() | ![]() ![]() |
![]() | ![]() ![]() |
已知的最好步长序列由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;
}
}
}
}
}
}
}