下面介绍三种线性时间复杂度的排序算法:计数排序、基数排序和桶排序。
这些算法是 用运算 而不是 用比较 来确定排序顺序的

image.png
image.png
image.png

这部分的描述,比《算法导论》好理解太多了。。 良心推荐

08 counting sort 计数排序 稳定

image.png
image.png
image.png
计数排序的核心在于将输入的数据值转化为键存储在额外开辟的数组空间中。作为一种线性时间复杂度的排序,计数排序要求输入的数据必须是有确定范围的整数。

应用:
如果一个问题的值域非常有限,应该考虑考虑会不会是计数排序

当输入的元素是 n 个 0 到 k 之间的整数时,它的运行时间是image.png。计数排序不是比较排序,因此不被 image.png的下界限制。

由于用来计数的数组 C 的长度取决于待排序数组中数据的范围(等于待排序数组的最大值与最小值的差加上1),这使得计数排序对于数据范围很大的数组,需要大量时间和内存。

例如:计数排序是用来排序 0 到 100 之间的数字的最好的算法,但是它不适合按字母顺序排序人名。但是,计数排序可以用在基数排序算法中,能够更有效的排序数据范围很大的数组。

通俗地理解,例如,有 10 个年龄不同的人,统计出有 8 个人的年龄比 A 小,那 A 的年龄就排在第 9 位,用这个方法可以得到其他每个人的位置,也就排好了序。当然,年龄有重复时需要特殊处理(保证稳定性),这就是为什么最后要反向填充目标数组,以及将每个数字的统计减去1。

算法的步骤如下:

  1. 找出待排序的数组中最大和最小的元素
  2. 统计数组中每个值为 i 的元素出现的次数,存入数组 C 的第 i 项
  3. 对所有的计数累加(从 C 中的第一个元素开始,每一项和前一项相加)(做前缀和)
  4. 反向填充目标数组:将每个元素 i 放在新数组的第 C[i] 项,每放一个元素就 C[i] 减去 1

时间复杂度:非比较排序 - 图9k代表待排序数据的值域大小
空间复杂度:非比较排序 - 图10
countingSort.gif

  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. const int N = 1e5 + 10, M = 110;
  4. int cnt[N], maxn = -1;
  5. int a[M], b[M], n;
  6. void counting_sort(){
  7. //cnt[x] 记录x出现的次数
  8. for (int i = 0; i < n; i++) cnt[a[i]]++;
  9. //维护出 比x小的有多少个
  10. for (int i = 1; i <= maxn; i++) cnt[i] += cnt[i - 1];
  11. //b[排第几个] = 原数组中的数值
  12. //从后往前拿,cnt从大到小,保证了稳定性
  13. for (int i = n - 1; i >= 0; i--) b[--cnt[a[i]]] = a[i];
  14. }
  15. int main(){
  16. cin >> n;
  17. for (int i = 0; i < n; i++){
  18. cin >> a[i];
  19. maxn = max(maxn, a[i]);
  20. }
  21. counting_sort();
  22. for (int i = 0; i < n; i++) cout << b[i] << ' ';
  23. puts("");
  24. return 0;
  25. }
  26. /*
  27. 5
  28. 3 4 2 1 5
  29. 1 2 3 4 5
  30. */

例题,CSP2019 入门组第一轮 计数排序

https://zhuanlan.zhihu.com/p/397169636

例题,P7072 [CSP-J2020] 直播获奖

  1. // 注意观察题目的数据范围,判断时间复杂度
  2. // 快排如何超时,怎么才能降低复杂度呢?

09 radix sort 基数排序 稳定

https://zh.wikipedia.org/zh-cn/%E5%9F%BA%E6%95%B0%E6%8E%92%E5%BA%8F
image.png
image.png

如果考虑和比较排序进行对照,基数排序的形式复杂度虽然不一定更小,但由于不进行比较,因此其基本操作的代价较小,而且在适当选择的B之下,k一般不大于logn,所以基数排序一般要快过基于比较的排序,比如快速排序。

基数排序是稳定的。

