概率论处理随机现象。给每个合理结果配上一个数,用来反映它发生的可能性,这套框架就是概率测度。从赌博到物理、生物、经济,定量地说“有多可能”,都靠它。要度量一个事件,先得知道实验有多少种结果,其中又有多少种落在这个事件里。计数就是这块有限基础。

计数原理

基本计数原理

先看最简单的情形。一枚均匀硬币,正面反面各占一半,一共两种结果。掷两次,用 TT 记正面、FF 记反面,四种排列是 (T,T)(T,T)、(T,F)(T,F)、(F,T)(F,T)、(F,F)(F,F)。这引出乘法原理。

定理乘法原理

两个实验相继进行。第一个有 mm 种结果,固定第一种结果后,第二个有 nn 种结果。则两个实验合在一起有 mnmn 种结果。

证明

按第一步结果把完整结果分组,得到 mm 个互不相交的组,每组恰有 nn 个结果。因此有限加法给出 n+⋯+n=mnn+\cdots+n=mn。各组第二步的结果不必具有相同标签,只需个数相同。这个计数命题不要求概率意义上的独立。

例年度母子

一个小社区有 1010 位母亲,每人有 33 个孩子。要选一对“年度母子”,有多少种选法?

解年度母子

选母亲有 m=10m=10 种,再选她的一个孩子有 n=3n=3 种,共 10×3=3010\times 3=30 种。

掷硬币 nn 次,和 nn 个布尔变量的真值表是同一件事:每个位置两种取值,总共 2n2^n 种。把两个实验推广到 rr 个,就是一般的乘法原理。

定理一般乘法原理

rr 个实验相继进行。第一个有 n1n_1 种结果;对已经出现的前一种结果,第二个有 n2n_2 种;对已经出现的前两种结果,第三个有 n3n_3 种;依此类推。则总共有

∏k=1rnk=n1n2⋯nr\prod_{k=1}^r n_k = n_1 n_2 \cdots n_r

种结果。

证明

对阶段数 rr 归纳。r=1r=1 时为 n1n_1。若前 rr 步有 n1⋯nrn_1\cdots n_r 个完整前缀,每个前缀又有 nr+1n_{r+1} 种延续,由两阶段乘法原理得到 n1⋯nrnr+1n_1\cdots n_r n_{r+1},归纳完成。空序列有一种实现,对应空积为一。

这就是计数里的乘法法则,后面算概率还会再用。

定理加法原理

一件事可以走 mm 条路完成,也可以走另外 nn 条路完成,两类路互不重叠,则总共 m+nm+n 种做法。

证明

设 AA、BB 有限且 A∩B=∅A\cap B=\emptyset。AA 有 nn 个元素,BB 有 mm 个。并集把两边的元素都收进来,没有重复,所以 ∣A∪B∣=n+m|A\cup B|=n+m。

例三份项目清单

学生从三份互不重复的项目清单里任选一个,三份分别有 2323、1515、1919 个项目。一共有多少种选法?

解三份项目清单

三类做法互斥,由加法原理:23+15+19=5723+15+19=57。

例车牌

车牌由三个大写英文字母后接三个数字组成,字母组合不加限制。能做出多少种不同车牌?

解车牌

每个字母位置 2626 种,每个数字位置 1010 种。由乘法原理

263×103=17576×1000=17,576,000.26^3\times 10^3=17576\times 1000=17{,}576{,}000.

两类做法若有重叠,直接相加会把公共部分数两遍,必须减回来。这就是减法,也是两个集合的容斥。

全集 UU 里丢掉子集 AA,剩下 ∣U∣−∣A∣|U|-|A| 个元素。两个集合时,

∣A∪B∣=∣A∣+∣B∣−∣A∩B∣.|A\cup B|=|A|+|B|-|A\cap B|.
定理容斥(两个做法)

一件事有 n1n_1 种做法,也有 n2n_2 种做法,则真正的种数是 n1+n2n_1+n_2 减去两种做法里重复的那些。

证明

把并集拆成互不相交的 A∖BA\setminus B、A∩BA\cap B、B∖AB\setminus A。在 ∣A∣+∣B∣|A|+|B| 中,中间部分算了两次,其余各一次;减去 ∣A∩B∣|A\cap B| 后,每个元素恰算一次。

例非虚构书

图书馆有 10001000 本书,其中 300300 本是小说。非小说有多少本?

∣U∣−∣A∣=1000−300=700.|U|-|A|=1000-300=700.
例茶与咖啡

“调查 200200 人,喜欢茶的 120120 人,喜欢咖啡的 150150 人,两者都喜欢的 5050 人”这组数据不可能成立:容斥给出 120+150−50=220>200120+150-50=220>200。交集至少应为 120+150−200=70120+150-200=70 人。若把交集改成 8080 人,并集就是 190190 人,另有 1010 人两者都不喜欢。

除法原理则用来抹掉“不在乎的差别”,和等价类是同一件事。

定理除法原理

一件事有 nn 种做法,而且每种真正在乎的结果 ww 恰好对应其中 dd 种做法,则不同的结果有 n/dn/d 种。

证明

按最终结果把过程结果划分。若最终有 kk 种结果,每个原像恰好含 d>0d>0 个过程结果,互不相交的加法给出 n=kdn=kd,故 k=n/dk=n/d。若各原像大小不同,就不能这样除。

例等大小的等价类

nn 元集被等价关系分成若干类,每类恰好 mm 个元素。不区分类内次序时,不同的块有 n/mn/m 个。

例图的边数

图的度数总和是 5858,有多少条边?

解图的边数

每条边被度数数了两次,所以边数是 58/2=2958/2=29。

鸽巢原理

88 只鸽子进 77 个巢,必有一巢至少两只。

定理鸽巢原理

把不少于 k+1k+1 个物体放进 kk 个盒子,至少有一个盒子装了两个或更多。

证明

用反证。若每个盒子最多一个物体,物体总数至多 kk,与“至少 k+1k+1 个”矛盾。

定理一般鸽巢原理

把 NN 个物体放进 kk 个盒子,至少有一个盒子不少于 ⌈N/k⌉\lceil N/k\rceil 个物体。

证明

仍用反证。若每个盒子都少于 ⌈N/k⌉\lceil N/k\rceil 个,则每个至多 ⌈N/k⌉−1\lceil N/k\rceil-1 个。总数至多

k(⌈N/k⌉−1).k\bigl(\lceil N/k\rceil-1\bigr).
注天花板函数

对任意实数 xx 有 ⌈x⌉≤x+1\lceil x\rceil\le x+1。

代入 x=N/kx=N/k 得 ⌈N/k⌉≤N/k+1\lceil N/k\rceil\le N/k+1,从而

k(⌈N/k⌉−1)<k(N/k+1−1)=N.k\bigl(\lceil N/k\rceil-1\bigr) < k\bigl(N/k+1-1\bigr) =N.

物体总数会严格小于 NN,矛盾。

例生日与首字母
  • 367367 人中必有两人生日相同,因为一年最多 366366 天。
  • 2727 个英文单词中必有两个同一字母开头,因为字母只有 2626 个。
例同一分数至少六人

离散数学课有 A,B,C,D,FA,B,C,D,F 五种成绩。要保证至少六人分数相同,最少需要多少学生?

证明

要某个盒子至少 m=6m=6 个,盒子数 k=5k=5。最坏情形是每种成绩恰好 55 人,共 2525 人;再来一人就有一种成绩达到 66。所以最少 2626 人。写成不等式就是 N>k(m−1)=25N>k(m-1)=25。

习题

练习一般乘法原理

证明一般乘法原理。

证明

对实验个数 rr 归纳。r=1r=1 时就是 n1n_1 种。假设 r=kr=k 时有 ∏i=1kni\prod_{i=1}^k n_i 种。再添第 k+1k+1 个实验,对每种已有结果又有 nk+1n_{k+1} 种,于是总数再乘 nk+1n_{k+1},得到 ∏i=1k+1ni\prod_{i=1}^{k+1} n_i。

