OFFSET
0,17
FORMULA
a(n) = [x^n] Product_{k=1..n} (x^(k^3) + 1/x^(k^3)).
PROG
(Python)
from functools import lru_cache
@lru_cache(maxsize=None)
def b(n, i):
if n > (i*(i+1)//2)**2: return 0
if i == 0: return 1
return b(n+i**3, i-1) + b(abs(n-i**3), i-1)
def a(n): return b(n, n)
print([a(n) for n in range(54)]) # Michael S. Branicky, Jan 28 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 28 2022
STATUS
approved