UA MATH564 概率论 QE练习题 信封问题

  • 2015年1月的第二题
  • 2015年5月的第一题

这一篇介绍QE理论中出现了两个信封问题相关的题目:2015年1月的第二题和2015年5月的第一题。信封问题最原始的设定是有nnn个信封和nnn个信件,把nnn个信件随机装入nnn个信封中。关于这个问题有个非常重要的结论,用PmP_mPm​表示恰好有mmm个信件被放入了正确的信封的概率,则
Pm=1m!∑j=0n−m(−1)jj!P_m = \frac{1}{m!}\sum_{j=0}^{n-m} \frac{(-1)^j}{j!}Pm​=m!1​j=0∑n−m​j!(−1)j​

这个式子的推导也非常容易,根据Waring公式(参考UA MATH564 概率论 计算至少有一个发生的概率:Waring公式)
Pm=∑k=mn(−1)k−mCkmSkP_m= \sum_{k=m}^n (-1)^{k-m}C_k^mS_kPm​=k=m∑n​(−1)k−mCkm​Sk​

假设事件AiA_iAi​表示第iii个信件被装在正确的信封中,则记号SkS_kSk​表示
Sk=∑1≤i1<i2<⋯<ik≤nP(Ai1∩Ai2∩⋯∩Aik)=Cnk(n−k)!n!=1k!S_k = \sum_{1 \le i_1 < i_2 < \cdots < i_k \le n} P(A_{i_1} \cap A_{i_2} \cap \cdots \cap A_{i_k}) = C_n^k \frac{(n-k)!}{n!} = \frac{1}{k!}Sk​=1≤i1​<i2​<⋯<ik​≤n∑​P(Ai1​​∩Ai2​​∩⋯∩Aik​​)=Cnk​n!(n−k)!​=k!1​

这个概率的计算比较简单,从nnn个信件中抽kkk个,这kkk个的位置是确定的,所以基本事件数目等于另外(n−k)(n-k)(n−k)个的全排列,基本事件总数等于nnn的全排列。因此
Pm=∑k=mn(−1)k−mCkmSk=∑k=mn(−1)k−mk!m!(k−m)!1k!=1m!∑j=0n−m(−1)jj!P_m= \sum_{k=m}^n (-1)^{k-m}C_k^mS_k = \sum_{k=m}^n (-1)^{k-m}\frac{k!}{m!(k-m)!}\frac{1}{k!} = \frac{1}{m!}\sum_{j=0}^{n-m} \frac{(-1)^j}{j!}Pm​=k=m∑n​(−1)k−mCkm​Sk​=k=m∑n​(−1)k−mm!(k−m)!k!​k!1​=m!1​j=0∑n−m​j!(−1)j​

2015年1月的第二题


Part a
这一问比较简单,不限制每个袋子容纳的球的数目,则基本事件总数为nnn^nnn,每个球都不在对应编号的袋子里的基本事件数目为(n−1)n(n-1)^n(n−1)n,因此概率为
nn−(n−1)nnn\frac{n^n - (n-1)^n}{n^n}nnnn−(n−1)n​

Part b
这一问就是典型的信封问题,每个袋子只能容纳一个球,每个球都不在对应编号的袋子里的概率,根据信封问题的公式,就是
P0=∑j=0n(−1)jj!P_0 = \sum_{j=0}^{n} \frac{(-1)^j}{j!}P0​=j=0∑n​j!(−1)j​

但!考试的时候是不能这样写的。这里提供两种推导的思路,分别从概率论与组合学的角度入手,概率论方法的核心是用示性函数的期望表示事件的概率,把复杂概率化归为简单概率计算;组合学方法的核心是容斥原理。

概率论方法
引入Bernoulli随机变量XiX_iXi​,Xi=1X_i = 1Xi​=1表示编号为iii的球被放在对应编号的袋子中,则
P0=E[∏i=1n(1−Xi)]=E[1−∑i=1nXi+(−1)m∑1≤i1<⋯<im≤nXi1⋯Xim⋯+(−1)nX1⋯Xn]P_0 = E \left[ \prod_{i=1}^n (1-X_i) \right] \\= E \left[ 1- \sum_{i=1}^n X_i + (-1)^{m}\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 \right]P0​=E[i=1∏n​(1−Xi​)]=E[1−i=1∑n​Xi​+(−1)m1≤i1​<⋯<im​≤n∑​Xi1​​⋯Xim​​⋯+(−1)nX1​⋯Xn​]