练习函数的个数

mm 元集到 nn 元集有多少个函数?

解函数的个数

定义域每个元素可独立选陪域里 nn 个象之一,共 nmn^m 个函数。

练习单射的个数

mm 元集到 nn 元集有多少个单射?

解单射的个数

m>nm>n 时没有单射。m≤nm\le n 时,第一个原象有 nn 种象,第二个 n−1n-1 种,第 kk 个 n−k+1n-k+1 种,共 n(n−1)⋯(n−m+1)n(n-1)\cdots(n-m+1) 个。

练习笛卡尔积的基数

证明:有限集 A1,…,AmA_1,\ldots,A_m 的笛卡尔积,其元素个数等于各集基数之积。

证明

构造一个 mm 元组,就是依次从 AiA_i 里各取一个元素。由乘法原理,

∣A1×⋯×Am∣=∣A1∣⋯∣Am∣.|A_1\times\cdots\times A_m|=|A_1|\cdots|A_m|.
练习没有单射

∣A∣=k+1|A|=k+1,∣B∣=k|B|=k。证明不存在单射 A→BA\to B。

解没有单射

k+1k+1 个原象进 kk 个象,鸽巢原理迫使某个象至少被两个原象用到,所以不可能单射。

练习扑克牌

一副 5252 张牌,要保证:

  1. 至少三张同花色,最少抽几张?
  2. 至少三张红心,最少抽几张?
解扑克牌

(a) 四个花色当四个盒子。抽 88 张时可能每种两张;第 99 张迫使某一花色达到三张。即 ⌈N/4⌉≥3\lceil N/4\rceil\ge 3 的最小 NN 是 99。

(b) 最坏情形先把另外三门各 1313 张共 3939 张抽完,再抽三张红心。要保证三张红心,最多需要 4242 张。

练习有限集上的关系计数

X={a,b,c,d}X=\{a,b,c,d\}。

  1. XX 上有多少个关系?
  2. 其中多少个自反?
  3. 其中多少个既自反又对称?
  4. 其中多少个等价关系?
  5. 若 ∣X∣=n|X|=n,(1)(2)(3) 各是多少?
解有限集上的关系计数

(a) X×XX\times X 有 1616 个有序对,每个可在关系里或不在,共 2162^{16} 个关系。

(b) 对角线四个位置必须是 11,其余 1212 个自由,共 2122^{12}。也就是 216/24=2122^{16}/2^4=2^{12}。矩阵形如

R=[1abcd1efgh1ijkl1],R=\begin{bmatrix} 1 & a & b & c \\ d & 1 & e & f \\ g & h & 1 & i \\ j & k & l & 1 \end{bmatrix},
注矩阵里的字母

这里的 a,…,la,\ldots,l 表示 00 或 11,不是 XX 的元素。

(c) 自反且对称:对角线已定,上三角(不含对角)有 66 个独立选择,共 262^6。

(d) 等价关系与划分一一对应。44 元集的划分数是 Bell 数 B4=15B_4=15。

(e) 一般 nn 元集:关系 2n22^{n^2} 个,自反 2n(n−1)2^{n(n-1)} 个,自反且对称 2n(n−1)/22^{n(n-1)/2} 个。

练习互素的一对

从 11 到 2n2n 中任取 n+1n+1 个互不相同的整数,证明其中必有两个互素。这个结论对只取 nn 个数还成立吗?

解互素的一对

把 {1,2},{3,4},…,{2n−1,2n}\{1,2\},\{3,4\},\ldots,\{2n-1,2n\} 看成 nn 个盒子。取 n+1n+1 个数,必有两个落在同一盒子,它们是相邻整数,最大公因数是 11。

取 nn 个数时结论不必成立:把 2,4,…,2n2,4,\ldots,2n 全取来,任意两个的最大公因数至少是 22。

排列与组合

上一节的四条原理已经能算不少问题。排列、组合把“有没有次序、能不能重复”这两条再抽象一层。

计数:顺序与放回。先判断顺序是否重要、是否允许重复,再选择对应的计数规则。

计数:顺序与放回。先判断顺序是否重要、是否允许重复,再选择对应的计数规则。

先判断顺序是否重要、是否允许重复,再选择对应的计数规则。

排列

字符串 abcabc 有 66 种排法:abc,acb,bac,bca,cab,cbaabc,acb,bac,bca,cab,cba。每一种都叫这些字母的一个排列。nn 个彼此可区分的对象,第一个位置 nn 种选法,第二个 n−1n-1 种,直到最后一个,总共

n(n−1)⋯2⋅1=n!.n(n-1)\cdots 2\cdot 1=n!.
定义全排列

nn 个对象的排列数为 n!n!。

记号阶乘

非负整数 nn 的阶乘是

n!=n×(n−1)×⋯×2×1.n!=n\times(n-1)\times\cdots\times 2\times 1.

约定 0!=10!=1。后面有一道习题专门问为什么不能是 00。

例成绩排名

概率课有 66 名男生、44 名女生,考试分数两两不同。一共有多少种排名?若男生女生分开排名呢?

解成绩排名

全体一起排是 (6+4)!=3,628,800(6+4)! = 3{,}628{,}800。男女分开排,是 6! 4!=720×24=17,2806!\,4!=720\times 24=17{,}280。

注阶乘的不等式

这个例子也说明:对正整数 aa、bb,总有 (a+b)!≥a! b!(a+b)!\ge a!\,b!。

例PEPPER

用字母 PEPPER\mathrm{PEPPER} 能排出多少个不同的字符串?

解PEPPER

若三个 P\mathrm{P}、两个 E\mathrm{E} 都加上标变成可区分的,有 6!6! 种排法。任意一种,例如 P1P2E1P3E2RP_1P_2E_1P_3E_2R,把三个 P\mathrm{P} 内部重排(3!3! 种)、两个 E\mathrm{E} 内部重排(2!2! 种),看起来仍是 PPEPER\mathrm{PPEPER}。由除法原理,不同字符串有

6!3! 2!=60\frac{6!}{3!\,2!}=60

种。

这就是有重复的排列。

定义有重复的排列

nn 个位置里,第 ii 种对象出现 nin_i 次(∑ni=n\sum n_i=n)。不区分同种对象时,排列数为

n!n1! n2!⋯nk!.\frac{n!}{n_1!\,n_2!\cdots n_k!}.

再看从一堆里取出若干个来排。

例可重复的三位数

用 1,2,3,41,2,3,4 做三位数,数字用过仍可再用,能做出多少个?

解可重复的三位数

每个位置 44 种,共 43=644^3=64。

例不可重复的三位数

同样四个数字,用过就不再用,能做出多少个三位数?

解不可重复的三位数

4×3×2=244\times 3\times 2=24。

第二种叫做 rr-排列。

定理rr-排列

nn 元集的 rr-排列个数(1≤r≤n1\le r\le n)是

P(n,r)=n(n−1)⋯(n−r+1).P(n,r)=n(n-1)\cdots(n-r+1).

也常写成 nPr{}^n P_r。

证明

依次填有序位置。已选 jj 个元素后,恰剩 n−jn-j 个未用元素,与具体选了哪段互异前缀无关。把 j=0,…,r−1j=0,\ldots,r-1 的个数相乘,便得乘积。再从 n!n! 中约去 1,…,n−r1,\ldots,n-r,得到阶乘商。r=0r=0 时恰有一个空排列,商式仍给出一。

写成闭形式:因为 P(n,0)=1P(n,0)=1,与 n!/(n−0)!n!/(n-0)! 一致,于是

P(n,r)=n!(n−r)!.P(n,r)=\frac{n!}{(n-r)!}.
推论rr-排列的闭形式

对 0≤r≤n0\le r\le n,

P(n,r)=n!(n−r)!.P(n,r)=\frac{n!}{(n-r)!}.

把开形式 n(n−1)⋯(n−r+1)n(n-1)\cdots(n-r+1) 上下同乘 (n−r)!(n-r)! 就是这个商。

