A344402 a(n) = denominator(R(n,3)), where R(n,d) = (Product_{j prime to d} Pochhammer(j/d, n)) / n!.

1, 9, 81, 2187, 19683, 177147, 4782969, 43046721, 387420489, 31381059609, 282429536481, 2541865828329, 68630377364883, 617673396283947, 5559060566555523, 150094635296999121, 1350851717672992089, 12157665459056928801, 984770902183611232881, 8862938119652501095929

Offset: 0

Author

Peter Luschny, May 20 2021

Keywords

• nonn

Crossrefs

R(n, 1) = A000012 / A000012.

R(n, 2) = A001790 / A046161.

R(n, 3) = A273194 / (this sequence).

Programs

Maple

coprimes := n -> select(j -> igcd(j, n) = 1, {$1..n}):
R := (n, d) -> mul(pochhammer(j/d, n), j in coprimes(d)) / n!:
seq(denom(R(n, 3)), n = 0..16);