cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A103970 Expansion of (1 - sqrt(1 - 4*x - 12*x^2))/(2*x).

Original entry on oeis.org

1, 4, 8, 32, 128, 576, 2688, 13056, 65024, 330752, 1710080, 8962048, 47497216, 254132224, 1370849280, 7447117824, 40707293184, 223731253248, 1235630948352, 6853893292032, 38166664839168, 213288826699776, 1195775593807872, 6723691157127168, 37908469021409280, 214260335517892608, 1213784937073737728, 6890689428042285056
Offset: 0

Views

Author

Paul Barry, Feb 23 2005

Keywords

Comments

Image of c(x), the g.f. of the Catalan numbers A000108 under the mapping g(x) -> (1+3x)g(x(1+3x)). In general, the image of the Catalan numbers under the mapping g(x) -> (1+i*x)g(x(1+i*x)) is given by a(n) = Sum_{k=0..n} i^(n-k)*C(k)*C(k+1,n-k).
Hankel transform is 4^C(n+1,2)*A128018(n). [Paul Barry, Nov 20 2009]
By following L. Comtet [Analyse Combinatoire Tomes 1 et 2, PUF, Paris 1970], we also obtain (n+1)*C(n) - 2*a*(2*n-1)*C(n-1) + 4*(n-2)*(a^2-b)*C(n-2) = 0. In the present case, we also have the asymptotic result: a(n) ~ sqrt(4/3)*2^(n-1)*3^(n+1)/sqrt(Pi*n^3) for large n. - Richard Choulet, Dec 17 2009

Links

Crossrefs

Cf. A025227, A025229, A103971, A103972.
Cf. A000108, A025228, A025230, A025231, A025232. [Richard Choulet, Dec 17 2009]

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 35); Coefficients(R!( (1-Sqrt(1-4*x-12*x^2))/(2*x) )); // G. C. Greubel, Mar 16 2019
  • Maple
    n:=30:a(0):=1:a(1):=4: k:=1: for k from 1 to n do a(k+1):=sum('a(p)*a(k-p)','p'=0..k):od:seq(a(k),k=0..n); # Richard Choulet, Dec 17 2009
taylor(((1-(1-4*z-12*z^2)^0.5)/(2*z)),z=0,32); # Richard Choulet, Dec 17 2009
  • Mathematica
    CoefficientList[Series[(1 - Sqrt[1-4x-12x^2])/(2x), {x, 0, 33}], x] (* Vincenzo Librandi, Aug 18 2017 *)
  • PARI
    my(x='x+O('x^35)); Vec((1-sqrt(1-4*x-12*x^2))/(2*x)) \\ G. C. Greubel, Mar 16 2019
  • Sage
    ((1-sqrt(1-4*x-12*x^2))/(2*x)).series(x, 35).coefficients(x, sparse=False) # G. C. Greubel, Mar 16 2019

Formula

G.f.: (1 - sqrt(1-4*x*(1+3*x)))/(2*x).
a(n) = Sum_{k=0..n} 3^(n-k)*C(k)*C(k+1, n-k).
D-finite with recurrence: (n+1)*a(n) = 2*(2*n-1)*a(n-1) + 12*(n-2)*a(n-2). - Richard Choulet, Dec 17 2009

A103971 Expansion of (1 - sqrt(1 - 4*x - 16*x^2))/(2*x).

Original entry on oeis.org

1, 5, 10, 45, 190, 930, 4660, 24445, 131190, 719830, 4013260, 22684370, 129661740, 748252580, 4353379560, 25508284445, 150392391590, 891549228430, 5310994644060, 31775749689670, 190860711108740, 1150473009844380
Offset: 0

Views

Author

Paul Barry, Feb 23 2005

Keywords

Comments

Image of c(x), the g.f. of the Catalan numbers A000108 under the mapping g(x) -> (1+4x)g(x(1+4x)). In general, the image of the Catalan numbers under the mapping g(x)->(1+i*x)g(x(1+i*x)) is given by a(n) = Sum_{k=0..n} i^(n-k)C(k)C(k+1,n-k).
More generally, the sequence C for which C(0)=a, C(1)=b and C(n+1) = sum(C(k)*C(n-k),k=0..n) has the following g.f. f: f(z) = (1-sqrt(1-4*z*(a-(a^2-b)*z)))/(2*z). We obtain: C(n)=(sum(-1)^(p-1)*2^{n-p}a^{n-2*p-1}*(a^2-b)^p*((2*n-2*p-1)*...*5*3*1/(p!*(n-2*p+1)!)),p=0..floor((n+1)/2)). By following Comtet [Analyse Combinatoire Tomes 1 et 2, PUF, Paris 1970], we obtain also: (n+1)*C(n) - 2*a*(2*n-1)*C(n-1) + 4*(n-2)*(a^2-b)*C(n-2) = 0. - Richard Choulet, Dec 17 2009

