当前位置:   article > 正文

希尔排序算法

希尔排序

本章介绍排序算法中的希尔排序。内容包括:
1. 希尔排序介绍
2. 希尔排序图文说明
3. 希尔排序的时间复杂度和稳定性
4. 希尔排序实现
4.1 希尔排序C实现
4.2 希尔排序C++实现
4.3 希尔排序Java实现

转载请注明出处:希尔排序 - 如果天空不死 - 博客园


更多内容:数据结构与算法系列 目录

希尔排序介绍

希尔排序(Shell Sort)是插入排序的一种,它是针对直接插入排序算法的改进。该方法又称缩小增量排序,因DL.Shell于1959年提出而得名。

希尔排序实质上是一种分组插入方法。它的基本思想是:对于n个待排序的数列,取一个小于n的整数gap(gap被称为步长)将待排序元素分成若干个组子序列,所有距离为gap的倍数的记录放在同一个组中;然后,对各组内的元素进行直接插入排序。 这一趟排序完成之后,每一个组的元素都是有序的。然后减小gap的值,并重复执行上述的分组和排序。重复这样的操作,当gap=1时,整个数列就是有序的。

希尔排序图文说明

希尔排序代码(一)

  1. /*
  2. * 希尔排序
  3. *
  4. * 参数说明:
  5. * a -- 待排序的数组
  6. * n -- 数组的长度
  7. */
  8. void shell_sort1(int a[], int n)
  9. {
  10. int i,j,gap;
  11. // gap为步长,每次减为原来的一半。
  12. for (gap = n / 2; gap > 0; gap /= 2)
  13. {
  14. // 共gap个组,对每一组都执行直接插入排序
  15. for (i = 0 ;i < gap; i++)
  16. {
  17. for (j = i + gap; j < n; j += gap)
  18. {
  19. // 如果a[j] < a[j-gap],则寻找a[j]位置,并将后面数据的位置都后移。
  20. if (a[j] < a[j - gap])
  21. {
  22. int tmp = a[j];
  23. int k = j - gap;
  24. while (k >= 0 && a[k] > tmp)
  25. {
  26. a[k + gap] = a[k];
  27. k -= gap;
  28. }
  29. a[k + gap] = tmp;
  30. }
  31. }
  32. }
  33. }
  34. }

在上面的希尔排序中,首先要选取步长gap的值。选取了gap之后,就将数列分成了gap个组,对于每一个组都执行直接插入排序。在排序完所有的组之后,将gap的值减半;继续对数列进行分组,然后进行排序。重复这样的操作,直到gap<0为止。此时,数列也就是有序的了。

为了便于观察,我们将希尔排序中的直接插入排序独立出来,得到代码(二)。

希尔排序代码(二)

  1. /*
  2. * 对希尔排序中的单个组进行排序
  3. *
  4. * 参数说明:
  5. * a -- 待排序的数组
  6. * n -- 数组总的长度
  7. * i -- 组的起始位置
  8. * gap -- 组的步长
  9. *
  10. * 组是"从i开始,将相隔gap长度的数都取出"所组成的!
  11. */
  12. void group_sort(int a[], int n, int i,int gap)
  13. {
  14. int j;
  15. for (j = i + gap; j < n; j += gap)
  16. {
  17. // 如果a[j] < a[j-gap],则寻找a[j]位置,并将后面数据的位置都后移。
  18. if (a[j] < a[j - gap])
  19. {
  20. int tmp = a[j];
  21. int k = j - gap;
  22. while (k >= 0 && a[k] > tmp)
  23. {
  24. a[k + gap] = a[k];
  25. k -= gap;
  26. }
  27. a[k + gap] = tmp;
  28. }
  29. }
  30. }
  31. /*
  32. * 希尔排序
  33. *
  34. * 参数说明:
  35. * a -- 待排序的数组
  36. * n -- 数组的长度
  37. */
  38. void shell_sort2(int a[], int n)
  39. {
  40. int i,gap;
  41. // gap为步长,每次减为原来的一半。
  42. for (gap = n / 2; gap > 0; gap /= 2)
  43. {
  44. // 共gap个组,对每一组都执行直接插入排序
  45. for (i = 0 ;i < gap; i++)
  46. group_sort(a, n, i, gap);
  47. }
  48. }

