20/09/29 21:29:02.84 a1+Kawn1.net
>>815
やってみた。
> # 総当たり
> dice4 <- function(num, N=6, digit = 4){
+ r=num%%N
+ q=num%/%N
+ while(q > 0 | digit > 1){
+ r=append(q%%N,r)
+ q=q%/%N
+ digit=digit-1
+ }
+ return(r+1)
+ }
>
> re=0
> for(i in 0:(6^4-1)){
+ d=dice4(i)
+ re <- re + (sum(d==6)>=2) # 2個以上6の目があるか？
+ }
> re
[1] 171
> 6^4
[1] 1296
> 171/1296
[1] 0.1319444
>
> # 1億回シミュレーション
> mean(rbinom(1e8,4,1/6)>=2)
[1] 0.1319338