A344402 - OEIS

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

(list; graph; refs; listen; history; text; internal format)

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.