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!)
A107096 G.f. A(x) satisfies: A(x)^2 = x*G(x) where G(G(x)) = (1+x)^2*G(x) and G(x) is the g.f. of A107095. 5
0, 1, 1, -2, 11, -88, 869, -9876, 124473, -1701630, 24870695, -384795184, 6257294780, -106377162620, 1882982975521, -34593496243070, 657935674477431, -12927331575084846, 261951066040220637, -5466177185459699916, 117315664923801661485, -2586804284853871362408 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

LINKS

Table of n, a(n) for n=0..21.

FORMULA

G.f. satisfies: A(x)^2 = x*A( A(x)^2/x )/(1+x).

G.f. satisfies: A(x) = x + x*Series_Reversion( A(x)^2/x ). - Paul D. Hanna, Mar 15 2010

PROG

(PARI) {a(n)=local(A, B, F); if(n<1, 0, F=x+2*x^2+x*O(x^n); A=F; for(j=0, n, for(i=0, j, B=serreverse(A); A=(A+subst(B, x, A*(1+x)^2))/2); A=round(A)); polcoeff(sqrt(x*A), n, x))}

(PARI) {a(n)=local(A=1+x+x*O(x^n)); for(i=1, n, A=1+serreverse(x*subst(A^2, x, x+x*O(x^n)))); polcoeff(x*A, n)} \\ Paul D. Hanna, Mar 15 2010

CROSSREFS

Cf. A107095.

Sequence in context: A305537 A036076 A047797 * A138739 A216831 A221864

Adjacent sequences:  A107093 A107094 A107095 * A107097 A107098 A107099

KEYWORD

sign

AUTHOR

Paul D. Hanna, May 12 2005

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 June 25 14:54 EDT 2022. Contains 354851 sequences. (Running on oeis4.)