20/03/15 04:25:33.20 OTl1KJku.net
>>635
ｎ を 1～２０でシミュレーションしてみた。
sim <- function(n,k=1e5){ # n:サイコロを振る回数 k:シミュレーション回数
sub <- function(n){
prod(sample(6,n,replace=TRUE))%%6==0 # n回の目の積の6で除算した剰余が0か？
}
mean(replicate(k,sub(n))) # 0となる割合を返す
}
p=sapply(1:20,function(n) sim(n))
data.frame(p)
p
1 0.16560
2 0.41369
3 0.61413
4 0.75417
5 0.84134
6 0.89669
7 0.93356
8 0.95765
9 0.97254
10 0.98093
11 0.98873
12 0.99175
13 0.99458
14 0.99611
15 0.99737
16 0.99849
17 0.99899
18 0.99924
19 0.99963
20 0.99978