22/04/22 09:02:07 qBIcbTPi.net
>>14
Σを使うことになるのでfor loop でプログラムしてみた(R言語ver4.1)
calc = \(n){
f=\(k){
re=0
for(i in 0:((k-n)%/%6)){
re=re+(-1)^i*choose(n,i)*choose(k-6*i-1,n-1)
}
re
}
j=n:(6*n)
p=sapply(j,f)/6^n
m=j/n
data.frame(j,m,p)
g=\(x) sum(1:6 >=x)
q=sapply(m,g)/6
sum(p*q)
}
n=1-30で実行
> sapply(1:30,calc)|> MASS::fractions()
[1] 7/12 13/24 19/36
[4] 3458/6641 12053/23328 37/72
[7] 286601/559872 64297/125970 8367/16430
[10] 8867/17444 529953/1044173 17845/35206
[13] 151219/298672 866532/1713169 24017/47524
[16] 3707/7341 139924/277287 64189/127284
[19] 641834/1273465 23774/47195 13626251/27063282
[22] 127525/253392 5137913/10213219 408529711/812390026
[25] 965461/1920556 4971/9892 35195/70053
[28] 4655713/9270612 232051/462151 112803/224545