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-10 of 30 results. Next

A026375 a(n) = Sum_{k=0..n} binomial(n,k)*binomial(2*k,k).

Original entry on oeis.org

1, 3, 11, 45, 195, 873, 3989, 18483, 86515, 408105, 1936881, 9238023, 44241261, 212601015, 1024642875, 4950790605, 23973456915, 116312293305, 565280386625, 2751474553575, 13411044301945, 65448142561035, 319756851757695
Offset: 0

Views

Author

Keywords

Comments

a(n) is the number of integer strings s(0),...,s(n) counted by array T in A026374 that have s(n)=0; also a(n)=T(2n,n).
Partial sums of A085362. Number of bilateral Schroeder paths (i.e., lattice paths consisting of steps U=(1,1), D=(1,-1) and H=(2,0)) from (0,0) to (2n,0) and with no H-steps at odd (positive or negative) levels. Example: a(2)=11 because we have HUD, UDH, UDUD, UUDD, UDDU, their reflections in the x-axis and HH. - Emeric Deutsch, Jan 30 2004
Largest coefficient of (1+3*x+x^2)^n; row sums of triangle in A124733. - Philippe Deléham, Oct 02 2007
Also number of paths from (0,0) to (n,0) using steps U=(1,1), H=(1,0) and D=(1,-1), the H steps come in three colors. - N-E. Fahssi, Feb 05 2008
Equals INVERT transform of A109033: (1, 2, 6, 22, 88, ...), INVERTi transform of A111966, binomial transform of A000984, and inverse Binomial transform of A081671. Convolved with A002212: (1, 3, 10, 36, ...) = A026376: (1, 6, 30, 144, ...). Equals convolution square root of A003463: (1, 6, 31, 156, 781, 3906, ...). - Gary W. Adamson, May 17 2009
Diagonal of array with rational generating function 1/(1 - (x^2 + 3*x*y + y^2)). - Gheorghe Coserea, Jul 29 2018
a(n) == 0 (mod 3) if and only if n is in A081606. - Fabio Visonà, Aug 03 2023

Examples

			G.f. = 1 + 3*x + 11*x^2 + 45*x^3 + 195*x^4 + 873*x^5 + 3989*x^6 + ...
		

Crossrefs

Column 3 of A292627. Column 1 of A110165. Central column of A272866.
First differences are in A085362. Bisection of A026380.
m-th binomial transforms of A000984: A126869 (m = -2), A002426 (m = -1 and m = -3 for signed version), A000984 (m = 0 and m = -4 for signed version), A026375 (m = 1 and m = -5 for signed version), A081671 (m = 2 and m = -6 for signed version), A098409 (m = 3 and m = -7 for signed version), A098410 (m = 4 and m = -8 for signed version), A104454 (m = 5 and m = -9 for signed version).

