OFFSET
0,9
LINKS
Eric Weisstein's World of Mathematics, Octagonal Number
FORMULA
G.f.: 1/(1 - Sum_{k>=1} x^(k*(3*k-2))).
MAPLE
h:= proc(n) option remember; `if`(n<1, 0, (t->
`if`(t*(3*t-2)>n, t-1, t))(1+h(n-1)))
end:
a:= proc(n) option remember; `if`(n=0, 1,
add(a(n-i*(3*i-2)), i=1..h(n)))
end:
seq(a(n), n=0..60); # Alois P. Heinz, Dec 28 2018
MATHEMATICA
nmax = 54; CoefficientList[Series[1/(1 - Sum[x^(k (3 k - 2)), {k, 1, nmax}]), {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Dec 26 2018
STATUS
approved