OFFSET
5,2
FORMULA
G.f.: ( x + Sum_{k>=1} x^prime(k) )^5.
MAPLE
b:= proc(n) option remember; series(`if`(n=0, 1, x*add(
`if`(j=1 or isprime(j), b(n-j), 0), j=1..n)), x, 6)
end:
a:= n-> coeff(b(n), x, 5):
seq(a(n), n=5..56); # Alois P. Heinz, Feb 24 2021
MATHEMATICA
nmax = 56; CoefficientList[Series[(x + Sum[x^Prime[k], {k, 1, nmax}])^5, {x, 0, nmax}], x] // Drop[#, 5] &
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Feb 24 2021
STATUS
approved