OFFSET
1,4
COMMENTS
Inverse Moebius transform of A033461.
MATHEMATICA
nmax = 90; CoefficientList[Series[Sum[-1 + Product[(1 + x^(k j^2)), {j, 1, Floor[nmax^(1/2)] + 1}], {k, 1, nmax}], {x, 0, nmax}], x] // Rest
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 13 2019
STATUS
approved