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,2

Links

Table of n, a(n) for n=0..19.

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);

Crossrefs

R(n, 1) = A000012 / A000012.

R(n, 2) = A001790 / A046161.

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

Sequence in context: A368446 A110853 A371640 * A167723 A203107 A168493

Adjacent sequences: A344399 A344400 A344401 * A344403 A344404 A344405

Keyword

nonn

Author

Peter Luschny, May 20 2021

Status

approved