OFFSET
1,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..400
MAPLE
s:= proc(n) s(n):= `if`(n<1, 0, ithprime(n)+s(n-1)) end:
b:= proc(n, i) option remember; `if`(n>s(i), 0, `if`(i=0, 1,
b(n, i-1)+b(n+ithprime(i), i-1)+b(abs(n-ithprime(i)), i-1)))
end:
a:= n-> b(ithprime(n), n):
seq(a(n), n=1..40); # Alois P. Heinz, Jan 22 2024
MATHEMATICA
Table[Coefficient[Product[(x^Prime[k] + 1 + 1/x^Prime[k]), {k, 1, n}], x, Prime[n]], {n, 1, 32}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 22 2024
STATUS
approved