The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A219537 G.f. satisfies: A(x) = 1 + x*(A(x)^2 - A(x)^3 + A(x)^4). 4
 1, 1, 3, 13, 66, 366, 2148, 13115, 82449, 530095, 3469401, 23037642, 154820262, 1050999343, 7196493255, 49644745965, 344704716018, 2407157839593, 16895247295947, 119121868831235, 843306880720218, 5992060655349521, 42718501097385207, 305476181765843358 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,3 LINKS FORMULA G.f. A(x) satisfies [from Paul D. Hanna, Mar 21 2016]: (1) A(x)^2 = 1 + x*(A(x)^2 + A(x)^5). (2) A(x)^3 = 1 + x*(A(x)^2 + A(x)^4 + A(x)^6). Let F(x) = (1+x - sqrt(1 - 2*x - 3*x^2)) / (2*x), then g.f. A(x) satisfies: (3) A(x) = sqrt( (1/x)*Series_Reversion(x/F(x)^2) ), (4) A(x) = F(x*A(x)^2) and F(x) = A(x/F(x)^2), where F(x) = 1 + x*M(x) such that M(x) = 1 + x*M(x) + x^2*M(x)^2 is the g.f. of the Motzkin numbers (A001006). Let G(x) = 1 + x*G(x)/(1 - x*G(x)^2), then g.f. A(x) satisfies: (5) A(x) = (1/x)*Series_Reversion(x/G(x)), (6) A(x) = G(x*A(x)) and G(x) = A(x/G(x)). where G(x) is the g.f. of A106228. Recurrence: 3*n*(3*n-1)*(3*n+1)*(5*n-11)*(5*n-8)*(5*n-6)*a(n) = 6*(5*n-11)*(900*n^5 - 3870*n^4 + 6033*n^3 - 4165*n^2 + 1238*n - 120)*a(n-1) - 2*(n-2)*(5*n-1)*(950*n^4 - 5510*n^3 + 11199*n^2 - 9207*n + 2430)*a(n-2) + 6*(n-3)*(n-2)*(2*n-5)*(5*n-6)*(5*n-3)*(5*n-1)*a(n-3). - Vaclav Kotesovec, Aug 19 2013 a(n) ~ sqrt(300+75*10^(2/3)+30*10^(1/3))/90 * (5/9*10^(2/3)+10/9*10^(1/3)+8/3)^n / (sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Aug 19 2013 Recurrence: 18*n*(2*n+1)*(n+2)*(n+1)*a(n)-(6*(n+1))*(n+2)*(46*n^2+120*n+81)*a(n+1)+(2*(n+2))*(470*n^3+2605*n^2+4748*n+2829)*a(n+2)-(n+3)*(845*n^3+6465*n^2+14764*n+8712)*a(n+3)-(2*(n+4))*(130*n^3+2080*n^2+10407*n+16641)*a(n+4)+(3*(n+5))*(153*n^3+2493*n^2+13520*n+24404)*a(n+5)-(6*(n+5))*(3*n+17)*(3*n+19)*(n+6)*a(n+6) = 0. - Robert Israel, Feb 25 2018 EXAMPLE G.f.: A(x) = 1 + x + 3*x^2 + 13*x^3 + 66*x^4 + 366*x^5 + 2148*x^6 +... Related expansions: A(x)^2 = 1 + 2*x + 7*x^2 + 32*x^3 + 167*x^4 + 942*x^5 + 5593*x^6 +... A(x)^3 = 1 + 3*x + 12*x^2 + 58*x^3 + 312*x^4 + 1794*x^5 + 10794*x^6 +... A(x)^4 = 1 + 4*x + 18*x^2 + 92*x^3 + 511*x^4 + 3000*x^5 + 18316*x^6 +... A(x)^5 = 1 + 5*x + 25*x^2 + 135*x^3 + 775*x^4 + 4651*x^5 + 28845*x^6 +... A(x)^6 = 1 + 6*x + 33*x^2 + 188*x^3 + 1116*x^4 + 6852*x^5 + 43204*x^6 +... where A(x) = 1 + x*(A(x)^2 - A(x)^3 + A(x)^4), and A(x)^2 = 1 + x*(A(x)^2 + A(x)^5), and A(x)^3 = 1 + x*(A(x)^2 + A(x)^4 + A(x)^6), and A(x)^4 = 1 + x*(A(x)^2 + A(x)^4 + A(x)^5 + A(x)^7), and A(x)^5 = 1 + x*(A(x)^2 + A(x)^4 + A(x)^5 + A(x)^6 + A(x)^8), etc. The g.f. satisfies A(x) = F(x*A(x)^2) and F(x) = A(x/F(x)^2) where F(x) = 1 + x + x^2 + 2*x^3 + 4*x^4 + 9*x^5 + 21*x^6 + 51*x^7 +...+ A001006(n-1)*x^n +... is a g.f. of the Motzkin numbers (A001006, shifted right 1 place). The g.f. satisfies A(x) = G(x*A(x)) and G(x) = A(x/G(x)) where G(x) = 1 + x + 2*x^2 + 6*x^3 + 21*x^4 + 80*x^5 + 322*x^6 +...+ A106228(n)*x^n +... satisfies G(x) = 1 + x*G(x)/(1 - x*G(x)^2). MAPLE rec := {(36*n^4+126*n^3+126*n^2+36*n)*a(n)+(-276*n^4-1548*n^3-3198*n^2-2898*n-972)*a(n+1)+(940*n^4+7090*n^3+19916*n^2+24650*n+11316)*a(n+2)+(-845*n^4-9000*n^3-34159*n^2-53004*n-26136)*a(n+3)+(-260*n^4-5200*n^3-37454*n^2-116538*n-133128)*a(n+4)+(459*n^4+9774*n^3+77955*n^2+276012*n+366060)*a(n+5)+(-54*n^4-1242*n^3-10686*n^2-40758*n-58140)*a(n+6), a(0) = 1, a(1) = 1, a(2) = 3, a(3) = 13, a(4) = 66, a(5) = 366}: f:= gfun:-rectoproc(rec, a(n), remember): map(f, [\$0..50]); # Robert Israel, Feb 25 2018 MATHEMATICA nmax = 23; sol = {a[0] -> 1}; Do[A[x_] = Sum[a[k] x^k, {k, 0, n}] /. sol; eq = CoefficientList[A[x] - (1 + x (A[x]^2 - A[x]^3 + A[x]^4)) + O[x]^(n + 1), x] == 0 /. sol; sol = sol ~Join~ Solve[eq][[1]], {n, 1, nmax}]; sol /. Rule -> Set; a /@ Range[0, nmax] (* Jean-François Alcover, Nov 01 2019 *) PROG (PARI) /* Formula A(x) = 1 + x*(A(x)^2 - A(x)^3 + A(x)^4): */ {a(n)=local(A=1); for(i=1, n, A=1+x*(A^2-A^3+A^4) +x*O(x^n)); polcoeff(A, n)} for(n=0, 25, print1(a(n), ", ")) (PARI) /* Formula using Series Reversion involving Motzkin numbers: */ {a(n)=local(A=1); A=(1+x-sqrt(1-2*x-3*x^2+x^3*O(x^n)))/(2*x); polcoeff(sqrt(1/x*serreverse(x/A^2)), n)} for(n=0, 25, print1(a(n), ", ")) CROSSREFS Cf. A001006, A106228, A219534, A219535, A219536, A219538. Sequence in context: A260783 A228987 A112807 * A045743 A110530 A142979 Adjacent sequences:  A219534 A219535 A219536 * A219538 A219539 A219540 KEYWORD nonn AUTHOR Paul D. Hanna, Nov 21 2012 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 January 25 01:10 EST 2021. Contains 340414 sequences. (Running on oeis4.)