OFFSET
10,2
FORMULA
G.f.: ( x + Sum_{k>=1} x^prime(k) )^10.
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, 11)
end:
a:= n-> coeff(b(n), x, 10):
seq(a(n), n=10..45); # Alois P. Heinz, Feb 25 2021
MATHEMATICA
nmax = 45; CoefficientList[Series[(x + Sum[x^Prime[k], {k, 1, nmax}])^10, {x, 0, nmax}], x] // Drop[#, 10] &
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Feb 24 2021
STATUS
approved