组合

例无序取三个字母

从 a,b,c,d,ea,b,c,d,e 里不放回地取 33 个字母,不讲究次序(abab 与 baba 算同一种)。有多少种?

解无序取三个字母

有次序时是 5×4×3=605\times 4\times 3=60。每种无序组合对应 3!3! 个有序串,除掉之后得

5×4×33!=10.\frac{5\times 4\times 3}{3!}=10.

也就是 P(5,3)/3!P(5,3)/3!。一般地

P(n,r)r!=n!(n−r)! r!.\frac{P(n,r)}{r!}=\frac{n!}{(n-r)!\,r!}.
定理组合数

nn 个不同对象里取 rr 个的组合数记作 C(n,r)C(n,r) 或 (nr)\binom{n}{r},叫做二项式系数:

(nr)=n!(n−r)! r!(r≤n).\binom{n}{r}=\frac{n!}{(n-r)!\,r!}\qquad(r\le n).

读作“nn 选 rr”。

证明

假设 0≤r≤n0\le r\le n 都为整数。每个无序 rr 元子集恰好产生 r!r! 个不同的有序选择,所以“忘掉次序”映射的每个原像大小均为 r!r!。由除法原理,(nr)=P(n,r)/r!\binom nr=P(n,r)/r!。r=0r=0 时使用空选择,结论仍成立。

例委员会

55 名女性和 77 名男性里,选 22 女 33 男组成委员会,有多少种?若其中两名男性不和,拒绝同席呢?

解委员会
(52)(73)=5⋅42⋅1⋅7⋅6⋅53⋅2⋅1=350.\binom{5}{2}\binom{7}{3}=\frac{5\cdot 4}{2\cdot 1}\cdot\frac{7\cdot 6\cdot 5}{3\cdot 2\cdot 1}=350.

含那对不和男性的三人组有 (22)(51)=5\binom{2}{2}\binom{5}{1}=5 种,所以不含他们的三人组有 35−5=3035-5=30 种,再配 22 名女性:

30×(52)=300.30\times\binom{5}{2}=300.
推论对称性

0≤r≤n0\le r\le n 时,C(n,r)=C(n,n−r)C(n,r)=C(n,n-r)。

证明
C(n,r)=n!r! (n−r)!=C(n,n−r).C(n,r)=\frac{n!}{r!\,(n-r)!}=C(n,n-r).
证明

组合证明:选 rr 个带上,等于选 n−rn-r 个留下。

例好香蕉与坏香蕉

店里 2020 根香蕉,其中 1919 根好的、11 根坏的。买 66 根:要么避开那根坏的,要么带着它再配 55 根好的。于是

(206)=(196)+(195)=19!6! 13!+19!5! 14!=27132+11628=38760.\begin{aligned} \binom{20}{6} &= \binom{19}{6}+\binom{19}{5} \\ &= \frac{19!}{6!\,13!}+\frac{19!}{5!\,14!} = 27132+11628 = 38760. \end{aligned}
推论Pascal 恒等式

若 0<k<n0<k<n,则

(nk)=(n−1k−1)+(n−1k).\binom{n}{k}=\binom{n-1}{k-1}+\binom{n-1}{k}.
证明
(n−1k)+(n−1k−1)=(n−1)!k! (n−k−1)!+(n−1)!(k−1)! (n−k)!=(n−k)(n−1)!k! (n−k)!+k(n−1)!k! (n−k)!=n!k! (n−k)!=(nk).\begin{aligned} \binom{n-1}{k}+\binom{n-1}{k-1} &= \frac{(n-1)!}{k!\,(n-k-1)!}+\frac{(n-1)!}{(k-1)!\,(n-k)!} \\ &= \frac{(n-k)(n-1)!}{k!\,(n-k)!}+\frac{k(n-1)!}{k!\,(n-k)!} \\ &= \frac{n!}{k!\,(n-k)!} = \binom{n}{k}. \end{aligned}

这是递归关系。组合解释:从 nn 个里选 kk 个,要么含某个指定元素(再从剩下 n−1n-1 个里选 k−1k-1 个),要么不含(从剩下 n−1n-1 个里选 kk 个)。

注Pascal 恒等式

这个恒等式叫做 Pascal 恒等式。链接里有组合证明。

用集合语言再看计数

四条原理、排列、组合,骨子里都是在操作集合。判断用哪一种,就看两件事:

  • 要不要次序(排列还是组合)
  • 能不能重复

nn 元集上的 kk-序列

定义kk-序列

若序列 SS 的陪域是 CC,称 SS 是 CC 上的序列。正整数 kk、nn 给定后,nn 元集上的 kk-序列是从 {1,…,k}\{1,\ldots,k\} 到某个 nn 元集 X={x1,…,xn}X=\{x_1,\ldots,x_n\} 的函数,写成

S=(s1,s2,…,sk),sj∈X.S=(s_1,s_2,\ldots,s_k),\qquad s_j\in X.

允许重复的 kk-排列,就是 nn 元集上的 kk-序列,共 nkn^k 个。例如 {1,2,3,4}\{1,2,3,4\} 上的 33-序列有 434^3 个。

不允许重复时,就是 P(n,k)P(n,k):

n(n−1)⋯(n−k+1).n(n-1)\cdots(n-k+1).

例如 66 元集上的 44-排列有 6×5×4×3=3606\times 5\times 4\times 3=360 个。44 元集上的 66-排列形式上会出现 00 因子,个数是 00,因为取不够。

nn 元集的子集个数是 2n2^n:每个元素可以在或不在特征序列里,位置上填 00 或 11,共 2n2^n 条特征序列,也就是 2n2^n 个子集。这是乘法原理的推论。

nn 元集的 kk-子集

不允许重复、不讲究次序,就是 (nk)\binom{n}{k}。

允许重复的组合对应多重集。从 nn 种里取 kk 个、同种可再取,种数是 (n+k−1k)\binom{n+k-1}{k}。

定义可重复组合

从 X={x1,…,xn}X=\{x_1,\ldots,x_n\} 中取 kk 个、允许重复,种数为

(n+k−1k)=(n+k−1)!k! (n−1)!.\binom{n+k-1}{k}=\frac{(n+k-1)!}{k!\,(n-1)!}.
证明

用星与杠。kk 颗星表示取到的对象,n−1n-1 根杠把它们分成 nn 组。一共 n+k−1n+k-1 个位置,从中选 kk 个放星(其余放杠),正是 (n+k−1k)\binom{n+k-1}{k}。

例两种水果

水果集 {苹果,香蕉,樱桃}\{\text{苹果},\text{香蕉},\text{樱桃}\},可重复地取 22 个:两苹果、苹果香蕉、苹果樱桃、两香蕉、香蕉樱桃、两樱桃,共

(3+2−12)=(42)=6\binom{3+2-1}{2}=\binom{4}{2}=6

种。

例三色球

红、蓝、绿三色球,可重复取 33 个:RRR,RRB,RRG,RBB,RBG,RGG,BBB,BBG,BGG,GGG\mathrm{RRR},\mathrm{RRB},\mathrm{RRG},\mathrm{RBB},\mathrm{RBG},\mathrm{RGG},\mathrm{BBB},\mathrm{BBG},\mathrm{BGG},\mathrm{GGG},共

(3+3−13)=(53)=10\binom{3+3-1}{3}=\binom{5}{3}=10

种。

对照:

  • nn 元集的 kk-子集:不重复、无序,(nk)\binom{n}{k}。
  • nn 种的 kk-多重子集:可重复、无序,(n+k−1k)\binom{n+k-1}{k}。
例书与水果

55 本不同的书选 33 本上架(书不能重复):(53)=10\binom{5}{3}=10。55 种水果各不限量,篮子装 33 个(可同种):(5+3−13)=(73)=35\binom{5+3-1}{3}=\binom{7}{3}=35。

二项式定理

