login
A294842
Expansion of Product_{k>=1} (1 + x^k)^(k^2*(k+1)/2).
8
1, 1, 6, 24, 73, 238, 722, 2175, 6343, 18177, 50982, 140671, 382227, 1023623, 2706184, 7067324, 18250671, 46635309, 117997008, 295794098, 735030985, 1811435607, 4429226677, 10749552338, 25903858181, 62000039513, 147435739522, 348431110651, 818549931526, 1912010876019, 4441687009798
OFFSET
0,3
COMMENTS
Weigh transform of the pentagonal pyramidal numbers (A002411).
LINKS
M. Bernstein and N. J. A. Sloane, Some canonical sequences of integers, Linear Alg. Applications, 226-228 (1995), 57-72; erratum 320 (2000), 210. [Link to arXiv version]
M. Bernstein and N. J. A. Sloane, Some canonical sequences of integers, Linear Alg. Applications, 226-228 (1995), 57-72; erratum 320 (2000), 210. [Link to Lin. Alg. Applic. version together with omitted figures]
N. J. A. Sloane, Transforms
Eric Weisstein's World of Mathematics, Pentagonal Pyramidal Number
FORMULA
G.f.: Product_{k>=1} (1 + x^k)^A002411(k).
a(n) ~ exp(-2401 * Pi^16 / (2^12 * 3^11 * 5^8 * Zeta(5)^3) + (343 * Pi^12 / (2^(38/5) * 3^(37/5) * 5^(36/5) * Zeta(5)^(11/5))) * n^(1/5) - (49*Pi^8 / (2^(31/5) * 3^(24/5) * 5^(22/5) * Zeta(5)^(7/5))) * n^(2/5) + (7*Pi^4 / (2^(14/5) * 3^(16/5) * 5^(8/5) * Zeta(5)^(3/5))) * n^(3/5) + (5 * 3^(2/5) * (5*Zeta(5))^(1/5) / 2^(12/5)) * n^(4/5)) * 3^(1/5) * Zeta(5)^(1/10) / (2^(167/240) * 5^(2/5) * sqrt(Pi) * n^(3/5)). - Vaclav Kotesovec, Nov 10 2017
MATHEMATICA
nmax = 30; CoefficientList[Series[Product[(1 + x^k)^(k^2 (k + 1)/2), {k, 1, nmax}], {x, 0, nmax}], x]
a[n_] := a[n] = If[n == 0, 1, Sum[Sum[(-1)^(k/d + 1) d^3 (d + 1)/2, {d, Divisors[k]}] a[n - k], {k, 1, n}]/n]; Table[a[n], {n, 0, 30}]
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 09 2017
STATUS
approved