Links

Crossrefs

Cf. A000108, A025227, A025228, A025229, A025230, A025231, A025232. - Richard Choulet, Dec 17 2009

Programs

  • Maple
    n:=30:a(0):=1:a(1):=5: for k from 1 to n do a(k+1):=sum('a(p)*a(k-p)','p'=0..k):od:seq(a(k),k=0..n); # Richard Choulet, Dec 17 2009
  • Mathematica
    CoefficientList[Series[(1-Sqrt[1-4x-16x^2])/(2x),{x,0,30}],x] (* Harvey P. Dale, Apr 02 2012 *)

Formula

G.f.: (1-sqrt(1-4*x*(1+4*x)))/(2*x).
a(n) = Sum_{k=0..n} 4^(n-k)*C(k)*C(k+1, n-k).
Another recurrence formula: (n+1)*a(n) = 2*(2*n-1)*a(n-1) + 16*(n-2)*a(n-2). - Richard Choulet, Dec 17 2009
a(n) ~ sqrt(10 + 2*sqrt(5))*(2 + 2*sqrt(5))^n/(2*sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Oct 17 2012
Equivalently, a(n) ~ 5^(1/4) * 2^(2*n) * phi^(n + 1/2) / (sqrt(Pi) * n^(3/2)), where phi = A001622 is the golden ratio. - Vaclav Kotesovec, Dec 08 2021

A103972 Expansion of (1-sqrt(1-4*x-20*x^2))/(2*x).

Original entry on oeis.org

1, 6, 12, 60, 264, 1392, 7392, 41424, 236640, 1384512, 8224896, 49554816, 301884672, 1856878080, 11514915840, 71915838720, 451938731520, 2855705994240, 18132621772800, 115637702461440, 740356410961920, 4756888756101120, 30662391191715840, 198229520200704000, 1285001080928845824
Offset: 0

Views

Author

Paul Barry, Feb 23 2005

Keywords

Comments

Image of c(x), the g.f. of the Catalan numbers A000108 under the mapping g(x)->(1+5x)g(x(1+5x)). In general, the image of the Catalan numbers under the mapping g(x)->(1+i*x)g(x(1+i*x)) is given by a(n)=sum{k=0..n, i^(n-k)C(k)C(k+1,n-k)}.
More generally, the sequence C for which C(0)=a, C(1)=b and C(n+1)=sum(C(k)*C(n-k),k=0..n) has the following G.f f: f(z)= (1-sqrt(1-4*z*(a-(a^2-b)*z)))/(2*z). We obtain: C(n)=(sum(-1)^(p-1)*2^{n-p}a^{n-2*p-1}*(a^2-b)^p*((2*n-2*p-1)*...*5*3*1/(p!*(n-2*p+1)!)),p=0..floor((n+1)/2)). By following L. Comtet [Analyse Combinatoire Tomes 1 et 2, PUF, Paris 1970], we obtain also: (n+1)*C(n)-2*a*(2*n-1)*C(n-1)+4*(n-2)*(a^2-b)*C(n-2)=0. - Richard Choulet, Dec 17 2009

Links

Crossrefs

Cf. A000108, A025227, A025228, A025229, A025230, A025231, A025232. [From Richard Choulet, Dec 17 2009]

Programs

  • Maple
    n:=30:a(0):=1:a(1):=6 :for k from 1 to n do a(k+1):=sum('a(p)*a(k-p)','p'=0..k):od:seq(a(k),k=0..n); # Richard Choulet, Dec 17 2009
  • Mathematica
    CoefficientList[Series[(1-Sqrt[1-4*x-20*x^2])/(2*x), {x, 0, 20}], x] (* Vaclav Kotesovec, Oct 17 2012 *)
  • PARI
    x='x+O('x^66); Vec((1-sqrt(1-4*x-20*x^2))/(2*x)) \\ Joerg Arndt, May 13 2013

Formula

G.f.: (1-sqrt(1-4*x*(1+5*x)))/(2*x).
a(n) = Sum_{k=0..n} 5^(n-k)*C(k)*C(k+1, n-k).
Another recurrence formula: (n+1)*a(n)=2*(2n-1)*a(n-1)+20*(n-2)*a(n-2). - Richard Choulet, Dec 17 2009
a(n) ~ sqrt(12+2*sqrt(6))*(2+2*sqrt(6))^n/(2*sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Oct 17 2012
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.