login
A368478
a(n) = [x^n] Product_{k=1..n} (x^(k^4) + 1 + 1/x^(k^4)).
1
1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 22, 88, 236, 568, 1147, 1999, 3030, 4333, 4793, 5620, 4628, 9353, 40048, 312436, 1638550, 7894061, 30866741, 105342660, 323671449, 863529149, 2152428536, 4750103818, 9796453345, 18572505806, 32181800057, 55712628236
OFFSET
0,13
MAPLE
b:= proc(n, i) option remember; (m-> `if`(n>m, 0, `if`(n=m, 1, b(n, i-1)+
b(abs(n-i^4), i-1)+b(n+i^4, i-1))))(i*(i+1)*(2*i+1)*(3*i^2+3*i-1)/30)
end:
a:= n-> b(n$2):
seq(a(n), n=0..33); # Alois P. Heinz, Jan 25 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 25 2024
EXTENSIONS
a(34)-a(38) from Alois P. Heinz, Jan 25 2024
STATUS
approved