OFFSET
1,8
COMMENTS
FORMULA
G.f.: Sum_{k>=1} x^(prime(k)^2)/(1 - x^(prime(k)^2)) / Product_{k>=1} (1 - x^(prime(k)^2)).
EXAMPLE
a(25) = 6 because we have [25], [9, 4, 4, 4, 4] and 1 + 5 = 6.
MATHEMATICA
nmax = 88; Rest[CoefficientList[Series[Sum[x^Prime[k]^2/(1 - x^Prime[k]^2), {k, 1, nmax}]/Product[1 - x^Prime[k]^2, {k, 1, nmax}], {x, 0, nmax}], x]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 27 2017
STATUS
approved