A291558 O.g.f.: exp( Sum_{n>=1} A251658(n) * x^n/n ), where A251658(n) = Sum_{k=1..n} C(n,k) * sigma(k)^(n-k).

1, 1, 2, 6, 25, 123, 695, 4710, 38573, 363290, 3762248, 42712427, 544108640, 7858522102, 125629094752, 2132933782538, 37180929264503, 653233840242274, 11489698066482812, 202366696882463928, 3583045708114956712, 64151951314880747470, 1170351942994473380196, 21993337390128064959103, 433693240259690293897796, 9276431005281628486906887, 226051070328728855263836057

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

Cf. A251658.

Sequence in context: A358499 A357949 A321186 * A357902 A370510 A370511

Adjacent sequences: A291555 A291556 A291557 * A291559 A291560 A291561

Keyword

nonn

Author

Paul D. Hanna, Aug 26 2017

Status

approved