下面以数列{80,30,60,40,20,10,50,70}为例,演示它的希尔排序过程。

第1趟:(gap=4)

当gap=4时,意味着将数列分为4个组: {80,20},{30,10},{60,50},{40,70}。 对应数列: {80,30,60,40,20,10,50,70}
对这4个组分别进行排序,排序结果: {20,80},{10,30},{50,60},{40,70}。 对应数列: {20,10,50,40,80,30,60,70}


第2趟:(gap=2)

当gap=2时,意味着将数列分为2个组:{20,50,80,60}, {10,40,30,70}。 对应数列: {20,10,50,40,80,30,60,70}
注意:{20,50,80,60}实际上有两个有序的数列{20,80}和{50,60}组成。
          {10,40,30,70}实际上有两个有序的数列{10,30}和{40,70}组成
对这2个组分别进行排序,排序结果:{20,50,60,80}, {10,30,40,70}。 对应数列: {20,10,50,30,60,40,80,70}

第3趟:(gap=1)

当gap=1时,意味着将数列分为1个组:{20,10,50,30,60,40,80,70}
注意:{20,10,50,30,60,40,80,70}实际上有两个有序的数列{20,50,60,80}和{10,30,40,70}组成。
对这1个组分别进行排序,排序结果:{10,20,30,40,50,60,70,80}

希尔排序的时间复杂度和稳定性

希尔排序时间复杂度
希尔排序的时间复杂度与增量(即,步长gap)的选取有关。例如,当增量为1时,希尔排序退化成了直接插入排序,此时的时间复杂度为O(N²),而Hibbard增量的希尔排序的时间复杂度为O(N3/2)。

希尔排序稳定性
希尔排序是不稳定的算法,它满足稳定算法的定义。对于相同的两个数,可能由于分在不同的组中而导致它们的顺序发生变化。
算法稳定性 -- 假设在数列中存在a[i]=a[j],若在排序之前,a[i]在a[j]前面;并且排序之后,a[i]仍然在a[j]前面。则这个排序算法是稳定的!

希尔排序实现

希尔排序C实现
实现代码(shell_sort.c)

  1. /**
  2. * 希尔排序:C 语言
  3. *
  4. * @author skywang
  5. * @date 2014/03/11
  6. */
  7. #include <stdio.h>
  8. // 数组长度
  9. #define LENGTH(array) ( (sizeof(array)) / (sizeof(array[0])) )
  10. /*
  11. * 希尔排序
  12. *
  13. * 参数说明:
  14. * a -- 待排序的数组
  15. * n -- 数组的长度
  16. */
  17. void shell_sort1(int a[], int n)
  18. {
  19. int i,j,gap;
  20. // gap为步长,每次减为原来的一半。
  21. for (gap = n / 2; gap > 0; gap /= 2)
  22. {
  23. // 共gap个组,对每一组都执行直接插入排序
  24. for (i = 0 ;i < gap; i++)
  25. {
  26. for (j = i + gap; j < n; j += gap)
  27. {
  28. // 如果a[j] < a[j-gap],则寻找a[j]位置,并将后面数据的位置都后移。
  29. if (a[j] < a[j - gap])
  30. {
  31. int tmp = a[j];
  32. int k = j - gap;
  33. while (k >= 0 && a[k] > tmp)
  34. {
  35. a[k + gap] = a[k];
  36. k -= gap;
  37. }
  38. a[k + gap] = tmp;
  39. }
  40. }
  41. }
  42. }
  43. }
  44. /*
  45. * 对希尔排序中的单个组进行排序
  46. *
  47. * 参数说明:
  48. * a -- 待排序的数组
  49. * n -- 数组总的长度
  50. * i -- 组的起始位置
  51. * gap -- 组的步长
  52. *
  53. * 组是"从i开始,将相隔gap长度的数都取出"所组成的!
  54. */
  55. void group_sort(int a[], int n, int i,int gap)
  56. {
  57. int j;
  58. for (j = i + gap; j < n; j += gap)
  59. {
  60. // 如果a[j] < a[j-gap],则寻找a[j]位置,并将后面数据的位置都后移。
  61. if (a[j] < a[j - gap])
  62. {
  63. int tmp = a[j];
  64. int k = j - gap;
  65. while (k >= 0 && a[k] > tmp)
  66. {
  67. a[k + gap] = a[k];
  68. k -= gap;
  69. }
  70. a[k + gap] = tmp;
  71. }
  72. }
  73. }
  74. /*
  75. * 希尔排序
  76. *
  77. * 参数说明:
  78. * a -- 待排序的数组
  79. * n -- 数组的长度
  80. */
  81. void shell_sort2(int a[], int n)
  82. {
  83. int i,gap;
  84. // gap为步长,每次减为原来的一半。
  85. for (gap = n / 2; gap > 0; gap /= 2)
  86. {
  87. // 共gap个组,对每一组都执行直接插入排序
  88. for (i = 0 ;i < gap; i++)
  89. group_sort(a, n, i, gap);
  90. }
  91. }
  92. void main()
  93. {
  94. int i;
  95. int a[] = {80,30,60,40,20,10,50,70};
  96. int ilen = LENGTH(a);
  97. printf("before sort:");
  98. for (i=0; i<ilen; i++)
  99. printf("%d ", a[i]);
  100. printf("\n");
  101. shell_sort1(a, ilen);
  102. //shell_sort2(a, ilen);
  103. printf("after sort:");
  104. for (i=0; i<ilen; i++)
  105. printf("%d ", a[i]);
  106. printf("\n");
  107. }

