OFFSET
1,4
FORMULA
G.f.: Sum_{k>=1} prime(k)^3 * x^(prime(k)^2) / (1 - x^prime(k)).
MATHEMATICA
Table[DivisorSum[n, #^3 &, # <= Sqrt[n] && PrimeQ[#] &], {n, 1, 75}]
nmax = 75; CoefficientList[Series[Sum[Prime[k]^3 x^(Prime[k]^2)/(1 - x^Prime[k]), {k, 1, nmax}], {x, 0, nmax}], x] // Rest
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Aug 20 2021
STATUS
approved