OFFSET
0,19
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..20000
Eric Weisstein's World of Mathematics, Heptagonal Number
FORMULA
G.f.: Product_{k>=1} (1 + x^(k*(5*k-3)/2))*(1 + x^(k*(5*k+3)/2)).
a(n) ~ zeta(3/2)^(1/3) * exp(3*Pi^(1/3) * zeta(3/2)^(2/3) * n^(1/3) / (2^(2/3) * (1 + sqrt(2))^(2/3) * 5^(1/3))) / (2^(4/3) * sqrt(3) * 5^(1/6) * (1 + sqrt(2))^(1/3) * Pi^(1/3) * n^(5/6)). - Vaclav Kotesovec, Mar 11 2026
EXAMPLE
a(18) = 2 because we have [18] and [13, 4, 1].
MATHEMATICA
nmax = 100; CoefficientList[Series[Product[(1 + x^(k (5 k - 3)/2)) (1 + x^(k (5 k + 3)/2)), {k, 1, Floor[Sqrt[9 + 40*nmax]/10] + 1}], {x, 0, nmax}], x] (* tuned for efficiency by Vaclav Kotesovec, Mar 11 2026 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 05 2017
STATUS
approved
