很無聊用Mathematica寫了一下：
(*真正的機率*)
f[n_]:=1-Times@@Range[365-n+1,365]/365^n//N
(*實驗5000次得到的機率*)
g[n_]:=N@Total[If[SortBy[Tally@RandomInteger[{1,365},n],Last][[-1,2]]>=2,1,0]&/@Range[5000]]/5000
(*將模擬資料與真實機率輸出*)
TableForm[{f[#],g[#]}&/@Range[2,40],TableHeadings->{Range[2,40],{f,g}}]
讀者回應 ( 0 意見 )
訂閱發佈留言 (Atom)
發佈留言
Please leave your name and tell me what you thought about this site. Comments, suggestions and views are welcomed.
如果這篇文章對你有幫助，那請留個訊息給我～