如何渐近分析冒泡排序的轮数的期望值？

对 ${\sum }_{k=1}^n\frac{k!k^{n-k}}{n!}$ 的渐进分析

已经有了很好的回答,在这里我只想谈谈

$n+1-{\sum }_{k=1}^n\frac{k!k^{n-k}}{n!}$ 是怎么来的

根据计算机程序设计艺术卷3:排序与查找(第二版) Donald E. Knuth著 贾洪峰译 人民邮电出版社 第84页的解说，我们实际上只需要理解为什么当 $1\le k\le n$ 时所有分量均小于 $k$ 的反序表数目是

$k!k^{n-k}$ 其实非常简单，由上面提到的著作第8页给出的不等式 $0\le b_1\le n-1$

$0\le b_2\le n-2-----0\le b_{n-1}\le 1b_n=0$

$b_1$ 有 $0---k-1$ 这 $k$ 种选择, $b_2$ 同样有 $0---k-1这k$ 种选择,-------

$b_{n-k}$ 同样有 $0---k-1$ 这 $k$ 种选择,

再往下 $b_{n-k+1}$ 有 $0---k-1$ 这 $k$ 种选择, $b_{n-k+2}$ 有 $0---k-2$ 这 $k-1$ 种选择

$b_{n-k+3}$ 有 $0---k-3$ 这 $k-2$ 种选择—以此类推 $b_{n-1}$ 有2种选择, $b_n$ 有1种选择,从而由乘法原理总共有 $k!k^{n-k}$ 种选择，故当 $1\le k\le n$ 时所有分量均小于 k 的反序表数目是

$k!k^{n-k}$有了这一结论其他也就不难理解了

另一种证明思路

设待排序序列为 $A,A$ 的长度 $n\ge 2$ , $A$ 中元素按从小到大的顺序排序为 $S=[S1,S2,---,Sn]S1<S2<---<Sn$

现考虑按从小到大的顺序排序元素的冒泡排序算法。我们用 $1,2,3,---,n$ 分别索引被排序元素所处的 $n$ 个槽位，初始时， $A$ 的 $n$ 个元素分别位于槽位 $1,2,3,---,n$ 上，记 $P_k(i,j)$ 为第 k 轮冒泡排序后 $S_i$ 位于槽位 $j$ 上的概率

由于排序的初始序列 $A$ 可以被认为是随机给定的,所以 $Si1<=i<=n$ 恰位于槽位 $j(1<=j<=n)$上的概率 $P_0(i,j)$ 为 $\frac{1}{n}$ (1)

考虑冒泡排序第 k 轮排序前( $1<=k<=n-1$ )的状态，此时槽位 $n-k+2$ 到 $n$ 上的元素从小到大排序，恰为 $A$ 的前 $k-1$ 大的数，而槽位 $1$ 到 $n-k+1$ 上的元素则可能处于无序状态，由数学归纳法和(1)，现在在 $P_{k-1}(i,j)(1\le i\le n-k+11\le j\le n-k+1)$ 已经求出的情况下计算 $P_k(i,j)(1\le i\le n-k1\le j\le n-k)$

事实上

$P_k(i,j)=(1-{\sum }_{s=j+2}^{n-k+1}{P}_{k-1}(i,s))[P_{k-1}(i,j)\ast {\prod }_{m=1}^{j-1}{\sum }_{q=1}^{i-1}{P}_{k-1}(q,m)\ast {\sum }_{s=i+1}^{n-k+1}{P}_{k-1}(s,j+1)+({\sum }_{s=1}^{j-1}(P_{k-1}(i,s)\ast {\prod }_{m=1,m\ne s}^j{\sum }_{q=1}^{i-1}{P}_{k-1}(q,m)))\ast {\sum }_{s=i+1}^{n-k+1}{P}_{k-1}(s,j+1)+P_{k-1}(i,j+1)(1-{\prod }_{m=1}^j{\sum }_{q=1}^{i-1}{P}_{k-1}(q,m))]$

(2)

注意, $j+2>n-k+1$ 时 ${\sum }_{s=j+2}^{n-k+1}{P}_{k-1}(i,s)=0$

$j=1$ 时 ${\prod }_{m=1}^{j-1}{\sum }_{q=1}^{i-1}{P}_{k-1}(q,m)=1$

$j=1$ 时 $({\sum }_{s=1}^{j-1}(P_{k-1}(i,s)\ast {\prod }_{m=1,m\ne s}^j{\sum }_{q=1}^{i-1}{P}_{k-1}(q,m)))=0$