假设事件AiA_iAi​表示编号为iii的球被装在对应编号的袋子中
E[Xi1⋯Xim]=P(Ai1∩Ai2∩⋯∩Aim)=(n−m)!n!E[X_{i_1}\cdots X_{i_m}] = P(A_{i_1} \cap A_{i_2} \cap \cdots \cap A_{i_m}) = \frac{(n-m)!}{n!}E[Xi1​​⋯Xim​​]=P(Ai1​​∩Ai2​​∩⋯∩Aim​​)=n!(n−m)!​

带入可得
P0=∑j=0n(−1)jj!P_0 = \sum_{j=0}^{n} \frac{(-1)^j}{j!}P0​=j=0∑n​j!(−1)j​

组合学方法
沿用上述记号,根据容斥原理,至少有一个球在对应编号的袋子里的基本事件总数为
#⋃i=1nAi=∑i=1n#Ai+(−1)m+1∑1≤i1<⋯<im≤n#⋂j=1mAij+(−1)n+1#⋂i=1nAi\# \bigcup_{i=1}^n A_i = \sum_{i=1}^n \#A_i +(-1)^{m+1} \sum_{1 \le i_1 < \cdots < i_m \le n}\# \bigcap_{j=1}^m A_{i_j} + (-1)^{n+1}\#\bigcap_{i=1}^n A_{i} #i=1⋃n​Ai​=i=1∑n​#Ai​+(−1)m+11≤i1​<⋯<im​≤n∑​#j=1⋂m​Aij​​+(−1)n+1#i=1⋂n​Ai​

其中
#⋂j=1mAij=(n−m)!\# \bigcap_{j=1}^m A_{i_j} = (n-m)!#j=1⋂m​Aij​​=(n−m)!

基本事件总数为n!n!n!,因此概率为
P0=n!−#⋃i=1nAin!=∑j=0n(−1)jj!P_0 = \frac{n! - \# \bigcup_{i=1}^n A_i}{n!} = \sum_{j=0}^{n} \frac{(-1)^j}{j!}P0​=n!n!−#⋃i=1n​Ai​​=j=0∑n​j!(−1)j​

2015年5月的第一题


这道题其实就是对上一题中概率论方法定义的Bernoulli随机变量的深入讨论。注意到Xi∼Ber(1/n)X_i \sim Ber(1/n)Xi​∼Ber(1/n),因此
E[Xi2]=12×1n+02×n−1n=1nE[X_i^2] = 1^2 \times \frac{1}{n} + 0^2 \times \frac{n-1}{n} = \frac{1}{n}E[Xi2​]=12×n1​+02×nn−1​=n1​

XiXjX_iX_jXi​Xj​也服从Bernoulli分布,
P(XiXj=1)=P(Xi=1,Xj=1)=(n−2)!n!=1n(n−1)E[X1Xj]=1n(n−1)P(X_iX_j = 1) = P(X_i = 1,X_j = 1) = \frac{(n-2)!}{n!} = \frac{1}{n(n-1)} \\ E[X_1X_j] = \frac{1}{n(n-1)}P(Xi​Xj​=1)=P(Xi​=1,Xj​=1)=n!(n−2)!​=n(n−1)1​E[X1​Xj​]=n(n−1)1​

有了这两个结论,剩下两个计算就十分容易了
E[Sn2]=E(∑i=1nXi)2=∑i=1nEXi2+2∑1≤i<j≤nEXiXj=1+2×Cn21n(n−1)=2E[Sn]=∑i=1nEXi=1Var(Sn)=E[Sn2]−(E[Sn])2=1E[S_n^2] = E \left( \sum_{i=1}^n X_i \right)^2 = \sum_{i=1}^n EX_i^2 + 2\sum_{1 \le i < j \le n} EX_iX_j = 1 + 2 \times C_n^2 \frac{1}{n(n-1)} = 2 \\ E[S_n] = \sum_{i=1}^n EX_i = 1 \\ Var(S_n) = E[S_n^2] - (E[S_n])^2 = 1E[Sn2​]=E(i=1∑n​Xi​)2=i=1∑n​EXi2​+21≤i<j≤n∑​EXi​Xj​=1+2×Cn2​n(n−1)1​=2E[Sn​]=i=1∑n​EXi​=1Var(Sn​)=E[Sn2​]−(E[Sn​])2=1

