数学建模社区-数学中国

标题: 常见的十种排序算法C++实现(附时空复杂度,稳定性分析) [打印本页]

作者: 杨利霞    时间: 2022-9-8 09:57
标题: 常见的十种排序算法C++实现(附时空复杂度,稳定性分析)
常见的十种排序算法C++实现(附时空复杂度,稳定性分析)6 D+ o$ j# O9 z4 l( o

# a" c8 o: Y+ p4 ]一、冒泡排序
$ {8 A+ V3 [4 F3 x/ B思路:每一轮排序使得当前无序的部分的最大值冒到当前无序数组的末尾,下一轮中无序部分的长度等于上一轮减一,原数组的末端不断的变成有序的部分。
- a9 {% |9 r* t5 m3 U* D. {7 Z# s# Q/ G7 L% |* e
void bubble_sort(vector<int> &nums)
- g5 Z8 k7 E: Q8 g{4 V/ ~+ N7 G* K) L0 r) V3 L
    // 用于标记当前轮是否有交换发生
) K1 G; n3 Q9 g8 T4 w6 D( w2 J    bool swaped;, J! y1 B3 P! R  o9 B8 d6 k3 ^
    // 从第一个数到倒数第二个数(n个数需要冒泡n-1次)7 x! E3 ]2 d* [8 s, R/ O
    for(int i = 0; i < nums.size() - 1; ++i)
# C& G6 e/ @! U+ v6 F2 G$ C+ ~    {
) b3 O8 n; e" B8 @2 Q5 g        swaped = false;
( ^1 N2 C: z5 Z9 Z" E8 `        // 因为是当前的和后面的比,所以结束条件是一直到冒泡完毕的元素的前一个
; N: D$ Y: ?4 P2 m& D; A; c' Q4 m        for (int j = 0; j < nums.size()-1-i; ++j)9 |: z" M3 y; k
        {    ! s& o# c; b: A! s8 L! p% g+ b( p& F5 m
            // 通过交换两个元素,使得最大的冒到数组的末尾去,并且标记交换过
( Y/ \8 P& l4 a7 P            // 若要从大到小排序,只需要把这个 '>' 改成 '<' 即可. M" V0 T$ I% X0 h& P
            if (nums[j] > nums[j+1])! |, X. J0 e5 v6 U
            {7 _9 A( g. g( ~% N
                int tmp = nums[j];. r5 G- U: b$ e8 I
                nums[j] = nums[j+1];
4 m) p  w( n7 @/ G& E                nums[j+1] = tmp;  Q$ Y' M$ e% {& M# j
                swaped = true;! P. W. R0 H1 A
            }( H; y  F- T! U( M8 T/ H' ?
        }, `! o& |/ z9 U0 a3 L# }# H
        // 若该轮没有交换发生,说明数组已经有序,跳出循环
9 v/ k) I9 f  o" f) j        if (!swaped)# S/ k" K! Z4 l2 ^' F- O3 _- U  u; C
            break;$ J/ C/ W7 I5 S* ^) S, C0 @+ ?
    }/ y9 D4 D7 Z: `6 R' ?: p
}* [1 ^4 ^) s1 h' P  m! j

, x/ l2 o$ {: Q9 ~) z/ X8 b  Q分析:时间复杂度最坏和平均情况是O(),最好的情况是已经有序所以是O(n),空间复杂度O(1),且因为两个相同的元素是不会交换位置的,因此是一种稳定的排序。
/ E+ v3 Y0 X: C# X& X9 U( x
1 Z% t' z: g7 i( s# c二、选择排序. |5 ]. o+ |$ ]6 K6 |
思路:每一轮都将未排序数组中最小的/最大的选择出来,然后与未排序数组的第一个交换。接着未排序数组的长度就应该减一,然后继续做这个操作,直到未排序数组的长度为1* k. N7 [' p0 M5 b; b7 ]/ }
  _/ \' c) C9 C
void select_sort(vector<int> &nums)& B# e; d4 n4 F( Q
{   
$ x1 \" E. u" b' n% A    // 用于存储每一轮中最小的数的下标# L; W( R  a4 \2 c) o2 Y
    int min_idx;
8 U0 ~0 E+ u1 I* Q* k, |& v# W    // 循环从第一个数开始到倒数第二个数$ U& o, x+ {/ z  k  @- \
    for(int i = 0; i < nums.size() - 1; ++i)
4 ?& \; f9 M  c7 ?8 Y# ~    {
5 y; N. y  d, `" @# ~# P- u        min_idx = i;8 C5 P/ Q8 ^; m& a1 F# l6 F
        // 在未排序部分寻找最小的数
$ C1 O! ~! z; v: c9 \# d9 [        for (int j = i+1; j < nums.size(); ++j), T  s  X; }2 f3 q* [5 _( P% b
        {" x# H& L( h1 L8 v4 m
            // 这里将 '<' 换成 '>' 就可以从大到小排序
9 {, f/ t. {6 g4 j, j            if (nums[j] < nums[min_idx])
% E) t+ T2 W6 V/ t6 `                min_idx = j;
# ]$ X" F, c" w$ i        }
) E$ C' }0 p' S6 [& C        // 将未排序部分最小的数与未排序部分的第一数交换位置0 R0 K. U" E& B& g2 ?
        int tmp = nums[min_idx];
7 J' Z9 J5 o( H        nums[min_idx] = nums;
7 w+ x3 K: M" M        nums = tmp;
& X" Q0 G( i! I8 n% M    }( a8 U2 [; X* \( J' u. }; @; r
}( \# l8 ]/ E) s' ~0 V7 e6 q; }" C

  ^9 x. W0 X. {; j- W6 z0 g7 Y  B4 G分析:时间复杂度最坏、最好和平均情况都是O(),空间复杂度O(1)。因为每次选择最小元素要与未排序部分的第一个元素进行交换,因此可能会出现两个相同元素相对位置发生改变的情况,因此者是一种不稳定的排序。
7 Y& W+ t& j6 y( m! z
' i3 ^: N* A* ?# {2 b, B9 ]三、插入排序/ L4 {; H  B: ?2 @! U
思路:从第二个元素作为当前位置,往前看,判断前面的元素是否大于当前位置,是的话就把前面的元素往后移动一个,依次类推,直到到了数组头或者碰到一个元素比当前元素小则循环终止,此时的位置就应该是当前位置元素的排序合适的位置。接着向后移动,把第三个元素作为当前位置,一直到最后一个元素。5 z- Y0 z2 ]9 m
: k! U- [; |* W* x. }
void insert_sort(vector<int> &nums)
: e( V6 Y  u. f2 F0 _{: N4 k( L$ ?* U1 }; ^1 x
    // 从第二个元素到最后一个元素4 J* X0 D+ f7 J: G7 U6 Y! q1 Y
    for(int i = 1; i < nums.size(); ++i)
$ _3 g7 O, O/ x$ m; V  }& n; {    {
" K3 t! N% X6 P+ Z9 n2 F        // tmp存下当前元素的值,将i赋给j
0 G6 v7 A( {- W: s        int tmp = nums;8 g7 g; V0 b& `
        int j = i;/ B) K& n9 D$ w+ p) W1 t
        // j从当前元素往前找,一直找到第二个元素位置
8 k# e7 v1 O, e        for (; j > 0; j--), M2 X* C" J+ c- h6 ?1 f- W, s* ?
        {$ N6 b7 W9 X: O3 l# W: a
            // 判断tmp是否小于j-1位置的元素,是的话j-1的元素往后移动0 ?6 |) v7 |! }% f, \' |& b
            // 要从大到小排序,只需要把 '<' 改为 '>' 即可
) ?; h% e. \' M7 q0 z; `            if (tmp < nums[j-1])# I( c. c" U% y% ]* G9 Y4 Y" q
                nums[j] = nums[j-1];9 V) h% j. }, ~7 j( O3 q
            // 否的话,j位置就应该是tmp的最终合适的位置
; y: y4 T& m  t            else. p& ~/ Y8 v. ]! `3 W
                break;$ V- T2 S5 r4 t% q- C8 s
        }- d& I+ ^4 E0 n- p
        // 将tmp赋给j位置的元素8 S9 I3 z; o) \
        nums[j] = tmp;* e, t8 f3 c( `. [
    }5 b( I  _  |# k5 S
}
1 @* K: _  r; @; L6 F/ ~8 O& S0 i" n! D; o3 J
分析: 时间复杂度最坏和平均情况都是O(),最好的情况是O(n)(最好情况就是已经有序,因此内层循环只需要判断一次就可以break),空间复杂度O(1)。因为是依次往前找的,元素移动也是相邻位置的移动,因此这个方法是稳定的。
: g6 k* _- h4 \' C- o! V' d/ s- d( z- w2 A1 c
四、希尔排序
2 R; s. H& _6 t" Q思路:希尔排序是插入排序的一种扩展,它不是一个的找然后插入。而是将与原数组根据间隔分成不同的组,比如十个数{1,2,3,4,5,6,7,8,9,10}可以设置gap为2,这样就分成了{1, 3, 5, 7, 9}和{2, 4, 6, 8, 10}两个组,分的组数就等于gap的数量,每组的包含的元素数量就等于数组的长度除以gap。通过这样操作,再最后gap为1的时候,数组就是基本有序的了,因此这时只需要很少的操作就可以将整个数组变为有序的。希尔排序gap的设置很关键,设置一组好的gap序列往往可以很大程度上提高希尔排序的效率) D% s9 N( U3 l! d  H

! y7 h. H  D* b# |3 Z1 Svoid shell_sort(vector<int> &nums)
# U  m  B4 b5 b{" Z2 W3 F: M9 G; W: S4 {
    // 初始将gap设置为数组长度的一半,每次gap等于上次gap的一半,一直循环一直到gap等于1# j8 _* _6 N! ?4 l! C
    for(int gap = nums.size()/2; gap > 0; gap/=2)
4 {) W- k3 D8 ^1 l* n    {8 E5 n, f9 a1 N8 W  n
        // 每一组的元素数量
/ G' o. Y" C4 C        int ele_sz = nums.size() / gap; & F$ k9 l2 i. U
        // 遍历每一组0 \& _# {# a5 T9 y  W9 ~
        for (int i = 0; i < gap; ++i)
- z6 U3 e9 T% F6 Z- A: d        {1 ?! c* C6 @6 d
            // 对每一个组进行插入排序,下面的代码和插入排序那基本一样2 `# s9 B, ^1 s8 `) p' r- e, c
            for (int j = 1; j < ele_sz; ++j)6 K; s" _* W- u5 z
            {  Y9 F7 V+ b$ j& H3 ?! c
                int tmp = nums[j*gap];
7 g3 W; C% \& @1 {) N                int k = j;
" {" X, d5 h* ?$ ~) G# x4 j4 n                for (; k > 0; --k)
" X% ?. O. v/ O% a                {
* \" ]9 A, D7 o                    // 注意这里的下标要符合希尔排序间隔\1 i, r3 t+ w* E: `7 @, c* L
                    // 要从大到小排序,只需要把 '<' 改为 '>' 即可$ V$ X9 ?  M) T/ B
                    if (tmp < nums[(k-1)*gap]); K  L. K6 a: ^: t
                        nums[k*gap] = nums[(k-1)*gap];
) P( s% y4 B/ c; l0 Y                    else
" X' C9 D% l& F" o" t2 r7 W                        break;
: o; B& j  n/ e+ C; a                }4 I& Z2 {0 n4 y, J8 O
                // 找到了最终插入的位置
  f' t; a5 H/ R- I                nums[k*gap] = tmp;
4 ]* v! l9 `0 N2 p% `            }# y. `. L, K* M' X8 k- ?
        }5 p9 ?; q* [2 k0 I: \
        // 因为每次gap等于上次一半,当gap为1的时候就陷入死循环,所以这里设置一个判断用于跳出
' A0 D) U. `1 d        if (gap == 1)$ x. S+ I( B: b9 H1 D0 A
            break;  b2 }. K0 n% Z- [
    }" `* r5 m' p0 K" s6 n" ?9 ]
}
& ?( B% t6 q$ q* a8 x
  |" i, q+ j& B8 ~* M$ l8 F$ `分析:希尔排序最好情况的时间复杂度为O(nlogn),而最坏情况和平均情况要根据它的增量序列来确定,好的增量序列可以把最坏情况的时间复杂度控制为O()。希尔排序是原地排序因此空间复杂度为O(1)。由于增量序列的存在,使得比较的元素存在间隔,因此就可能导致相同元素的相对位置发生改变,因此希尔排序是不稳定的。
2 S2 [" Y% o* y. N1 I) r
- Q( x2 @$ r6 H8 T) u, v9 j五、归并排序
/ }" a1 F' B* ]- B思路:归并排序是分治法的典型应用之一。首先先把大的数组拆分成两半,再继续拆,直到只剩一个数位置,拆完之后要合并,合并的时候就需要对两边的数进行排序,合并完成的时候就变成了有序的数组。以此类推,对左边的已经有序和右边已经有序的数组进行合并,最终得到了整体有序的数组。/ O/ d5 W0 s: W+ m5 a) n; E2 O1 {
8 E  i; F4 f' }( j0 d* _+ ]6 m" Z' F6 g
// 合的过程" }. }1 @# H* p7 i  `
vector<int> merge_vec(vector<int> ivec1, vector<int> ivec2)
3 I( c- }7 F0 X, s5 K{/ \5 _. e. o& c6 G7 K9 i' B* r8 ?
    // 定义结果数组res, i和j分别用于遍历ivec1和ivec29 _3 I% l( q0 {5 |
    vector<int> res;4 o, m0 V) w  e0 g& ?1 T) v: ?
    int i = 0, j = 0;
6 h0 f# G4 N6 Y8 M1 y' L) r5 {    // 当两个数组都没遍历到底/ ^) N( y0 o( V6 |# r1 R
    while(i < ivec1.size() && j < ivec2.size())" d5 I5 @7 [0 w
    {
' a* d5 P/ }! E, ~9 t% L6 u        // 谁的当前数字小,谁就加入结果数组,并且下标后移
$ H+ Z9 O7 d2 c& d# E/ I        // 这里将 '<' 改为 '>' 就可以从大到小排序
1 Q5 {# P$ G% ~% ?        if (ivec1 < ivec2[j])
# V& J& S6 U) n. G* s            res.push_back(ivec1[i++]);! T$ Q" k6 x. E* ^
        else$ u* v( R3 e7 D
            res.push_back(ivec2[j++]);
* H0 ]% U8 c$ t- M) @' Y# Q* k    }) w8 g3 r8 E0 [* U5 N  M0 z- L
    // 处理未遍历完的数组,直接全部加入最后面3 I  U2 F4 y: W2 A. _3 m4 M  R
    while (j < ivec2.size())
& f% w9 r6 Z  H# q" A        res.push_back(ivec2[j++]);7 F! p$ x/ M' Z4 R( X: \
    while (i < ivec1.size())5 q9 K1 @/ v  e4 t! M7 k4 Z
        res.push_back(ivec1[i++]);
3 w; p0 J/ q0 `. N+ m3 X    // 返回合并好的数组
6 y! {# v8 ]- D$ ^    return res;4 n& p% P5 U. ?/ |, H: t
}
3 N2 T7 A. d# o& B) B; z% A2 E2 P
, s& L0 P6 T- j: C$ t& a& {// 分的过程1 [9 P! C2 o  Z
vector<int> merge_sort(vector<int> nums)5 B* l! U& ?% i& f. g
{5 n. V: q3 D( j& o8 r5 B
    // 递归结束条件是只剩一个元素,则无序再分,直接返回$ R: H' t9 f( g) _
    if (nums.size() == 1)
' y/ D4 T/ g8 C        return nums;. p  C& n5 s! V5 b
6 H2 T0 i3 g! h6 M1 {$ @9 y
    // 将大数组分治为左边的和右边的1 d- B9 i# d1 F; D* e
    vector<int> ivec1 = merge_sort(vector<int> (nums.begin(), nums.begin()+nums.size()/2));
7 I; ]6 S4 C1 N5 M    vector<int> ivec2 = merge_sort(vector<int> (nums.begin()+nums.size()/2, nums.end()));! l" w# l; S. G% X
    // 返回合并后的大数组, c$ \$ h( b5 U5 h- P6 r
    return merge_vec(ivec1, ivec2);
6 y; g5 n0 A; t, v}
: w  h3 B* b% d: ]: b7 B# c2 Z! \5 n2 y5 Q& V/ S8 {$ ~; ~
分析:拆的过程时间复杂度是O(logn),合的过程时间复杂度是O(n),因此总的最好、最坏和平均时间复杂度是O(nlogn)。这种排序方法并非是原地排序,所以其空间复杂度是O(n)。这种方法并不会破坏相等元素之间的相对位置因此这是一种稳定的排序算法。
$ t* @& ^/ f- e& R, I( C6 B: w7 l% o* V, W, K; `$ M1 n* ^
六、快速排序
4 u/ O2 w1 t4 w# F, Q/ m+ U2 }思路:快速排序同样也是利用分治法的思想,但是快速排序可以实现原地排序因此无序额外的空间。快速排序和归并排序的分的过程不同,快速排序分是根据小于和大于某个元素来分成两部分,而归并排序总是以中间的地方来划分。
4 l3 N$ D0 {' H+ I* Y5 e4 }0 M' L" ^  X
注意:算法中每一轮quicksort会确定一个元素的最终位置,我们递归处理这个元素左边的剩余数组和右边的剩余数组。lp和rp在走的时候,需要注意的只要rp指向的元素大于等于pivot就一直往左,lp指向的元素小于等于pivot就一直往右,若不加入等于这个判断条件可能就会在中途停下来而陷入死循环。而且pivot的选取也很关键,我们是以pivot大小来划分左右的,常见的方法可以以最左边或者最右边的元素作为pivot(但是这样在已经有序的数组中效率很差),或者随机选择pivot。
$ W3 w' F1 r8 A' _9 U
) U- \) I8 r) g5 Ovoid quicksort(vector<int> &nums, vector<int>::iterator lp, vector<int>::iterator rp)  A7 J7 Q% W- P  a! Q. g5 A
{' ]7 v$ d$ W; V& h2 X
    // 若左右相等,说明只有一个元素,就不需要排序4 V$ m, @. C8 H; Q! M# m
    if (rp - lp > 0)7 O# h7 N  o, Y) E  f8 J! M9 ?# m
    {
2 s, {$ Y" \6 a        // old和orp存取起始位置和结束为止" m0 G: S" u2 W' F1 |% z
        auto olp = lp, orp = rp;( Y9 q( {  z% d

# C) p& o# j0 T7 c0 M        // 下面两行代码可以要也可以去掉,其作用是随机选定一个范围在lp和rp之间的数作为pivot
- U+ j5 s1 ]" g3 @) r! o        // 将这个数与最左边的数调换位置就可以通过*lp得到这个数
4 k; G9 x" z) E        int random = rand() % (rp - lp + 1);- s# \- h: @9 e! G+ A3 J
        swap(*lp, *(lp + random));6 K& @& H6 C2 `/ B6 _8 W' H
+ V+ s( Y2 ]) |7 B# P3 O
        int pivot = *lp;
; C: d3 H6 [3 D- J% d; X, {- b1 S9 C2 Q3 i/ s) @2 H+ k) U
        // 当lp和rp没遇到的时候2 [  A2 `+ [7 I7 A( {
        while (lp != rp)9 ?4 Y- \% r$ H6 `3 t- u7 k
        {
( }6 E; u* A( f7 f( T9 g9 _8 v            // 只要rp指向的元素大于等于pivot而且lp和rp没遇到,rp自减
/ q" q0 Q( [4 I% ~* p5 m            // 将这里的 '>=' 换成 '<=' 并且把下面的 '<=' 换成 '>=' 就是从大到小
- R! x  [6 n" U. y- O            while (*rp >= pivot && lp != rp)
9 W7 v4 h5 J  H& ^% s                --rp;0 H# Q8 T+ A, o& A0 y- w
            // 只要lp指向的元素小于等于pivot而且lp和rp没遇到,lp自增5 N" ^0 H; m  a! z( B
            while (*lp <= pivot && lp != rp)
2 c4 M9 g( _1 A  l& E* ]# ?                ++lp;
; ], g/ e' z/ i3 V4 L            // 停下来的时候交换值
5 A- b  m8 K/ y1 B            swap(*lp, *rp);% _: m& q" ?: u7 G1 T
        }5 K* v, G# {  B' B
6 a: l3 l3 C, U% M, Q
        // lp遇到rp说明比pivot小的已经放在左边,比pivot大的已经放在右边了. k) c# N; L2 z  x# L" }! y, E5 B3 ^/ V# S
        // 因为while循环中rp先走,说明rp所指的元素一定是小于pivot或者正好是pivot( v# y  {: d# ?5 q- f) ~0 u
        // 因此需要交换二者的值,同时pivot就被放在了正确的排序位置上
# @$ U- O# ~" n8 F5 r- W2 s        swap(*olp, *rp);
0 }2 H# }. A& n        // 递归处理pivot的左边和右边* F& M+ M+ ]; S' ~
        quicksort(nums, olp, lp-1);5 U1 Y  _; `7 H4 b# t& C
        quicksort(nums, rp+1, orp);- f1 {! h$ W6 Q
    }
  F$ i2 h5 g% i' f  I, X}- R( H1 a; H5 A; V" k* s' T" l
( F* |- R% \, O& V! a
; X9 W9 n) x" T6 G
void quick_sort(vector<int> &nums)4 g6 C) Y% B. F7 I; N- y; ]
{: X1 H4 D9 o1 {- s% c! d$ K" A
    quicksort(nums, nums.begin(), nums.end()-1);1 T* |: q3 Q2 G( g6 Y3 [& g. a
}
2 y) z7 l" B: g8 ]- h+ Q# j5 _1 p7 r' Z' e# i7 C: d
分析:快速排序算法的时间复杂度分是O(logn),治是O(n),因此其最好和平均的复杂度是O(nlogn)。考虑一种情况,当数组有序而且我们选取pivot的方法是选最左边的或者最右边的元素时复杂度退化为O()。原地排序因此不需要额外空间,但是由于是递归函数,因此需要栈来保存状态,最好和平均情况下需要O(logn),最坏的情况下需要O(n)。由于快速排序涉及到随机位置值的交换,因此快速排序是一种不稳定的排序方法。+ ?6 i  a( e# P; M% x% [# d2 |

. b9 d2 o, C$ ?) g, x七、堆排序/ |! l! B, l1 H$ u
思路:堆排序顾名思义需要一个堆,堆是一个完全二叉树,可以分为大顶堆和小顶堆,其中大顶堆是满足子节点的值小于父节点的堆,小顶堆是满足子节点的值大于父节点的堆。建好了堆之后我们就可以通过取堆顶元素再不断调整堆的方法来得到排序结果。因为堆是完全二叉树,因此其父节点和子节点的关系容易得到,所以我们可以用数组来模拟一棵树(一个堆),下面这个视频将堆排序讲的很清晰。, D+ q! ]" F5 q6 u2 g7 p
9 I3 l1 N! B" U  K8 G
【从堆的定义到优先队列、堆排序】 10分钟看懂必考的数据结构——堆_哔哩哔哩_bilibili
6 n; k- X, v0 P+ g6 i$ ^7 X9 g$ E# B8 {" [4 s
建堆:建立堆的方式可以分为从下往上和从上往下。以大顶堆为例,从下往上的话,对于每个父节点,判断其子节点值是否存在大于父节点的情况,是的话就交换这父节点和子节点的值,然后继续调整子节点(若子节点也存在节点的话则继续往下调整)。如果是从上往下建堆,就相当于是堆从从空的开始,不断插入元素再进行调整,判断当前插入的选择值是否大于其父节点的值,若是的话则交换当前节点与父节点,并继续向上调整直到无需调整或者到根节点。可以判断,从下往上的方式其时间复杂度为O(n),而从上往下调整的方式其时间复杂度为O(nlogn)。同时需要注意的是,对于同一个数组,采用从下往上和从上往下的方式建立可能会导致节点的排列方式不一样,但是都满足大顶堆/小顶堆的性质。# N+ E* H7 U4 o0 g, z2 k1 K- ^
. L3 @; _& |, u
下面是从下往上建大/小顶堆的代码:
. X; o) G1 j( k3 k) ^3 U" ~% f+ Y$ ~, M( R: c/ C
// 接收参数为数组nums,当前调整节点下标i,当前堆的大小sz
7 U" k8 G" [" }4 jvoid adjust_from_buttom(vector<int> &nums, int i, int sz)
3 j5 k8 W$ v) w& X{  m7 @4 w$ L1 ^& \# o/ u
    // lpos和rpos分别代表左孩子和右孩子,largest用来存储左右孩子中较大的那个的下标. u/ B" H, @0 O' G
    // 建小顶堆需将largest改为minest,此时这个变量用于存储较小的元素( x# D9 A$ a: Z+ X  L: y( j
    int lpos = 2*i+1, rpos = 2*i+2, largest = i;
) ^' }6 T7 }# D" H- |) b1 K6 e$ _( Z, i    // 若有左孩子(下标小于sz)而且左孩子数值大于当前节点数值,就记录下左孩子下标2 s! _2 q0 A6 ?& w5 ?( p- B% \2 w
    // 建小顶堆需将 '>' 改为 '<', Z% X) G5 a6 \0 H8 n4 {
    if (lpos < sz && nums[lpos] > nums)
6 w; V! Y% v* T        largest = lpos;/ j% s7 I* D7 V/ X
    // 若有右孩子(下标小于sz)而且右孩子数值大于最大的节点值,就记录下右孩子下标
) @/ y/ `  H5 W/ \( c3 z- F* z    // 建小顶堆需将 '>' 改为 '<'" z  s' b- |: l8 u' k8 u
    if (rpos < sz && nums[rpos] > nums[largest])
2 O' ]6 F) b' g: f+ x1 f- q4 E; u        largest = rpos;' h5 i9 N+ c; x9 \2 Y. k
    // 若存在某个孩子的值比当前值大
/ D" C# w: K5 n. D- X  g    if (largest != i)9 @" V; _/ x" q8 }
    {2 ?" X. \  x( o% T# r) z/ Z. X
        // 交换当前节点和较大的那个孩子的值
' J/ T' F, M+ r3 X8 A: H        swap(nums[largest], nums);% x* u1 b. w& M* ]  [& y
        // 交换完之后继续往上调整树
7 g; q% q, J8 Q' b2 J        adjust_from_buttom(nums, largest, sz);
6 K) X" O5 L9 Z; G* L    }; }: o7 ]/ U4 C$ ?0 N
}
$ f- ~2 t& ^8 |
$ r3 C( T) R( l5 hvoid build_heap_from_buttom(vector<int> &nums)
( m- x# X+ _3 r$ M$ \: V{! i; p* c" }9 ~( R2 F) |2 k
    // 若数组只有一个元素则无需建立堆4 d# d( e4 K% J) z
    if (nums.size() != 1)+ G- x5 z5 F7 a7 d2 S  y' F
        // 从最后一个父节点开始循环,一直到根节点% d7 s$ l. z2 Q$ O+ j7 W
        for (int i = nums.size()/2 -1; i >= 0; --i)* C9 h5 i7 \( g# Z# a( {: \9 p
            // 不断对堆进行调整  ^- f) N, S4 W- S: r) K
            adjust(nums, i, nums.size());
0 X2 ]. H. b3 x9 E2 Y3 c}
) D; r8 u- W  b: f" P7 B% X% e! l& O
下面是从上往下建大/小顶堆的代码:/ t8 S1 b2 m9 v& U
. b) v% y& x+ E# P: u0 P( @; O
void adjust_from_top(vector<int> &nums, int i, int sz); D  o( C5 ]- \* f
{
& @5 `. T1 e) n5 U. n    // 如果i不等于08 O$ P5 Y: R8 G( c8 z1 P6 G) n
    if (i)+ I2 Z' S/ v* C" d8 k; S
    {+ V( }' x4 Z7 }7 d$ ~
        // 若i为偶数,其父节点下标为i/2-1。若满足当前节点的值比父节点的值大0 l% ?1 `: [- @9 h+ u2 z
        // 建小顶堆需将 '>' 改为 '<'
' f5 G8 g3 O1 {0 m* V% x9 l7 @        if (!(i%2) && nums > nums[i/2-1])3 h# ~4 ]7 \4 M, L9 N* q+ Q
        {   / O% e2 H  ]% O; J- X
            // 交换当前节点与其父节点,再检查父节点需不需要调整
0 f) t1 f8 u* a) I# ~, y% ]% J$ h            swap(nums, nums[i/2-1]);% L; h% @% E! i; `8 G2 S* G# Y
            adjust_from_top(nums, i/2-1, sz);! U1 z6 g2 j. n; }: {6 F8 d
        }  w9 ]) O6 y; r
        // 若i为奇数,其父节点下标为i/2,则当前节点比父节点的值大! w$ ^  C; ~) N9 M
        // 建小顶堆需将 '>' 改为 '<'
7 H" Z: M- U1 X7 M0 j' Q0 w        if (i%2 && nums > nums[i/2])
7 A* K2 d' E" I1 c$ n5 ]9 G        {: T' {- P1 O9 e- [$ f. D0 {
            // 交换当前节点与其父节点,再检查父节点需不需要调整
( Q: x' c! W7 }4 v7 x1 y2 J( s            swap(nums, nums[i/2]);! c& v/ p8 L2 f, n! a8 w2 o
            adjust_from_top(nums, i/2, sz);
) O6 C. a. _3 d' F9 U( {% b        }6 e( V; |" z; @2 W! L7 q3 K7 S/ g
    }5 X7 W- A1 _% @0 o$ d
}
: {0 y  _. v  c) C# C! U7 Q
/ ?' t* K8 q  Q2 u* uvoid build_heap_from_top(vector<int> &nums)9 a( D; r1 ~4 b' M! r- K
{   // 从第一个元素开始到最后一个元素结束
3 P, K' L5 z; U! a, F6 j0 s: Z    for (int i = 0; i < nums.size(); ++i)
3 F7 j4 M" i2 U  P+ y6 N4 V5 I        // 插入新的元素并调整
3 w% ~5 l/ {  Z( a; g4 }6 H  X8 R        adjust_from_top(nums, i, i+1);
0 m" e! H( u6 R}
( X  Z% e$ r0 o1 B7 G8 X
& J5 k6 F8 H8 e; n1 c建立好的堆之后就可以通过取出堆的根节点,再不断调整堆来进行堆排序,我们可以把取出的节点与堆的最后一个节点进行交换,然后把堆的大小减一,这样就可以利用堆原有的空间来存储排序的结果,而无需额外开辟空间来存储。堆排序代码如下:+ w' d& w2 [: n1 ^" e3 ^
5 `$ V0 I" l. t! p: l( J
void heap_sort(vector<int> &nums)5 W$ h0 _' o! P7 B# x3 ~
{3 O4 d- N; Z; S; l  R) q$ }
    build_heap_from_buttom(nums); // 或者是build_heap_from_buttom(nums);4 S" w+ J; Q- M
    // 从最后一个节点开始往上走+ E$ Q; ~9 ]9 p2 a: U( \$ v
    for(int i = nums.size()-1; i > 0; --i)
- z' S0 B5 s/ Q5 J/ P  ^    {) s+ ^- ^9 R7 P0 l
        // 将最后一个节点与堆的根节点交换4 G4 I0 _( n: R  l( z3 F- B
        swap(nums, nums[0]);' \, l0 s# b) ?3 j' M1 M( S
        // 调整堆,堆的大小为原大小减一
, a; M# w0 Z( y        adjust_from_buttom(nums, 0, i);# ]5 f* ]7 u' |8 V  @! [, Y/ V, s6 L
    }$ S, I2 o6 S' R1 t* Q$ Z; @
}
1 _4 g2 W# {: k, O$ s我们可以注意到堆排序的过程交换了节点之后,只有根节点发生了改变,因此只有可能是根节点不满足堆的性质,所以在调整的时候我们调用了adjust_from_buttom()方法来堆根节点进行调整。6 Q  I4 t$ Z7 U2 Q9 ?$ u" B# o- b+ n

8 o  S& j0 \% }) h! I+ }% K0 {/ @分析:堆排序的时间复杂度包括建堆和调整,从下往上建堆时间复杂度为O(n),调整堆的时间复杂度为(logn),因此堆排序最好、最坏和平均情况的时间复杂度为O(nlogn)。堆排序不需要额外的空间因此其空间复杂度为O(1),但是由于建堆的过程中可能出现位置调整,因此堆排序是一种不稳定的排序方法。3 e+ j0 K5 U3 A3 E8 c+ J. P

. V/ M/ {9 J3 [5 c2 _八、计数排序
4 ~1 N& p* x8 D- y& P* m- _思路:计数排序主要是针对一定范围内的整数来进行了,先确定序列中的最大值和最小值,从而开辟max-min+1大小的空间,空间的每一个位置代表一个数,若遍历到了某个数就将该空间的计数值自增,那么遍历完原数组。排序结果可以通过遍历这个空间来得到。在实现上C++中可以通过map来实现这个操作,因为map本身是保证有序的。当然因为map内部要进去排序,所以如果要实现极致的空间换时间的话,需要自己手动来实现这样的map。1 y* {$ B% W$ X; e+ r- h; B

9 L! s0 A8 O) v' {1 Jvoid count_sort(vector<int> &nums)
9 z- r1 y5 |' g# ^- @! Z0 z$ p{   
. _$ Y3 v! R6 Z! i    // 用于存储数与出现次数的对应关系% v3 p+ }. i! [0 v4 M' O5 T& h& K) u4 ?
    map<int, int> cnt_map;
1 P2 I, P9 ~+ C6 ]4 G    // 遍历数组以将元素及出现次数添加进map2 i* [4 G" E, ^! ~! M1 x+ U( }
    for (auto i : nums)9 C' D, n  c. b& y, G4 G
    {, I- X' s! {; |  i: A: [: ?
        //
4 o! B0 Q- v  G  H2 k/ W5 ?        if (cnt_map.find(i) == cnt_map.end())' O/ ?: R. g* l* e. t4 O0 o
            cnt_map = 1;6 \  E: @' `$ ~( m
        else& l; A+ {& \  e  R& U! G' `4 K
            ++cnt_map;
  C! A+ y4 Q9 i- D$ L    }( [1 ?5 _3 v) W" ]1 `, G! ~* T
    // idx代表原数组下标
5 t7 x0 q3 N. h( _* A: m/ p    int idx = 0;8 H6 u" `5 R" G& D
    // 遍历map以将map中的元素放回原数组* R' v* ]& {' O; i
    for (auto it = cnt_map.begin(); it != cnt_map.end(); ++it)
" g" J% t$ v% \( {  q4 t    {7 `' U6 _- I" \9 `
        while (it->second != 0)- s8 d8 D& u( V1 @9 c
        {
5 j2 M) h; q9 R% j            nums[idx++] = it->first;8 g% c: E8 y+ D: C
            --it->second;
% l$ {9 E; U4 d. M  A; C* Q        }
0 b* }# o$ a' z    }2 Y) Y* j: k2 T& s5 {
}! ?) T9 o. D, o! M
) i# O! ?0 f5 u$ Y
分析:计数排序最好、最坏和平均情况下的时间复杂度都是O(n+k),n是表示遍历原数组的时间,k(数据范围)是表示遍历辅助空间的时间,其空间复杂度为O(k)。计数排序可以看做是一种稳定的排序,比如我们把计数排序存储每个元素的位置看成队列,那么就符合先进先出,所以是稳定的。; P4 ~) b. ?- W$ k* E
3 h- ]# {2 q' v+ y( W2 |6 l
九、基数排序 ( V3 R  ?. |1 u9 E2 v
基数排序通常可以用在整数或字符串上,比如整数上就可以分0-9十个桶,通过比较个位,十位一直到最高位,分别把数装进相应的桶中,再依次取出,比较更高位再装进桶中,再取出,直到所有的数都在0这个桶中说明所有的数都排好序了。  P4 O; E9 }. d1 k5 r6 \! ]
4 z) L" d- C5 j% T. ]
但是需要注意的是,一般的基数排序无法处理有负数的情况。针对负数,一般有两种解决方法:
4 g& x3 Y, @7 a: z. D& |2 `# Q: d* K3 S: D0 ^$ K
1、另外开辟10个桶用来存储负数的情况(需要占用更多的空间)' I/ \) Y8 u6 Z( j

: |) {( M% K' Y' a+ F1 a2、将所有元素变为正数再处理(可以将所有元素都加上最小负数的绝对值,但存在溢出风险),下面代码实现这种方式
$ u& Z" f2 P, `9 Z) y" R
; }2 }  j. x' A4 \3 S) ~8 h! T2 Gvoid radix_sort(vector<int> &nums)1 J  {* c0 a; w" @
{5 j- {. y0 w! c" z4 q# B8 k
    // 开辟二维数组,第一维10代表十个桶,第二维当前为空
: q9 ^1 `4 t& ?: y" _0 X6 M    vector<vector<int>> ivec(10, vector<int>(0));, H4 \, h3 e7 m! Q; O
    // 用于取出最小元素,将所有元素都变为正数- y( I. H. u  C2 f5 }
    int min_num = *min_element(nums.begin(), nums.end());
) V; c/ u( ~7 h! V    if (min_num < 0)
5 i% M8 h8 V5 [, P+ G& u9 E        for (auto &i : nums), e( ?: Q3 K% F+ L$ X  O1 C* i
            i -= min_num;6 h" {2 y3 O' ^" f2 W& Z) a
    // 用来控制比较的位数,初始为1代表比较个位数
% C: d4 N( |4 K1 ~7 R. D, h: A    int digit = 1;
: j) F( P: h( z, g7 S0 s    // 循环一直进行: R1 v' y& i( y7 K- a. w  m+ r
    while (true)7 Q( ]/ r, I# h0 _5 H. `: s4 X' ~
    {$ d1 y) H# L$ [5 _/ x* R: n. D. e$ D
        // 对nums中所有元素遍历,将每个元素放入对应桶中
8 H4 t3 ]8 ^8 l) v2 b  V        for(auto i : nums)
( h, J- ~  w( U            ivec[i%(digit*10)/digit].push_back(i);
( Y9 Q" b/ \( L) J& P: T+ i5 V5 L# n* }$ j
        // 这是循环终止条件,即所有元素都在0号桶中
3 Z  Y* z5 ?; e# G        if (ivec[0].size() == nums.size())
; _( \/ p, ^# f% N& C/ G* ?: g            break;
- z- V. i  l7 X3 o, ^        // 记录nums中的下标用于放回nums中
5 |% _# ?! D: n& e% }        int idx = 0;
; H+ b5 T$ Z5 z% Q3 X        // 依次从每个桶中拿出元素放回nums,取完之后清空桶
# R" ?2 P4 A' L) {/ V2 f        for (int i = 0; i < 10; ++i)- c6 l- c2 M  R+ |6 f
        {
# W+ V- C5 ~& M0 m% u" G            for(int j = 0; j < ivec.size(); ++j)
  n. g- V# P' `                nums[idx++] = ivec[j];
8 l: z: E9 t8 S0 R, _            ivec.clear();
$ U, A+ g0 O; q9 i& m  n        }
9 D6 o2 ]5 l8 `7 ~7 \, q1 [        // 下一次循环比较更高位
8 T  j7 V: h5 f6 x: {        digit *= 10;
" X/ x2 F9 t; o( V0 d    }
5 H, X7 w9 g2 A) [; t+ N! @    // 排序完毕后将元素恢复为初始值! ?* `, |3 g% L7 S! g
    if (min_num < 0)
3 k$ i- s- F2 C' G        for (auto &i : nums)  M( {" R/ ]4 t6 v% S6 `/ p
            i += min_num;$ y9 Z9 t) x% |9 o* {
}
' W0 Z8 i7 H4 p7 Y4 w7 q8 {
; t/ ]! `4 \8 v4 O分析:基数排序的时间复杂度取决于待排序的元素中的最大位数k,如最大元素为1023,那么最大位数k=4。对于每一位我们需要遍历两遍数组,但是这里的两遍可以省略,因此最好、最坏和平均的时间复杂度都是O(k*n)。由于我们需要开辟一个额外的空间来把每个数放入桶中,因此空间复杂度为O(n)。又因为我们放入桶和从桶中拿出的操作都保证了先进先出,因此基数排序是一种稳定的排序方法。* t3 U9 ?9 }' [* q

; m1 M; i% o( P7 d十、桶排序 & J5 X9 T* \, P! I( n% ?, M4 E
前面计数排序和基数排序都是桶排序的一种特例。桶排序是这样一种排序,它将原来的数组分到不同的桶里,再分别对每个桶内部采用排序算法进行排序(可以用快排、可以用选择、冒泡等都没问题),最后将每个桶里的元素依次取出就得到了排好序的结果。因此桶排序如何分桶是一个关键的问题,如果分的好的话,每个桶都能均匀地承载一部分数据,最后时间复杂度就会比较理想。若果分的不好,很多数据都集中在一个桶里,那么最后的时间复杂度就很差(最差情况下会略大于桶内排序算法的效率),此外桶排序还需要额外的空间。
/ L0 i- A% m# M! A# X* C1 L
: ~' z# A8 r8 K; k4 }# Mvoid bucket_sort(vector<int> &nums, int bucket_cnt)
' l+ ?0 d! _* e0 k% F{* }* g7 F7 W0 U6 C( @3 G
    // 建立bucket_cnt个桶,每个桶初始为空/ A# v4 M# O- Y3 O. @
    vector<vector<int>> ivec(bucket_cnt, vector<int>(0));
& E; K: a" J9 P1 {1 q    // 拿出数组中最大和最小的元素来确定每个桶之间的间隔
9 c" n% P5 _) E" s! V: S    int min_ele = *min_element(nums.begin(), nums.end());
0 Z( d# L  b* G& T# H" n; e* O! G9 @' s    int max_ele = *max_element(nums.begin(), nums.end());
0 N1 W$ F2 A) l' r4 b% N# o    int gap = (max_ele - min_ele) / bucket_cnt;* m, F) }, ~! G. q+ [4 `) R
    // 如果间隔比桶个数都小,那么直接排序可能更快- t4 X9 S3 l% x/ Z
    if (gap > bucket_cnt)7 F! ^9 o# o) p! n; |
    {+ x4 C8 S8 [# o9 q
        // 遍历数组中每个元素
. f, q2 B. ^; R! S3 Y" t        for (auto i : nums)
0 e/ K3 `; ]* A        {4 W, E' P) j; j; o; I3 i1 y- U
            // 确定元素所在的桶序号
# C7 V& U6 ^; i2 e. U1 j            int bucket_num = (i-min_ele)/gap;9 u4 i" s! M; M, l
            // 若算出的序号等于桶个数,那么说明越界了(比如数组中最大元素就会出现这样的情况)6 r( R3 ^4 G2 F! w
            if (bucket_num == bucket_cnt)
% u* o; V& `5 n+ Z+ d                --bucket_num;
: F( T' O0 C7 i6 S6 h# q            // 将该元素装入桶中
- }' q6 R8 Z" u# c            ivec[bucket_num].push_back(i);
7 c' g- s6 V# h9 X9 Q: {        }
; U' d: }+ L9 Z3 @- p# J        // idx记录原数组nums中的下标用于装回元素
5 F5 Y- a+ z  l3 r        int idx = 0;
8 M2 g- o3 ]0 h# r- m" n$ H% I        // 对每个桶先排序再往nums中装入元素,最后清空桶. i  b/ Z- F5 M- M8 a8 U$ _8 L
        for (auto &vec : ivec)
, s; b, s8 }, {6 J. I8 q9 v        {0 a8 F1 C! X3 w* N, k6 X% N
            quick_sort(vec);$ B5 k" D# a& X. S5 ?" `1 N) W
            for (const auto ele : vec)
7 ~3 U& q% ?, l  t! [                nums[idx++] = ele;% X- @7 K, X- F) U3 w0 Y' G6 c
            vec.clear();
( \+ i) W& s. s' S        }
: ?8 |* S! h& n- x# N/ _. o    } 0 M" ^; ^* q0 e
    // 直接排序
* }2 s) [% z9 H+ ~  T, O0 J6 Z    else- C; n+ `7 l2 I# Q/ G" @# w
        quick_sort(nums);. ?8 L! t1 a8 R( J# x& k
}! J) c* \. b5 z! v/ E, V( U" r% b7 {
% i& I, v* l3 |/ B4 P* l( ~
分析:桶排序将原数组分为k个桶来处理,假设数组总数为n,那么每个桶装了n/k个数据,假设我们在桶的内部使用快速排序,其时间复杂度为O(klogk),那么对于所有桶来说,总的时间复杂度就是k*O((n/k)*log(n/k)) = O(n*log(n/k))。最好的情况是每个数据都装进一个桶,这样就无需排序了,因此最好的时间复杂度是O(n),平均情况的时间复杂度是O(n*log(n/k)),最坏的情况的可以接近O(nlogn),但是比O(nlogn)要更差(因为存在分桶等操作)。桶排序的空间复杂度是O(n+k),因为桶排序除了要装下原始数组长度的元素,还需要额外的k个桶,这k个桶系统会默认分配内存。桶排序是否稳定依赖于其内部的排序算法,若内部使用快排则不稳定,若内部使用归并则是稳定的。0 M& w5 i5 I6 Y& E# t- h  W
————————————————
  ^/ V" E$ y& e  x版权声明:本文为CSDN博主「Xaiver_97」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。: p" D( B5 L" y7 R" C, o% u, Y+ X
原文链接:https://blog.csdn.net/Xavier_97/article/details/126722423
" {4 x8 G9 I& d& s7 p* H4 ]0 ~: v% Q  b
7 A0 k, x5 w5 F% M+ I5 d





欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5