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

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A107086 G.f. A(x) satisfies: A(x)^4 = A(x^2)^2 + 4*x. 9
 1, 1, -1, 2, -5, 13, -35, 99, -289, 857, -2578, 7864, -24252, 75430, -236348, 745431, -2364399, 7536482, -24127482, 77544613, -250098478, 809169322, -2625483810, 8541037140, -27851360659, 91018956200, -298052119611, 977825373366, -3213513271929, 10577811289462, -34870732260397 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,4 COMMENTS Self-convolution is A107087. Self-convolution 4th power is A107088. LINKS Paul D. Hanna, Table of n, a(n) for n = 0..1000 EXAMPLE A(x)^4 = 1 + 4*x + 2*x^2 - x^4 + 2*x^6 - 5*x^8 + 12*x^10 - 30*x^12 +... A(x^2)^2 = 1 + 2*x^2 - x^4 + 2*x^6 - 5*x^8 + 12*x^10 - 30*x^12 +... MATHEMATICA nmin = 0; nmax = 30; sol = {a[0] -> 1}; Do[A[x_] = Sum[a[k] x^k, {k, 0, n}] /. sol; eq = CoefficientList[A[x]^4 - A[x^2]^2 - 4x + O[x]^(n+1), x][[2;; ]] == 0 /. sol; sol = sol ~Join~ Solve[eq][[1]], {n, 2, nmax}]; a /@ Range[nmin, nmax] /. sol (* Jean-François Alcover, Nov 07 2019 *) PROG (PARI) {a(n)=local(A=1+x); for(i=1, n, A=(subst(A, x, x^2)^2+4*x+x*O(x^n))^(1/4)); polcoeff(A, n, x)} for(n=0, 40, print1(a(n), ", ")) CROSSREFS Cf. A107087, A107088, A187814, A228711. Sequence in context: A037247 A308273 A126221 * A294790 A234643 A089846 Adjacent sequences: A107083 A107084 A107085 * A107087 A107088 A107089 KEYWORD sign AUTHOR Paul D. Hanna, May 11 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.

Last modified February 25 04:52 EST 2024. Contains 370310 sequences. (Running on oeis4.)