时间复杂度:非比较排序 - 图14
空间复杂度:非比较排序 - 图15
image.png
radixSort.gif基数排序.gif

  1. // n个数,每个数有k个关键字
  2. // 对这n个数,进行k个关键字排序
  3. #include <bits/stdc++.h>
  4. using namespace std;
  5. const int N = 1e5 + 10;
  6. int n, k; //n个数,每个数有k个关键字
  7. int maxn = 100; //每个关键字最大是100
  8. int cnt[N];
  9. struct node{
  10. int key[110]; //模拟100位数
  11. bool operator< (const node& W)const{
  12. for (int i = 1; i <= k; i++){
  13. if (key[i] == W.key[i]) continue;
  14. return key[i] < W.key[i];
  15. }
  16. return false;
  17. }
  18. }a[N], b[N];
  19. void counting_sort(int p){
  20. memset(cnt, 0, sizeof cnt);
  21. for (int i = 0; i < n; i++) cnt[a[i].key[p]]++;
  22. for (int i = 1; i <= maxn; i++) cnt[i] += cnt[i - 1];
  23. // 为保证排序的稳定性,此处循环i应从n到1
  24. // 即当两元素关键字的值相同时,原先排在后面的元素在排序后仍应排在后面
  25. for (int i = n - 1; i >= 0; i--) b[--cnt[a[i].key[p]]] = a[i];
  26. memcpy(a, b, sizeof a);
  27. }
  28. void radix_sort(){
  29. for (int i = k; i >= 1; i--)
  30. counting_sort(i);
  31. }
  32. int main(){
  33. cin >> n >> k;
  34. for (int i = 0; i < n; i++)
  35. for (int j = 1; j <= k; j++) cin >> a[i].key[j];
  36. //sort(a, a + n); //自定义一个快排,模拟数据进行验证
  37. radix_sort();
  38. for (int i = 0; i < n; i++){
  39. for (int j = 1; j <= k; j++) cout << a[i].key[j] << ' ';
  40. puts("");
  41. }
  42. return 0;
  43. }
  44. /*
  45. 3 3
  46. 33 44 55
  47. 11 22 33
  48. 33 55 66
  49. 3 3
  50. 33 11 22
  51. 11 22 33
  52. 88 22 10
  53. */
  1. // 下面这份代码,a[i]是1-index
  2. // 这在从桶里拿数的时候,是b[cnt[a[i].key[p]]--] = a[i];
  3. // 而不是b[--cnt[a[i].key[p]]] = a[i];
  4. // 这和排n, n-1是有关系的
  5. #include <bits/stdc++.h>
  6. using namespace std;
  7. const int N = 1e5 + 10;
  8. int n, k; //n个数,每个数有k个关键字
  9. int maxn = 100; //每个关键字最大是100
  10. int cnt[N];
  11. struct node{
  12. int key[110]; //模拟100位数
  13. bool operator< (const node& W)const{
  14. for (int i = 1; i <= k; i++){
  15. if (key[i] == W.key[i]) continue;
  16. return key[i] < W.key[i];
  17. }
  18. return false;
  19. }
  20. }a[N], b[N];
  21. void counting_sort(int p){
  22. memset(cnt, 0, sizeof cnt);
  23. for (int i = 1; i <= n; i++) cnt[a[i].key[p]]++;
  24. for (int i = 1; i <= maxn; i++) cnt[i] += cnt[i - 1];
  25. // 为保证排序的稳定性,此处循环i应从n到1
  26. // 即当两元素关键字的值相同时,原先排在后面的元素在排序后仍应排在后面
  27. for (int i = n; i >= 1; i--) b[cnt[a[i].key[p]]--] = a[i];
  28. memcpy(a, b, sizeof a);
  29. }
  30. void radix_sort(){
  31. for (int i = k; i >= 1; i--)
  32. counting_sort(i);
  33. }
  34. int main(){
  35. cin >> n >> k;
  36. for (int i = 1; i <= n; i++)
  37. for (int j = 1; j <= k; j++) cin >> a[i].key[j];
  38. //sort(a + 1, a + 1 + n); //自定义一个快排,模拟数据进行验证
  39. radix_sort();
  40. for (int i = 1; i <= n; i++){
  41. for (int j = 1; j <= k; j++) cout << a[i].key[j] << ' ';
  42. puts("");
  43. }
  44. return 0;
  45. }
  46. /*
  47. 3 3
  48. 33 44 55
  49. 11 22 33
  50. 33 55 66
  51. 3 3
  52. 33 11 22
  53. 11 22 33
  54. 88 22 10
  55. */

10 bucket sort(或 bin sort) 桶排 稳定

桶排序(Bucket sort)或所谓的箱排序,是一个排序算法,工作的原理是将数组分到有限数量的桶里。每个桶再个别排序(有可能再使用别的排序算法或是以递归方式继续使用桶排序进行排序)。

