login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A141366 G.f.: A(x) = x - A(-A(x)^2). 1
1, 1, 2, 4, 8, 16, 32, 65, 138, 316, 792, 2142, 6052, 17316, 49160, 137108, 374632, 1004668, 2656872, 6974510, 18310220, 48416528, 129637152, 352339898, 971336204, 2707410220, 7596357720, 21363062480, 60019857952, 168132717792 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,3
COMMENTS
What is Limit a(n+1)/a(n) = ? ; a(501)/a(500) = 2.91525180...
The limit is equal to 2.9240038... - Vaclav Kotesovec, Jan 06 2023
LINKS
FORMULA
G.f.: A(x) = Series_Reversion( x + A(-x^2) ).
a(n) ~ c * d^n / n^(3/2), where d = 2.924003847... and c = 0.28352589... - Vaclav Kotesovec, Jan 06 2023
EXAMPLE
G.f.: A(x) = x + x^2 + 2*x^3 + 4*x^4 + 8*x^5 + 16*x^6 + 32*x^7 + 65*x^8 +...
A(x)^2 = x^2 + 2*x^3 + 5*x^4 + 12*x^5 + 28*x^6 + 64*x^7 + 144*x^8 +...
Series_Reversion( A(x) ) = x + A(-x^2), thus:
x = A(x - x^2 + x^4 - 2*x^6 + 4*x^8 - 8*x^10 + 16*x^12 -+...).
PROG
(PARI) {a(n)=local(A=x+x^2); for(i=0, n, A=x-subst(A, x, -(A+x*O(x^n))^2)); polcoeff(A, n)}
CROSSREFS
Sequence in context: A100137 A325917 A210542 * A049142 A367661 A303023
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jul 18 2008
STATUS
approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 25 01:41 EDT 2024. Contains 374585 sequences. (Running on oeis4.)