OFFSET
0,7
FORMULA
G.f.: Product_{p prime, k>=1} (1 + x^(2*p^k) / (1 - x^(p^k))).
EXAMPLE
a(12) = 5 because we have [4, 4, 4], [4, 4, 2, 2], [3, 3, 3, 3], [3, 3, 2, 2, 2] and [2, 2, 2, 2, 2, 2].
MATHEMATICA
nmax = 65; CoefficientList[Series[Product[1 + Boole[PrimePowerQ[k]] x^(2 k)/(1 - x^k), {k, 1, nmax}], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 28 2020
STATUS
approved