桶排序假设输入数据服从均匀分布,平均情况下它的时间代价为非比较排序 - 图19。与计数排序类似,因为对输入数据做了某种假设,桶排序的速度也很快。具体来说,计数排序假设输入数据都属于一个小区间内的整数,而桶排序则假设输入是由一个过程产生,该过程将元素均匀、独立地分布在[0, 1)区间上。

桶排序将[0, 1)区间划分为 n 个相同大小的子区间,或称为。然后,将 n 个输入数分别放到各个桶中。因为输入数据是均匀、独立地分布在[0, 1)上,所以一般不会出现很多数落在同一个桶中的情况。为了得到输出结果,我们先对每个桶中的数进行排序,然后遍历每个桶,按照次序把各个桶中的元素取出来即可。

桶排序是鸽巢排序的一种归纳结果。当要被排序的数组内的数值是均匀分配的时候,桶排序使用线性时间非比较排序 - 图20。但桶排序并不是比较排序,他不受到非比较排序 - 图21下限的影响。

算法过程:

  1. 设置一个定量的数组当作空桶子。
  2. 寻访序列,并且把项目一个一个放到对应的桶子去。
  3. 对每个不是空的桶子进行排序。
  4. 从不是空的桶子里把项目再放回原来的序列中。

桶排序是计数排序的升级版。它利用了函数的映射关系,高效与否的关键就在于这个映射函数的确定。为了使桶排序更加高效,我们需要做到这两点:

  1. 在额外空间充足的情况下,尽量增大桶的数量
  2. 使用的映射函数能够将输入的 N 个数据均匀的分配到 K 个桶中

同时,对于桶中元素的排序,选择何种比较排序算法对于性能的影响至关重要。

1. 什么时候最快
当输入的数据可以均匀的分配到每一个桶中
2. 什么时候最慢
当输入的数据被分配到了同一个桶中

时间复杂度:非比较排序 - 图22,k表示数字大小的范围
空间复杂度:非比较排序 - 图23

如果每个桶只存放一种数字,则不需要对桶内数字进行排序,直接从小到大输出桶内的数字即可。这种排序算法的总时间复杂度为非比较排序 - 图24。这种特殊的桶排序算法被称为计数排序,它只适用于 k 不是特别大的情况。

元素分配到桶中
image.png
对桶中元素排序
image.png

  1. // n个数,取值上限是w
  2. // 每个桶中的排序,用的是插入排序
  3. #include <bits/stdc++.h>
  4. using namespace std;
  5. const int N = 1010;
  6. vector<int> bucket[N];
  7. int n; // n个桶
  8. int w; // max{a[i]}
  9. int a[N];
  10. void insertion_sort(vector<int>& A) {
  11. for (int i = 1; i < A.size(); ++i) {
  12. int key = A[i];
  13. int j = i - 1;
  14. while (j >= 0 && A[j] > key) {
  15. A[j + 1] = A[j];
  16. --j;
  17. }
  18. A[j + 1] = key;
  19. }
  20. }
  21. void bucket_sort() {
  22. int bucket_size = w / n + 1;
  23. for (int i = 0; i < n; ++i) {
  24. bucket[i].clear();
  25. }
  26. for (int i = 1; i <= n; ++i) {
  27. bucket[a[i] / bucket_size].push_back(a[i]); //按规则落入桶中
  28. }
  29. for (int i = 0; i < n; i++) insertion_sort(bucket[i]);
  30. int p = 0;
  31. for (int i = 0; i < n; i++){
  32. for (int j = 0; j < bucket[i].size(); ++j) {
  33. a[++p] = bucket[i][j];
  34. }
  35. }
  36. }
  37. int main(){
  38. cin >> n >> w;
  39. for (int i = 1; i <= n; i++) cin >> a[i];
  40. bucket_sort();
  41. for (int i = 1; i <= n; i++) cout << a[i] << ' ';
  42. puts("");
  43. return 0;
  44. }
  1. // 特殊的桶排序,计数排序
  2. // 10 100
  3. // 1 2 1 2 1 2 3 4 5 4
  4. // 从小到大排序
  5. #include <bits/stdc++.h>
  6. using namespace std;
  7. const int N = 1010;
  8. int n, m;
  9. int a[N], book[N];
  10. int main(){
  11. cin >> n >> m;
  12. for (int i = 1; i <= n; i++){
  13. cin >> a[i];
  14. book[a[i]]++;
  15. }
  16. for (int j = 0; j <= m; j++)
  17. while (book[j]--) cout << j << ' ';
  18. puts("");
  19. return 0;
  20. }