UA MATH564 概率论 计算至少有一个发生的概率:Waring公式

    技术2022-07-12  82

    UA MATH564 概率论 计算至少有一个发生的概率:Waring公式

    基于组合学的证明方法基于概率论的证明方法

    在推出Poincare公式后,计算 n n n个事件中至少有一个事件发生的概率的问题就彻底解决了。但数学家们又开始纠结了,Poincare公式只能够计算至少有一个发生的概率,那有没有能计算恰好有 m m m个事件同时发生的概率的方法呢?

    假设用 B m B_m Bm表示 A 1 , ⋯   , A n A_1,\cdots,A_n A1,,An中正好有 m m m个发生的事件,则 P ( B m ) = ∑ k = m n ( − 1 ) k − m C k m S k P(B_m) = \sum_{k=m}^n (-1)^{k-m}C_k^mS_k P(Bm)=k=mn(1)kmCkmSk

    这个公式叫做Waring公式。其中 S k S_k Sk就是上一讲定义的记号 S m = ∑ 1 ≤ i 1 < i 2 < ⋯ < i m ≤ n P ( A i 1 ∩ A i 2 ∩ ⋯ ∩ A i m ) S_m = \sum_{1 \le i_1 < i_2 < \cdots < i_m \le n} P(A_{i_1} \cap A_{i_2} \cap \cdots \cap A_{i_m}) Sm=1i1<i2<<imnP(Ai1Ai2Aim)

    基于组合学的证明方法

    Feller的Introduction to Probability and its application中给出了基于组合学的证明方法。用反证法。我们要证明的公式等价于:假设 E E E是一个恰好包含在 l l l个事件中的子事件,当且仅当 l = m l=m l=m时, P ( E ) P(E) P(E)才是 P ( B m ) P(B_m) P(Bm)的一部分。显然这时 P ( E ) P(E) P(E)也是 S 1 , S 2 , ⋯   , S l S_1,S_2,\cdots,S_l S1,S2,,Sl的一部分,并且 P ( E ) P(E) P(E)只会出现在 S l S_l Sl的一项中。如果 l < m l <m l<m,则 P ( E ) P(E) P(E)不是 P ( B m ) P(B_m) P(Bm)的一部分。如果 l > m l>m l>m,则 P ( E ) P(E) P(E)会出现在 S k , k ≥ m S_k,k\ge m Sk,km C l k C_l^k Clk项中,因此对 S k , k ≥ m S_k,k\ge m Sk,km的贡献是 C l k C k m C_l^k C_k^m ClkCkm根据容斥原理, P ( E ) P(E) P(E) P ( B m ) P(B_m) P(Bm)的贡献为 C l m C m m − C l m + 1 C m + 1 m + C l m + 2 C m + 2 m + 1 − ⋯ = C l m ( C l − m 0 − C l − m 1 + ⋯   ) = C l m ( 1 − 1 ) n − m = 0 C_l^mC_m^m - C_l^{m+1}C_{m+1}^m + C_l^{m+2}C_{m+2}^{m+1} - \cdots \\ = C_l^m(C_{l-m}^0 - C_{l-m}^1 + \cdots) = C_l^m(1-1)^{n-m} = 0 ClmCmmClm+1Cm+1m+Clm+2Cm+2m+1=Clm(Clm0Clm1+)=Clm(11)nm=0

    这就说明了当且仅当 l = m l=m l=m时, P ( E ) P(E) P(E)才是 P ( B m ) P(B_m) P(Bm)的一部分,因此上述公式成立。

    基于概率论的证明方法

    施利亚耶夫的概率论习题集提到了基于事件的Indicator function与期望的证明方法。我们记恰好发生的 m m m个事件为 i 1 , i 2 , ⋯   , i m i_1,i_2,\cdots,i_m i1,i2,,im,记没有发生的事件为 j 1 , j 2 , ⋯   , j n − m j_1,j_2,\cdots,j_{n-m} j1,j2,,jnm,记随机变量 X = I A X = I_A X=IA,则恰好有 m m m个事件发生的随机变量是 ∑ X i 1 X i 2 ⋯ X i m ( 1 − X j 1 ) ⋯ ( 1 − X j n − m ) \sum X_{i_1}X_{i_2}\cdots X_{i_m} (1-X_{j_1})\cdots (1-X_{j_{n-m}}) Xi1Xi2Xim(1Xj1)(1Xjnm)

    上式是对任意可能的 i 1 , ⋯   , i m i_1,\cdots,i_m i1,,im的选取的求和。事件的Indicator function的期望等于事件发生的概率,因此 P ( B m ) = E ∑ X i 1 X i 2 ⋯ X i m ( 1 − X j 1 ) ⋯ ( 1 − X j n − m ) P(B_m) = E\sum X_{i_1}X_{i_2}\cdots X_{i_m} (1-X_{j_1})\cdots (1-X_{j_{n-m}}) P(Bm)=EXi1Xi2Xim(1Xj1)(1Xjnm)

    注意 ( 1 − X j 1 ) ⋯ ( 1 − X j n − m ) (1-X_{j_1})\cdots (1-X_{j_{n-m}}) (1Xj1)(1Xjnm)的展开式中一共有 2 n − m 2^{n-m} 2nm项,其中一项为 1 1 1,因此概率中有一项为 E ∑ X i 1 X i 2 ⋯ X i m = ∑ P ( A i 1 ∩ A i 2 ∩ ⋯ ∩ A i m ) E\sum X_{i_1}X_{i_2}\cdots X_{i_m} = \sum P(A_{i_1} \cap A_{i_2} \cap \cdots \cap A_{i_m}) EXi1Xi2Xim=P(Ai1Ai2Aim)

    其余展开式期望的计算与此类似,事件的Indicator function乘积的期望等于事件的交发生的概率。下面研究 ( 1 − X j 1 ) ⋯ ( 1 − X j n − m ) (1-X_{j_1})\cdots (1-X_{j_{n-m}}) (1Xj1)(1Xjnm)的展开式。

    引理1 用示性函数表示的容斥原理为: X i = I A i X_i = I_{A_i} Xi=IAi I ⋃ i = 1 n A i = 1 − ∏ i = 1 n ( 1 − X i ) = ∑ i = 1 n X i − ∑ 1 ≤ i 1 < i 2 ≤ n X i 1 X i 2 + ⋯ + ( − 1 ) m + 1 ∑ 1 ≤ i 1 < ⋯ < i m ≤ n X i 1 ⋯ X i m ⋯ + ( − 1 ) n X 1 ⋯ X n I_{\bigcup_{i=1}^n A_i} =1 - \prod_{i=1}^n (1-X_i) = \sum_{i=1}^n X_i - \sum_{1 \le i_1 < i_2 \le n} X_{i_1}X_{i_2} + \cdots \\+ (-1)^{m+1}\sum_{1 \le i_1 < \cdots < i_m \le n} X_{i_1}\cdots X_{i_m} \cdots +(-1)^n X_1\cdots X_n Ii=1nAi=1i=1n(1Xi)=i=1nXi1i1<i2nXi1Xi2++(1)m+11i1<<imnXi1Xim+(1)nX1Xn 基于这个引理把 ( 1 − X j 1 ) ⋯ ( 1 − X j n − m ) (1-X_{j_1})\cdots (1-X_{j_{n-m}}) (1Xj1)(1Xjnm)展开,然后把期望改写成概率就可以得到Waring公式了。

    Processed: 0.023, SQL: 9