OFFSET
0,9
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000 (first 501 terms from T. D. Noe)
FORMULA
G.f.: 1 / (1 - Sum_{n>=1} x^(n^3) ). - Joerg Arndt, Mar 30 2014
a(n) ~ c * d^n, where d = 1.2338881403372741887535479..., c = 0.418031200641837887398653... - Vaclav Kotesovec, May 01 2014
MAPLE
a:= proc(n) option remember; `if`(n=0, 1,
`if`(n<0, 0, add(a(n-i^3), i=1..iroot(n, 3))))
end:
seq(a(n), n=0..80); # Alois P. Heinz, Sep 08 2014
MATHEMATICA
a[n_] := a[n] = If[n==0, 1, If[n<0, 0, Sum[a[n-i^3], {i, 1, Floor[n^(1/3)]}]]]; Table[a[n], {n, 0, 80}] (* Jean-François Alcover, Apr 08 2015, after Alois P. Heinz *)
PROG
(PARI) E=6; N=E^3-1; q='q+O('q^N);
gf=1/(1 - sum(n=1, E, q^(n^3) ) ); \\ test, several similar seqs.
v=Vec(gf) \\ Joerg Arndt, Mar 30 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved