login
This site is supported by donations to The OEIS Foundation.

 

Logo


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

Table of n, a(n) for n=1..20.

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.

License Agreements, Terms of Use, Privacy Policy. .

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