OFFSET
0,3
COMMENTS
Weigh transform of the hexagonal pyramidal numbers (A002412).
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, Hexagonal Pyramidal Number
FORMULA
G.f.: Product_{k>=1} (1 + x^k)^A002412(k).
a(n) ~ exp(-2401 * Pi^16 / (671846400000000 * Zeta(5)^3) - 49*Pi^8 * Zeta(3) / (518400000 * Zeta(5)^2) - Zeta(3)^2 / (2400*Zeta(5)) + (343 * Pi^12 / (77760000000 * 15^(1/5) * Zeta(5)^(11/5)) + 7*Pi^4*Zeta(3) / (72000 * 15^(1/5) * Zeta(5)^(6/5))) * n^(1/5) - (49*Pi^8 / (8640000 * 15^(2/5) * Zeta(5)^(7/5)) + Zeta(3) / (8 * (15*Zeta(5))^(2/5))) * n^(2/5) + (7*Pi^4 / (720 * (15*Zeta(5))^(3/5))) * n^(3/5) + (5*(15*Zeta(5))^(1/5)/4) * n^(4/5)) * (3*Zeta(5))^(1/10) / (2^(173/360) * 5^(2/5) * sqrt(Pi) * n^(3/5)). - Vaclav Kotesovec, Nov 10 2017
MATHEMATICA
nmax = 29; CoefficientList[Series[Product[(1 + x^k)^(k (k + 1)(4 k - 1)/6), {k, 1, nmax}], {x, 0, nmax}], x]
a[n_] := a[n] = If[n == 0, 1, Sum[Sum[(-1)^(k/d + 1) d^2 (d + 1)(4 d - 1)/6, {d, Divisors[k]}] a[n - k], {k, 1, n}]/n]; Table[a[n], {n, 0, 29}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 09 2017
STATUS
approved