【算法学习】线性时间排序-计数排序、基数排序和桶排序详解与编程实现

计数排序


计数排序假设n个输入元素中的每一个都是介于0到k之间的整数。此处k为某个整数(输入数据在一个小范围内)。

算法思想

计数排序的基本思想是对每一个输入元素x,确定出小于x的元素的个数。然后再将x直接放置在它在最终输出数组中的位置上。


由于数组中可能有相等的数,在处理时需要注意。


时间复杂度和空间复杂度分析


算法总时间Θ(k + n)。当k=O(n)时,计数排序的运行时间是Θ(n)。

空间复杂度是O(n+k)。需要两个辅助数组:存放排序结果的数组B[n],存放临时结果的C[k]。


计数排序是稳定的排序算法。


编程实现(CPP)

  1. //计数排序-《算法导论(第二版)》P98 8.2计数排序  
  2. //Author:江南烟雨  
  3. //E-Mail:xiajunhust@gmail.com  
  4.   
  5. #include <iostream>  
  6. #include <cstdlib>  
  7.   
  8. using namespace std;  
  9.   
  10. void CountSort(int *a,const int num,int *result)  
  11. {  
  12.     int MaxVal = -99999;  
  13.     for(int i = 0;i < num;++i)  
  14.     {  
  15.         if(MaxVal < *(a + i))  
  16.             MaxVal = *(a + i);  
  17.     }  
  18.     int *tempResult = new int[MaxVal + 5];//记录中间结果  
  19.     for(int i = 0;i < MaxVal  + 5;++i)  
  20.         *(tempResult + i) = 0;  
  21.     //result[i]记录数组中值等于i的元素的个数  
  22.     for(int i = 0;i < num;++i)  
  23.         *(tempResult + *(a + i)) = *(tempResult + *(a + i)) + 1;  
  24.     //result[i]记录数组中值小于等于i的元素的个数  
  25.     for(int i = 1;i < MaxVal + 5;++i)  
  26.         *(tempResult + i) = *(tempResult + i) + *(tempResult + i - 1);  
  27.     //注意,数组中可能存在相等的元素  
  28.     //将数组中各元素直接放入正确的位置  
  29.     for (int i = num - 1;i >= 0;--i)  
  30.     {  
  31.         *(result + *(tempResult + *(a + i))) = *(a + i);  
  32.         *(tempResult + *(a + i)) = *(tempResult + *(a + i)) - 1;  
  33.     }  
  34.   
  35.     delete[] tempResult;  
  36. }  
  37.   
  38. int main()  
  39. {  
  40.      int num = 7;  
  41.     int *a = new int[num];  
  42.     for(int i = 0;i < num;++i)  
  43.         *(a + i) = rand();  
  44.   
  45.     cout << "Before sort: " << endl;  
  46.     for(int i = 0;i < num;++i)  
  47.         cout << *(a + i) << " ";  
  48.     cout << endl;  
  49.   
  50.     int *result = new int[num + 5];  
  51.   
  52.     CountSort(a,num,result);  
  53.   
  54.     cout << "After sort: " << endl;  
  55.     for(int i = 1;i <= num;++i)  
  56.         cout << *(result + i) << " ";  
  57.     cout << endl;  
  58.   
  59.     delete[] a;  
  60.     delete[] result;  
  61. }  


基数排序


算法思想

基数排序是从低位到高位依次对所有的数进行排序。如果所有的数最高位数是d,那么先按最低有效位数字进行排序,得到一个结果。然后往高位重复这个过程。

需要注意的是,按位排序必须是稳定的排序算法。经常采用的是计数排序。


