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-2 of 2 results.

A003518 a(n) = 8*binomial(2*n+1,n-3)/(n+5).

Original entry on oeis.org

1, 8, 44, 208, 910, 3808, 15504, 62016, 245157, 961400, 3749460, 14567280, 56448210, 218349120, 843621600, 3257112960, 12570420330, 48507033744, 187187399448, 722477682080, 2789279908316, 10772391370048, 41620603020640, 160878516023680, 622147386185325
Offset: 3

Views

Author

Keywords

Comments

a(n-6) is the number of n-th generation nodes in the tree of sequences with unit increase labeled by 7 (cf. Zoran Sunic reference). - Benoit Cloitre, Oct 07 2003
Number of standard tableaux of shape (n+4,n-3). - Emeric Deutsch, May 30 2004

Examples

			G.f. = x^3 + 8*x^4 + 44*x^5 + 208*x^6 + 910*x^7 + 3808*x^8 + 15504*x^9 + ...
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A002057.
First differences are in A026018.
A diagonal of any of the essentially equivalent arrays A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

  • Magma
    [8*Binomial(2*n+1,n-3)/(n+5): n in [3..30]]; // Vincenzo Librandi, Jan 23 2017
  • Mathematica
    Table[8 Binomial[2 n + 1, n - 3]/(n + 5), {n, 3, 25}] (* Michael De Vlieger, Oct 26 2016 *)
    CoefficientList[Series[((1 - Sqrt[1 - 4 x])/(2 x))^8, {x, 0, 30}], x] (* Vincenzo Librandi, Jan 23 2017 *)
  • PARI
    {a(n) = if( n<3, 0, 8 * binomial(2*n + 1, n-3) / (n + 5))}; /* Michael Somos, Mar 14 2011 */
    
  • PARI
    my(x='x+O('x^50)); Vec(x^3*((1-(1-4*x)^(1/2))/(2*x))^8) \\ Altug Alkan, Nov 01 2015
    

Formula

G.f.: x^3*C(x)^8, where C(x)=(1-sqrt(1-4*x))/(2*x) is g.f. for the Catalan numbers (A000108). - Emeric Deutsch, May 30 2004
The convolution of A002057 with itself. - Gerald McGarvey, Nov 08 2007
Let A be the Toeplitz matrix of order n defined by: A[i,i-1]=-1, A[i,j]=Catalan(j-i), (i<=j), and A[i,j]=0, otherwise. Then, for n>=7, a(n-4)=(-1)^(n-7)*coeff(charpoly(A,x),x^7). - Milan Janjic, Jul 08 2010
a(n) = A214292(2*n,n-4) for n > 3. - Reinhard Zumkeller, Jul 12 2012
Integral representation as the n-th moment of the signed weight function W(x) on (0,4), i.e.: a(n+3) = Integral_{x=0..4} x^n*W(x) dx, n >= 0, with W(x) = (1/2)*x^(7/2)*(x-2)*(x^2-4*x+2)*sqrt(4-x)/Pi. - Karol A. Penson, Oct 26 2016
From Ilya Gutkovskiy, Jan 22 2017: (Start)
E.g.f.: 4*BesselI(4,2*x)*exp(2*x)/x.
a(n) ~ 4^(n+2)/(sqrt(Pi)*n^(3/2)). (End)
D-finite with recurrence: -(n+5)*(n-3)*a(n) +2*n*(2*n+1)*a(n-1)=0. - R. J. Mathar, Feb 20 2020
From Amiram Eldar, Jan 02 2022: (Start)
Sum_{n>=3} 1/a(n) = 43*Pi/(36*sqrt(3)) - 81/80.
Sum_{n>=3} (-1)^(n+1)/a(n) = 6213*log(phi)/(50*sqrt(5)) - 10339/400, where phi is the golden ratio (A001622). (End)

Extensions

More terms from Jon E. Schoenfield, May 06 2010

A026009 Triangular array T read by rows: T(n,0) = 1 for n >= 0; T(1,1) = 1; and for n >= 2, T(n,k) = T(n-1,k-1) + T(n-1,k) for k = 1,2,...,[(n+1)/2]; T(n,n/2 + 1) = T(n-1,n/2) if n is even.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 4, 6, 3, 1, 5, 10, 9, 1, 6, 15, 19, 9, 1, 7, 21, 34, 28, 1, 8, 28, 55, 62, 28, 1, 9, 36, 83, 117, 90, 1, 10, 45, 119, 200, 207, 90, 1, 11, 55, 164, 319, 407, 297, 1, 12, 66, 219, 483, 726, 704, 297, 1, 13, 78, 285, 702, 1209, 1430, 1001, 1, 14, 91, 363, 987, 1911, 2639, 2431, 1001
Offset: 0

Views

Author

Keywords

Examples

			From _Jonathon Kirkpatrick_, Jul 01 2016: (Start)
Triangle begins:
  1;
  1,  1;
  1,  2,  1;
  1,  3,  3;
  1,  4,  6,   3;
  1,  5, 10,   9;
  1,  6, 15,  19,   9;
  1,  7, 21,  34,  28;
  1,  8, 28,  55,  62,   28;
  1,  9, 36,  83, 117,   90;
  1, 10, 45, 119, 200,  207,   90;
  1, 11, 55, 164, 319,  407,  297;
  1, 12, 66, 219, 483,  726,  704,  297;
  1, 13, 78, 285, 702, 1209, 1430, 1001;
  ... (End)
		

Crossrefs

Sums involving this sequence: A026010, A027287, A027288, A027289, A027290, A027291, A027292.

Programs

  • Magma
    [1] cat [Binomial(n,k) - Binomial(n,k-3): k in [0..Floor((n+2)/2)], n in [1..15]]; // G. C. Greubel, Mar 18 2021
  • Mathematica
    T[n_, k_]:= Binomial[n, k] - Binomial[n, k-3];
    Join[{1}, Table[T[n, k], {n,14}, {k,0,Floor[(n+2)/2]}]//Flatten] (* G. C. Greubel, Mar 18 2021 *)
  • Sage
    [1]+flatten([[binomial(n,k) - binomial(n,k-3) for k in (0..(n+2)//2)] for n in (1..15)]) # G. C. Greubel, Mar 18 2021
    

Formula

T(n, k) = binomial(n, k) - binomial(n, k-3). - Darko Marinov (marinov(AT)lcs.mit.edu), May 17 2001
Sum_{k=0..floor((n+2)/2)} T(n, k) = A026010(n). - G. C. Greubel, Mar 18 2021
Showing 1-2 of 2 results.