UA MATH564 概率论 QE练习题 信封问题相关推荐

  1. UA MATH564 概率论 QE练习题6

    UA MATH564 概率论 QE练习题6 第二题 第三题 这是2016年5月的2-3题. 第二题 Part a Denote O as center of the circle. Randomly ...

  2. UA MATH564 概率论 QE练习题 概率极限理论

    UA MATH564 概率论 QE练习题 概率极限理论 2015/5/3 2016/1/3 这是2015年5月的3题.2016年1月的3题 2015/5/3 这个题干有点意思,有一列随机变量但并不是互 ...

  3. UA MATH564 概率论 QE练习题3

    UA MATH564 概率论 QE练习题3 第一题 第二题 第三题 这是2015年1月的1-3题. 第一题 Part a Obvioulsy, EU=EW=0EU = EW = 0EU=EW=0, V ...

  4. UA MATH564 概率论 QE练习题2

    UA MATH564 概率论 QE练习题2 第一题 第二题 第三题 这是2014年5月理论的1-3题. 第一题 For X2X_2X2​, P(X2=0∣X1=0)=b+1b+r+1P(X2=1∣X1 ...

  5. UA MATH564 概率论 QE练习题1

    UA MATH564 概率论 QE练习题1 第一题 第二题 第三题 2014年1月理论的1-3题. 第一题 Part (a) Marginal density of YYY is fY(y)=∫−∞∞ ...

  6. UA MATH564 概率论 QE练习 Glivenko–Cantelli定理

    UA MATH564 概率论 QE练习 Glivenko–Cantelli定理 Glivenko–Cantelli定理 Part a Notice P(Xi≤t)=F(t),i=1,⋯,nP(X_i ...

  7. UA MATH564 概率论VI 数理统计基础3 卡方分布的正态近似

    UA MATH564 概率论VI 数理统计基础3 卡方分布的正态近似 卡方分布的正态近似 Gamma分布的正态近似 基于卡方分布的近似 基于指数分布的近似 在做UA MATH566 统计理论 QE练习 ...

  8. UA MATH564 概率论 多元随机变量的变换 理论与应用2

    UA MATH564 概率论 多元随机变量的变换 几个例题 例5 X1,X2,X3∼iidEXP(λ)X_1,X_2,X_3 \sim_{iid} EXP(\lambda)X1​,X2​,X3​∼ii ...

  9. UA MATH564 概率论 依概率收敛的一个例题

    UA MATH564 概率论 依概率收敛的一个例题 Part (a) Let Yn∼U(−1/n,1/n)Y_n \sim U(-1/n,1/n)Yn​∼U(−1/n,1/n), Zn∼U(n,n+1 ...

最新文章

  1. 12岁AI开发者现身DuerOS发布会:得开发者得天下
  2. Mrtg搭建流量监控服务器
  3. 程序员面试题精选100题(39)-颠倒栈[数据结构]
  4. 实时数仓入门训练营:Hologres性能调优实践
  5. 吵架后女生和男生的夜晚!所有男生都这样吗?
  6. P1445-[Violet]樱花【数学】
  7. 微软dns能做cname吗_为什么域的根不能是CNAME以及有关DNS的其他花絮
  8. SVG 教程 (二)矩形
  9. C++基础部分_C++文件操作_二进制文件的写操作---C++语言工作笔记078
  10. EnterCriticalSection 多线程操作相同数据遇到的问题(线程锁)
  11. 学会拒绝别人的6个技巧_4个秘诀教你怎么拒绝别人,还能不伤和气
  12. thinkphp5 mysql加1_ThinkPHP5.1的数据库链接和增删改查
  13. Java for Web学习笔记(十六):JSP(6)jspx
  14. windows驱动备份还原-(驱动软件无法更新)
  15. PetaLinux 添加启动后自动执行脚本
  16. C#使用libVLC制作视频播放器时自定义鼠标事件的捕获或者忽略
  17. PDF转图片怎样转换?有哪些好用PDF转换的方法?
  18. 时空跳跃者的追捕行动模拟赛
  19. 大数据高级开发工程师——Spark学习笔记(10)
  20. 最好的远程桌面破解软件frdpb_v2

热门文章

  1. 算法导论之计算几何学
  2. Hive查询结果输出文件
  3. 虚拟光驱安装服务器无法运行,windows7虚拟光驱无法正常打开怎么办
  4. Rest 微服务工程搭建03——工程重构
  5. Move from SQL to SOQL
  6. Manage Jenkins管理界面提示“依赖错误: 部分插件由于缺少依赖无法加载...“问题解决办法
  7. C语言实现测量数据处理,求出实际值的范围
  8. 模拟电路技术之基础知识(三)
  9. CTFshow 命令执行 web46
  10. 第三周实践项目8 以第一个元素为分界线,将所有小于它的元素移到该元素的前面,将所有大于它的元素移到该元素的后面。