20/11/16 19:01:49.96 wldb2lUp.net
>>484
プログラムを一般化して、サイコロをふる回数を１～６回にすると
fn <- function(n=3){
pm=expand.grid(replicate(n,1:6,simplify = FALSE))
f <- function(x){
P=numeric(n)
for(i in 1:n) P[i]=sum(x[1:i])%%5
length(unique(P))==n
}
num=sum(apply(pm,1,f))
den=nrow(pm)
MASS::fractions(num/den)
}
lapply(1:6,fn)
> lapply(1:6,fn)
[[1]]
[1] 1
[[2]]
[1] 5/6
[[3]]
[1] 19/36
[[4]]
[1] 49/216
[[5]]
[1] 65/1296
[[6]]
[1] 0