OFFSET
0,11
COMMENTS
Convolution inverse of A010057.
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..10000
FORMULA
a(0) = 1; a(n) = -Sum_{k=1..n} A010057(k) * a(n-k). - Seiichi Manyama, Mar 19 2022
MAPLE
a:=series(1/add(x^(k^3), k=0..100), x=0, 59): seq(coeff(a, x, n), n=0..58); # Paolo P. Lava, Apr 02 2019
MATHEMATICA
nmax = 58; CoefficientList[Series[1/Sum[x^k^3, {k, 0, nmax}], {x, 0, nmax}], x]
a[0] = 1; a[n_] := a[n] = -Sum[Boole[IntegerQ[k^(1/3)]] a[n - k], {k, 1, n}]; Table[a[n], {n, 0, 58}]
PROG
(PARI) my(N=66, x='x+O('x^N)); Vec(1/sum(k=0, N^(1/3), x^k^3)) \\ Seiichi Manyama, Mar 19 2022
(PARI) a(n) = if(n==0, 1, -sum(k=1, n, ispower(k, 3)*a(n-k))); \\ Seiichi Manyama, Mar 19 2022
CROSSREFS
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Jan 21 2019
STATUS
approved