OFFSET
1,1
COMMENTS
Inverse Möbius transform of A015128.
LINKS
N. J. A. Sloane, Transforms
FORMULA
MAPLE
a:=series(add(-1+mul((1+x^(k*j))/(1-x^(k*j)), j=1..100), k=1..100), x=0, 44): seq(coeff(a, x, n), n=1..43); # Paolo P. Lava, Apr 02 2019
MATHEMATICA
nmax = 43; Rest[CoefficientList[Series[Sum[-1 + Product[(1 + x^(k j))/(1 - x^(k j)), {j, 1, nmax}], {k, 1, nmax}], {x, 0, nmax}], x]]
nmax = 43; Rest[CoefficientList[Series[Sum[1/EllipticTheta[4, 0, x^k] - 1, {k, 1, nmax}], {x, 0, nmax}], x]]
Table[Sum[Sum[PartitionsP[d - k] PartitionsQ[k], {k, 0, d}], {d, Divisors[n]}], {n, 43}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Oct 24 2018
STATUS
approved