OFFSET
0,3
COMMENTS
a(n)/n! is the Euler transform of [1, 3/2, 4/3, 7/4, 6/5, ... = sums of reciprocals of divisors of 1, 2, 3, 4, 5, ...].
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..438
Lida Ahmadi, Ricardo Gómez Aíza, and Mark Daniel Ward, A unified treatment of families of partition functions, La Matematica (2024). Preprint available as arXiv:2303.02240 [math.CO], 2023.
N. J. A. Sloane, Transforms
FORMULA
E.g.f.: exp(Sum_{k>=1} Sum_{j>=1} x^(j*k)/(j*k*(1 - x^(j*k)))).
log(a(n)/n!) ~ sqrt(n) * Pi^2 / 3. - Vaclav Kotesovec, Sep 04 2018
MAPLE
with(numtheory): a := proc(n) option remember; `if`(n = 0, 1, add(add(sigma(d), d = divisors(j))*a(n-j), j = 1..n)/n) end proc; seq(n!*a(n), n = 0..20); # Vaclav Kotesovec, Sep 04 2018
MATHEMATICA
nmax = 21; CoefficientList[Series[Product[1/(1 - x^k)^(DivisorSigma[1, k]/k), {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]!
nmax = 21; CoefficientList[Series[Exp[Sum[Sum[x^(j k)/(j k (1 - x^(j k))), {j, 1, nmax}], {k, 1, nmax}]], {x, 0, nmax}], x] Range[0, nmax]!
a[n_] := a[n] = If[n == 0, 1, Sum[Sum[d DivisorSigma[-1, d], {d, Divisors[k]}] a[n - k], {k, 1, n}]/n]; Table[n! a[n], {n, 0, 21}]
nmax = 21; s = 1 - x; Do[s *= Sum[Binomial[DivisorSigma[1, k]/k, j]*(-1)^j*x^(j*k), {j, 0, nmax/k}]; s = Expand[s]; s = Take[s, Min[nmax + 1, Exponent[s, x] + 1, Length[s]]]; , {k, 2, nmax}]; CoefficientList[Series[1/s, {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, Sep 03 2018 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, May 26 2018
STATUS
approved