Programs

  • GAP
    List([0..25],n->Sum([0..n],k->Binomial(n,k)*Binomial(2*k,k))); # Muniru A Asiru, Jul 29 2018
  • Haskell
    a026375 n = a026374 (2 * n) n  -- Reinhard Zumkeller, Feb 22 2014
    
  • Maple
    seq( add(binomial(n,k)*binomial(2*k,k), k=0..n), n=0..30 ); # Detlef Pauly (dettodet(AT)yahoo.de), Nov 08 2001
    a := n -> simplify(GegenbauerC(n, -n, -3/2)):
    seq(a(n), n=0..22); # Peter Luschny, May 09 2016
  • Mathematica
    Table[SeriesCoefficient[1/Sqrt[1-6*x+5*x^2],{x,0,n}],{n,0,20}] (* Vaclav Kotesovec, Oct 08 2012 *)
    (* From Michael Somos, May 11 2014: (Start) *)
    a[ n_] := Sum[ Binomial[n, k] Binomial[2 k, k], {k, 0, n}];
    a[ n_] := If[ n < 0, 0, Hypergeometric2F1[-n, 1/2, 1, -4]];
    a[ n_] := If[ n < 0, 0, Coefficient[(1 + 3 x + x^2)^n, x, n]];
    a[ n_] := If[ n < 0, 0, n! SeriesCoefficient[Exp[3 x] BesselI[0,2 x], {x, 0, n}]];
    (* (End) *)
  • Maxima
    A026375(n):=coeff(expand((1+3*x+x^2)^n),x,n);
    makelist(A026375(n),n,0,12); /* Emanuele Munarini, Mar 02 2011 */
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( (1 + 3*x + x^2)^n, n))}; /* Michael Somos, Sep 09 2002 */
    
  • PARI
    a(n)={my(v=Vec((1-x-x^2)^n)); sum(k=1,#v, v[k]^2);} \\ Joerg Arndt, Jul 06 2011
    
  • PARI
    {a(n) = sum(k=0, n, 5^(n-k)*(-1)^k*binomial(n, k)*binomial(2*k, k))} \\ Seiichi Manyama, Apr 22 2019
    
  • PARI
    {a(n) = sum(k=0, n\2, 3^(n-2*k)*binomial(n, 2*k)*binomial(2*k, k))} \\ Seiichi Manyama, May 04 2019
    

Formula

Representation by Gauss's hypergeometric function, in Maple notation: a(n)=hypergeom([ -n, 1/2 ], [ 1 ], -4). - Karol A. Penson, Apr 20 2001
This sequence is the binomial transform of A000984. - John W. Layman, Aug 11 2000; proved by Emeric Deutsch, Oct 26 2002
E.g.f.: exp(3*x)*I_0(2x), where I_0 is Bessel function. - Michael Somos, Sep 17 2002
G.f.: 1/sqrt(1-6*x+5*x^2). - Emeric Deutsch, Oct 26 2002
D-finite with recurrence: n*a(n)-3*(2*n-1)*a(n-1)+5*(n-1)*a(n-2)=0 for n > 1. - Emeric Deutsch, Jan 24 2004
From Emeric Deutsch, Jan 30 2004: (Start)
a(n) = [t^n](1+3*t+t^2)^n;
a(n) = Sum_{j=ceiling(n/2)..n} 3^(2*j-n)*binomial(n, j)*binomial(j, n-j). (End)
a(n) = A026380(2*n-1) (n>0). - Emeric Deutsch, Feb 18 2004
G.f.: 1/(1-x-2*x/(1-x/(1-x-x/(1-x/(1-x-x/(1-x/(1-x-x/(1-x... (continued fraction). - Paul Barry, Jan 06 2009
a(n) = sum of squared coefficients of (1+x-x^2)^n - see triangle A084610. - Paul D. Hanna, Jul 18 2009
a(n) = sum of squares of coefficients of (1-x-x^2)^n. - Joerg Arndt, Jul 06 2011
a(n) = (1/Pi)*Integral_{x=-2..2} ((3+x)^n/sqrt((2-x)*(2+x))) dx. - Peter Luschny, Sep 12 2011
a(n) ~ 5^(n+1/2)/(2*sqrt(Pi*n)). - Vaclav Kotesovec, Oct 08 2012
G.f.: G(0)/(1-x), where G(k) = 1 + 4*x*(4*k+1)/( (4*k+2)*(1-x) - 2*x*(1-x)*(2*k+1)*(4*k+3)/(x*(4*k+3) + (1-x)*(k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 24 2013
0 = a(n)*(+25*a(n+1) - 45*a(n+2) + 10*a(n+3)) + a(n+1)*(-15*a(n+1) + 36*a(n+2) - 9*a(n+3)) + a(n+2)*(-3*a(n+2) + a(n+3)) for all n in Z. - Michael Somos, May 11 2014
a(n) = GegenbauerC(n, -n, -3/2). - Peter Luschny, May 09 2016
a(n) = Sum_{k=0..n} 5^(n-k) * (-1)^k * binomial(n,k) * binomial(2*k,k). - Seiichi Manyama, Apr 22 2019
a(n) = Sum_{k=0..floor(n/2)} 3^(n-2*k) * binomial(n,2*k) * binomial(2*k,k). - Seiichi Manyama, May 04 2019
a(n) = (1/Pi) * Integral_{x = -1..1} (1 + 4*x^2)^n/sqrt(1 - x^2) dx = (1/Pi) * Integral_{x = -1..1} (5 - 4*x^2)^n/sqrt(1 - x^2) dx. - Peter Bala, Jan 27 2020
From Peter Bala, Jan 10 2022: (Start)
1 + x*exp(Sum_{n >= 1} a(n)*x^n/n) = 1 + x + 3*x^2 + 10*x^3 + 36*x^4 + ... is the o.g.f. of A002212.
The Gauss congruences a(n*p^k) == a(n*p^(k-1)) (mod p^k) hold for prime p and positive integers n and k. (End)
a(n) = (1/4)^n * Sum_{k=0..n} 5^k * binomial(2*k,k) * binomial(2*(n-k),n-k). - Seiichi Manyama, Aug 18 2025

Extensions

Definition simplified by N. J. A. Sloane, Feb 16 2012

A110170 First differences of the central Delannoy numbers (A001850).

Original entry on oeis.org

1, 2, 10, 50, 258, 1362, 7306, 39650, 217090, 1196834, 6634890, 36949266, 206549250, 1158337650, 6513914634, 36718533570, 207412854786, 1173779487810, 6653482333450, 37770112857074, 214694383882498, 1221832400430482, 6961037946938250, 39697830840765090, 226596964146630658
Offset: 0

Views

Author

Emeric Deutsch, Jul 14 2005

Keywords

Comments

Number of Delannoy paths of length n that do not start with a (1, 1) step (a Delannoy path of length n is a path from (0, 0) to (n, n), consisting of steps E = (1, 0), N = (0, 1) and D = (1, 1)). Example: a(1) = 2 because we have NE and EN. Column 0 of A110169 (also nonzero entries in each column of A110169).
For n > 0: a(n) = A128966(2*n,n). - Reinhard Zumkeller, Jul 20 2013

Crossrefs

Programs

  • Haskell
    a110170 0 = 1
    a110170 n = a128966 (2 * n) n  -- Reinhard Zumkeller, Jul 20 2013
  • Maple
    with(orthopoly): a:=proc(n) if n=0 then 1 else P(n,3)-P(n-1,3) fi end: seq(a(n),n=0..25);
    a := n -> `if`(n=0, 1, 2*hypergeom([1 - n, -n], [1], 2)):
    seq(simplify(a(n)), n=0..24); # Peter Luschny, May 22 2017
  • Mathematica
    CoefficientList[Series[(1 - x)/Sqrt[1 - 6 * x + x^2], {x, 0, 20}], x] (* Vaclav Kotesovec, Oct 18 2012 *)
  • PARI
    x='x+O('x^66); Vec((1-x)/sqrt(1-6*x+x^2)) \\ Joerg Arndt, May 16 2013
    

Formula

G.f.: (1-z)/sqrt(1-6*z+z^2).
a(n) = P_n(3) - P_{n-1}(3) (n >= 1), where P_j is j-th Legendre polynomial.
From Paul Barry, Oct 18 2009: (Start)
G.f.: (1-x)/(1-x-2x/(1-x-x/(1-x-x/(1-x-x/(1-... (continued fraction);
G.f.: 1/(1-2x/((1-x)^2-x/(1-x/((1-x)^2-x/(1-x/((1-x)^2-x/(1-... (continued fraction);
a(n) = Sum_{k = 0..n} (0^(n + k) + C(n + k - 1, 2k - 1)) * C(2k, k) = 0^n + Sum_{k = 0..n} C(n + k - 1, 2k - 1) * C(2k, k). (End)
D-finite with recurrence: n*(2*n-3)*a(n) = 2*(6*n^2-12*n+5)*a(n-1) - (n-2)*(2*n-1)*a(n-2). - Vaclav Kotesovec, Oct 18 2012
a(n) ~ 2^(-1/4)*(3+2*sqrt(2))^n/sqrt(Pi*n). - Vaclav Kotesovec, Oct 18 2012
a(n) = A277919(2n,n). - John P. McSorley, Nov 23 2016
a(n) = 2*hypergeom([1 - n, -n], [1], 2) for n>0. - Peter Luschny, May 22 2017
D-finite with recurrence: n*a(n) +(-7*n+5)*a(n-1) +(7*n-16)*a(n-2) +(-n+3)*a(n-3)=0. - R. J. Mathar, Jan 15 2020
a(0) = 1; a(n) = (2/n) * Sum_{k=0..n-1} (n^2-k^2) * a(k). - Seiichi Manyama, Mar 28 2023
G.f.: Sum_{n >= 0} binomial(2*n, n)*x^n/(1 - x)^(2*n) = 1 + 2*x + 10*x^2 + 50*x^3 + .... - Peter Bala, Apr 17 2024

A361375 Expansion of 1/(1 - 9*x/(1 - x))^(1/3).

Original entry on oeis.org

1, 3, 21, 165, 1380, 11982, 106626, 965442, 8854725, 82022115, 765787773, 7195638909, 67973370618, 644991134880, 6143707229880, 58714212503784, 562741793028282, 5407273475087934, 52074626299010130, 502513862912425650, 4857975310180620720
Offset: 0

Views

Author

Seiichi Manyama, Mar 28 2023

Keywords

Crossrefs

Programs

  • Maple
    a := n -> if n = 0 then 1 else 3*hypergeom([1 - n, 4/3], [2], -9) fi:
    seq(simplify(a(n)), n = 0..20); # Peter Luschny, Mar 30 2023
  • PARI
    my(N=30, x='x+O('x^N)); Vec(1/(1-9*x/(1-x))^(1/3))

Formula

a(n) = Sum_{k=0..n} (-9)^k * binomial(-1/3,k) * binomial(n-1,n-k).
a(0) = 1; a(n) = (3/n) * Sum_{k=0..n-1} (n+2*k) * a(k).
n*a(n) = (11*n-8)*a(n-1) - 10*(n-2)*a(n-2) for n > 1.
a(n) ~ 3^(2/3) * 10^(n - 1/3) / (Gamma(1/3) * n^(2/3)). - Vaclav Kotesovec, Mar 28 2023
a(n) = 3*hypergeom([1 - n, 4/3], [2], -9) for n >= 1. - Peter Luschny, Mar 30 2023

A085363 a(0)=1, for n>0: a(n) = 4*9^(n-1) - (1/2)*Sum_{i=1..n-1} a(i)*a(n-i).

Original entry on oeis.org

1, 4, 28, 212, 1676, 13604, 112380, 940020, 7936620, 67494980, 577309148, 4961187092, 42801458764, 370478720356, 3215827927228, 27982214082612, 244004165618220, 2131710838837380, 18654504783815580, 163488269572628820
Offset: 0

Views

Author

Mario Catalani (mario.catalani(AT)unito.it), Jun 25 2003

Keywords

Comments

Apparently, the number of 2-D directed walks of semilength n starting at (0,0) and ending on the X-axis using steps NE, SE, NW and SW avoiding adjacent NW/SE and adjacent NE/SW. - David Scambler, Jun 20 2013
Form an array with m(0,n) = m(n,0) = 2^n; m(i,j) equals the sum of the terms to the left of m(i,j) and the sum of the terms above m(i,j), which is m(i,j) = Sum_{k-0..j-1} m(i,k) + Sum_{k=0..i-1} m(k,j). m(n,n) = a(n). - J. M. Bergot, Jul 10 2013
From G. C. Greubel, May 22 2020: (Start)
This sequence is part of a class of sequences, for m >= 0, with the properties:
a(n) = 2*m*(4*m+1)^(n-1) - (1/2)*Sum_{k=1..n-1} a(k)*a(n-k).
a(n) = Sum_{k=0..n} m^k * binomial(n-1, n-k) * binomial(2*k, k).
a(n) = (2*m) * Hypergeometric2F1(-n+1, 3/2; 2; -4*m), for n > 0.
n*a(n) = 2*((2*m+1)*n - (m+1))*a(n-1) - (4*m+1)*(n-2)*a(n-2).
(4*m + 1)^n = Sum_{k=0..n} Sum_{j=0..k} a(j)*a(k-j).
G.f.: sqrt( (1 - t)/(1 - (4*m+1)*t) ).
This sequence is the case of m=2. (End)
The number of elements in the free group on two generators of length 2n that are zero exponent sum. - Tey Berendschot, Aug 09 2021

Crossrefs

Cf. A001019 (9^n), A084771, A085362, A085364, diagonal of A348595.

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 30); Coefficients(R!( Sqrt((1-x)/(1-9*x)) )); // G. C. Greubel, May 23 2020
    
  • Maple
    seq(coeff(series(sqrt((1-x)/(1-9*x)), x, n+1), x, n), n = 0..30); # G. C. Greubel, May 23 2020
  • Mathematica
    CoefficientList[Series[Sqrt[(1-x)/(1-9x)], {x, 0, 25}], x]
  • PARI
    my(x='x+O('x^66)); Vec(sqrt((1-x)/(1-9*x)) ) \\ Joerg Arndt, May 10 2013
    
  • Sage
    def A085363_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( sqrt((1-x)/(1-9*x)) ).list()
    A085363_list(30) # G. C. Greubel, May 23 2020

Formula

G.f.: sqrt((1-x)/(1-9*x)).
Sum_{i=0..n} Sum_{j=0..i} a(j)*a(i-j) = 9^n.
From Vladeta Jovovic, Oct 10 2003: (Start)
First differences of A084771.
a(n) = Sum_{k=1..n} 2^k * binomial(n-1, k-1) * binomial(2*k, k). (End)
D-finite with recurrence n*a(n) = (10*n-6)*a(n-1) - (9*n-18)*a(n-2). - Vladeta Jovovic, Jul 16 2004
a(n) ~ 2*sqrt(2)*3^(2*n-1)/sqrt(Pi*n). - Vaclav Kotesovec, Oct 14 2012
a(0) = 1; a(n) = (4/n) * Sum_{k=0..n-1} (n+k) * a(k). - Seiichi Manyama, Mar 28 2023
From Seiichi Manyama, Aug 22 2025: (Start)
a(n) = (1/4)^n * Sum_{k=0..n} 9^k * binomial(2*k,k) * binomial(2*(n-k),n-k)/(1-2*(n-k)).
a(n) = Sum_{k=0..n} (-2)^k * 9^(n-k) * binomial(2*k,k)/(1-2*k) * binomial(n-1,n-k). (End)

A359758 Expansion of 1/sqrt(1 - 4*x/(1-x)^5).

Original entry on oeis.org

1, 2, 16, 110, 770, 5512, 40066, 294484, 2182850, 16288430, 122198926, 920820578, 6964483628, 52840433000, 401990254180, 3065365241440, 23422905551018, 179302895759782, 1374785979255880, 10556280995419090, 81161958814162700, 624750086745027388
Offset: 0

Views

Author

Seiichi Manyama, Mar 24 2023

Keywords

Crossrefs

Programs

  • PARI
    my(N=30, x='x+O('x^N)); Vec(1/sqrt(1-4*x/(1-x)^5))
    
  • PARI
    a(n)=sum(k=0, n, binomial(2*k,k) * binomial(n+4*k-1,n-k)) \\ Winston de Greef, Mar 24 2023

Formula

a(n) = Sum_{k=0..n} binomial(2*k,k) * binomial(n+4*k-1,n-k).
n*a(n) = (10*n-8)*a(n-1) - (19*n-46)*a(n-2) + 20*(n-3)*a(n-3) - 15*(n-4)*a(n-4) + 6*(n-5)*a(n-5) - (n-6)*a(n-6) for n > 5.
a(0) = 1; a(n) = (2/n) * Sum_{k=0..n-1} (n+k) * binomial(n+3-k,4) * a(k).

A360132 Expansion of 1/sqrt(1 - 4*x/(1-x)^6).

Original entry on oeis.org

1, 2, 18, 134, 1010, 7788, 60978, 482708, 3853338, 30964238, 250150176, 2029781310, 16530857930, 135051216620, 1106287906140, 9083459084364, 74734798117570, 615998603183550, 5085522355488150, 42045309424052250, 348067638153560040, 2884832348569699340
Offset: 0

Views

Author

Seiichi Manyama, Mar 24 2023

Keywords

Crossrefs

Programs

  • PARI
    my(N=30, x='x+O('x^N)); Vec(1/sqrt(1-4*x/(1-x)^6))

Formula

a(n) = Sum_{k=0..n} binomial(2*k,k) * binomial(n+5*k-1,n-k).
n*a(n) = (11*n-9)*a(n-1) - (25*n-60)*a(n-2) + 35*(n-3)*a(n-3) - 35*(n-4)*a(n-4) + 21*(n-5)*a(n-5) - 7*(n-6)*a(n-6) + (n-7)*a(n-7) for n > 6.
a(0) = 1; a(n) = (2/n) * Sum_{k=0..n-1} (n+k) * binomial(n+4-k,5) * a(k).

A360318 a(n) = Sum_{k=0..n} 3^(n-k) * binomial(n-1,n-k) * binomial(2*k,k).

Original entry on oeis.org

1, 2, 12, 74, 466, 2982, 19320, 126390, 833220, 5527190, 36852052, 246751854, 1658106394, 11176100138, 75528743352, 511600414554, 3472363279170, 23609924743590, 160788499672020, 1096566516149790, 7488135911236806, 51193972101241362, 350368409215623192
Offset: 0

Views

Author

Seiichi Manyama, Feb 03 2023

Keywords

Crossrefs

Programs

  • PARI
    a(n) = sum(k=0, n, 3^(n-k)*binomial(n-1, n-k)*binomial(2*k, k));
    
  • PARI
    my(N=30, x='x+O('x^N)); Vec(sqrt((1-3*x)/(1-7*x)))

Formula

G.f.: sqrt( (1-3*x)/(1-7*x) ).
n*a(n) = 2*(5*n-4)*a(n-1) - 21*(n-2)*a(n-2).
Sum_{i=0..n} Sum_{j=0..i} (1/3)^i * a(j) * a(i-j) = (7/3)^n.
a(n) = 2 * A122898(n-1) for n > 0.
a(n) ~ 2 * 7^(n - 1/2) / sqrt(Pi*n). - Vaclav Kotesovec, Feb 04 2023
From Seiichi Manyama, Aug 22 2025: (Start)
a(n) = (1/4)^n * Sum_{k=0..n} 7^k * 3^(n-k) * binomial(2*k,k) * binomial(2*(n-k),n-k)/(1-2*(n-k)).
a(n) = Sum_{k=0..n} (-1)^k * 7^(n-k) * binomial(2*k,k)/(1-2*k) * binomial(n-1,n-k). (End)

A360317 a(n) = Sum_{k=0..n} 2^(n-k) * binomial(n-1,n-k) * binomial(2*k,k).

Original entry on oeis.org

1, 2, 10, 52, 278, 1516, 8388, 46920, 264678, 1503052, 8581676, 49215256, 283297660, 1635904376, 9472214344, 54975423504, 319729353606, 1862896455180, 10871759717916, 63539265366264, 371837338366740, 2178604586281128, 12778264475444280, 75022726995053808
Offset: 0

Views

Author

Seiichi Manyama, Feb 03 2023

Keywords

Crossrefs

Programs

  • PARI
    a(n) = sum(k=0, n, 2^(n-k)*binomial(n-1, n-k)*binomial(2*k, k));
    
  • PARI
    my(N=30, x='x+O('x^N)); Vec(sqrt((1-2*x)/(1-6*x)))

Formula

G.f.: sqrt( (1-2*x)/(1-6*x) ).
n*a(n) = 2*(4*n-3)*a(n-1) - 12*(n-2)*a(n-2).
Sum_{i=0..n} Sum_{j=0..i} (1/2)^i * a(j) * a(i-j) = 3^n.
a(n) = 2 * A005573(n-1) for n > 0.
a(n) ~ 2^(n + 1/2) * 3^(n - 1/2) / sqrt(Pi*n). - Vaclav Kotesovec, Feb 04 2023
From Seiichi Manyama, Aug 22 2025: (Start)
a(n) = (1/2)^n * Sum_{k=0..n} 3^k * binomial(2*k,k) * binomial(2*(n-k),n-k)/(1-2*(n-k)).
a(n) = Sum_{k=0..n} (-1)^k * 6^(n-k) * binomial(2*k,k)/(1-2*k) * binomial(n-1,n-k). (End)

A377199 Expansion of 1/(1 - 4*x/(1-x))^(5/2).

Original entry on oeis.org

1, 10, 80, 570, 3790, 24062, 147780, 885190, 5199560, 30065870, 171623328, 969151710, 5422863630, 30105497970, 165993714540, 909770119914, 4959840748350, 26912374137150, 145411035749600, 782681600883950, 4198276264607290, 22448626776903450, 119690255236279100
Offset: 0

Views

Author

Seiichi Manyama, Oct 19 2024

Keywords

Crossrefs

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 35); Coefficients(R!( 1/(1 - 4*x/(1-x))^(5/2))); // Vincenzo Librandi, May 10 2025
  • Mathematica
    Table[Sum[(-4)^k*Binomial[-5/2,k]*Binomial[n-1,n-k],{k,0,n}],{n,0,35}] (* Vincenzo Librandi, May 10 2025 *)
  • PARI
    a(n) = sum(k=0, n, (-4)^k*binomial(-5/2, k)*binomial(n-1, n-k));
    

Formula

a(0) = 1; a(n) = 2 * Sum_{k=0..n-1} (5-3*k/n) * a(k).
a(n) = (2*(3*n+2)*a(n-1) - 5*(n-2)*a(n-2))/n for n > 1.
a(n) = Sum_{k=0..n} (-4)^k * binomial(-5/2,k) * binomial(n-1,n-k).
a(n) ~ 128 * n^(3/2) * 5^(n - 5/2) / (3*sqrt(Pi)). - Vaclav Kotesovec, Oct 26 2024
a(n) = 10*hypergeom([7/2, 1-n], [2], -4) for n > 0. - Stefano Spezia, May 08 2025

A085364 a(0)=1, for n>0: a(n) = 6*13^(n-1) - (1/2)*Sum_{i=1..n-1} a(i)*a(n-i).

Original entry on oeis.org

1, 6, 60, 654, 7458, 87378, 1042152, 12587730, 153479508, 1885010946, 23285957604, 289018502682, 3601315495050, 45023019250398, 564465885846216, 7094214579174558, 89351097367355826, 1127492973620753010
Offset: 0

Views

Author

Mario Catalani (mario.catalani(AT)unito.it), Jun 25 2003

Keywords

Comments

From G. C. Greubel, May 23 2020: (Start)
This sequence is part of a class of sequences, for m >= 0, with the properties:
a(n) = 2*m*(4*m+1)^(n-1) - (1/2)*Sum_{k=1..n-1} a(k)*a(n-k).
a(n) = Sum_{k=0..n} m^k * binomial(n-1, n-k) * binomial(2*k, k).
n*a(n) = 2*((2*m+1)*n - (m+1))*a(n-1) - (4*m+1)*(n-2)*a(n-2).
a(n) = (2*m) * Hypergeometric2F1(-n+1, 3/2; 2; -4*m), for n > 0.
(4*m + 1)^n = Sum_{k=0..n} Sum_{j=0..k} a(j)*a(k-j).
G.f.: sqrt( (1 - t)/(1 - (4*m+1)*t) ).
This sequence is the case of m=3. (End)

Crossrefs

Cf. A001022 (13^n), A085362, A085363.

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 30); Coefficients(R!( Sqrt((1-x)/(1-13*x)) )); // G. C. Greubel, May 23 2020
    
  • Maple
    seq(coeff(series( sqrt((1-x)/(1-13*x)) , x, n+1), x, n), n = 0..30); # G. C. Greubel, May 23 2020
  • Mathematica
    CoefficientList[Series[Sqrt[(1-x)/(1-13x)], {x, 0, 25}], x]
  • PARI
    my(x='x+O('x^66)); Vec(sqrt((1-x)/(1-13*x))) \\ Joerg Arndt, May 10 2013
    
  • Sage
    def A085362_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( sqrt((1-x)/(1-13*x)) ).list()
    A085362_list(30) # G. C. Greubel, May 23 2020

Formula

G.f.: sqrt((1-x)/(1-13*x))
Sum_{i=0..n} Sum_{j=0..i} a(j)*a(i-j) = 13^n.
D-finite with recurrence: n*a(n) = 2*(7*n-4)*a(n-1) - 13*(n-2)*a(n-2). - Vaclav Kotesovec, Oct 14 2012
a(n) ~ 2*sqrt(3)*13^(n-1/2)/sqrt(Pi*n). - Vaclav Kotesovec, Oct 14 2012
a(0) = 1; a(n) = (6/n) * Sum_{k=0..n-1} (n+k) * a(k). - Seiichi Manyama, Mar 28 2023
From Seiichi Manyama, Aug 22 2025: (Start)
a(n) = (1/4)^n * Sum_{k=0..n} 13^k * binomial(2*k,k) * binomial(2*(n-k),n-k)/(1-2*(n-k)).
a(n) = Sum_{k=0..n} (-3)^k * 13^(n-k) * binomial(2*k,k)/(1-2*k) * binomial(n-1,n-k). (End)
Showing 1-10 of 30 results. Next