"生日悖论"补充

从排列与组合的python实现到”生日问题”的解释

首先定义我们的记号（notation），用整数 $i=1,2,\ldots,k$ 来对房间里的人进行编号，其中 $k$是房间里的总人数。不考虑闰年，所有的年份都是 $n=365$，令 $b_i$表示第 $i$个人的生日。我们假设每个人的生日是独立的（这是一个重要假设），则 $i$ 和 $j$ 生日落在同一天 $r$ 的概率为：

$$\mathrm {Pr}\{b_i=r,b_j=r\}=\mathrm {Pr}\{b_i=r\}\mathrm {Pr}\{b_j=r\}=\frac1{n^2}$$

他们的生日落在同一天的概率为：

$$\sum_{r=1}^n\mathrm{Pr}\{b_i=r,b_j=r\}=\sum_{r=1}^n\mathrm {Pr}\{b_i=r\}\mathrm {Pr}\{b_j=r\}=n\cdot\frac1{n^2}=\frac1n$$

所以这有一个结论，如果两个人生日彼此独立的话，一个人生日在某一天的概率为 $\frac1n$，两个人生日在同一天的概率也为 $\frac1n$。

再来看，至少两人生日在同一天的概率：

$$\begin{split} 1-\frac{A_n^k}{n^k}=&1-\frac{\frac{n!}{(n-k)!}}{n^k}\\ =&1-\frac{n(n-1)\cdots(n-k+1)}{n^k}\\ =&1-1\cdot\frac{n-1}n\frac{n-2}n\cdots\frac{n-k+1}n \end{split}$$

又由 $e^x$的泰勒展开为：

$$e^x=\sum_{n=0}^\infty\frac{x^n}{n!}=1+x+\frac{x^2}2+\cdots$$
也即
$$1+x\leq e^x$$

所以如果要求至少有两人生日在同一天的概率 $\geq 1$，也即 $1-\frac{A_n^k}{n^k}\geq \frac12$，也即：

$$1\cdot\frac{n-1}n\frac{n-2}n\cdots\frac{n-k+1}n\leq \frac12$$
又由于 $1+x\leq e^x$，所以，
$$1\cdot\frac{n-1}n\frac{n-2}n\cdots\frac{n-k+1}n\leq e^{-1/n}e^{-2/n}\cdots e^{-(k-1)/n}=e^{-k(k-1)/2n}\leq 1/2$$
求解得 $k\geq (1+\sqrt{1+(8\ln 2)n})/2$，当 $n=365$时，必须有 $k\geq 23$。

>>> from math import sqrt, log
>>> (1+sqrt(1+8*log(2)*365))/2

22.99994315123396

也即，如果至少有23个人主宰一个房间里，那么至少有两人生日相同的概率至少为$1/2$。

在火星上，一年有669个火星日，也即 $n=669$，

>>> (1+sqrt(1+8*log(2)*669))/2
30.957854940707936

则在火星上要达到同样的效果，必须有31个火星人。

利用指示器随机变量进行近似分析

利用指示器随机变量（Indicator Random Variable），可以给生日悖论一个简单而近似的分析。对房间中 $k$个人中的每一对 $(i,j)$，$1\leq i<j\leq k$，定义指示器随机变量 $X_{ij}$如下：

$$X_{ij}=\mathbf{1}_{i,j生日相同}=\left \{ \begin{array}{ll} 1,&如果\;i,j\; 生日相同\\ 0,&否则 \end{array} \right.$$
又知两个人生日在同一天的概率为 $1/n$，所以有：
$$E[X_{ij}]=\mathrm{Pr}\{i,j\;生日相同\}=\frac1n$$
令 $X$表示计数具有相同生日的两人对出现的数目的随机变量：
$$X=\sum_{i=1}^k\sum_{j=i+1}^kX_{ij}$$
对两边取期望并应用 期望的线性性质，得到：
$$E(X)=\sum_{i=1}^k\sum_{j=i+1}^kE(X_{ij})=\binom k2\frac1n$$

因此当 $k(k-1)\geq 2n$时，有相同数目的两人对的期望数目大于1。当 $n=365$时，$k=28$。