OFFSET
0,3
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..200
FORMULA
O.g.f.: exp( Sum_{n>=1} x^n/(1 - sigma(n)*x)^n / n ). - Paul D. Hanna, Oct 26 2017
EXAMPLE
O.g.f: A(x) = 1 + x + 2*x^2 + 6*x^3 + 25*x^4 + 123*x^5 + 695*x^6 + 4710*x^7 + 38573*x^8 + 363290*x^9 + 3762248*x^10 + 42712427*x^11 + 544108640*x^12 + 7858522102*x^13 + 125629094752*x^14 + 2132933782538*x^15 +...
such that
log(A(x)) = x + 3*x^2/2 + 13*x^3/3 + 75*x^4/4 + 471*x^5/5 + 3273*x^6/6 + 26916*x^7/7 + 262027*x^8/8 + 2851375*x^9/9 +...+ A251658(n)*x^n/n +...
Also,
log(A(x)) = x/(1-x) + x^2/(1-3*x)^2/2 + x^3/(1-4*x)^3/3 + x^4/(1-7*x)^4/4 + x^5/(1-6*x)^5/5 + x^6/(1-12*x)^6/6 + x^7/(1-8*x)^7/7 + x^8/(1-15*x)^8/8 + x^9/(1-13*x)^9/9 + x^10/(1-18*x)^10/10 +...+ x^n/(1-sigma(n)*x)^n/n +...
PROG
(PARI) {A251658(n)=sum(k=1, n, binomial(n, k)*sigma(k)^(n-k))}
{a(n) = polcoeff( exp( sum(m=1, n+1, A251658(m)*x^m/m) +x*O(x^n) ), n)}
for(n=0, 30, print1(a(n), ", "))
(PARI) {a(n) = polcoeff( exp( sum(m=1, n+1, x^m/(1-sigma(m)*x+x*O(x^n))^m/m) ), n)}
for(n=0, 30, print1(a(n), ", ")) \\ - Paul D. Hanna, Oct 26 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Aug 26 2017
STATUS
approved