编程实现(CPP)

  1. //基数排序  
  2. //《算法导论(第二版)》P100 8.3 基数排序  
  3. //Author:江南烟雨  
  4. //E-Mail:xiajunhust@gmail.com  
  5.   
  6. #include <iostream>  
  7. #include <cstdlib>  
  8. #include <ctime>  
  9.   
  10. using namespace std;  
  11.   
  12. //得到某个整数第i位的数值  
  13. int getDigitNun(int a,int digit);  
  14. //按位排序  
  15. void DigitSort(int *a,int n,int digit,int *result);  
  16. //基数排序算法  
  17. void RadixSort(int *a,int n,int d);  
  18.   
  19. int main()  
  20. {  
  21.     int n = 7,i;  
  22.     int *a = new int[n];  
  23.     srand(time(NULL));  
  24.     for(i = 0;i < n;++i)  
  25.         *(a + i) = rand();  
  26.     //判断最大的数的位数  
  27.     int MaxVal = -1,d = 0;  
  28.     cout << "Before sort : " << endl;  
  29.     for(i = 0;i < n;++i)  
  30.     {  
  31.         cout << *(a + i) << " ";  
  32.         MaxVal = MaxVal < *(a + i) ? *(a + i) : MaxVal;  
  33.     }  
  34.     cout << endl;  
  35.     while(MaxVal > 0)  
  36.     {  
  37.         ++d;  
  38.         MaxVal /= 10;  
  39.     }  
  40.   
  41.     RadixSort(a,n,d);  
  42.   
  43.     cout << "After sort : " << endl;  
  44.     for(i = 0;i < n;++i)  
  45.         cout << *(a + i) << " ";  
  46.     cout << endl;  
  47. }  
  48.   
  49. //基数排序算法  
  50. void RadixSort(int *a,int n,int d)  
  51. {  
  52.     int *result = new int[n + 5];  
  53.     //循环执行按位排序操作  
  54.     for (int i =1;i <= d;++i)  
  55.     {  
  56.         DigitSort(a,n,i,result);  
  57.         for (int j = 0;j < n;++j)  
  58.         {  
  59.             *(a + j) = *(result + j + 1);  
  60.         }  
  61.     }  
  62.   
  63.     delete[] result;  
  64. }  
  65.   
  66. //得到某个整数第i位的数值  
  67. int getDigitNun(int a,int digit)  
  68. {  
  69.     while(--digit)  
  70.     {  
  71.         a /= 10;  
  72.     }  
  73.   
  74.     return a % 10;  
  75. }  
  76.   
  77. //按位排序  
  78. //这里采用选择排序  
  79. void DigitSort(int *a,int n,int digit,int *result)  
  80. {  
  81.     //记录中间结果  
  82.     const int num = 15;  
  83.     int *tempResult = new int[num];  
  84.     for(int i = 0;i < num;++i)  
  85.         *(tempResult + i) = 0;//初始化  
  86.   
  87.     //tempResult[i]记录数组中等于i的数的个数  
  88.     for(int i = 0;i < n;++i)  
  89.         *(tempResult + getDigitNun(*(a + i),digit)) = *(tempResult + getDigitNun(*(a + i),digit)) + 1;  
  90.   
  91.     //tempResult[i]记录数组中小于等于i的数的个数  
  92.     for(int i = 1;i < num;++i)  
  93.         *(tempResult + i) = *(tempResult + i) + *(tempResult + i - 1);  
  94.   
  95.     //将个元素直接放入正确的位置  
  96.     for(int i = n - 1;i >= 0;--i)  
  97.     {  
  98.         *(result + *(tempResult + getDigitNun(*(a + i),digit))) = *(a + i);  
  99.         *(tempResult + getDigitNun(*(a + i),digit)) = *(tempResult + getDigitNun(*(a + i),digit)) - 1;  
  100.     }  
  101.   
  102.     delete[] tempResult;  
  103. }  

时间复杂度和空间复杂度分析

给定n个d位数,每一个数位可能取值中数是k,如果所用的稳定的按位排序时间复杂度是Θ(n+k),基数排序时间复杂度是Θ(d(n+k))。空间复杂度O(n+k)。

当d为常数,k=O(n)时,基数排序有线性时间复杂度。


关于如何将每个关键字分解成若干数位方面,有另外一个定理:

给定n个b维数和任何正整数r<=b,基数排序能在Θ((b/r)(n+2^r))时间内对这些数进行排序。

这里,对一个值r<=b,将每个关键字看做是有d = floor(b/r)个数字,每个数字含有r位,再进行计数排序。

上述式子可以推导得到Θ(n)复杂度。


但是这并不意味着基数排序比基于比较的排序算法比如快排更好!因为隐含在记号中的常数因子是不同的。哪一个排序算法更好取决于底层机器的实现特性,比如快排同==排通常可以更有效地利用硬件缓存。同时还取决于输入数据。而且利用计数排序作为中间稳定排序不是原地排序。


桶排序

当输入数据符合均匀分布时,即可以以线性期望时间运行。即使输入不满足线性关系,桶排序也仍然可以以线性时间运行。只要输入满足这样一个性质,即各个桶尺寸的平方和与总的元素数呈线性关系。


