■幾何分布と誕生日の問題(その42)

 k=365として,

  pn,k=Σ(1,k)(−1)^k-jkCjj^n/k^n

  pn,k=Σ(1,k)(−1)^k-jkCj(j/k)^n

  p>0.5,p>0.9

となるnを計算してみたところ,誤差が累積して計算不能であった.

 k=10の場合,

  p>0.5→n=27

  p>0.9→n=44

 kが小さいときは計算できるのであるが,何かうまい近似法はあるのだろうか?

===================================