This site is supported by donations to The OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A227457 E.g.f. equals the series reversion of x - x*log(1+x). 0
 1, 2, 9, 68, 720, 9804, 163184, 3210192, 72870120, 1874721360, 53905894152, 1713195438624, 59633476003920, 2256257009704320, 92196226214092800, 4046446853549201664, 189845257963376620800, 9481546020840245199360, 502242773970728703225600, 28124368575613839072714240 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,2 LINKS FORMULA a(n) = Sum_{k=0..n-1} k! * Stirling1(n-1,k) * binomial(n+k-1,n-1). [From a formula in A052819 due to Vladimir Kruchinin] E.g.f. A(x) satisfies: (1) A(x - x*log(1+x)) = x. (2) A(x) = x/(1 - log(1+A(x))). (3) A(x) = x + Sum_{n>=1} d^(n-1)/dx^(n-1) x^n * log(1+x)^n / n!. (4) A(x) = x*exp( Sum_{n>=1} d^(n-1)/dx^(n-1) x^(n-1) * log(1+x)^n / n! ). a(n) ~ n^(n-1) * (1-c) / (c*sqrt(1+c) * exp(n) * (c-2+1/c)^n), where c = LambertW(1) = 0.5671432904... (see A030178). - Vaclav Kotesovec, Jan 10 2014 EXAMPLE E.g.f.: A(x) = x + 2*x^2/2! + 9*x^3/3! + 68*x^4/4! + 720*x^5/5! +... where A(x) = x/(1 - log(1+A(x))). The e.g.f. satisfies: (3) A(x) = x + x*log(1+x) + d/dx x^2*log(1+x)^2/2! + d^2/dx^2 x^3*log(1+x)^3/3! + d^3/dx^3 x^4*log(1+x)^4/4! +... (4) log(A(x)/x) = log(1+x) + d/dx x*log(1+x)^2/2! + d^2/dx^2 x^2*log(1+x)^3/3! + d^3/dx^3 x^3*log(1+x)^4/4! +... MATHEMATICA Rest[CoefficientList[InverseSeries[Series[x - x*Log[1+x], {x, 0, 20}], x], x] * Range[0, 20]!] (* Vaclav Kotesovec, Jan 10 2014 *) PROG (PARI) {a(n)=n!*polcoeff(serreverse(x-x*log(1+x +x*O(x^n))), n)} for(n=1, 25, print1(a(n), ", ")) (PARI) {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D} {a(n)=local(A=x); A=x+sum(m=1, n, Dx(m-1, x^m*log(1+x+x*O(x^n))^m/m!)); n!*polcoeff(A, n)} for(n=1, 25, print1(a(n), ", ")) (PARI) {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D} {a(n)=local(A=x+x^2+x*O(x^n)); A=x*exp(sum(m=1, n, Dx(m-1, x^(m-1)*log(1+x+x*O(x^n))^m/m!)+x*O(x^n))); n!*polcoeff(A, n)} for(n=1, 25, print1(a(n), ", ")) (PARI) {Stirling1(n, k)=n!*polcoeff(binomial(x, n), k)} {a(n)=sum(k=0, n-1, k!*Stirling1(n-1, k)*binomial(n+k-1, n-1))} for(n=1, 25, print1(a(n), ", ")) CROSSREFS Cf. A052819, A030178. Sequence in context: A120020 A200248 A180747 * A134200 A217066 A038037 Adjacent sequences:  A227454 A227455 A227456 * A227458 A227459 A227460 KEYWORD nonn AUTHOR Paul D. Hanna, Jul 12 2013 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified April 19 06:30 EDT 2019. Contains 322237 sequences. (Running on oeis4.)