OFFSET
0,5
COMMENTS
Euler transform of A048050.
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..2000
FORMULA
G.f.: Product_{k>=1} 1 / (1 - x^k)^A048050(k).
a(n) ~ exp(3^(2/3) * ((Pi^2 - 6)*Zeta(3))^(1/3) * n^(2/3)/2 - Pi^2 * (3/((Pi^2 - 6)*Zeta(3)))^(1/3) * n^(1/3)/4 - Pi^4 / (32*(Pi^2 - 6)*Zeta(3)) - 1/8) * A^(3/2)* (2*Pi)^(1/24) / (3^(1/8) * ((Pi^2 - 6)*Zeta(3))^(3/8) * n^(1/8)), where A is the Glaisher-Kinkelin constant A074962. - Vaclav Kotesovec, Oct 26 2019
MAPLE
with(numtheory):
b:= proc(n) option remember; `if`(n<4, 0, sigma(n)-1-n) end:
a:= proc(n) option remember; `if`(n=0, 1, add(add(
d*b(d), d=divisors(j))*a(n-j), j=1..n)/n)
end:
seq(a(n), n=0..50); # Alois P. Heinz, Oct 20 2019
MATHEMATICA
nmax = 47; CoefficientList[Series[Product[1/(1 - x^k)^(DivisorSigma[1, k] - k - 1), {k, 2, nmax}], {x, 0, nmax}], x]
a[n_] := a[n] = If[n == 0, 1, Sum[Sum[If[d == 1, 0, d (DivisorSigma[1, d] - d - 1)], {d, Divisors[k]}] a[n - k], {k, 1, n}]/n]; Table[a[n], {n, 0, 47}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Oct 20 2019
STATUS
approved