login
A281612
Expansion of Sum_{i = p*q, p prime, q prime} x^i/(1 - x^i) / Product_{j>=1} (1 - x^j).
0
0, 0, 0, 1, 1, 3, 4, 8, 12, 20, 28, 45, 62, 92, 127, 181, 244, 340, 452, 614, 809, 1077, 1401, 1841, 2371, 3071, 3923, 5026, 6363, 8078, 10149, 12769, 15939, 19899, 24676, 30604, 37726, 46489, 57007, 69849, 85211, 103871, 126119, 152987, 184955, 223349, 268898, 323384, 387830, 464587, 555168, 662619, 789084
OFFSET
1,6
COMMENTS
Total number of semiprime parts (A001358) in all partitions of n.
Convolution of A000041 and A086971.
FORMULA
G.f.: Sum_{i = p*q, p prime, q prime} x^i/(1 - x^i) / Product_{j>=1} (1 - x^j).
EXAMPLE
a(6) = 3 because we have [6], [5, 1], [4, 2], [4, 1, 1], [3, 3], [3, 2, 1], [3, 1, 1, 1], [2, 2, 2], [2, 1, 1, 1, 1], [1, 1, 1, 1, 1, 1] and 1 + 0 + 1 + 1 + 0 + 0 + 0 + 0 + 0 + 0 + 0 = 3.
MATHEMATICA
nmax = 53; Rest[CoefficientList[Series[Sum[Floor[PrimeOmega[i]/2] Floor[2/PrimeOmega[i]] x^i/(1 - x^i), {i, 2, nmax}]/Product[1 - x^j, {j, 1, nmax}], {x, 0, nmax}], x]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jan 25 2017
STATUS
approved