Pascal 恒等式是二项式定理的根。中学里的 (a±b)2=a2±2ab+b2(a\pm b)^2=a^2\pm 2ab+b^2 只是 n=2n=2 的特例。

定理二项式定理

对实数 xx、yy 和非负整数 nn,

(x+y)n=∑k=0n(nk)xkyn−k.(x+y)^n=\sum_{k=0}^{n}\binom{n}{k}x^k y^{n-k}.
证明

对 nn 归纳。n=1n=1 时两边都是 x+yx+y。假设 n−1n-1 成立,则

(x+y)n=(x+y)∑k=0n−1(n−1k)xkyn−1−k.(x+y)^n=(x+y)\sum_{k=0}^{n-1}\binom{n-1}{k}x^k y^{n-1-k}.

拆成两段求和,第一段令 i=k+1i=k+1,第二段令 i=ki=k,再用 Pascal 恒等式把中间系数合成 (ni)\binom{n}{i},得到

(x+y)n=∑i=0n(ni)xiyn−i.(x+y)^n=\sum_{i=0}^{n}\binom{n}{i}x^i y^{n-i}.

更干净的组合证明见 AoPS:展开时每个因子贡献 xx 或 yy,选 kk 个因子贡献 xx 的方式有 (nk)\binom{n}{k} 种。

例低次展开

另一种常见写法是

(a+b)n=∑k=0n(nk)an−kbk.(a+b)^n=\sum_{k=0}^{n}\binom{n}{k}a^{n-k}b^k.
注两种写法

只是把 xx、yy 的角色对调,不是另一条定理。

(a+b)3=a3+3a2b+3ab2+b3,(a+b)4=a4+4a3b+6a2b2+4ab3+b4,(a+b)5=a5+5a4b+10a3b2+10a2b3+5ab4+b5.\begin{aligned} (a+b)^3&=a^3+3a^2b+3ab^2+b^3, \\ (a+b)^4&=a^4+4a^3b+6a^2b^2+4ab^3+b^4, \\ (a+b)^5&=a^5+5a^4b+10a^3b^2+10a^2b^3+5ab^4+b^5. \end{aligned}

系数排成 Pascal 三角:每一位等于肩上两个之和。例如 (a+b)5(a+b)^5 里 a3b2a^3b^2 的系数 1010,来自 (a+b)4(a+b)^4 里 a4ba^4b 的 44 加上 a3b2a^3b^2 的 66。更多见 Pascal 三角。

子集个数也可以用二项式定理再解释一遍。

证明

大小为 kk 的子集有 (nk)\binom{n}{k} 个,所以

∑k=0n(nk)=(1+1)n=2n.\sum_{k=0}^{n}\binom{n}{k}=(1+1)^n=2^n.

多项式定理

三个以上的加项,(a+b+c)n(a+b+c)^n 怎么展开?先看多项式系数。

例SUCCESS

单词 SUCCESS\mathrm{SUCCESS} 有 77 个字母:33 个 S\mathrm{S}、22 个 C\mathrm{C},其余各 11。不同排法

(73,2,1,1)=7!3! 2! 1! 1!=420.\binom{7}{3,2,1,1}=\frac{7!}{3!\,2!\,1!\,1!}=420.
例分组

nn 个不同对象分成 rr 个有标号的组,大小分别为 n1,…,nrn_1,\ldots,n_r(∑ni=n\sum n_i=n),种数是

(nn1,…,nr)=n!n1!⋯nr!=(nn1)(n−n1n2)⋯(nrnr).\binom{n}{n_1,\ldots,n_r} = \frac{n!}{n_1!\cdots n_r!} = \binom{n}{n_1}\binom{n-n_1}{n_2}\cdots\binom{n_r}{n_r}.
记号多项式系数

若 n1+⋯+nr=nn_1+\cdots+n_r=n,定义

(nn1,…,nr)=n!n1!⋯nr!.\binom{n}{n_1,\ldots,n_r}=\frac{n!}{n_1!\cdots n_r!}.
定义多项式定理

对正整数 nn 以及满足 n1+⋯+nk=nn_1+\cdots+n_k=n 的非负整数,

(a1+⋯+ak)n=∑(nn1,…,nk)a1n1⋯aknk,(a_1+\cdots+a_k)^n = \sum\binom{n}{n_1,\ldots,n_k}a_1^{n_1}\cdots a_k^{n_k},

求和跑遍所有这样的 (n1,…,nk)(n_1,\ldots,n_k)。

归纳证明留给习题。这里只写组合证明。

证明

多项式系数 (nn1,…,nk)\binom{n}{n_1,\ldots,n_k} 统计把 nn 个位置分成大小为 nin_i 的 kk 组的方法。展开 (x1+⋯+xk)n(x_1+\cdots+x_k)^n 时,每个因子贡献一个 xix_i,得到单项 x1n1⋯xknkx_1^{n_1}\cdots x_k^{n_k} 的方式数正是这个系数。

这与有重复的排列是同一套数。

例(a+b+c)4(a+b+c)^4
(a+b+c)4=a4+4a3b+4a3c+6a2b2+12a2bc+6a2c2+4ab3+12ab2c+12abc2+4ac3+b4+4b3c+6b2c2+4bc3+c4.\begin{aligned} (a+b+c)^4 &= a^4+4a^3b+4a^3c+6a^2b^2+12a^2bc \\ &\quad+6a^2c^2+4ab^3+12ab^2c+12abc^2+4ac^3 \\ &\quad+b^4+4b^3c+6b^2c^2+4bc^3+c^4. \end{aligned}

Catalan 数

本节待撰写 / 占位标记

本小节关于 Catalan 数(Cn=1n+1(2nn)C_n = \frac{1}{n+1}\binom{2n}{n})、Dyck 路径与合法括号序列计数的内容已列入后续撰写队列。

字符串计数

本节待撰写 / 占位标记

本小节关于字符串排列、字符重数与组合计数的内容已列入后续撰写队列。

习题

练习为什么 0!=10!=1

课文把 0!0! 规定成 11。为什么不能是 00?任选一种说得通的理由即可。

解为什么 0!=10!=1
  • 空排列只有一种:什么都不排。
  • 递推 n!=n(n−1)!n!=n(n-1)! 在 n=1n=1 时要求 1=1⋅0!1=1\cdot 0!,所以 0!=10!=1。
练习一个奇怪的递推

定义 f:N0→N0f:\mathbb{N}_0\to\mathbb{N}_0 为

