>384の神コードをRに移植。

f <- \(n, m, stack = list()) {
key = paste0(n, "_", m)
if (!is.null(stack[[key]])) {
return(stack[[key]])
}
if (n == 1 && m > 0 && m < 7) {
re = m
} else {
if (m < n || m > 6 * n) {
re = 0
} else {
re = 0
for (k in 1:6) {
tmp = k * f(n - 1, m - k, stack)
re = re + tmp
}
}
}
stack[[key]] = re
return(re)
}
calc <-\(k) f_r(10, k) / 21^10

> sapply(43:45,calc)
[1] 0.08256769 0.08356575 0.08100966
>384の　N[]の方の数値とほぼ同じ。

In[4]:= Table[{k,N[f[10,k]/21^10,6]},{k,43,45}]
Out[4]= {{43, 0.0825677}, {44, 0.0835657}, {45, 0.0810097}}
自分のスキルではRで分数解を正しく算出するのは無理でした。