$i=1$ 时${\sum }_{q=1}^{i-1}{P}_{k-1}(q,m)=0$

于是第 k 轮排序后$Si$ $1<=i<=n-k$落到槽位 j ( $1<=j<=n-k$ )上的概率为 $P_k(i,j)$

现考虑第k轮排序,记事件 $B_i^k$ 表示第k轮排序完成后 $Si$ 恰位于槽位 i 上( $1<=i<=n-k$ )

记事件 $C^k$ 表示第k轮排序完成后 $1,2,---,n-k$ 槽位上的元素完全有序(从小到大)

则 $P(C^k)=P(B_1^k{B}_2^k---B_{n-k}^k)$

$=P(B_1^k)P(B_2^k|B_1^k)---P(B_{n-k}^k|B_1^k---B_{n-k-1}^k)$

注意到 $P(B_{n-k}^k|B_1^k---B_{n-k-1}^k)=1$ 故 $P(C^k)=P(B_1^k)P(B_2^k|B_1^k)---P(B_{n-k-1}^k|B_1^k---B_{n-k-2}^k)$

现在我们来计算 $P(B_m^k|B_1^k---B_{m-1}^k)2<=m<=n-k-1$

如果事件 $B_1^k---B_{m-1}^k$ 发生,即第 k 轮排序结束后 $S1---Sm-1$ 分别位于槽位 $1,2,---,m-1$ 上 对冒泡排序而言,在第 k 轮排序中,当槽位 $m-1$ 和 $m$ 上的元素比较交换完成后,前 $m-1$ 个槽位上的元素就完全确定了,和第k轮排序完成后的结果一致，恰为前 $m-1$ 小的元素(由小到大排列),故此时要使 $B_m^k$ 发生,则 $Sm$ 必须位于槽位 $m$ 或$m+1$ 上

当 $Sm$ 位于槽位 $m$ 上时, $m+1$ 上的元素 $S_{m+1}^{{}^{\prime }}$ 必须大于 $Sm$ (3)

当 $Sm$ 位于 $m+1$ 上时, $m$ 上的元素 $S_m^{{}^{\prime }}$ 必须大于 $Sm$ (4)(以上论断的正确性请读者自己思考)

不难看出(3)自动成立，这是因为此时槽位 $1,---,m-1$ 由前 $m-1$ 小的元素占据，从而槽位 $m,m+1,---,n-k$上排列着第 $m,m+1,---,n-k$ 小的元素，故此时当 $Sm$ 位于槽位 $m$ 上时,槽位 $m+1$ 上的元素 $S_{m+1}^{{}^{\prime }}$ 只能是第 $m+1$ 小或第 $m+2$ 小—或第 $n-k$ 小元素,从而必有 $S_{m+1}^{{}^{\prime }}>Sm$ (5)

类似的不难看出(4)自动成立 (6)

另外此时由(2) $Sm$ 位于槽位 $m+1$ 位置的概率为 $P_{k-1}(m,m+1)$ , 而$Sm$ 位于 $m$ 位置的概率为 ${\sum }_{q=1}^m{P}_{k-1}(m,q)$

由(5)当 $Sm$ 位于槽位 $m$ 上时, $m+1$ 上的元素 $S_{m+1}^{{}^{\prime }}$ 大于 $Sm$

的概率为 1

而由(6)当 $Sm$ 位于 $m+1$ 上时,槽位 $m$ 上的 $S_m^{{}^{\prime }}>Sm$ 的概率为 1

故 $P(B_m^k|B_1^k---B_{m-1}^k)={\sum }_{q=1}^m{P}_{k-1}(m,q)\ast 1+P_{k-1}(m,m+1)\ast 1$ (7)

注意当 $m=1$ 时(7)依然成立

这样我们有 $P(C^k)={\prod }_{m=1}^{n-k-1}({\sum }_{q=1}^m{P}_{k-1}(m,q)+P_{k-1}(m,m+1))$ (8)

$P(排序轮数恰为i轮)={\prod }_{k=1}^{i-1}(1-P(C^k))P(C^i)$ (9)

冒泡排序排序轮数的期望次数 $E={\sum }_{i=1}^{n-1}iP(排序轮数恰为i轮)$ (10)

注意,(9)式中 $i=1$ 时 ${\prod }_{k=1}^{i-1}(1-P(C^k))=1$

(8)式中 $k=n-1$ 时(8)式即为 1

联立(8)(9)(10)式可得冒泡排序轮数期望值的展开式，但该展开式的渐进分析比较繁琐困难，也不是本文的重点，有兴趣的读者可以尝试以下分析该展开式的渐进表达式