f(n)={1,n=0,n⋅f(f(n−1)),n>0.f(n)= \begin{cases} 1,&n=0,\\ n\cdot f\bigl(f(n-1)\bigr),&n>0. \end{cases}
  1. f(n)f(n) 是否恒等于 n!n!?证明或反驳。
  2. 若不是,写出 f(n)f(n) 的显式,并举例说明它与 n!n! 如何分道。

提示:对照阶乘自己的递推 n!=n⋅(n−1)!n!=n\cdot(n-1)!。

解一个奇怪的递推

f(0)=1f(0)=1,f(1)=1⋅f(f(0))=f(1)f(1)=1\cdot f(f(0))=f(1) 会自指;若约定 f(1)=1f(1)=1,则 f(2)=2f(f(1))=2f(1)=2f(2)=2f(f(1))=2f(1)=2,f(3)=3f(f(2))=3f(2)=6f(3)=3f(f(2))=3f(2)=6,f(4)=4f(f(3))=4f(6)f(4)=4f(f(3))=4f(6)。问题在于求 f(4)f(4) 需要 f(6)f(6),而 f(6)f(6) 又依赖更大的值,递推并不能沿用阶乘那条链走下去。因此 ff 并没有被定义成处处等于阶乘的函数;阶乘的递推用的是 f(n−1)f(n-1),这里却是 f(f(n−1))f(f(n-1))。

练习证明 rr-排列公式

证明:若 1≤r≤n1\le r\le n,则 P(n,r)=n(n−1)⋯(n−r+1)P(n,r)=n(n-1)\cdots(n-r+1)。

证明

第一位 nn 种,第二位 n−1n-1 种,第 rr 位 n−r+1n-r+1 种。由乘法原理即得。

练习单词重排

下列单词各有多少种字母排法?

  1. Fluke
  2. Propose
  3. Mississippi
  4. Arrange
解单词重排

用 n!n1!⋯nk!\dfrac{n!}{n_1!\cdots n_k!}。

  1. Fluke 无重复:5!=1205!=120。
  2. Propose 中 P\mathrm{P} 两次:7!2!=2520\dfrac{7!}{2!}=2520。
  3. Mississippi 中 S\mathrm{S} 四次、I\mathrm{I} 四次、P\mathrm{P} 两次:11!4! 4! 2!=34650\dfrac{11!}{4!\,4!\,2!}=34650。
  4. Arrange 中 A\mathrm{A} 两次、R\mathrm{R} 两次:7!2! 2!=1260\dfrac{7!}{2!\,2!}=1260。
练习辨认计数模型

回答下列问题。

  1. 1010 人里选主席、财务、秘书,有多少种?
  2. 1010 人里选一个三人小组,有多少种?
  3. (1) 和 (2) 本质差别是什么?哪个更大?不做计算能看出来吗?
  4. 1010 种口味里选三勺冰淇淋(同口味可重复),有多少种?
  5. 五份不同奖品分给 Anastasia、Becky、Cadel(可以有人空手),有多少种?
  6. 六匹马跑完全程且没有并列,终点顺序有多少种?
解辨认计数模型
  1. 职位不同,用排列:P(10,3)=10×9×8=720P(10,3)=10\times 9\times 8=720。
  2. 无序组合:(103)=120\binom{10}{3}=120。
  3. 差别就是要不要次序。排列比组合大,因为每种三人组对应 3!3! 种职务分配。
  4. 若勺的次序也算,是 103=100010^3=1000。若只关心口味组合,是可重复组合 (10+3−13)=(123)=220\binom{10+3-1}{3}=\binom{12}{3}=220。
  5. 每份奖独立送给三人之一:35=2433^5=243。
  6. 6!=7206!=720。
练习用二项式定理归纳证明多项式定理

用二项式定理,对加项个数作归纳,证明多项式定理。

解用二项式定理归纳证明多项式定理

k=2k=2 就是二项式定理。假设对 kk 成立。写

(x1+⋯+xk+xk+1)n=((x1+⋯+xk)+xk+1)n=∑i=0n(ni)(x1+⋯+xk)n−ixk+1i.(x_1+\cdots+x_k+x_{k+1})^n = \bigl((x_1+\cdots+x_k)+x_{k+1}\bigr)^n = \sum_{i=0}^{n}\binom{n}{i}(x_1+\cdots+x_k)^{n-i}x_{k+1}^i.

内层用归纳假设展开成多项式系数,再与 (ni)\binom{n}{i} 合成 k+1k+1 项的多项式系数。

概率公理

有了“有多少种结果”,就可以谈某个事件占多大份额。

样本空间与事件

掷一枚骰子,结果写成 S={1,2,3,4,5,6}S=\{1,2,3,4,5,6\},掷出 33 的可能性是 1/∣S∣=1/61/|S|=1/6。SS 叫做样本空间,掷出 33 是一个事件,事件是样本空间的子集。

定义样本空间

实验或随机试验的样本空间 SS 是全部可能结果的集合。这些结果互斥,并且合在一起穷尽所有可能性。

定义事件

事件是 SS 的任意子集,表示若干可能结果的汇集。可以空(空事件 ∅\emptyset),也可以是整个 SS。

例抽一张牌

标准 5252 张牌抽一张。样本空间有 5252 个元素。例如:

  • DD:抽到 J、Q 或 K
  • EE:抽到红心
  • FF:抽到 A
例掷两次硬币

公平硬币掷两次。样本空间

S={HH,HT,TH,TT}.S=\{HH,HT,TH,TT\}.
  • GG:至少一次正面,G={HH,HT,TH}G=\{HH,HT,TH\}
  • HH:第二次是反面,H={HT,TT}H=\{HT,TT\}

事件也是集合,交、并照常做。上例里“至少一次正面,且第二次不是反面”就是

G∩H={HT}.G\cap H=\{HT\}.
注交的简写

概率里常把 G∩HG\cap H 写成 GHGH。

事件很多时,用 ⋃\bigcup、⋂\bigcap 写一串运算。

记号可数并与交

事件列 {En}n=1∞\{E_n\}_{n=1}^{\infty} 的并 ⋃n=1∞En\bigcup_{n=1}^{\infty}E_n 由至少属于某个 EnE_n 的结果组成;交 ⋂n=1∞En\bigcap_{n=1}^{\infty}E_n 由属于每一个 EnE_n 的结果组成。

补事件写 HcH^c。上例 Hc=S∖H={HH,TH}H^c=S\setminus H=\{HH,TH\}。

例De Morgan

布尔代数和集合论里的 De Morgan 律,对任意多个事件仍成立:

(⋃i=1nEi)c=⋂i=1nEic,(⋂i=1nEi)c=⋃i=1nEic.\begin{aligned} \Bigl(\bigcup_{i=1}^n E_i\Bigr)^c &= \bigcap_{i=1}^n E_i^c, \\ \Bigl(\bigcap_{i=1}^n E_i\Bigr)^c &= \bigcup_{i=1}^n E_i^c. \end{aligned}
注补的几种写法

ScS^c、S′S'、S‾\overline{S} 表示同一件事。

概率公理

口语里“概率就是机会大小”,数学上要先给定义。两种常见说法如下。

定义频率极限

若事件 EE 在 nn 次试验里出现 nEn_E 次,且极限存在,则

P(E)=lim⁡n→∞nEn.P(E)=\lim_{n\to\infty}\frac{n_E}{n}.
定义古典概型

若样本空间有限且每个结果等可能,则

P(E)=∣E∣∣S∣.P(E)=\frac{|E|}{|S|}.

Kolmogorov 在 19331933 年给出的三条公理,是现代概率的骨架。满足这些公理的系统都叫做概率空间。

公理非负性

对样本空间 SS 中任意事件 EE,

0≤P(E)≤1.0\le P(E)\le 1.
公理规范性
P(S)=1.P(S)=1.
公理可数可加性

对两两不交的事件列 {Ei}\{E_i\},

P(⋃i=1∞Ei)=∑i=1∞P(Ei).P\Bigl(\bigcup_{i=1}^{\infty}E_i\Bigr)=\sum_{i=1}^{\infty}P(E_i).

由公理立刻得到若干命题。

命题补事件
P(Ec)=1−P(E).P(E^c)=1-P(E).
证明

SS 剖成 EE 与 EcE^c 两块,不相交,所以 P(S)=P(E)+P(Ec)=1P(S)=P(E)+P(E^c)=1。

命题单调性

若 E⊆FE\subseteq F,则 P(E)≤P(F)P(E)\le P(F)。

证明

F=E∪(F∖E)F=E\cup(F\setminus E),两块不相交,所以 P(F)=P(E)+P(F∖E)P(F)=P(E)+P(F\setminus E)。后一项非负,故 P(E)≤P(F)P(E)\le P(F)。

命题两个事件的容斥
P(E∪F)=P(E)+P(F)−P(E∩F).P(E\cup F)=P(E)+P(F)-P(E\cap F).
证明

E∪F=E∪(F∖E)E\cup F=E\cup(F\setminus E),两块不相交。而 P(F∖E)=P(F)−P(E∩F)P(F\setminus E)=P(F)-P(E\cap F)。

集合论里的容斥可以写成对所有非空下标子集交替求和。

定理一般容斥(基数)

对有限集 A1,…,AnA_1,\ldots,A_n,

∣⋃i=1nAi∣=∑∅≠I⊆[n](−1)∣I∣+1∣⋂i∈IAi∣.\left|\bigcup_{i=1}^{n}A_i\right| = \sum_{\emptyset\neq I\subseteq[n]}(-1)^{|I|+1}\left|\bigcap_{i\in I}A_i\right|.
证明

在全集 XX 上令 fif_i 为 AiA_i 的特征函数,

F(x)=∏i=1n(1−fi(x)).F(x)=\prod_{i=1}^{n}(1-f_i(x)).

F(x)=1F(x)=1 当且仅当 xx 不属于任何 AiA_i。展开乘积:

F(x)=∑I⊆[n](−1)∣I∣∏i∈Ifi(x).F(x)=\sum_{I\subseteq[n]}(-1)^{|I|}\prod_{i\in I}f_i(x).

对 x∈Xx\in X 求和,左边是 ∣X∖⋃Ai∣|X\setminus\bigcup A_i|,右边是各重交的交替和。整理后即得容斥。

概率版本不必另证。

定理一般容斥(概率)
P(⋃i=1nAi)=∑∅≠I⊆[n](−1)∣I∣+1P(⋂i∈IAi).P\Bigl(\bigcup_{i=1}^{n}A_i\Bigr) = \sum_{\emptyset\neq I\subseteq[n]}(-1)^{|I|+1}P\Bigl(\bigcap_{i\in I}A_i\Bigr).
证明

示性函数逐点满足

1∪iAi=1−∏i(1−1Ai)=∑∅≠I⊆[n](−1)∣I∣+11∩i∈IAi.\mathbf1_{\cup_i A_i} =1-\prod_i(1-\mathbf1_{A_i}) =\sum_{\emptyset\ne I\subseteq[n]}(-1)^{|I|+1}\mathbf1_{\cap_{i\in I}A_i}.

两边取期望就是概率公式。这里只有有限项,不涉及交换无穷级数,单个样本点概率为零的空间也同样适用。

组合地看:某个结果若恰好落在 m>0m>0 个 AiA_i 里,它对并的概率贡献一次。在容斥展开里它被加 (m1)\binom{m}{1} 次、减 (m2)\binom{m}{2} 次,如此交替。因为

∑i=0m(mi)(−1)i=(1−1)m=0,\sum_{i=0}^{m}\binom{m}{i}(-1)^i=(1-1)^m=0,

所以从 i=1i=1 加到 mm 恰好等于 11。每个结果被正确地点了一次。

习题

练习Boole 不等式

证明 Boole 不等式:对有限样本空间中的事件 E1,…,EnE_1,\ldots,E_n,

P(⋃i=1nEi)≤∑i=1nP(Ei).P\Bigl(\bigcup_{i=1}^n E_i\Bigr)\le\sum_{i=1}^n P(E_i).
定理Boole 不等式
P(⋃i=1nEi)≤∑i=1nP(Ei).P\Bigl(\bigcup_{i=1}^n E_i\Bigr)\le\sum_{i=1}^n P(E_i).
证明

对 nn 归纳。n=1n=1 显然。假设对 nn 成立。则

P(⋃i=1n+1Ei)=P((⋃i=1nEi)∪En+1)≤P(⋃i=1nEi)+P(En+1),P\Bigl(\bigcup_{i=1}^{n+1}E_i\Bigr) = P\Bigl(\bigl(\bigcup_{i=1}^{n}E_i\bigr)\cup E_{n+1}\Bigr) \le P\Bigl(\bigcup_{i=1}^{n}E_i\Bigr)+P(E_{n+1}),

再用归纳假设。

证明

两个事件时,P(E1∪E2)=P(E1)+P(E2)−P(E1∩E2)≤P(E1)+P(E2)P(E_1\cup E_2)=P(E_1)+P(E_2)-P(E_1\cap E_2)\le P(E_1)+P(E_2)。再对个数归纳,或者直接看一般容斥里后面那些交都带着符号,整体不会超过第一项之和。

注何时取等

当且仅当诸 EiE_i 两两不交时等号成立,此时交的概率全是 00,可加性给出等式。

练习灌铅骰子

一枚六面骰,掷出 33 的可能性是掷出其他某个指定点数的两倍。求各面概率。

解灌铅骰子

设 p(3)=2p(1)p(3)=2p(1),且 p(1)=p(2)=p(4)=p(5)=p(6)p(1)=p(2)=p(4)=p(5)=p(6)。总和为 11:

5p(1)+2p(1)=1  ⟹  p(1)=17,p(3)=27.5p(1)+2p(1)=1\implies p(1)=\frac17,\qquad p(3)=\frac27.
练习并与交的下界

p(E)=0.8p(E)=0.8,p(F)=0.6p(F)=0.6。证明 p(E∪F)≥0.8p(E\cup F)\ge 0.8 且 p(E∩F)≥0.4p(E\cap F)\ge 0.4。

解并与交的下界

p(E∪F)=0.8+0.6−p(E∩F)≤1p(E\cup F)=0.8+0.6-p(E\cap F)\le 1,故 p(E∩F)≥0.4p(E\cap F)\ge 0.4。另一方面 E⊆E∪FE\subseteq E\cup F,所以 p(E∪F)≥p(E)=0.8p(E\cup F)\ge p(E)=0.8。

练习Bonferroni 不等式

证明:对事件 EE、FF,

p(E∩F)≥p(E)+p(F)−1.p(E\cap F)\ge p(E)+p(F)-1.
定理Bonferroni 不等式
p(E∩F)≥p(E)+p(F)−1.p(E\cap F)\ge p(E)+p(F)-1.
证明

p(E∪F)=p(E)+p(F)−p(E∩F)≤1p(E\cup F)=p(E)+p(F)-p(E\cap F)\le 1,移项即得。

练习一般 Bonferroni 不等式

用归纳把 Bonferroni 不等式推到 nn 个事件:

P(⋂i=1nEi)≥∑i=1nP(Ei)−(n−1).P\Bigl(\bigcap_{i=1}^n E_i\Bigr)\ge\sum_{i=1}^n P(E_i)-(n-1).

也可用容斥另证。

定理一般 Bonferroni 不等式
P(⋂i=1nEi)≥∑i=1nP(Ei)−(n−1).P\Bigl(\bigcap_{i=1}^n E_i\Bigr)\ge\sum_{i=1}^n P(E_i)-(n-1).
证明

弱归纳。n=2n=2 就是 Bonferroni。假设对 kk 成立。把前 kk 个的交与 Ek+1E_{k+1} 再用二元 Bonferroni:

P((⋂i=1kEi)∩Ek+1)≥P(⋂i=1kEi)+P(Ek+1)−1,P\Bigl(\bigl(\bigcap_{i=1}^k E_i\bigr)\cap E_{k+1}\Bigr) \ge P\Bigl(\bigcap_{i=1}^k E_i\Bigr)+P(E_{k+1})-1,

代入归纳假设即得 k+1k+1 的情形。

证明

强归纳可以从 n=1n=1(此时右边是 P(E1)P(E_1))起,步骤相同。

注弱归纳与强归纳

弱归纳只假设 n=kn=k,强归纳假设所有更小的情形。问题只依赖前一步时,弱归纳更干净;需要整段历史时,用强归纳。

证明

用容斥。

P(⋂i=1nEi)=1−P(⋃i=1nEic)≥1−∑i=1nP(Eic)=1−∑i=1n(1−P(Ei))=∑i=1nP(Ei)−(n−1).P\Bigl(\bigcap_{i=1}^n E_i\Bigr) = 1-P\Bigl(\bigcup_{i=1}^n E_i^c\Bigr) \ge 1-\sum_{i=1}^n P(E_i^c) = 1-\sum_{i=1}^n\bigl(1-P(E_i)\bigr) = \sum_{i=1}^n P(E_i)-(n-1).
练习可数可加性作为极限

设 E1,E2,…E_1,E_2,\ldots 两两不交。通过取极限证明

p(⋃i=1∞Ei)=∑i=1∞p(Ei).p\Bigl(\bigcup_{i=1}^\infty E_i\Bigr)=\sum_{i=1}^\infty p(E_i).
解可数可加性作为极限

有限可加性给出 p(⋃i=1nEi)=∑i=1np(Ei)p(\bigcup_{i=1}^n E_i)=\sum_{i=1}^n p(E_i)。并随着 nn 递增,由概率的下连续性,

p(⋃i=1∞Ei)=lim⁡n→∞p(⋃i=1nEi)=∑i=1∞p(Ei).p\Bigl(\bigcup_{i=1}^\infty E_i\Bigr) = \lim_{n\to\infty}p\Bigl(\bigcup_{i=1}^n E_i\Bigr) = \sum_{i=1}^\infty p(E_i).
练习两个骰子上的事件

掷两枚骰子。EE 为点数之和为奇数,FF 为至少一枚是 11,GG 为和为 55。描述 EFEF、E∪FE\cup F、FGFG、EFcEF^c、EFGEFG。

解两个骰子上的事件
  • EFEF:和为奇数且至少一枚是 11,即 {(1,2),(1,4),(1,6),(2,1),(4,1),(6,1)}\{(1,2),(1,4),(1,6),(2,1),(4,1),(6,1)\}。
  • E∪FE\cup F:和为奇数,或至少一枚是 11。
  • FGFG:至少一枚是 11 且和为 55,即 {(1,4),(4,1)}\{(1,4),(4,1)\}。
  • EFcEF^c:和为奇数且没有一枚是 11。
  • EFGEFG:与 FGFG 相同,因为和为 55 已经是奇数。
练习三家报纸

某镇 100,000100{,}000 人,三家报纸 I、II、III 的阅读比例为:I 10%10\%,II 30%30\%,III 5%5\%;I 与 II 8%8\%,I 与 III 2%2\%,II 与 III 4%4\%;三家都看 1%1\%。

  1. 只看一家的人数
  2. 至少看两家的人数
  3. 若 I、III 是早报,II 是晚报,至少看一份早报再加晚报的人数
  4. 谁也不看的人数
  5. 只看一份早报加一份晚报的人数
解三家报纸
∣I∣=10,000,∣II∣=30,000,∣III∣=5,000,∣I∩II∣=8,000,∣I∩III∣=2,000,∣II∩III∣=4,000,∣I∩II∩III∣=1,000.\begin{aligned} |I|&=10{,}000,& |II|&=30{,}000,& |III|&=5{,}000, \\ |I\cap II|&=8{,}000,& |I\cap III|&=2{,}000,& |II\cap III|&=4{,}000, \\ |I\cap II\cap III|&=1{,}000. \end{aligned}

只看 I:10,000−(8,000+2,000−1,000)=1,00010{,}000-(8{,}000+2{,}000-1{,}000)=1{,}000。只看 II:19,00019{,}000。只看 III:00。

  1. 1,000+19,000+0=20,0001{,}000+19{,}000+0=20{,}000。
  2. 8,000+2,000+4,000−2×1,000=12,0008{,}000+2{,}000+4{,}000-2\times 1{,}000=12{,}000。
  3. 至少一份早报加晚报对应 (I∩II)∪(III∩II)(I\cap II)\cup(III\cap II)。三家都看的读者只能计一次,因此人数为 8,000+4,000−1,000=11,0008{,}000+4{,}000-1{,}000=11{,}000。
  4. 看报总人数用容斥:10,000+30,000+5,000−8,000−2,000−4,000+1,000=32,00010{,}000+30{,}000+5{,}000-8{,}000-2{,}000-4{,}000+1{,}000=32{,}000,谁也不看的人数为 100,000−32,000=68,000100{,}000-32{,}000=68{,}000。
  5. 恰好一份早报加晚报必须排除三家都看的读者:只看 I、II 的有 7,0007{,}000 人,只看 III、II 的有 3,0003{,}000 人,合计 10,00010{,}000。

第 3 问允许两份早报都看,第 5 问要求恰好一份早报;两问都统计人数,所以每位读者只能计一次。

练习队员的职业与党派

1515 名队员,每人 22 种职业、33 种政治归属。

  • 样本空间有多少个结果?
  • “至少一名蓝领”有多少个结果?
  • “没有人自认为独立人士”有多少个结果?
解队员的职业与党派

每人 2×3=62\times 3=6 种组合,全体 6156^{15}。

全是白领时每人只剩 33 种党派,故至少一名蓝领:615−3156^{15}-3^{15}。

没有独立人士时每人 22 种职业 ×\times 22 种党派,共 4154^{15}。

练习两点数之和

掷两枚骰子,点数之和等于 ii 的概率,i=2,3,…,12i=2,3,\ldots,12。

解两点数之和

总共 3636 种等可能结果。和为 22 到 1212 的方式数依次是 1,2,3,4,5,6,5,4,3,2,11,2,3,4,5,6,5,4,3,2,1。例如和为 77 有 66 种,概率 6/36=1/66/36=1/6。

练习5 先于 7

反复掷一对骰子,直到出现和为 55 或和为 77。求先出现 55 的概率。

解5 先于 7

和为 55 有 44 种,和为 77 有 66 种,所以单次 P(5)=4/36=1/9P(5)=4/36=1/9,P(7)=6/36=1/6P(7)=6/36=1/6。

EnE_n 的概率

记 EnE_n 为:前 n−1n-1 次既没有 55 也没有 77,第 nn 次出现 55。单次既不是 55 也不是 77 的概率是 26/36=13/1826/36=13/18,于是

P(En)=(2636)n−1⋅436.P(E_n)=\Bigl(\frac{26}{36}\Bigr)^{n-1}\cdot\frac{4}{36}.

对 P(En)P(E_n) 求和

∑n=1∞P(En)=4/361−26/36=410=25.\sum_{n=1}^{\infty}P(E_n) = \frac{4/36}{1-26/36} = \frac{4}{10} = \frac{2}{5}.
练习Bell 数

集合 SS 的划分是一族互不相交的非空子集,并起来等于 SS。记 TnT_n 为 {1,…,n}\{1,\ldots,n\} 的划分个数。已知 T1=1T_1=1,T2=2T_2=2。

  1. 列出全部划分,验证 T3=5T_3=5,T4=15T_4=15。
  2. 证明 Tn+1=1+∑k=1n(nk)TkT_{n+1}=1+\sum_{k=1}^{n}\binom{n}{k}T_k,并由此算 T10T_{10}。
注Bell 数

TnT_n 就是 Bell 数,清点 nn 元集的划分。

解Bell 数

(a) T3T_3:整块一块;11 单独、另两个一起(三种);三个都单独。共 55。

T4T_4:整块;1+31+3 型四种;2+22+2 型三种;2+1+12+1+1 型六种;四个单点。共 1515。

(b) 把 n+1n+1 当作特殊元素。它若独自成块,剩下 TnT_n 种。它若与某 kk 个元素同块(1≤k≤n1\le k\le n),先从 nn 个里选这 kk 个,剩下 n−kn-k 个有 Tn−kT_{n-k} 种划分。约定 T0=1T_0=1,可改写成

Tn+1=1+∑k=1n(nk)Tk.T_{n+1}=1+\sum_{k=1}^{n}\binom{n}{k}T_k.

递推得 T5=52T_5=52,T6=203T_6=203,T7=877T_7=877,T8=4140T_8=4140,T9=21147T_9=21147,T10=115975T_{10}=115975。

用计数求概率

许多问题假定样本空间里每个结果等可能。此时

P(E)=∣E∣∣S∣,P(E)=\frac{|E|}{|S|},

分子分母都可以用前面的计数方法来算。套公式之前仍要先想清楚:什么算不同的结果。

一些基本问题

例和为 7

掷两枚骰子,点数之和为 77 的概率?

解和为 7

3636 种等可能结果,有利的是 {(1,6),(2,5),(3,4),(4,3),(5,2),(6,1)}\{(1,6),(2,5),(3,4),(4,3),(5,2),(6,1)\},共 66 种,概率 6/36=1/66/36=1/6。

例委员会里的男女

66 男 99 女,随机选 55 人委员会。恰好 33 男 22 女的概率?

解委员会里的男女
(63)=20,(92)=36,(155)=3003,\binom{6}{3}=20, \qquad \binom{9}{2}=36, \qquad \binom{15}{5}=3003,

所求

P=20×363003=7203003=2401001≈0.2398.P=\frac{20\times 36}{3003}=\frac{720}{3003}=\frac{240}{1001}\approx 0.2398.
例特殊球

罐中 nn 个球,其中一个特别。不放回地取 kk 个,每次剩下的球等可能。特殊球被取到的概率?

解特殊球
P=(n−1k−1)(nk)=kn.P=\frac{\binom{n-1}{k-1}}{\binom{n}{k}}=\frac{k}{n}.

另一种算法:记 AiA_i 为特殊球在第 ii 次被抽到,P(Ai)=1/nP(A_i)=1/n,诸 AiA_i 不交,于是 P(⋃i=1kAi)=k/nP(\bigcup_{i=1}^k A_i)=k/n。

例三项运动

俱乐部里打网球 3636 人,壁球 2828 人,羽毛球 1818 人;网球且壁球 2222 人,网球且羽毛球 1212 人,壁球且羽毛球 99 人,三项都会 44 人。至少参加一项的有多少人?

解三项运动

随机抽一名会员,令 P(C)P(C) 为抽到的人属于集合 CC 的概率。容斥给出

P(T∪S∪B)=36+28+18−22−12−9+4N=43N,P(T\cup S\cup B) = \frac{36+28+18-22-12-9+4}{N} = \frac{43}{N},

所以有 4343 人至少参加一项。

再深入一些的问题

下面几题更绕,适合当作选做。

例桥牌

5252 张牌分给 44 人。求:

  1. 有一人拿齐 1313 张黑桃的概率
  2. 每人恰好一张 A 的概率
解桥牌

(a) 记 EiE_i 为第 ii 家拿齐黑桃。四家互斥,且

P(Ei)=1(5213),P(⋃i=14Ei)=4(5213)≈6.3×10−12.P(E_i)=\frac{1}{\binom{52}{13}}, \qquad P\Bigl(\bigcup_{i=1}^{4}E_i\Bigr)=\frac{4}{\binom{52}{13}}\approx 6.3\times 10^{-12}.

(b) 先把 44 张 A 放在一边,其余 4848 张按每人 1212 张来分,有 (4812,12,12,12)\binom{48}{12,12,12,12} 种;A 再按 4!4! 种方式每人一张。总发牌方式 (5213,13,13,13)\binom{52}{13,13,13,13}。于是

4! (4812,12,12,12)(5213,13,13,13)≈0.105.\frac{4!\,\binom{48}{12,12,12,12}}{\binom{52}{13,13,13,13}}\approx 0.105.
例室友

球队 2020 名进攻队员、2020 名防守队员,随机两两配对住。没有“进攻-防守”混合室友对的概率?

解室友

4040 人分成 2020 个无序对,种数是 40!220 20!\dfrac{40!}{2^{20}\,20!}。进攻内部互配、防守内部互配的种数是 (20!210 10!)2\bigl(\dfrac{20!}{2^{10}\,10!}\bigr)^2。于是

P0=(20!/(210 10!))240!/(220 20!)=(20!)3(10!)2 (40!).P_0 = \frac{\bigl(20!/(2^{10}\,10!)\bigr)^2}{40!/(2^{20}\,20!)} = \frac{(20!)^3}{(10!)^2\,(40!)}.
例生日问题

房间里 nn 个人,没有两人生日相同的概率?nn 要多大,这个概率才小于 1/21/2?

解生日问题

忽略闰日,总共 365n365^n 种生日指派。没有重复的概率是

P=365×364×⋯×(365−n+1)365n.P=\frac{365\times 364\times\cdots\times(365-n+1)}{365^n}.

n≥23n\ge 23 时 P<1/2P<1/2。2323 远小于 365365,但一对对来看已有 (232)=253\binom{23}{2}=253 对。n=50n=50 时至少两人同日的概率约 0.9700.970;n=100n=100 时优势超过三百万比一。

例圆桌夫妇

1010 对夫妇随机坐圆桌,没有一对夫妇女邻座的概率?

解圆桌夫妇

记 EiE_i 为第 ii 对邻座。所求是 1−P(⋃i=110Ei)1-P(\bigcup_{i=1}^{10}E_i)。圆桌 2020 人有 19!19! 种排法。指定 nn 对必须邻座时,把每对捆成一块有 2n2^n 种朝向,剩下 20−n20-n 块(把每对看成一人)绕圆排有 (19−n)!(19-n)! 种,于是

P(Ei1⋯Ein)=2n(19−n)!19!.P(E_{i_1}\cdots E_{i_n})=\frac{2^n(19-n)!}{19!}.

代入容斥,没有一对邻座的概率约 0.33950.3395。

习题

练习(i,j)(i,j)-可分等差数列

设 mm 为正整数,a1,…,a4m+2a_1,\ldots,a_{4m+2} 是公差非零的等差数列。若任意去掉两项 ai,aja_i,a_j(i<ji<j)后,剩下 4m4m 项能分成 mm 组、每组四个数仍成等差,则称原数列为 (i,j)(i,j)-可分数列。

  1. 写出使 a1,…,a6a_1,\ldots,a_6 成为 (i,j)(i,j)-可分数列的全部 (i,j)(i,j),1≤i<j≤61\le i<j\le 6。
  2. 证明 m≥3m\ge 3 时,a1,…,a4m+2a_1,\ldots,a_{4m+2} 是 (2,13)(2,13)-可分数列。
  3. 从 1,2,…,4m+21,2,\ldots,4m+2 中随机取两个整数 i<ji<j,记原数列为 (i,j)(i,j)-可分的概率为 pmp_m。证明 pm>1/8p_m>1/8。
解(i,j)(i,j)-可分等差数列

(1) 去掉开头两项、末尾两项,或首尾各一项,剩下四个数仍成等差。于是 (1,2)(1,2)、(5,6)(5,6)、(1,6)(1,6)。

(2) 公差为 dd 时 an=a1+(n−1)da_n=a_1+(n-1)d。m>3m>3 时,下标 1515 以后的连续四项本身已成等差,只需检查前 1414 项。取公差 3d3d 的三组

{a1,a4,a7,a10},{a3,a6,a9,a12},{a5,a8,a11,a14}\{a_1,a_4,a_7,a_{10}\},\quad \{a_3,a_6,a_9,a_{12}\},\quad \{a_5,a_8,a_{11},a_{14}\}

覆盖去掉 a2a_2、a13a_{13} 之后的前 1212 项,后面按四项一组切即可。公差 2d2d 组不起来,但 3d3d 已经够用。

(3) 可分的 (i,j)(i,j) 由两组底样生成:(1,2)(1,2) 与 (2,9)(2,9)(后者在 m=2m=2 时出现),再沿下标平移 44 的倍数。表中打勾的位置给出计数 (m+1)2−m(m+1)^2-m。总选法 (4m+22)\binom{4m+2}{2},于是

pm=m2+m+18m2+6m+1.p_m=\frac{m^2+m+1}{8m^2+6m+1}.

令 g(m)=pm−1/8g(m)=p_m-1/8,则

g(m)=2m+78(8m2+6m+1)>0g(m)=\frac{2m+7}{8(8m^2+6m+1)}>0

对一切正整数 mm 成立。

m≤3m\le 3 与 m≤5m\le 5 的 (i,j)(i,j) 表如下(列为较小下标 1,5,9,…1,5,9,\ldots,行为 2,6,10,…2,6,10,\ldots)。

m≤3m\le 3:

15913
2✓✓✓
6✓✓✓
10✓✓✓
14✓✓✓✓

m≤5m\le 5:

159131721
2✓✓✓✓✓
6✓✓✓✓✓
10✓✓✓✓✓
14✓✓✓✓✓
18✓✓✓✓✓
22✓✓✓✓✓✓