login
A370752
a(n) = 3^n * [x^n] Product_{k>=1} ((1 + 3*x^k)/(1 - 3*x^k))^(1/3).
2
1, 6, 36, 360, 1998, 18792, 121176, 1123632, 7537860, 72078174, 510702408, 4896308088, 35923749480, 345406994280, 2600934294816, 24985346997888, 191735328374478, 1838307293836560, 14317601666954364, 136953233511162840, 1079293961918593800, 10299943344889922832
OFFSET
0,2
COMMENTS
In general, if d > 1, m >= 1 and g.f. = Product_{k>=1} ((1 + d*x^k)/(1 - d*x^k))^(1/m), then a(n) ~ QPochhammer(-1, 1/d)^(1/m) * d^n / (Gamma(1/m) * QPochhammer(1/d)^(1/m) * n^(1 - 1/m)).
FORMULA
G.f.: Product_{k>=1} ((1 + 3*(3*x)^k)/(1 - 3*(3*x)^k))^(1/3).
a(n) ~ QPochhammer(-1, 1/3)^(1/3) * 9^n / (Gamma(1/3) * QPochhammer(1/3)^(1/3) * n^(2/3)).
MATHEMATICA
nmax = 30; CoefficientList[Series[Product[(1 + 3*x^k)/(1 - 3*x^k), {k, 1, nmax}]^(1/3), {x, 0, nmax}], x] * 3^Range[0, nmax]
nmax = 30; CoefficientList[Series[Product[(1 + 3*(3*x)^k)/(1 - 3*(3*x)^k), {k, 1, nmax}]^(1/3), {x, 0, nmax}], x]
CROSSREFS
Cf. A303390 (d=3,m=1), A370751 (d=3,m=2), A370752 (d=3,m=3).
Cf. A261584 (d=2,m=1), A303346 (d=2,m=2), A370750 (d=2,m=3), A370749 (d=2,m=4).
Cf. A015128 (d=1,m=1), A303307 (d=1,m=2), A303342 (d=1,m=3).
Sequence in context: A367490 A053335 A222925 * A108733 A350525 A229530
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Feb 29 2024
STATUS
approved