希尔排序C++实现
实现代码(ShellSort.cpp)

  1. /**
  2. * 希尔排序:C++
  3. *
  4. * @author skywang
  5. * @date 2014/03/11
  6. */
  7. #include <iostream>
  8. using namespace std;
  9. /*
  10. * 希尔排序
  11. *
  12. * 参数说明:
  13. * a -- 待排序的数组
  14. * n -- 数组的长度
  15. */
  16. void shellSort1(int* a, int n)
  17. {
  18. int i,j,gap;
  19. // gap为步长,每次减为原来的一半。
  20. for (gap = n / 2; gap > 0; gap /= 2)
  21. {
  22. // 共gap个组,对每一组都执行直接插入排序
  23. for (i = 0 ;i < gap; i++)
  24. {
  25. for (j = i + gap; j < n; j += gap)
  26. {
  27. // 如果a[j] < a[j-gap],则寻找a[j]位置,并将后面数据的位置都后移。
  28. if (a[j] < a[j - gap])
  29. {
  30. int tmp = a[j];
  31. int k = j - gap;
  32. while (k >= 0 && a[k] > tmp)
  33. {
  34. a[k + gap] = a[k];
  35. k -= gap;
  36. }
  37. a[k + gap] = tmp;
  38. }
  39. }
  40. }
  41. }
  42. }
  43. /*
  44. * 对希尔排序中的单个组进行排序
  45. *
  46. * 参数说明:
  47. * a -- 待排序的数组
  48. * n -- 数组总的长度
  49. * i -- 组的起始位置
  50. * gap -- 组的步长
  51. *
  52. * 组是"从i开始,将相隔gap长度的数都取出"所组成的!
  53. */
  54. void groupSort(int* a, int n, int i,int gap)
  55. {
  56. int j;
  57. for (j = i + gap; j < n; j += gap)
  58. {
  59. // 如果a[j] < a[j-gap],则寻找a[j]位置,并将后面数据的位置都后移。
  60. if (a[j] < a[j - gap])
  61. {
  62. int tmp = a[j];
  63. int k = j - gap;
  64. while (k >= 0 && a[k] > tmp)
  65. {
  66. a[k + gap] = a[k];
  67. k -= gap;
  68. }
  69. a[k + gap] = tmp;
  70. }
  71. }
  72. }
  73. /*
  74. * 希尔排序
  75. *
  76. * 参数说明:
  77. * a -- 待排序的数组
  78. * n -- 数组的长度
  79. */
  80. void shellSort2(int* a, int n)
  81. {
  82. int i,gap;
  83. // gap为步长,每次减为原来的一半。
  84. for (gap = n / 2; gap > 0; gap /= 2)
  85. {
  86. // 共gap个组,对每一组都执行直接插入排序
  87. for (i = 0 ;i < gap; i++)
  88. groupSort(a, n, i, gap);
  89. }
  90. }
  91. int main()
  92. {
  93. int i;
  94. int a[] = {80,30,60,40,20,10,50,70};
  95. int ilen = (sizeof(a)) / (sizeof(a[0]));
  96. cout << "before sort:";
  97. for (i=0; i<ilen; i++)
  98. cout << a[i] << " ";
  99. cout << endl;
  100. shellSort1(a, ilen);
  101. //shellSort2(a, ilen);
  102. cout << "after sort:";
  103. for (i=0; i<ilen; i++)
  104. cout << a[i] << " ";
  105. cout << endl;
  106. return 0;
  107. }

