23/03/04 15:42:16.60 ekLrbL79H.net
> calc=\(n){
+ paste0("expand.grid(",paste0(rep("1:6",n),collapse=","),')') |> str2lang() |> eval() -> gr
+ apply(gr,1,\(x) prod(x)%%6==0) |> sum() -> nu
+ MASS::fractions(nu/6^n) |> print()
+ invisible(nu/6^n)
+ }
> calc(3)
[1] 133/216
> n=1:10
> pn=sapply(n,calc)
[1] 1/6
[1] 5/12
[1] 133/216
[1] 325/432
[1] 6541/7776
[1] 4655/5184
[1] 35120/37597
[1] 2417/2525
[1] 592/609
[1] 4935/5027
>