OFFSET
0,5
FORMULA
G.f.: (1 + x^2 / (1 - x)) * Product_{p prime, k>=1} (1 + x^(2*p^k) / (1 - x^(p^k))).
EXAMPLE
a(6) = 4 because we have [3, 3], [2, 2, 2], [2, 2, 1, 1] and [1, 1, 1, 1, 1, 1].
MATHEMATICA
nmax = 56; CoefficientList[Series[(1 + x^2/(1 - x)) 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