OFFSET
0,6
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..10000
Eric Weisstein's World of Mathematics, Octagonal Number
FORMULA
G.f.: Product_{k>=1} 1/((1 - x^(k*(3*k-2)))*(1 - x^(k*(3*k+2)))).
a(n) ~ zeta(3/2) * exp(Pi^(1/3) * ((3/2)*zeta(3/2))^(2/3) * n^(1/3)) / (2*Pi*3^(3/2)*n^(3/2)). - Vaclav Kotesovec, Mar 11 2026
EXAMPLE
a(8) = 3 because we have [8], [5, 1, 1, 1] and [1, 1, 1, 1, 1, 1, 1, 1].
MATHEMATICA
nmax = 74; CoefficientList[Series[Product[1/((1 - x^(k (3 k - 2))) (1 - x^(k (3 k + 2)))), {k, 1, Floor[Sqrt[1 + 3*nmax]/3] + 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
