login
A368348
a(n) = [x^(n^4)] Product_{k=1..n} (x^(k^4) + 1/x^(k^4)).
0
1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 3, 0, 0, 19, 26, 0, 0, 40, 129, 0, 0, 616, 785, 0, 0, 4080, 9309, 0, 0, 44775, 72659, 0, 0, 430297, 781505, 0, 0, 3934457, 7765047, 0, 0, 44740433, 78818429, 0, 0, 463089552, 900950811, 0, 0, 5344766190, 9806206864, 0, 0
OFFSET
0,18
MAPLE
b:= proc(n, i) option remember; (m-> `if`(n>m, 0, `if`(n=m, 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-> `if`(irem(n, 4)>1, 0, b(n^4, n)):
seq(a(n), n=0..43); # Alois P. Heinz, Jan 25 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 25 2024
EXTENSIONS
a(46)-a(59) from Alois P. Heinz, Jan 25 2024
STATUS
approved