18/07/01 19:08:12.35 lSk1MYzX.net
３人以上で誕生月日が一致する確率が50％を超えるのは何人のときか？
シミュレーションしてみた。
date=1:366
p=c(97/400,rep(1,365)) # leap year
Birth <- function(n,lwr=2,upr=1e6,k=1e4){
same.birth <- function(n,lwr=2,upr=1e6){
x=sample(date,n,replace=TRUE,prob=p)
di=max(table(x))
lwr<=di & di<=upr
}
mean(replicate(k,same.birth(n,lwr,upr)))
}
vB=Vectorize(Birth)
x=1:100
y3=vB(x,3,k=1e4)
plot(x,y3,bty='l',pch=19,ylab='probability',xlab='subjects')
abline(h=0.5,lty=3,col=4)
min(x[which(y3 > 0.5)])
vB(87:88,3,k=1e5)
結果をグラフにすると
URLﾘﾝｸ(i.imgur.com)