OFFSET
1,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 1..450
FORMULA
E.g.f.: Sum_{k>=1} Sum_{j>=1} x^(j*k)/(j*k).
E.g.f.: -log(Product_{k>=1} (1 - x^k)^(1/k)).
E.g.f.: A(x) = log(B(x)), where B(x) = e.g.f. of A028342.
a(p^k) = (k + 1)*(p^k - 1)!, where p is a prime.
MATHEMATICA
Table[(n - 1)! DivisorSigma[0, n], {n, 1, 24}]
nmax = 24; Rest[CoefficientList[Series[Sum[Sum[x^(j k)/(j k), {j, 1, nmax}], {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]!]
nmax = 24; Rest[CoefficientList[Series[-Log[Product[(1 - x^k)^(1/k), {k, 1, nmax}]], {x, 0, nmax}], x] Range[0, nmax]!]
PROG
(PARI) a(n) = (n-1)!*numdiv(n); \\ Michel Marcus, Aug 22 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Aug 22 2018
STATUS
approved