桶排序的思想:


将区间[0,1)分成n个相同大小的子区间,或称为桶。然后将n个输入元素分布到各个桶中去。每个桶中的元素用一个链表来存储。


编程实现(CPP)

  1. //桶排序  
  2. //《算法导论(第二版)》P102 8.4 桶排序  
  3. //Author:江南烟雨(2013-03027)  
  4. //E-Mail:xiajunhust@gmail.com  
  5.   
  6. #include <iostream>  
  7. #include <cstdlib>  
  8. #include <ctime>  
  9. #include <cmath>  
  10.   
  11. using namespace std;  
  12.   
  13. //桶中链表节点数据结构  
  14. typedef struct StructLinkNode{  
  15.     double elem;  
  16.     struct StructLinkNode *next;  
  17. }LinkNode,*LinkNodePtr;  
  18.   
  19. //桶排序  
  20. void BucketSort(double *a,int n);  
  21. //删除一条链表  
  22. void deleteLinkList(LinkNodePtr head);  
  23.   
  24. int main()  
  25. {  
  26.     srand(time(NULL));  
  27.     int n = 8;  
  28.     double *a = new double[n];  
  29.     for(int i = 0;i < n;++i)  
  30.         *(a + i) = rand() * 1.0 / RAND_MAX;  
  31.   
  32.     cout << "Before sort : " << endl;  
  33.     for(int i = 0;i < n;++i)  
  34.         cout << *(a + i) << "  ";  
  35.     cout << endl;  
  36.   
  37.     BucketSort(a,n);  
  38.   
  39.     cout << "After sort : " << endl;  
  40.     for(int i = 0;i < n;++i)  
  41.         cout << *(a + i) << "  ";  
  42.     cout << endl;  
  43. }  
  44.   
  45. //桶排序  
  46. void BucketSort(double *a,int n)  
  47. {  
  48.     //存放链表的数组  
  49.     LinkNodePtr *linkListArr = new LinkNodePtr[n];  
  50.     //初始化  
  51.     for (int i = 0;i < n;++i)  
  52.     {  
  53.         linkListArr[i] = new LinkNode;  
  54.         linkListArr[i]->elem = -1;  
  55.         linkListArr[i]->next = NULL;  
  56.     }  
  57.   
  58.     //将n个输入元素依次放入n个桶中  
  59.     for (int i = 0;i < n;++i)  
  60.     {  
  61.         LinkNodePtr newNode = new LinkNode;  
  62.         newNode->elem = *(a + i);  
  63.         newNode->next = NULL;  
  64.   
  65.         //将新元素插入对应桶的链表的正确位置  
  66.         int index = floor(n * *(a + i));  
  67.         LinkNodePtr loopPtr = linkListArr[index]->next;  
  68.         LinkNodePtr prevPtr = linkListArr[index];  
  69.         while(loopPtr != NULL && *(a + i) > loopPtr->elem)  
  70.         {  
  71.             prevPtr = loopPtr;  
  72.             loopPtr = loopPtr->next;  
  73.         }  
  74.         newNode->next = loopPtr;  
  75.         prevPtr->next = newNode;  
  76.     }  
  77.   
  78.     int count = 0;  
  79.     for (int i = 0;i < n;++i)  
  80.     {  
  81.         LinkNodePtr loopPtr = linkListArr[i]->next;  
  82.         while(loopPtr != NULL)  
  83.         {  
  84.             *(a + count) = loopPtr->elem;  
  85.             ++count;  
  86.             loopPtr = loopPtr->next;  
  87.         }  
  88.     }  
  89.   
  90.     for (int i = 0;i < n;++i)  
  91.         deleteLinkList(linkListArr[i]);  
  92. }  
  93.   
  94. //删除一条链表  
  95. void deleteLinkList(LinkNodePtr head)  
  96. {  
  97.     if (NULL == head)  
  98.     {  
  99.         return;  
  100.     }  
  101.     deleteLinkList(head->next);  
  102.     delete head;  
  103. }  

时间和空间复杂度分析

时间复杂度是O(n)。

空间复杂度是O(n)。需要一个辅助数组来存放桶(链表)。


即使输入不满足均匀分布,桶排序也仍然可以以线性时间运行,只要输入满足这样一个条件:各个桶尺寸的平方和与总的元素呈线性关系。


桶排序是稳定排序算法。


转自:http://blog.csdn.net/xiajun07061225/article/details/8731179


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值