希尔排序Java实现
实现代码(ShellSort.java)

  1. /**
  2. * 希尔排序:Java
  3. *
  4. * @author skywang
  5. * @date 2014/03/11
  6. */
  7. public class ShellSort {
  8. /**
  9. * 希尔排序
  10. *
  11. * 参数说明:
  12. * a -- 待排序的数组
  13. * n -- 数组的长度
  14. */
  15. public static void shellSort1(int[] a, int n) {
  16. // gap为步长,每次减为原来的一半。
  17. for (int gap = n / 2; gap > 0; gap /= 2) {
  18. // 共gap个组,对每一组都执行直接插入排序
  19. for (int i = 0 ;i < gap; i++) {
  20. for (int j = i + gap; j < n; j += gap) {
  21. // 如果a[j] < a[j-gap],则寻找a[j]位置,并将后面数据的位置都后移。
  22. if (a[j] < a[j - gap]) {
  23. int tmp = a[j];
  24. int k = j - gap;
  25. while (k >= 0 && a[k] > tmp) {
  26. a[k + gap] = a[k];
  27. k -= gap;
  28. }
  29. a[k + gap] = tmp;
  30. }
  31. }
  32. }
  33. }
  34. }
  35. /**
  36. * 对希尔排序中的单个组进行排序
  37. *
  38. * 参数说明:
  39. * a -- 待排序的数组
  40. * n -- 数组总的长度
  41. * i -- 组的起始位置
  42. * gap -- 组的步长
  43. *
  44. * 组是"从i开始,将相隔gap长度的数都取出"所组成的!
  45. */
  46. public static void groupSort(int[] a, int n, int i,int gap) {
  47. for (int j = i + gap; j < n; j += gap) {
  48. // 如果a[j] < a[j-gap],则寻找a[j]位置,并将后面数据的位置都后移。
  49. if (a[j] < a[j - gap]) {
  50. int tmp = a[j];
  51. int k = j - gap;
  52. while (k >= 0 && a[k] > tmp) {
  53. a[k + gap] = a[k];
  54. k -= gap;
  55. }
  56. a[k + gap] = tmp;
  57. }
  58. }
  59. }
  60. /**
  61. * 希尔排序
  62. *
  63. * 参数说明:
  64. * a -- 待排序的数组
  65. * n -- 数组的长度
  66. */
  67. public static void shellSort2(int[] a, int n) {
  68. // gap为步长,每次减为原来的一半。
  69. for (int gap = n / 2; gap > 0; gap /= 2) {
  70. // 共gap个组,对每一组都执行直接插入排序
  71. for (int i = 0 ;i < gap; i++)
  72. groupSort(a, n, i, gap);
  73. }
  74. }
  75. public static void main(String[] args) {
  76. int i;
  77. int a[] = {80,30,60,40,20,10,50,70};
  78. System.out.printf("before sort:");
  79. for (i=0; i<a.length; i++)
  80. System.out.printf("%d ", a[i]);
  81. System.out.printf("\n");
  82. shellSort1(a, a.length);
  83. //shellSort2(a, a.length);
  84. System.out.printf("after sort:");
  85. for (i=0; i<a.length; i++)
  86. System.out.printf("%d ", a[i]);
  87. System.out.printf("\n");
  88. }
  89. }

上面3种实现的原理和输出结果都是一样的。下面是它们的输出结果:

before sort:80 30 60 40 20 10 50 70 
after  sort:10 20 30 40 50 60 70 80 
声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/知新_RL/article/detail/886006
推荐阅读
相关标签
  

闽ICP备14008679号