OFFSET
1,4
FORMULA
a(n) = Sum_{p|n, p prime} binomial(n/p+p-2, p-1).
MATHEMATICA
nmax = 65; CoefficientList[Series[Sum[x^Prime[k]/(1 - x^Prime[k])^Prime[k], {k, 1, nmax}], {x, 0, nmax}], x] // Rest
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Mar 30 2025
STATUS
approved
