OFFSET
0,22
LINKS
Eric Weisstein's World of Mathematics, Octagonal Number
FORMULA
G.f.: Product_{k>=1} (1 + x^(k*(3*k-2)))*(1 + x^(k*(3*k+2))).
EXAMPLE
a(21) = 2 because we have [21] and [16, 5].
MATHEMATICA
nmax = 100; CoefficientList[Series[Product[(1 + x^(k (3 k - 2))) (1 + x^(k (3 k + 2))), {k, 1, nmax}], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 05 2017
STATUS
approved