総当り解 R言語ver4.21

> pm=as.matrix(expand.grid(1:6,1:6,1:6))
> apply(pm,1,prod) |> table() -> p
> max(p)/6^3 |> MASS::fractions()
[1] 5/72
> as.numeric(names(p[p==7]))
[1] 8