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.

Previous Showing 11-20 of 146 results. Next

A056040 Swinging factorial, a(n) = 2^(n-(n mod 2))*Product_{k=1..n} k^((-1)^(k+1)).

Original entry on oeis.org

1, 1, 2, 6, 6, 30, 20, 140, 70, 630, 252, 2772, 924, 12012, 3432, 51480, 12870, 218790, 48620, 923780, 184756, 3879876, 705432, 16224936, 2704156, 67603900, 10400600, 280816200, 40116600, 1163381400, 155117520, 4808643120, 601080390, 19835652870, 2333606220
Offset: 0

Views

Author

Labos Elemer, Jul 25 2000

Keywords

Comments

a(n) is the number of 'swinging orbitals' which are enumerated by the trinomial n over [floor(n/2), n mod 2, floor(n/2)].
Similar to but different from A001405(n) = binomial(n, floor(n/2)), a(n) = lcm(A001405(n-1), A001405(n)) (for n>0).
A055773(n) divides a(n), A001316(floor(n/2)) divides a(n).
Exactly p consecutive multiples of p follow the least positive multiple of p if p is an odd prime. Compare with the similar property of A100071. - Peter Luschny, Aug 27 2012
a(n) is the number of vertices of the polytope resulting from the intersection of an n-hypercube with the hyperplane perpendicular to and bisecting one of its long diagonals. - Didier Guillet, Jun 11 2018 [Edited by Peter Munn, Dec 06 2022]

Examples

			a(10) = 10!/5!^2 = trinomial(10,[5,0,5]);
a(11) = 11!/5!^2 = trinomial(11,[5,1,5]).
		

Crossrefs

Programs

  • Magma
    [(Factorial(n)/(Factorial(Floor(n/2)))^2): n in [0..40]]; // Vincenzo Librandi, Sep 11 2011
    
  • Maple
    SeriesCoeff := proc(s,n) series(s(w,n),w,n+2);
    convert(%,polynom); coeff(%,w,n) end;
    a1 := proc(n) local k;
    2^(n-(n mod 2))*mul(k^((-1)^(k+1)),k=1..n) end:
    a2 := proc(n) option remember;
    `if`(n=0,1,n^irem(n,2)*(4/n)^irem(n+1,2)*a2(n-1)) end;
    a3 := n -> n!/iquo(n,2)!^2;
    g4 := z -> BesselI(0,2*z)*(1+z);
    a4 := n -> n!*SeriesCoeff(g4,n);
    g5 := z -> (1+z/(1-4*z^2))/sqrt(1-4*z^2);
    a5 := n -> SeriesCoeff(g5,n);
    g6 := (z,n) -> (1+z^2)^n+n*z*(1+z^2)^(n-1);
    a6 := n -> SeriesCoeff(g6,n);
    a7 := n -> combinat[multinomial](n,floor(n/2),n mod 2,floor(n/2));
    h := n -> binomial(n,floor(n/2)); # A001405
    a8 := n -> ilcm(h(n-1),h(n));
    F := [a1, a2, a3, a4, a5, a6, a7, a8];
    for a in F do seq(a(i), i=0..32) od;
  • Mathematica
    f[n_] := 2^(n - Mod[n, 2])*Product[k^((-1)^(k + 1)), {k, n}]; Array[f, 33, 0] (* Robert G. Wilson v, Aug 02 2010 *)
    f[n_] := If[OddQ@n, n*Binomial[n - 1, (n - 1)/2], Binomial[n, n/2]]; Array[f, 33, 0] (* Robert G. Wilson v, Aug 10 2010 *)
    sf[n_] := With[{f = Floor[n/2]}, Pochhammer[f+1, n-f]/f!]; (* or, twice faster: *) sf[n_] := n!/Quotient[n, 2]!^2; Table[sf[n], {n, 0, 32}] (* Jean-François Alcover, Jul 26 2013, updated Feb 11 2015 *)
  • PARI
    a(n)=n!/(n\2)!^2 \\ Charles R Greathouse IV, May 02 2011
    
  • Sage
    def A056040():
        r, n = 1, 0
        while True:
            yield r
            n += 1
            r *= 4/n if is_even(n) else n
    a = A056040(); [next(a) for i in range(36)]  # Peter Luschny, Oct 24 2013

Formula

a(n) = n!/floor(n/2)!^2. [Essentially the original name.]
a(0) = 1, a(n) = n^(n mod 2)*(4/n)^(n+1 mod 2)*a(n-1) for n>=1.
E.g.f.: (1+x)*BesselI(0, 2*x). - Vladeta Jovovic, Jan 19 2004
O.g.f.: a(n) = SeriesCoeff_{n}((1+z/(1-4*z^2))/sqrt(1-4*z^2)).
P.g.f.: a(n) = PolyCoeff_{n}((1+z^2)^n+n*z*(1+z^2)^(n-1)).
a(2n+1) = A046212(2n+1) = A100071(2n+1). - M. F. Hasler, Jan 25 2012
a(2*n) = binomial(2*n,n); a(2*n+1) = (2*n+1)*binomial(2*n,n). Central terms of triangle A211226. - Peter Bala, Apr 10 2012
D-finite with recurrence: n*a(n) + (n-2)*a(n-1) + 4*(-2*n+3)*a(n-2) + 4*(-n+1)*a(n-3) + 16*(n-3)*a(n-4) = 0. - Alexander R. Povolotsky, Aug 17 2012
Sum_{n>=0} 1/a(n) = 4/3 + 8*Pi/(9*sqrt(3)). - Alexander R. Povolotsky, Aug 18 2012
E.g.f.: U(0) where U(k)= 1 + x/(1 - x/(x + (k+1)*(k+1)/U(k+1))); (continued fraction, 3-step). - Sergei N. Gladkovskii, Oct 19 2012
Central column of the coefficients of the swinging polynomials A162246. - Peter Luschny, Oct 22 2013
a(n) = Sum_{k=0..n} A189231(n, 2*k). (Cf. A212303 for the odd case.) - Peter Luschny, Oct 30 2013
a(n) = hypergeometric([-n,-n-1,1/2],[-n-2,1],2)*2^(n-1)*(n+2). - Peter Luschny, Sep 22 2014
a(n) = 4^floor(n/2)*hypergeometric([-floor(n/2), (-1)^n/2], [1], 1). - Peter Luschny, May 19 2015
Sum_{n>=0} (-1)^n/a(n) = 4/3 - 4*Pi/(9*sqrt(3)). - Amiram Eldar, Mar 10 2022

Extensions

Extended and edited by Peter Luschny, Jun 28 2009

A039599 Triangle formed from even-numbered columns of triangle of expansions of powers of x in terms of Chebyshev polynomials U_n(x).

Original entry on oeis.org

1, 1, 1, 2, 3, 1, 5, 9, 5, 1, 14, 28, 20, 7, 1, 42, 90, 75, 35, 9, 1, 132, 297, 275, 154, 54, 11, 1, 429, 1001, 1001, 637, 273, 77, 13, 1, 1430, 3432, 3640, 2548, 1260, 440, 104, 15, 1, 4862, 11934, 13260, 9996, 5508, 2244, 663, 135, 17, 1
Offset: 0

Views

Author

Keywords

Comments

T(n,k) is the number of lattice paths from (0,0) to (n,n) with steps E = (1,0) and N = (0,1) which touch but do not cross the line x - y = k and only situated above this line; example: T(3,2) = 5 because we have EENNNE, EENNEN, EENENN, ENEENN, NEEENN. - Philippe Deléham, May 23 2005
The matrix inverse of this triangle is the triangular matrix T(n,k) = (-1)^(n+k)* A085478(n,k). - Philippe Deléham, May 26 2005
Essentially the same as A050155 except with a leading diagonal A000108 (Catalan numbers) 1, 1, 2, 5, 14, 42, 132, 429, .... - Philippe Deléham, May 31 2005
Number of Grand Dyck paths of semilength n and having k downward returns to the x-axis. (A Grand Dyck path of semilength n is a path in the half-plane x>=0, starting at (0,0), ending at (2n,0) and consisting of steps u=(1,1) and d=(1,-1)). Example: T(3,2)=5 because we have u(d)uud(d),uud(d)u(d),u(d)u(d)du,u(d)duu(d) and duu(d)u(d) (the downward returns to the x-axis are shown between parentheses). - Emeric Deutsch, May 06 2006
Riordan array (c(x),x*c(x)^2) where c(x) is the g.f. of A000108; inverse array is (1/(1+x),x/(1+x)^2). - Philippe Deléham, Feb 12 2007
The triangle may also be generated from M^n*[1,0,0,0,0,0,0,0,...], where M is the infinite tridiagonal matrix with all 1's in the super and subdiagonals and [1,2,2,2,2,2,2,...] in the main diagonal. - Philippe Deléham, Feb 26 2007
Inverse binomial matrix applied to A124733. Binomial matrix applied to A089942. - Philippe Deléham, Feb 26 2007
Number of standard tableaux of shape (n+k,n-k). - Philippe Deléham, Mar 22 2007
From Philippe Deléham, Mar 30 2007: (Start)
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=x*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+y*T(n-1,k)+T(n-1,k+1) for k>=1. Other triangles arise by choosing different values for (x,y):
(0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970
(1,0) -> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877;
(1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598;
(2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954;
(3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791;
(4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. (End)
The table U(n,k) = Sum_{j=0..n} T(n,j)*k^j is given in A098474. - Philippe Deléham, Mar 29 2007
Sequence read mod 2 gives A127872. - Philippe Deléham, Apr 12 2007
Number of 2n step walks from (0,0) to (2n,2k) and consisting of step u=(1,1) and d=(1,-1) and the path stays in the nonnegative quadrant. Example: T(3,0)=5 because we have uuuddd, uududd, ududud, uduudd, uuddud; T(3,1)=9 because we have uuuudd, uuuddu, uuudud, ududuu, uuduud, uduudu, uudduu, uduuud, uududu; T(3,2)=5 because we have uuuuud, uuuudu, uuuduu, uuduuu, uduuuu; T(3,3)=1 because we have uuuuuu. - Philippe Deléham, Apr 16 2007, Apr 17 2007, Apr 18 2007
Triangular matrix, read by rows, equal to the matrix inverse of triangle A129818. - Philippe Deléham, Jun 19 2007
Let Sum_{n>=0} a(n)*x^n = (1+x)/(1-mx+x^2) = o.g.f. of A_m, then Sum_{k=0..n} T(n,k)*a(k) = (m+2)^n. Related expansions of A_m are: A099493, A033999, A057078, A057077, A057079, A005408, A002878, A001834, A030221, A002315, A033890, A057080, A057081, A054320, A097783, A077416, A126866, A028230, A161591, for m=-3,-2,-1,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15, respectively. - Philippe Deléham, Nov 16 2009
The Kn11, Kn12, Fi1 and Fi2 triangle sums link the triangle given above with three sequences; see the crossrefs. For the definitions of these triangle sums, see A180662. - Johannes W. Meijer, Apr 20 2011
4^n = (n-th row terms) dot (first n+1 odd integer terms). Example: 4^4 = 256 = (14, 28, 20, 7, 1) dot (1, 3, 5, 7, 9) = (14 + 84 + 100 + 49 + 9) = 256. - Gary W. Adamson, Jun 13 2011
The linear system of n equations with coefficients defined by the first n rows solve for diagonal lengths of regular polygons with N= 2n+1 edges; the constants c^0, c^1, c^2, ... are on the right hand side, where c = 2 + 2*cos(2*Pi/N). Example: take the first 4 rows relating to the 9-gon (nonagon), N = 2*4 + 1; with c = 2 + 2*cos(2*Pi/9) = 3.5320888.... The equations are (1,0,0,0) = 1; (1,1,0,0) = c; (2,3,1,0) = c^2; (5,9,5,1) = c^3. The solutions are 1, 2.53208..., 2.87938..., and 1.87938...; the four distinct diagonal lengths of the 9-gon (nonagon) with edge = 1. (Cf. comment in A089942 which uses the analogous operations but with c = 1 + 2*cos(2*Pi/9).) - Gary W. Adamson, Sep 21 2011
Also called the Lobb numbers, after Andrew Lobb, are a natural generalization of the Catalan numbers, given by L(m,n)=(2m+1)*Binomial(2n,m+n)/(m+n+1), where n >= m >= 0. For m=0, we get the n-th Catalan number. See added reference. - Jayanta Basu, Apr 30 2013
From Wolfdieter Lang, Sep 20 2013: (Start)
T(n, k) = A053121(2*n, 2*k). T(n, k) appears in the formula for the (2*n)-th power of the algebraic number rho(N):= 2*cos(Pi/N) = R(N, 2) in terms of the odd-indexed diagonal/side length ratios R(N, 2*k+1) = S(2*k, rho(N)) in the regular N-gon inscribed in the unit circle (length unit 1). S(n, x) are Chebyshev's S polynomials (see A049310):
rho(N)^(2*n) = Sum_{k=0..n} T(n, k)*R(N, 2*k+1), n >= 0, identical in N > = 1. For a proof see the Sep 21 2013 comment under A053121. Note that this is the unreduced version if R(N, j) with j > delta(N), the degree of the algebraic number rho(N) (see A055034), appears.
For the odd powers of rho(n) see A039598. (End)
Unsigned coefficients of polynomial numerators of Eqn. 2.1 of the Chakravarty and Kodama paper, defining the polynomials of A067311. - Tom Copeland, May 26 2016
The triangle is the Riordan square of the Catalan numbers in the sense of A321620. - Peter Luschny, Feb 14 2023

Examples

			Triangle T(n, k) begins:
  n\k     0     1     2     3     4     5    6   7   8  9
  0:      1
  1:      1     1
  2:      2     3     1
  3:      5     9     5     1
  4:     14    28    20     7     1
  5:     42    90    75    35     9     1
  6:    132   297   275   154    54    11    1
  7:    429  1001  1001   637   273    77   13   1
  8:   1430  3432  3640  2548  1260   440  104  15   1
  9:   4862 11934 13260  9996  5508  2244  663 135  17  1
  ... Reformatted by _Wolfdieter Lang_, Dec 21 2015
From _Paul Barry_, Feb 17 2011: (Start)
Production matrix begins
  1, 1,
  1, 2, 1,
  0, 1, 2, 1,
  0, 0, 1, 2, 1,
  0, 0, 0, 1, 2, 1,
  0, 0, 0, 0, 1, 2, 1,
  0, 0, 0, 0, 0, 1, 2, 1 (End)
From _Wolfdieter Lang_, Sep 20 2013: (Start)
Example for rho(N) = 2*cos(Pi/N) powers:
n=2: rho(N)^4 = 2*R(N,1) + 3*R(N,3) + 1*R(N, 5) =
  2 + 3*S(2, rho(N)) + 1*S(4, rho(N)), identical in N >= 1. For N=4 (the square with only one distinct diagonal), the degree delta(4) = 2, hence R(4, 3) and R(4, 5) can be reduced, namely to R(4, 1) = 1 and R(4, 5) = -R(4,1) = -1, respectively. Therefore, rho(4)^4 =(2*cos(Pi/4))^4 = 2 + 3 -1 = 4. (End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 796.
  • T. Myers and L. Shapiro, Some applications of the sequence 1, 5, 22, 93, 386, ... to Dyck paths and ordered trees, Congressus Numerant., 204 (2010), 93-104.

Crossrefs

Row sums: A000984.
Triangle sums (see the comments): A000958 (Kn11), A001558 (Kn12), A088218 (Fi1, Fi2).

Programs

  • Magma
    /* As triangle */ [[Binomial(2*n, k+n)*(2*k+1)/(k+n+1): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Oct 16 2015
    
  • Maple
    T:=(n,k)->(2*k+1)*binomial(2*n,n-k)/(n+k+1): for n from 0 to 12 do seq(T(n,k),k=0..n) od; # yields sequence in triangular form # Emeric Deutsch, May 06 2006
    T := proc(n, k) option remember; if k = n then 1 elif k > n then 0 elif k = 0 then T(n-1, 0) + T(n-1,1) else T(n-1, k-1) + 2*T(n-1, k) + T(n-1, k+1) fi end:
    seq(seq(T(n, k), k = 0..n), n = 0..9) od; # Peter Luschny, Feb 14 2023
  • Mathematica
    Table[Abs[Differences[Table[Binomial[2 n, n + i], {i, 0, n + 1}]]], {n, 0,7}] // Flatten (* Geoffrey Critzer, Dec 18 2011 *)
    Join[{1},Flatten[Table[Binomial[2n-1,n-k]-Binomial[2n-1,n-k-2],{n,10},{k,0,n}]]] (* Harvey P. Dale, Dec 18 2011 *)
    Flatten[Table[Binomial[2*n,m+n]*(2*m+1)/(m+n+1),{n,0,9},{m,0,n}]] (* Jayanta Basu, Apr 30 2013 *)
  • PARI
    a(n, k) = (2*n+1)/(n+k+1)*binomial(2*k, n+k)
    trianglerows(n) = for(x=0, n-1, for(y=0, x, print1(a(y, x), ", ")); print(""))
    trianglerows(10) \\ Felix Fröhlich, Jun 24 2016
  • Sage
    # Algorithm of L. Seidel (1877)
    # Prints the first n rows of the triangle
    def A039599_triangle(n) :
        D = [0]*(n+2); D[1] = 1
        b = True ; h = 1
        for i in range(2*n-1) :
            if b :
                for k in range(h,0,-1) : D[k] += D[k-1]
                h += 1
            else :
                for k in range(1,h, 1) : D[k] += D[k+1]
            if b : print([D[z] for z in (1..h-1)])
            b = not b
    A039599_triangle(10)  # Peter Luschny, May 01 2012
    

Formula

T(n,k) = C(2*n-1, n-k) - C(2*n-1, n-k-2), n >= 1, T(0,0) = 1.
From Emeric Deutsch, May 06 2006: (Start)
T(n,k) = (2*k+1)*binomial(2*n,n-k)/(n+k+1).
G.f.: G(t,z)=1/(1-(1+t)*z*C), where C=(1-sqrt(1-4*z))/(2*z) is the Catalan function. (End)
The following formulas were added by Philippe Deléham during 2003 to 2009: (Start)
Triangle T(n, k) read by rows; given by A000012 DELTA A000007, where DELTA is Deléham's operator defined in A084938.
T(n, k) = C(2*n, n-k)*(2*k+1)/(n+k+1). Sum(k>=0; T(n, k)*T(m, k) = A000108(n+m)); A000108: numbers of Catalan.
T(n, 0) = A000108(n); T(n, k) = 0 if k>n; for k>0, T(n, k) = Sum_{j=1..n} T(n-j, k-1)*A000108(j).
T(n, k) = A009766(n+k, n-k) = A033184(n+k+1, 2k+1).
G.f. for column k: Sum_{n>=0} T(n, k)*x^n = x^k*C(x)^(2*k+1) where C(x) = Sum_{n>=0} A000108(n)*x^n is g.f. for Catalan numbers, A000108.
T(0, 0) = 1, T(n, k) = 0 if n<0 or n=1, T(n, k) = T(n-1, k-1) + 2*T(n-1, k) + T(n-1, k+1).
a(n) + a(n+1) = 1 + A000108(m+1) if n = m*(m+3)/2; a(n) + a(n+1) = A039598(n) otherwise.
T(n, k) = A050165(n, n-k).
Sum_{j>=0} T(n-k, j)*A039598(k, j) = A028364(n, k).
Matrix inverse of the triangle T(n, k) = (-1)^(n+k)*binomial(n+k, 2*k) = (-1)^(n+k)*A085478(n, k).
Sum_{k=0..n} T(n, k)*x^k = A000108(n), A000984(n), A007854(n), A076035(n), A076036(n) for x = 0, 1, 2, 3, 4.
Sum_{k=0..n} (2*k+1)*T(n, k) = 4^n.
T(n, k)*(-2)^(n-k) = A114193(n, k).
Sum_{k>=h} T(n,k) = binomial(2n,n-h).
Sum_{k=0..n} T(n,k)*5^k = A127628(n).
Sum_{k=0..n} T(n,k)*7^k = A115970(n).
T(n,k) = Sum_{j=0..n-k} A106566(n+k,2*k+j).
Sum_{k=0..n} T(n,k)*6^k = A126694(n).
Sum_{k=0..n} T(n,k)*A000108(k) = A007852(n+1).
Sum_{k=0..floor(n/2)} T(n-k,k) = A000958(n+1).
Sum_{k=0..n} T(n,k)*(-1)^k = A000007(n).
Sum_{k=0..n} T(n,k)*(-2)^k = (-1)^n*A064310(n).
T(2*n,n) = A126596(n).
Sum_{k=0..n} T(n,k)*(-x)^k = A000007(n), A126983(n), A126984(n), A126982(n), A126986(n), A126987(n), A127017(n), A127016(n), A126985(n), A127053(n) for x=1,2,3,4,5,6,7,8,9,10 respectively.
Sum_{j>=0} T(n,j)*binomial(j,k) = A116395(n,k).
T(n,k) = Sum_{j>=0} A106566(n,j)*binomial(j,k).
T(n,k) = Sum_{j>=0} A127543(n,j)*A038207(j,k).
Sum_{k=0..floor(n/2)} T(n-k,k)*A000108(k) = A101490(n+1).
T(n,k) = A053121(2*n,2*k).
Sum_{k=0..n} T(n,k)*sin((2*k+1)*x) = sin(x)*(2*cos(x))^(2*n).
T(n,n-k) = Sum_{j>=0} (-1)^(n-j)*A094385(n,j)*binomial(j,k).
Sum_{j>=0} A110506(n,j)*binomial(j,k) = Sum_{j>=0} A110510(n,j)*A038207(j,k) = T(n,k)*2^(n-k).
Sum_{j>=0} A110518(n,j)*A027465(j,k) = Sum_{j>=0} A110519(n,j)*A038207(j,k) = T(n,k)*3^(n-k).
Sum_{k=0..n} T(n,k)*A001045(k) = A049027(n), for n>=1.
Sum_{k=0..n} T(n,k)*a(k) = (m+2)^n if Sum_{k>=0} a(k)*x^k = (1+x)/(x^2-m*x+1).
Sum_{k=0..n} T(n,k)*A040000(k) = A001700(n).
Sum_{k=0..n} T(n,k)*A122553(k) = A051924(n+1).
Sum_{k=0..n} T(n,k)*A123932(k) = A051944(n).
Sum_{k=0..n} T(n,k)*k^2 = A000531(n), for n>=1.
Sum_{k=0..n} T(n,k)*A000217(k) = A002457(n-1), for n>=1.
Sum{j>=0} binomial(n,j)*T(j,k)= A124733(n,k).
Sum_{k=0..n} T(n,k)*x^(n-k) = A000012(n), A000984(n), A089022(n), A035610(n), A130976(n), A130977(n), A130978(n), A130979(n), A130980(n), A131521(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively.
Sum_{k=0..n} T(n,k)*A005043(k) = A127632(n).
Sum_{k=0..n} T(n,k)*A132262(k) = A089022(n).
T(n,k) + T(n,k+1) = A039598(n,k).
T(n,k) = A128899(n,k)+A128899(n,k+1).
Sum_{k=0..n} T(n,k)*A015518(k) = A076025(n), for n>=1. Also Sum_{k=0..n} T(n,k)*A015521(k) = A076026(n), for n>=1.
Sum_{k=0..n} T(n,k)*(-1)^k*x^(n-k) = A033999(n), A000007(n), A064062(n), A110520(n), A132863(n), A132864(n), A132865(n), A132866(n), A132867(n), A132869(n), A132897(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 respectively.
Sum_{k=0..n} T(n,k)*(-1)^(k+1)*A000045(k) = A109262(n), A000045:= Fibonacci numbers.
Sum_{k=0..n} T(n,k)*A000035(k)*A016116(k) = A143464(n).
Sum_{k=0..n} T(n,k)*A016116(k) = A101850(n).
Sum_{k=0..n} T(n,k)*A010684(k) = A100320(n).
Sum_{k=0..n} T(n,k)*A000034(k) = A029651(n).
Sum_{k=0..n} T(n,k)*A010686(k) = A144706(n).
Sum_{k=0..n} T(n,k)*A006130(k-1) = A143646(n), with A006130(-1)=0.
T(n,2*k)+T(n,2*k+1) = A118919(n,k).
Sum_{k=0..j} T(n,k) = A050157(n,j).
Sum_{k=0..2} T(n,k) = A026012(n); Sum_{k=0..3} T(n,k)=A026029(n).
Sum_{k=0..n} T(n,k)*A000045(k+2) = A026671(n).
Sum_{k=0..n} T(n,k)*A000045(k+1) = A026726(n).
Sum_{k=0..n} T(n,k)*A057078(k) = A000012(n).
Sum_{k=0..n} T(n,k)*A108411(k) = A155084(n).
Sum_{k=0..n} T(n,k)*A057077(k) = 2^n = A000079(n).
Sum_{k=0..n} T(n,k)*A057079(k) = 3^n = A000244(n).
Sum_{k=0..n} T(n,k)*(-1)^k*A011782(k) = A000957(n+1).
(End)
T(n,k) = Sum_{j=0..k} binomial(k+j,2j)*(-1)^(k-j)*A000108(n+j). - Paul Barry, Feb 17 2011
Sum_{k=0..n} T(n,k)*A071679(k+1) = A026674(n+1). - Philippe Deléham, Feb 01 2014
Sum_{k=0..n} T(n,k)*(2*k+1)^2 = (4*n+1)*binomial(2*n,n). - Werner Schulte, Jul 22 2015
Sum_{k=0..n} T(n,k)*(2*k+1)^3 = (6*n+1)*4^n. - Werner Schulte, Jul 22 2015
Sum_{k=0..n} (-1)^k*T(n,k)*(2*k+1)^(2*m) = 0 for 0 <= m < n (see also A160562). - Werner Schulte, Dec 03 2015
T(n,k) = GegenbauerC(n-k,-n+1,-1) - GegenbauerC(n-k-1,-n+1,-1). - Peter Luschny, May 13 2016
T(n,n-2) = A014107(n). - R. J. Mathar, Jan 30 2019
T(n,n-3) = n*(2*n-1)*(2*n-5)/3. - R. J. Mathar, Jan 30 2019
T(n,n-4) = n*(n-1)*(2*n-1)*(2*n-7)/6. - R. J. Mathar, Jan 30 2019
T(n,n-5) = n*(n-1)*(2*n-1)*(2*n-3)*(2*n-9)/30. - R. J. Mathar, Jan 30 2019

Extensions

Corrected by Philippe Deléham, Nov 26 2009, Dec 14 2009

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

Original entry on oeis.org

1, 3, 9, 29, 99, 351, 1275, 4707, 17577, 66197, 250953, 956385, 3660541, 14061141, 54177741, 209295261, 810375651, 3143981871, 12219117171, 47564380971, 185410909791, 723668784231, 2827767747951, 11061198475551, 43308802158651, 169719408596403, 665637941544507
Offset: 0

Views

Author

Keywords

Comments

The expression a(n) = B^n*Sum_{ k=0..n } binomial(2*k,k)/B^k gives A006134 for B=1, A082590 (B=2), A132310 (B=3), A002457 (B=4), A144635 (B=5). - N. J. A. Sloane, Jan 21 2009
T(n+1,1) from table A045912 of characteristic polynomial of negative Pascal matrix. - Michael Somos, Jul 24 2002
p divides a((p-3)/2) for p=11, 13, 23, 37, 47, 59, 61, 71, 73, 83, 97, 107, 109, 131, 157, 167, ...: A097933. Also primes congruent to {1, 2, 3, 11} mod 12 or primes p such that 3 is a square mod p (excluding 2 and 3) A038874. - Alexander Adamchuk, Jul 05 2006
Partial sums of the even central binomial coefficients. For p prime >=5, a(p-1) = 1 or -1 (mod p) according as p = 1 or -1 (mod 3) (see Pan and Sun link). - David Callan, Nov 29 2007
First column of triangle A187887. - Michel Marcus, Jun 23 2013
From Gus Wiseman, Apr 20 2023: (Start)
Also the number of nonempty subsets of {1,...,2n+1} with median n+1, where the median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length). The odd/even-length cases are A000984 and A006134(n-1). For example, the a(0) = 1 through a(2) = 9 subsets are:
{1} {2} {3}
{1,3} {1,5}
{1,2,3} {2,4}
{1,3,4}
{1,3,5}
{2,3,4}
{2,3,5}
{1,2,4,5}
{1,2,3,4,5}
Alternatively, a(n-1) is the number of nonempty subsets of {1,...,2n-1} with median n.
(End)

Examples

			1 + 3*x + 9*x^2 + 29*x^3 + 99*x^4 + 351*x^5 + 1275*x^6 + 4707*x^7 + 17577*x^8 + ...
		

References

  • Marko Petkovsek, Herbert Wilf and Doron Zeilberger, A=B, A K Peters, 1996, p. 22.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000984 (first differences), A097933, A038874, A132310.
Equals A066796 + 1.
Odd bisection of A100066.
Row sums of A361654 (also column k = 2).
A007318 counts subsets by length, A231147 by median, A013580 by integer median.
A359893 and A359901 count partitions by median.

Programs

  • MATLAB
    n=10; x=pascal(n); trace(x)
    
  • Magma
    &cat[ [&+[ Binomial(2*k, k): k in [0..n]]]: n in [0..30]]; // Vincenzo Librandi, Aug 13 2015
  • Maple
    A006134 := proc(n) sum(binomial(2*k,k),k=0..n); end;
    a := n -> -binomial(2*(n+1),n+1)*hypergeom([1,n+3/2],[n+2], 4) - I/sqrt(3):
    seq(simplify(a(n)), n=0..24); # Peter Luschny, Oct 29 2015
    # third program:
    A006134 := series(exp(2*x)*BesselI(0, 2*x) + exp(x)*int(BesselI(0, 2*x)*exp(x), x), x = 0, 25):
    seq(n!*coeff(A006134, x, n), n=0..24); # Mélika Tebni, Feb 27 2024
  • Mathematica
    Table[Sum[((2k)!/(k!)^2),{k,0,n}], {n,0,50}] (* Alexander Adamchuk, Jul 05 2006 *)
    a[ n_] := (4/3) Binomial[ 2 n, n] Hypergeometric2F1[ 1/2, 1, -n + 1/2, -1/3] (* Michael Somos, Jun 20 2012 *)
    Accumulate[Table[Binomial[2n,n],{n,0,30}]] (* Harvey P. Dale, Jan 11 2015 *)
    CoefficientList[Series[1/((1 - x) Sqrt[1 - 4 x]), {x, 0, 33}], x] (* Vincenzo Librandi, Aug 13 2015 *)
  • Maxima
    makelist(sum(binomial(2*k,k),k,0,n),n,0,12); /* Emanuele Munarini, Mar 15 2011 */
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( charpoly( matrix( n+1, n+1, i, j, -binomial( i+j-2, i-1))), 1))} \\ Michael Somos, Jul 10 2002
    
  • PARI
    {a(n)=binomial(2*n,n)*sum(k=0,2*n,(-1)^k*polcoeff((1+x+x^2)^n,k)/binomial(2*n,k))} \\ Paul D. Hanna, Aug 21 2007
    
  • PARI
    my(x='x+O('x^100)); Vec(1/((1-x)*sqrt(1-4*x))) \\ Altug Alkan, Oct 29 2015
    

Formula

From Alexander Adamchuk, Jul 05 2006: (Start)
a(n) = Sum_{k=0..n} (2k)!/(k!)^2.
a(n) = A066796(n) + 1, n>0. (End)
G.f.: 1/((1-x)*sqrt(1-4*x)).
D-finite with recurrence: (n+2)*a(n+2) - (5*n+8)*a(n+1) + 2*(2*n+3)*a(n) = 0. - Emanuele Munarini, Mar 15 2011
a(n) = C(2n,n) * Sum_{k=0..2n} (-1)^k*trinomial(n,k)/C(2n,k) where trinomial(n,k) = [x^k] (1 + x + x^2)^n. E.g. a(2) = C(4,2)*(1/1 - 2/4 + 3/6 - 2/4 + 1/1) = 6*(3/2) = 9 ; a(3) = C(6,3)*(1/1 - 3/6 + 6/15 - 7/20 + 6/15 - 3/6 + 1/1) = 20*(29/20) = 29. - Paul D. Hanna, Aug 21 2007
From Alzhekeyev Ascar M, Jan 19 2012: (Start)
a(n) = Sum_{ k=0..n } b(k)*binomial(n+k,k), where b(k)=0 for n-k == 2 (mod 3), b(k)=1 for n-k == 0 or 1 (mod 6), and b(k)=-1 for n-k== 3 or 4 (mod 6).
a(n) = Sum_{ k=0..n-1 } c(k)*binomial(2n,k) + binomial(2n,n), where c(k)=0 for n-k == 0 (mod 3), c(k)=1 for n-k== 1 (mod 3), and c(k)=-1 for n-k==2 (mod 3). (End)
a(n) ~ 2^(2*n+2)/(3*sqrt(Pi*n)). - Vaclav Kotesovec, Nov 06 2012
G.f.: G(0)/2/(1-x), where G(k)= 1 + 1/(1 - 2*x*(2*k+1)/(2*x*(2*k+1) + (k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 24 2013
G.f.: G(0)/(1-x), where G(k)= 1 + 4*x*(4*k+1)/( (4*k+2) - x*(4*k+2)*(4*k+3)/(x*(4*k+3) + (k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 26 2013
a(n) = Sum_{k = 0..n} binomial(n+1,k+1)*A002426(k). - Peter Bala, Oct 29 2015
a(n) = -binomial(2*(n+1),n+1)*hypergeom([1,n+3/2],[n+2], 4) - i/sqrt(3). - Peter Luschny, Oct 29 2015
a(n) = binomial(2*n, n)*hypergeom([1,-n], [1/2-n], 1/4). - Peter Luschny, Mar 16 2016
From Gus Wiseman, Apr 20 2023: (Start)
a(n+1) - a(n) = A000984(n).
a(n) = A013580(2n+1,n+1) (conjectured).
a(n) = 2*A024718(n) - 1.
a(n) = A100066(2n+1).
a(n) = A231147(2n+1,n+1) (conjectured). (End)
a(n) = Sum_{k=0..floor(n/3)} 3^(n-3*k) * binomial(n-k,2*k) * binomial(2*k,k) (Sawhney, 2017). - Amiram Eldar, Feb 24 2024
From Mélika Tebni, Feb 27 2024: (Start)
Limit_{n -> oo} a(n) / A281593(n) = 2.
E.g.f.: exp(2*x)*BesselI(0,2*x) + exp(x)*integral( BesselI(0,2*x)*exp(x) ) dx. (End)
a(n) = [(x*y)^n] 1/((1 - (x + y))*(1 - x*y)). - Stefano Spezia, Feb 16 2025
a(n) = Sum_{k=0..floor(n/2)} (-1)^k*binomial(2*n+1-k, n-2*k). - Michael Weselcouch, Jun 17 2025
a(n) = binomial(1+2*n, n)*hypergeom([1, (1-n)/2, -n/2], [-1-2*n, 2+n], 4). - Stefano Spezia, Jun 18 2025

Extensions

Simpler definition from Alexander Adamchuk, Jul 05 2006

A039598 Triangle formed from odd-numbered columns of triangle of expansions of powers of x in terms of Chebyshev polynomials U_n(x). Sometimes called Catalan's triangle.

Original entry on oeis.org

1, 2, 1, 5, 4, 1, 14, 14, 6, 1, 42, 48, 27, 8, 1, 132, 165, 110, 44, 10, 1, 429, 572, 429, 208, 65, 12, 1, 1430, 2002, 1638, 910, 350, 90, 14, 1, 4862, 7072, 6188, 3808, 1700, 544, 119, 16, 1, 16796, 25194, 23256, 15504, 7752, 2907, 798, 152, 18, 1
Offset: 0

Views

Author

Keywords

Comments

T(n,k) is the number of leaves at level k+1 in all ordered trees with n+1 edges. - Emeric Deutsch, Jan 15 2005
Riordan array ((1-2x-sqrt(1-4x))/(2x^2),(1-2x-sqrt(1-4x))/(2x)). Inverse array is A053122. - Paul Barry, Mar 17 2005
T(n,k) is the number of walks of n steps, each in direction N, S, W, or E, starting at the origin, remaining in the upper half-plane and ending at height k (see the R. K. Guy reference, p. 5). Example: T(3,2)=6 because we have ENN, WNN, NEN, NWN, NNE and NNW. - Emeric Deutsch, Apr 15 2005
Triangle T(n,k), 0<=k<=n, read by rows given by T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0) = 2*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + 2*T(n-1,k) + T(n-1,k+1) for k>=1. - Philippe Deléham, Mar 30 2007
Number of (2n+1)-step walks from (0,0) to (2n+1,2k+1) and consisting of steps u=(1,1) and d=(1,-1) in which the path stays in the nonnegative quadrant. Examples: T(2,0)=5 because we have uuudd, uudud, uuddu, uduud, ududu; T(2,1)=4 because we have uuuud, uuudu, uuduu, uduuu; T(2,2)=1 because we have uuuuu. - Philippe Deléham, Apr 16 2007, Apr 18 2007
Triangle read by rows: T(n,k)=number of lattice paths from (0,0) to (n,k) that do not go below the line y=0 and consist of steps U=(1,1), D=(1,-1) and two types of steps H=(1,0); example: T(3,1)=14 because we have UDU, UUD, 4 HHU paths, 4 HUH paths and 4 UHH paths. - Philippe Deléham, Sep 25 2007
This triangle belongs to the family of triangles defined by T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0) = x*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + y*T(n-1,k) + T(n-1,k+1) for k>=1. Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0) -> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
With offset [1,1] this is the (ordinary) convolution triangle a(n,m) with o.g.f. of column m given by (c(x)-1)^m, where c(x) is the o.g.f. for Catalan numbers A000108. See the Riordan comment by Paul Barry.
T(n, k) is also the number of order-preserving full transformations (of an n-chain) with exactly k fixed points. - Abdullahi Umar, Oct 02 2008
T(n,k)/2^(2n+1) = coefficients of the maximally flat lowpass digital differentiator of the order N=2n+3. - Pavel Holoborodko (pavel(AT)holoborodko.com), Dec 19 2008
The signed triangle S(n,k) := (-1)^(n-k)*T(n,k) provides the transformation matrix between f(n,l) := L(2*l)*5^n*F(2*l)^(2*n+1) (F=Fibonacci numbers A000045, L=Lucas numbers A000032) and F(4*l*(k+1)), k = 0, ..., n, for each l>=0: f(n,l) = Sum_{k=0..n} S(n,k)*F(4*l*(k+1)), n>=0, l>=0. Proof: the o.g.f. of the l.h.s., G(l;x) := Sum_{n>=0} f(n,l)*x^n = F(4*l)/(1 - 5*F(2*l)^2*x) is shown to match the o.g.f. of the r.h.s.: after an interchange of the n- and k-summation, the Riordan property of S = (C(x)/x,C(x)) (compare with the above comments by Paul Barry), with C(x) := 1 - c(-x), with the o.g.f. c(x) of A000108 (Catalan numbers), is used, to obtain, after an index shift, first Sum_{k>=0} F(4*l*(k))*GS(k;x), with the o.g.f of column k of triangle S which is GS(k;x) := Sum_{n>=k} S(n,k)*x^n = C(x)^(k+1)/x. The result is GF(l;C(x))/x with the o.g.f. GF(l,x) := Sum_{k>=0} F(4*l*k)*x^k = x*F(4*l)/(1-L(4*l)*x+x^2) (see a comment on A049670, and A028412). If one uses then the identity L(4*n) - 5*F(2*n)^2 = 2 (in Koshy's book [reference under A065563] this is No. 15, p. 88, attributed to Lucas, 1876), the proof that one recovers the o.g.f. of the l.h.s. from above boils down to a trivial identity on the Catalan o.g.f., namely 1/c^2(-x) = 1 + 2*x - (x*c(-x))^2. - Wolfdieter Lang, Aug 27 2012
O.g.f. for row polynomials R(x) := Sum_{k=0..n} a(n,k)*x^k:
((1+x) - C(z))/(x - (1+x)^2*z) with C the o.g.f. of A000108 (Catalan numbers). From Riordan ((C(x)-1)/x,C(x)-1), compare with a Paul Barry comment above. This coincides with the o.g.f. given by Emeric Deutsch in the formula section. - Wolfdieter Lang, Nov 13 2012
The A-sequence for this Riordan triangle is [1,2,1] and the Z-sequence is [2,1]. See a W. Lang link under A006232 with details and references. - Wolfdieter Lang, Nov 13 2012
From Wolfdieter Lang, Sep 20 2013: (Start)
T(n, k) = A053121(2*n+1, 2*k+1). T(n, k) appears in the formula for the (2*n+1)-th power of the algebraic number rho(N) := 2*cos(Pi/N) = R(N, 2) in terms of the even-indexed diagonal/side length ratios R(N, 2*(k+1)) = S(2*k+1, rho(N)) in the regular N-gon inscribed in the unit circle (length unit 1). S(n, x) are Chebyshev's S polynomials (see A049310): rho(N)^(2*n+1) = Sum_{k=0..n} T(n, k)*R(N, 2*(k+1)), n >= 0, identical in N >= 1. For a proof see the Sep 21 2013 comment under A053121. Note that this is the unreduced version if R(N, j) with j > delta(N), the degree of the algebraic number rho(N) (see A055034), appears. For the even powers of rho(n) see A039599. (End)
The tridiagonal Toeplitz production matrix P in the Example section corresponds to the unsigned Cartan matrix for the simple Lie algebra A_n as n tends to infinity (cf. Damianou ref. in A053122). - Tom Copeland, Dec 11 2015 (revised Dec 28 2015)
T(n,k) is the number of pairs of non-intersecting walks of n steps, each in direction N or E, starting at the origin, and such that the end points of the two paths are separated by a horizontal distance of k. See Shapiro 1976. - Peter Bala, Apr 12 2017
Also the convolution triangle of the Catalan numbers A000108. - Peter Luschny, Oct 07 2022

Examples

			Triangle T(n,k) starts:
n\k     0      1      2      3      4     5    6    7   8  9 10
0:      1
1:      2      1
2:      5      4      1
3:     14     14      6      1
4:     42     48     27      8      1
5:    132    165    110     44     10     1
6:    429    572    429    208     65    12    1
7:   1430   2002   1638    910    350    90   14    1
8:   4862   7072   6188   3808   1700   544  119   16   1
9:  16796  25194  23256  15504   7752  2907  798  152  18  1
10: 58786  90440  87210  62016  33915 14364 4655 1120 189 20  1
... Reformatted and extended by _Wolfdieter Lang_, Nov 13 2012.
Production matrix begins:
2, 1
1, 2, 1
0, 1, 2, 1
0, 0, 1, 2, 1
0, 0, 0, 1, 2, 1
0, 0, 0, 0, 1, 2, 1
0, 0, 0, 0, 0, 1, 2, 1
0, 0, 0, 0, 0, 0, 1, 2, 1
- _Philippe Deléham_, Nov 07 2011
From _Wolfdieter Lang_, Nov 13 2012: (Start)
Recurrence: T(5,1) = 165 = 1*42 + 2*48 +1*27. The Riordan A-sequence is [1,2,1].
Recurrence from Riordan Z-sequence [2,1]: T(5,0) = 132 = 2*42 + 1*48. (End)
From _Wolfdieter Lang_, Sep 20 2013: (Start)
  Example for rho(N) = 2*cos(Pi/N) powers:
  n=2: rho(N)^5 = 5*R(N, 2) + 4*R(N, 4) + 1*R(N, 6) = 5*S(1, rho(N)) + 4*S(3, rho(N)) + 1*S(5, rho(N)), identical in N >= 1. For N=5 (the pentagon with only one distinct diagonal) the degree delta(5) = 2, hence R(5, 4) and R(5, 6) can be reduced, namely to R(5, 1) = 1 and R(5, 6) = -R(5,1) = -1, respectively. Thus rho(5)^5 = 5*R(N, 2) + 4*1  + 1*(-1) = 3 + 5*R(N, 2) = 3 + 5*rho(5), with the golden section rho(5). (End)
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 796.
  • B. A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8.

Crossrefs

Mirror image of A050166. Row sums are A001700.

Programs

  • Magma
    /* As triangle: */ [[Binomial(2*n,n-k) - Binomial(2*n,n-k-2): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Jul 22 2015
    
  • Maple
    T:=(n,k)->binomial(2*n, n-k) - binomial(2*n, n-k-2); # N. J. A. Sloane, Aug 26 2013
    # Uses function PMatrix from A357368. Adds row and column above and to the left.
    PMatrix(10, n -> binomial(2*n, n) / (n + 1)); # Peter Luschny, Oct 07 2022
  • Mathematica
    Flatten[Table[Binomial[2n, n-k] - Binomial[2n, n-k-2], {n,0,9}, {k,0,n}]] (* Jean-François Alcover, May 03 2011 *)
  • PARI
    T(n,k)=binomial(2*n,n-k) - binomial(2*n,n-k-2) \\ Charles R Greathouse IV, Nov 07 2016
  • Sage
    # Algorithm of L. Seidel (1877)
    # Prints the first n rows of the triangle.
    def A039598_triangle(n) :
        D = [0]*(n+2); D[1] = 1
        b = True; h = 1
        for i in range(2*n) :
            if b :
                for k in range(h,0,-1) : D[k] += D[k-1]
                h += 1
            else :
                for k in range(1,h, 1) : D[k] += D[k+1]
            b = not b
            if b : print([D[z] for z in (1..h-1) ])
    A039598_triangle(10)  # Peter Luschny, May 01 2012
    

Formula

Row n: C(2n, n-k) - C(2n, n-k-2).
a(n, k) = C(2n+1, n-k)*2*(k+1)/(n+k+2) = A050166(n, n-k) = a(n-1, k-1) + 2*a(n-1, k)+ a (n-1, k+1) [with a(0, 0) = 1 and a(n, k) = 0 if n<0 or nHenry Bottomley, Sep 24 2001
From Philippe Deléham, Feb 14 2004: (Start)
T(n, 0) = A000108(n+1), T(n, k) = 0 if n0, T(n, k) = Sum_{j=1..n} T(n-j, k-1)*A000108(j).
G.f. for column k: Sum_{n>=0} T(n, k)*x^n = x^k*C(x)^(2*k+2) where C(x) = Sum_{n>=0} A000108(n)*x^n is g.f. for Catalan numbers, A000108.
Sum_{k>=0} T(m, k)*T(n, k) = A000108(m+n+1). (End)
T(n, k) = A009766(n+k+1, n-k) = A033184(n+k+2, 2k+2). - Philippe Deléham, Feb 14 2004
Sum_{j>=0} T(k, j)*A039599(n-k, j) = A028364(n, k). - Philippe Deléham, Mar 04 2004
Antidiagonal Sum_{k=0..n} T(n-k, k) = A000957(n+3). - Gerald McGarvey, Jun 05 2005
The triangle may also be generated from M^n * [1,0,0,0,...], where M = an infinite tridiagonal matrix with 1's in the super- and subdiagonals and [2,2,2,...] in the main diagonal. - Gary W. Adamson, Dec 17 2006
G.f.: G(t,x) = C^2/(1-txC^2), where C = (1-sqrt(1-4x))/(2x) is the Catalan function. From here G(-1,x)=C, i.e., the alternating row sums are the Catalan numbers (A000108). - Emeric Deutsch, Jan 20 2007
Sum_{k=0..n} T(n,k)*x^k = A000957(n+1), A000108(n), A000108(n+1), A001700(n), A049027(n+1), A076025(n+1), A076026(n+1) for x=-2,-1,0,1,2,3,4 respectively (see square array in A067345). - Philippe Deléham, Mar 21 2007, Nov 04 2011
Sum_{k=0..n} T(n,k)*(k+1) = 4^n. - Philippe Deléham, Mar 30 2007
Sum_{j>=0} T(n,j)*binomial(j,k) = A035324(n,k), A035324 with offset 0 (0 <= k <= n). - Philippe Deléham, Mar 30 2007
T(n,k) = A053121(2*n+1,2*k+1). - Philippe Deléham, Apr 16 2007, Apr 18 2007
T(n,k) = A039599(n,k) + A039599(n,k+1). - Philippe Deléham, Sep 11 2007
Sum_{k=0..n+1} T(n+1,k)*k^2 = A029760(n). - Philippe Deléham, Dec 16 2007
Sum_{k=0..n} T(n,k)*A059841(k) = A000984(n). - Philippe Deléham, Nov 12 2008
G.f.: 1/(1-xy-2x-x^2/(1-2x-x^2/(1-2x-x^2/(1-2x-x^2/(1-2x-x^2/(1-.... (continued fraction).
Sum_{k=0..n} T(n,k)*x^(n-k) = A000012(n), A001700(n), A194723(n+1), A194724(n+1), A194725(n+1), A194726(n+1), A194727(n+1), A194728(n+1), A194729(n+1), A194730(n+1) for x = 0,1,2,3,4,5,6,7,8,9 respectively. - Philippe Deléham, Nov 03 2011
From Peter Bala, Dec 21 2014: (Start)
This triangle factorizes in the Riordan group as ( C(x), x*C(x) ) * ( 1/(1 - x), x/(1 - x) ) = A033184 * A007318, where C(x) = (1 - sqrt(1 - 4*x))/(2*x) is the o.g.f. for the Catalan numbers A000108.
Let U denote the lower unit triangular array with 1's on or below the main diagonal and zeros elsewhere. For k = 0,1,2,... define U(k) to be the lower unit triangular block array
/I_k 0\
\ 0 U/ having the k X k identity matrix I_k as the upper left block; in particular, U(0) = U. Then this array equals the bi-infinite product (...*U(2)*U(1)*U(0))*(U(0)*U(1)*U(2)*...). (End)
From Peter Bala, Jul 21 2015: (Start)
O.g.f. G(x,t) = (1/x) * series reversion of ( x/f(x,t) ), where f(x,t) = ( 1 + (1 + t)*x )^2/( 1 + t*x ).
1 + x*d/dx(G(x,t))/G(x,t) = 1 + (2 + t)*x + (6 + 4*t + t^2)*x^2 + ... is the o.g.f for A094527. (End)
Conjecture: Sum_{k=0..n} T(n,k)/(k+1)^2 = H(n+1)*A000108(n)*(2*n+1)/(n+1), where H(n+1) = Sum_{k=0..n} 1/(k+1). - Werner Schulte, Jul 23 2015
From Werner Schulte, Jul 25 2015: (Start)
Sum_{k=0..n} T(n,k)*(k+1)^2 = (2*n+1)*binomial(2*n,n). (A002457)
Sum_{k=0..n} T(n,k)*(k+1)^3 = 4^n*(3*n+2)/2.
Sum_{k=0..n} T(n,k)*(k+1)^4 = (2*n+1)^2*binomial(2*n,n).
Sum_{k=0..n} T(n,k)*(k+1)^5 = 4^n*(15*n^2+15*n+4)/4. (End)
The o.g.f. G(x,t) is such that G(x,t+1) is the o.g.f. for A035324, but with an offset of 0, and G(x,t-1) is the o.g.f. for A033184, again with an offset of 0. - Peter Bala, Sep 20 2015
Denote this lower triangular array by L; then L * transpose(L) is the Cholesky factorization of the Hankel matrix ( 1/(i+j)*binomial(2*i+2*j-2, i+j-1) )A172417%20read%20as%20a%20square%20array.%20See%20Chamberland,%20p.%201669.%20-%20_Peter%20Bala">i,j >= 1 = A172417 read as a square array. See Chamberland, p. 1669. - _Peter Bala, Oct 15 2023

Extensions

Typo in one entry corrected by Philippe Deléham, Dec 16 2007

A003506 Triangle of denominators in Leibniz's Harmonic Triangle a(n,k), n >= 1, 1 <= k <= n.

Original entry on oeis.org

1, 2, 2, 3, 6, 3, 4, 12, 12, 4, 5, 20, 30, 20, 5, 6, 30, 60, 60, 30, 6, 7, 42, 105, 140, 105, 42, 7, 8, 56, 168, 280, 280, 168, 56, 8, 9, 72, 252, 504, 630, 504, 252, 72, 9, 10, 90, 360, 840, 1260, 1260, 840, 360, 90, 10, 11, 110, 495, 1320, 2310, 2772, 2310, 1320, 495, 110, 11
Offset: 1

Views

Author

Keywords

Comments

Array 1/Beta(n,m) read by antidiagonals. - Michael Somos, Feb 05 2004
a(n,3) = A027480(n-2); a(n,4) = A033488(n-3). - Ross La Haye, Feb 13 2004
a(n,k) = total size of all of the elements of the family of k-size subsets of an n-element set. For example, a 2-element set, say, {1,2}, has 3 families of k-size subsets: one with 1 0-size element, one with 2 1-size elements and one with 1 2-size element; respectively, {{}}, {{1},{2}}, {{1,2}}. - Ross La Haye, Dec 31 2006
Second slice along the 1-2-plane in the cube a(m,n,o) = a(m-1,n,o) + a(m,n-1,o) + a(m,n,o-1) with a(1,0,0)=1 and a(m<>1=0,n>=0,0>=o)=0, for which the first slice is Pascal's triangle (slice read by antidiagonals). - Thomas Wieder, Aug 06 2006
Triangle, read by rows, given by [2,-1/2,1/2,0,0,0,0,0,0,...] DELTA [2,-1/2,1/2,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Oct 07 2007
This sequence * [1/1, 1/2, 1/3, ...] = (1, 3, 7, 15, 31, ...). - Gary W. Adamson, Nov 14 2007
n-th row = coefficients of first derivative of corresponding Pascal's triangle row. Example: x^4 + 4x^3 + 6x^2 + 4x + 1 becomes (4, 12, 12, 4). - Gary W. Adamson, Dec 27 2007
From Paul Curtz, Jun 03 2011: (Start)
Consider
1 1/2 1/3 1/4 1/5
-1/2 -1/6 -1/12 -1/20 -1/30
1/3 1/12 1/30 1/60 1/105
-1/4 -1/20 -1/60 -1/140 -1/280
1/5 1/30 1/105 1/280 1/630
This is an autosequence (the inverse binomial transform is the sequence signed) of the second kind: the main diagonal is 2 times the first upper diagonal.
Note that 2, 12, 60, ... = A005430(n+1), Apery numbers = 2*A002457(n). (End)
From Louis Conover (for the 9th grade G1c mathematics class at the Chengdu Confucius International School), Mar 02 2015: (Start)
The i-th order differences of n^-1 appear in the (i+1)th row.
1, 1/2, 1/3, 1/4, 1/5, 1/6, 1/7, 1/8, ...
1/2, 1/6, 1/12, 1/20, 1/30, 1/42, 1/56, 1/72, ...
1/3, 1/12, 1/30, 1/60, 1/105, 1/168, 1/252, 1/360, ...
1/4, 1/20, 1/60, 1/140, 1/280, 1/504, 1/840, 1/1320, ...
1/5, 1/30, 1/105, 1/280, 1/630, 1/1260, 1/2310, 1/3960, ...
1/6, 1/42, 1/168, 1/504, 1/1260, 1/2772, 1/5544, 1/12012, ...
(End)
T(n,k) is the number of edges of distance k from a fixed vertex in the n-dimensional hypercube. - Simon Burton, Nov 04 2022

Examples

			The triangle begins:
  1;
  1/2, 1/2;
  1/3, 1/6, 1/3;
  1/4, 1/12, 1/12, 1/4;
  1/5, 1/20, 1/30, 1/20, 1/5;
  ...
The triangle of denominators begins:
   1
   2   2
   3   6   3
   4  12  12    4
   5  20  30   20    5
   6  30  60   60   30    6
   7  42 105  140  105   42    7
   8  56 168  280  280  168   56    8
   9  72 252  504  630  504  252   72   9
  10  90 360  840 1260 1260  840  360  90  10
  11 110 495 1320 2310 2772 2310 1320 495 110 11
		

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, see 130.
  • B. A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8. English translation published by Fibonacci Association, Santa Clara Univ., Santa Clara, CA, 1993; see p. 38.
  • G. Boole, A Treatise On The Calculus of Finite Differences, Dover, 1960, p. 26.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 83, Problem 25.
  • M. Elkadi and B. Mourrain, Symbolic-numeric methods for solving polynomial equations and applications, Chap 3. of A. Dickenstein and I. Z. Emiris, eds., Solving Polynomial Equations, Springer, 2005, pp. 126-168. See p. 152.
  • D. Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, 35.

Crossrefs

Row sums are in A001787. Central column is A002457. Half-diagonal is in A090816. A116071, A215652.
Denominators of i-th order differences of n^-1 are given in: (1st) A002378, (2nd) A027480, (3rd) A033488, (4th) A174002, (5th) A253946. - Louis Conover, Mar 02 2015
Columns k >= 1 (offset 1): A000027, A002378, A027480, A033488, A174002, A253946(n+4), ..., with sum of reciprocals: infinity, 1, 1/2, 1/3, 1/4, 1/5, ..., respectively. - Wolfdieter Lang, Jul 20 2022

Programs

  • Haskell
    a003506 n k = a003506_tabl !! (n-1) !! (n-1)
    a003506_row n = a003506_tabl !! (n-1)
    a003506_tabl = scanl1 (\xs ys ->
       zipWith (+) (zipWith (+) ([0] ++ xs) (xs ++ [0])) ys) a007318_tabl
    a003506_list = concat a003506_tabl
    -- Reinhard Zumkeller, Nov 14 2013, Nov 17 2011
    
  • Maple
    with(combstruct):for n from 0 to 11 do seq(m*count(Combination(n), size=m), m = 1 .. n) od; # Zerinvary Lajos, Apr 09 2008
    A003506 := (n,k) -> k*binomial(n,k):
    seq(print(seq(A003506(n,k),k=1..n)),n=1..7); # Peter Luschny, May 27 2011
  • Mathematica
    L[n_, 1] := 1/n; L[n_, m_] := L[n, m] = L[n - 1, m - 1] - L[n, m - 1]; Take[ Flatten[ Table[ 1 / L[n, m], {n, 1, 12}, {m, 1, n}]], 66]
    t[n_, m_] = Gamma[n]/(Gamma[n - m]*Gamma[m]); Table[Table[t[n, m], {m, 1, n - 1}], {n, 2, 12}]; Flatten[%] (* Roger L. Bagula and Gary W. Adamson, Sep 14 2008 *)
    Table[k*Binomial[n,k],{n,1,7},{k,1,n}] (* Peter Luschny, May 27 2011 *)
    t[n_, k_] := Denominator[n!*k!/(n+k+1)!]; Table[t[n-k, k] , {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Nov 28 2013 *)
  • PARI
    A(i,j)=if(i<1||j<1,0,1/subst(intformal(x^(i-1)*(1-x)^(j-1)),x,1))
    
  • PARI
    A(i,j)=if(i<1||j<1,0,1/sum(k=0,i-1,(-1)^k*binomial(i-1,k)/(j+k)))
    
  • PARI
    {T(n, k) = (n + 1 - k) * binomial( n, k - 1)} /* Michael Somos, Feb 06 2011 */
    
  • SageMath
    T_row = lambda n: (n*(x+1)^(n-1)).list()
    for n in (1..10): print(T_row(n)) # Peter Luschny, Feb 04 2017
    # Assuming offset 0:
    def A003506(n, k):
        return falling_factorial(n+1,n)//(factorial(k)*factorial(n-k))
    for n in range(9): print([A003506(n, k) for k in range(n+1)]) # Peter Luschny, Aug 13 2022

Formula

a(n, 1) = 1/n; a(n, k) = a(n-1, k-1) - a(n, k-1) for k > 1.
Considering the integer values (rather than unit fractions): a(n, k) = k*C(n, k) = n*C(n-1, k-1) = a(n, k-1)*a(n-1, k-1)/(a(n, k-1) - a(n-1, k-1)) = a(n-1, k) + a(n-1, k-1)*k/(k-1) = (a(n-1, k) + a(n-1, k-1))*n/(n-1) = k*A007318(n, k) = n*A007318(n-1, k-1). Row sums of integers are n*2^(n-1) = A001787(n); row sums of the unit fractions are A003149(n-1)/A000142(n). - Henry Bottomley, Jul 22 2002
From Vladeta Jovovic, Nov 01 2003: (Start)
G.f.: x*y/(1-x-y*x)^2.
E.g.f.: x*y*exp(x+x*y). (End)
T(n,k) = n*binomial(n-1,k-1) = n*A007318(n-1,k-1). - Philippe Deléham, Aug 04 2006
Binomial transform of A128064(unsigned). - Gary W. Adamson, Aug 29 2007
From Roger L. Bagula and Gary W. Adamson, Sep 14 2008: (Start)
t(n,m) = Gamma(n)/(Gamma(n - m)*Gamma(m)).
f(s,n) = Integral_{x=0..oo} exp(-s*x)*x^n dx = Gamma(n)/s^n; t(n,m) = f(s,n)/(f(s,n-m)*f(s,m)) = Gamma(n)/(Gamma(n - m)*Gamma(m)); the powers of s cancel out. (End)
From Reinhard Zumkeller, Mar 05 2010: (Start)
T(n,5) = T(n,n-4) = A174002(n-4) for n > 4.
T(2*n,n) = T(2*n,n+1) = A005430(n). (End)
T(n,k) = 2*T(n-1,k) + 2*T(n-1,k-1) - T(n-2,k) - 2*T(n-2,k-1) - T(n-2,k-2), T(1,1) = 1 and, for n > 1, T(n,k) = 0 if k <= 1 or if k > n. - Philippe Deléham, Mar 17 2012
T(n,k) = Sum_{i=1..k} i*binomial(k,i)*binomial(n+1-k,k+1-i). - Mircea Merca, Apr 11 2012
If we include a main diagonal of zeros so that the array is in the form
0
1 0
2 2 0
3 6 3 0
4 12 12 4 0
...
then we obtain the exponential Riordan array [x*exp(x),x], which factors as [x,x]*[exp(x),x] = A132440*A007318. This array is the infinitesimal generator for A116071. A signed version of the array is the infinitesimal generator for A215652. - Peter Bala, Sep 14 2012
a(n,k) = (n-1)!/((n-k)!(k-1)!) if k > n/2 and a(n,k) = (n-1)!/((n-k-1)!k!) otherwise. [Forms 'core' for Pascal's recurrence; gives common term of RHS of T(n,k) = T(n-1,k-1) + T(n-1,k)]. - Jon Perry, Oct 08 2013
Assuming offset 0: T(n, k) = FallingFactorial(n + 1, n) / (k! * (n - k)!). The counterpart using the rising factorial is A356546. - Peter Luschny, Aug 13 2022

Extensions

Edited by N. J. A. Sloane, Oct 07 2007

A063007 T(n,k) = binomial(n,k)*binomial(n+k,k), 0 <= k <= n, triangle read by rows.

Original entry on oeis.org

1, 1, 2, 1, 6, 6, 1, 12, 30, 20, 1, 20, 90, 140, 70, 1, 30, 210, 560, 630, 252, 1, 42, 420, 1680, 3150, 2772, 924, 1, 56, 756, 4200, 11550, 16632, 12012, 3432, 1, 72, 1260, 9240, 34650, 72072, 84084, 51480, 12870, 1, 90, 1980, 18480, 90090, 252252, 420420, 411840, 218790, 48620
Offset: 0

Views

Author

Henry Bottomley, Jul 02 2001

Keywords

Comments

T(n,k) is the number of compatible k-sets of cluster variables in Fomin and Zelevinsky's Cluster algebra of finite type B_n. Take a row of this triangle regarded as a polynomial in x and rewrite as a polynomial in y := x+1. The coefficients of the polynomial in y give a row of triangle A008459 (squares of binomial coefficients). For example, x^2+6*x+6 = y^2+4*y+1. - Paul Boddington, Mar 07 2003
T(n,k) is the number of lattice paths from (0,0) to (n,n) using steps E=(1,0), N=(0,1) and D=(1,1) (i.e., bilateral Schroeder paths), having k N=(0,1) steps. E.g. T(2,0)=1 because we have DD; T(2,1) = 6 because we have NED, NDE, EDN, END, DEN and DNE; T(2,2)=6 because we have NNEE, NENE, NEEN, EENN, ENEN and ENNE. - Emeric Deutsch, Apr 20 2004
Another version of [1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, ...] DELTA [0, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, ...] = 1; 1, 0; 1, 2, 0; 1, 6, 6, 0; 1, 12, 30, 20, 0; ..., where DELTA is the operator defined in A084938. - Philippe Deléham Apr 15 2005
Terms in row n are the coefficients of the Legendre polynomial P(n,2x+1) with increasing powers of x.
From Peter Bala, Oct 28 2008: (Start)
Row n of this triangle is the f-vector of the simplicial complex dual to an associahedron of type B_n (a cyclohedron) [Fomin & Reading, p.60]. See A008459 for the corresponding h-vectors for associahedra of type B_n and A001263 and A033282 respectively for the h-vectors and f-vectors for associahedra of type A_n.
An alternative description of this triangle in terms of f-vectors is as follows. Let A_n be the root lattice generated as a monoid by {e_i - e_j: 0 <= i,j <= n+1}. Let P(A_n) be the polytope formed by the convex hull of this generating set. Then the rows of this array are the f-vectors of a unimodular triangulation of P(A_n) [Ardila et al.]. A008459 is the corresponding array of h-vectors for these type A_n polytopes. See A127674 (without the signs) for the array of f-vectors for type C_n polytopes and A108556 for the array of f-vectors associated with type D_n polytopes.
The S-transform on the ring of polynomials is the linear transformation of polynomials that is defined on the basis monomials x^k by S(x^k) = binomial(x,k) = x(x-1)...(x-k+1)/k!. Let P_n(x) denote the S-transform of the n-th row polynomial of this array. In the notation of [Hetyei] these are the Stirling polynomials of the type B associahedra. The first few values are P_1(x) = 2*x + 1, P_2(x) = 3*x^2 + 3*x + 1 and P_3(x) = (10*x^3 + 15*x^2 + 11*x + 3)/3. These polynomials have their zeros on the vertical line Re x = -1/2 in the complex plane, that is, the polynomials P_n(-x) satisfy a Riemann hypothesis. See A142995 for further details. The sequence of values P_n(k) for k = 0,1,2,3, ... produces the n-th row of A108625. (End)
This is the row reversed version of triangle A104684. - Wolfdieter Lang, Sep 12 2016
T(n, k) is also the number of (n-k)-dimensional faces of a convex n-dimensional Lipschitz polytope of real functions f defined on the set X = {1, 2, ..., n+1} which satisfy the condition f(n+1) = 0 (see Gordon and Petrov). - Stefano Spezia, Sep 25 2021
The rows seem to give (up to sign) the coefficients in the expansion of the integer-valued polynomial ((x+1)*(x+2)*(x+3)*...*(x+n) / n!)^2 in the basis made of the binomial(x+i,i). - F. Chapoton, Oct 09 2022
Chapoton's observation above is correct: the precise expansion is ((x+1)*(x+2)*(x+3)*...*(x+n)/ n!)^2 = Sum_{k = 0..n} (-1)^k*T(n,n-k)*binomial(x+2*n-k, 2*n-k), as can be verified using the WZ algorithm. For example, n = 3 gives ((x+1)*(x+2)*(x+3)/3!)^2 = 20*binomial(x+6,6) - 30*binomial(x+5,5) + 12*binomial(x+4,4) - binomial(x+3,3). - Peter Bala, Jun 24 2023

Examples

			The triangle T(n, k) starts:
  n\k 0  1    2     3     4      5      6      7      8     9
  0:  1
  1:  1  2
  2:  1  6    6
  3:  1 12   30    20
  4:  1 20   90   140    70
  5:  1 30  210   560   630    252
  6:  1 42  420  1680  3150   2772    924
  7:  1 56  756  4200 11550  16632  12012   3432
  8:  1 72 1260  9240 34650  72072  84084  51480  12870
  9:  1 90 1980 18480 90090 252252 420420 411840 218790 48620
... reformatted by _Wolfdieter Lang_, Sep 12 2016
From _Petros Hadjicostas_, Jul 11 2020: (Start)
Its inverse (from Table II, p. 92, in Ser's book) is
   1;
  -1/2,  1/2;
   1/3, -1/2,    1/6;
  -1/4,  9/20,  -1/4,   1/20;
   1/5, -2/5,    2/7,  -1/10,  1/70;
  -1/6,  5/14, -25/84,  5/36, -1/28,  1/252;
   1/7, -9/28,  25/84, -1/6,   9/154, -1/84, 1/924;
   ... (End)
		

References

  • J. M. Borwein and P. B. Borwein, Pi and the AGM, Wiley, 1987, p. 366.
  • J. Ser, Les Calculs Formels des Séries de Factorielles. Gauthier-Villars, Paris, 1933, Table I, p. 92.
  • D. Zagier, Integral solutions of Apery-like recurrence equations, in: Groups and Symmetries: from Neolithic Scots to John McKay, CRM Proc. Lecture Notes 47, Amer. Math. Soc., Providence, RI, 2009, pp. 349-366.

Crossrefs

See A331430 for an essentially identical triangle, except with signed entries.
Columns include A000012, A002378, A033487 on the left and A000984, A002457, A002544 on the right.
Main diagonal is A006480.
Row sums are A001850. Alternating row sums are A033999.
Cf. A033282 (f-vectors type A associahedra), A108625, A080721 (f-vectors type D associahedra).
The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692,A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)

Programs

  • Haskell
    a063007 n k = a063007_tabl !! n !! k
    a063007_row n = a063007_tabl !! n
    a063007_tabl = zipWith (zipWith (*)) a007318_tabl a046899_tabl
    -- Reinhard Zumkeller, Nov 18 2014
    
  • Magma
    /* As triangle: */ [[Binomial(n,k)*Binomial(n+k,k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Sep 03 2015
  • Maple
    p := (n,x) -> orthopoly[P](n,1+2*x): seq(seq(coeff(p(n,x),x,k), k=0..n), n=0..9);
  • Mathematica
    Flatten[Table[Binomial[n, k]Binomial[n + k, k], {n, 0, 10}, {k, 0, n}]] (* Harvey P. Dale, Dec 24 2011 *)
    Table[CoefficientList[Hypergeometric2F1[-n, n + 1, 1, -x], x], {n, 0, 9}] // Flatten
    (* Peter Luschny, Mar 09 2018 *)
  • PARI
    {T(n, k) = local(t); if( n<0, 0, t = (x + x^2)^n; for( k=1, n, t=t'); polcoeff(t, k) / n!)} /* Michael Somos, Dec 19 2002 */
    
  • PARI
    {T(n, k) = binomial(n, k) * binomial(n+k, k)} /* Michael Somos, Sep 22 2013 */
    
  • PARI
    {T(n, k) = if( k<0 || k>n, 0, (n+k)! / (k!^2 * (n-k)!))} /* Michael Somos, Sep 22 2013 */
    

Formula

T(n, k) = (n+k)!/(k!^2*(n-k)!) = T(n-1, k)*(n+k)/(n-k) = T(n, k-1)*(n+k)*(n-k+1)/k^2 = T(n-1, k-1)*(n+k)*(n+k-1)/k^2.
binomial(x, n)^2 = Sum_{k>=0} T(n,k) * binomial(x, n+k). - Michael Somos, May 11 2012
T(n, k) = A109983(n, k+n). - Michael Somos, Sep 22 2013
G.f.: G(t, z) = 1/sqrt(1-2*z-4*t*z+z^2). Row generating polynomials = P_n(1+2z), i.e., T(n, k) = [z^k] P_n(1+2*z), where P_n are the Legendre polynomials. - Emeric Deutsch, Apr 20 2004
Sum_{k>=0} T(n, k)*A000172(k) = Sum_{k>=0} T(n, k)^2 = A005259(n). - Philippe Deléham, Jun 08 2005
1 + z*d/dz(log(G(t,z))) = 1 + (1 + 2*t)*z + (1 + 8*t + 8*t^2)*z^2 + ... is the o.g.f. for a signed version of A127674. - Peter Bala, Sep 02 2015
If R(n,t) denotes the n-th row polynomial then x^3 * exp( Sum_{n >= 1} R(n,t)*x^n/n ) = x^3 + (1 + 2*t)*x^4 + (1 + 5*t + 5*t^2)*x^5 + (1 + 9*t + 21*t^2 + 14*t^3)*x^6 + ... is an o.g.f for A033282. - Peter Bala, Oct 19 2015
P(n,x) := 1/(1 + x)*Integral_{t = 0..x} R(n,t) dt are (modulo differences of offset) the row polynomials of A033282. - Peter Bala, Jun 23 2016
From Peter Bala, Mar 09 2018: (Start)
R(n,x) = Sum_{k = 0..n} binomial(2*k,k)*binomial(n+k,n-k)*x^k.
R(n,x) = Sum_{k = 0..n} binomial(n,k)^2*x^k*(1 + x)^(n-k).
n*R(n,x) = (1 + 2*x)*(2*n - 1)*R(n-1,x) - (n - 1)*R(n-2,x).
R(n,x) = (-1)^n*R(n,-1 - x).
R(n,x) = 1/n! * (d/dx)^n ((x^2 + x)^n). (End)
The row polynomials are R(n,x) = hypergeom([-n, n + 1], [1], -x). - Peter Luschny, Mar 09 2018
T(n,k) = C(n+1,k)*A009766(n,k). - Bob Selcoe, Jan 18 2020 (Connects this triangle with the Catalan triangle. - N. J. A. Sloane, Jan 18 2020)
If we let A(n,k) = (-1)^(n+k)*(2*k+1)*(n*(n-1)*...*(n-(k-1)))/((n+1)*...*(n+(k+1))) for n >= 0 and k = 0..n, and we consider both T(n,k) and A(n,k) as infinite lower triangular arrays, then they are inverses of one another. (Empty products are by definition 1.) See the example below. The rational numbers |A(n,k)| appear in Table II on p. 92 in Ser's (1933) book. - Petros Hadjicostas, Jul 11 2020
From Peter Bala, Nov 28 2021: (Start)
Row polynomial R(n,x) = Sum_{k >= n} binomial(k,n)^2 * x^(k-n)/(1+x)^(k+1) for x > -1/2.
R(n,x) = 1/(1 + x)^(n+1) * hypergeom([n+1, n+1], [1], x/(1 + x)).
R(n,x) = (1 + x)^n * hypergeom([-n, -n], [1], x/(1 + x)).
R(n,x) = hypergeom([(n+1)/2, -n/2], [1], -4*x*(1 + x)).
If we set R(-1,x) = 1, we can run the recurrence n*R(n,x) = (1 + 2*x)*(2*n - 1)*R(n-1,x) - (n - 1)*R(n-2,x) backwards to give R(-n,x) = R(n-1,x).
R(n,x) = [t^n] ( (1 + t)*(1 + x*(1 + t)) )^n. (End)
n*T(n,k) = (2*n-1)*T(n-1,k) + (4*n-2)*T(n-1,k-1) - (n-1)*T(n-2,k). - Fabián Pereyra, Jun 30 2022
From Peter Bala, Oct 07 2024: (Start)
n-th row polynomial R(n,x) = Sum_{k = 0..n} binomial(n, k) * x^k o (1 + x)^(n-k), where o denotes the black diamond product of power series as defined by Dukes and White (see Bala, Section 4.4, exercise 3).
Denote this triangle by T. Then T * transpose(T) = A143007, the square array of crystal ball sequences for the A_n X A_n lattices.
Let S denote the triangle ((-1)^(n+k)*T(n, k))n,k >= 0, a signed version of this triangle. Then S^(-1) * T = A007318, Pascal's triangle; it appears that T * S^(-1) = A110098.
T = A007318 * A115951. (End)

A001803 Numerators in expansion of (1 - x)^(-3/2).

Original entry on oeis.org

1, 3, 15, 35, 315, 693, 3003, 6435, 109395, 230945, 969969, 2028117, 16900975, 35102025, 145422675, 300540195, 9917826435, 20419054425, 83945001525, 172308161025, 1412926920405, 2893136075115, 11835556670925
Offset: 0

Views

Author

Keywords

Comments

a(n) is the denominator of the integral from 0 to Pi of (sin(x))^(2*n+1). - James R. Buddenhagen, Aug 17 2008
a(n) is the denominator of (2n)!!/(2*n + 1)!! = 2^(2*n)*n!*n!/(2*n + 1)! (see Andersson). - N. J. A. Sloane, Jun 27 2011
a(n) = (2*n + 1)*A001790(n). A046161(n)/a(n) = 1, 2/3, 8/15, 16/35, 128/315, 256/693, ... is binomial transform of Madhava-Gregory-Leibniz series for Pi/4 (i.e., 1 - 1/3 + 1/5 - 1/7 + ... ). See A173384 and A173396. - Paul Curtz, Feb 21 2010
a(n) is the denominator of Integral_{x=-oo..oo} sech(x)^(2*n+2) dx. The corresponding numerator is A101926(n). - Mohammed Yaseen, Jul 25 2023

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 798.
  • G. Prévost, Tables de Fonctions Sphériques. Gauthier-Villars, Paris, 1933, pp. 156-157.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 6, equation 6:14:9 at page 51.

Crossrefs

The denominator is given in A046161.
Largest odd divisors of A001800, A002011, A002457, A005430, A033876, A086228.
Bisection of A004731, A004735, A086116.
Second column of triangle A100258.
Cf. A002596 (numerators in expansion of (1-x)^(1/2)).
Cf. A161198 (triangle related to the series expansions of (1-x)^((-1-2*n)/2)).
A163590 is the odd part of the swinging factorial, A001790 at even indices. - Peter Luschny, Aug 01 2009

Programs

  • Julia
    A001803(n) = sum(<<(A001790(k), A005187(n) - A005187(k)) for k in 0:n) # Peter Luschny, Oct 03 2019
    
  • Magma
    A001803:= func< n | Numerator(Binomial(n+2,2)*Catalan(n+1)/4^n) >;
    [A001803(n): n in [0..30]]; // G. C. Greubel, Apr 27 2025
    
  • Maple
    swing := proc(n) option remember; if n = 0 then 1 elif irem(n, 2) = 1 then swing(n-1)*n else 4*swing(n-1)/n fi end:
    sigma := n -> 2^(add(i,i= convert(iquo(n,2),base,2))):
    a := n -> swing(2*n+1)/sigma(2*n+1); # Peter Luschny, Aug 01 2009
    A001803 := proc(n) (2*n+1)*binomial(2*n,n)/4^n ; numer(%) ; end proc: # R. J. Mathar, Jul 06 2011
    a := n -> denom(Pi*binomial(n, -1/2)): seq(a(n), n = 0..22); # Peter Luschny, Dec 06 2024
  • Mathematica
    Numerator/@CoefficientList[Series[(1-x)^(-3/2),{x,0,25}],x]  (* Harvey P. Dale, Feb 19 2011 *)
    Table[Denominator[Beta[1, n + 1, 1/2]], {n, 0, 22}] (* Gerry Martens, Nov 13 2016 *)
  • PARI
    a(n) = numerator((2*n+1)*binomial(2*n,n)/(4^n)); \\ Altug Alkan, Sep 06 2018
    
  • SageMath
    def A001803(n): return numerator((n+1)*binomial(2*n+2,n+1)/2^(2*n+1))
    print([A001803(n) for n in range(31)]) # G. C. Greubel, Apr 27 2025

Formula

a(n) = (2*n + 1)! /(n!^2*2^A000120(n)) = (n + 1)*binomial(2*n+2,n+1)/2^(A000120(n)+1). - Ralf Stephan, Mar 10 2004
From Johannes W. Meijer, Jun 08 2009: (Start)
a(n) = numerator( (2*n+1)*binomial(2*n,n)/(4^n) ).
(1 - x)^(-3/2) = Sum_{n>=0} ((2*n+1)*binomial(2*n,n)/4^n)*x^n. (End)
Truncations of rational expressions like those given by the numerator or denominator operators are artifacts in integer formulas and have many disadvantages. A pure integer formula follows. Let n$ denote the swinging factorial and sigma(n) = number of '1's in the base-2 representation of floor(n/2). Then a(n) = (2*n+1)$ / sigma(2*n+1) = A056040(2*n+1) / A060632(2*n+2). Simply said: This sequence gives the odd part of the swinging factorial at odd indices. - Peter Luschny, Aug 01 2009
a(n) = denominator(Pi*binomial(n, -1/2)). - Peter Luschny, Dec 06 2024

A002697 a(n) = n*4^(n-1).

Original entry on oeis.org

0, 1, 8, 48, 256, 1280, 6144, 28672, 131072, 589824, 2621440, 11534336, 50331648, 218103808, 939524096, 4026531840, 17179869184, 73014444032, 309237645312, 1305670057984, 5497558138880, 23089744183296
Offset: 0

Views

Author

Keywords

Comments

Coefficient of x^(2n-2) in Chebyshev polynomial T(2n) is -a(n).
Let M_n be the n X n matrix m_(i,j) = 1 + 2*abs(i-j); then det(M_n) = (-1)^(n-1)*a(n-1). - Benoit Cloitre, May 28 2002
Number of subsequences 00 in all words of length n+1 on the alphabet {0,1,2,3}. Example: a(2)=8 because we have 000,001,002,003,100,200,300 (the other 57=A125145(3) words of length 3 have no subsequences 00). a(n) = Sum_{k=0..n} k*A128235(n+1, k). - Emeric Deutsch, Feb 27 2007
Let P(A) be the power set of an n-element set A. Then a(n) = the sum of the size of the symmetric difference of x and y for every subset {x,y} of P(A). - Ross La Haye, Dec 30 2007 (See the comment from Bernard Schott below.)
Let P(A) be the power set of an n-element set A and B be the Cartesian product of P(A) with itself. Then remove (y,x) from B when (x,y) is in B and x != y and call this R35. Then a(n) = the sum of the size of the symmetric difference of x and y for every (x,y) of R35. [proposed edit of comment just above; by Ross La Haye]
The numbers in this sequence are the Wiener indices of the graphs of n-cubes (Boolean hypercubes). For example, the 3-cube is the graph of the standard cube whose Wiener index is 48. - K.V.Iyer, Feb 26 2009
From Gary W. Adamson, Sep 06 2009: (Start)
Starting (1, 8, 48, ...) = 4th binomial transform of [1, 4, 0, 0, 0, ...].
Equals the sum of terms in 2^n X 2^n semi-magic square arrays in which each row and column is composed of a binomial frequency of terms in the set (1, 3, 5, 7, ...).
The first few such arrays = [1] [1,3; 3,1]; /Q.
[1, 3, 5, 3;
3, 1, 3, 5;
5, 3, 1, 3;
3, 5, 3, 1]
(sum of terms = 48, with a binomial frequency of (1, 2, 1) as to (1, 3, 5) in each row and column)
[1, 3, 5, 3, 5, 7, 5, 3;
3, 1, 3, 5, 7, 5, 3, 5;
5, 3, 1, 3, 5, 3, 5, 7;
3, 5, 3, 1, 3, 5, 7, 5;
5, 7, 5, 3, 1, 3, 5, 3;
7, 5, 3, 5, 3, 1, 3, 5;
5, 3, 5, 7, 5, 3, 1, 3;
3, 5, 7, 5, 3, 5, 3, 1]
(sum of terms = 256, with each row and column composed of one 1, three 3's, three 5's, and one 7)
... (End)
Let P(A) be the power set of an n-element set A and B be the Cartesian product of P(A) with itself. Then a(n) = the sum of the size of the intersection of x and y for every (x,y) of B. - Ross La Haye, Jan 05 2013
Following the last comment of Ross, A002699 is the similar sequence when "intersection" is replaced by "symmetric difference" and A212698 is the similar sequence when "intersection" is replaced by "union". - Bernard Schott, Jan 04 2013
Also, following the first comment of Ross, A082134 is the similar sequence when "symmetric difference" is replaced by "intersection" and A133224 is the similar sequence when "symmetric difference" is replaced by "union". - Bernard Schott, Jan 15 2013
Let [n] denote the set {1,2,3,...,n} and denote an n-permutation of the elements of [n] by p = p(1)p(2)p(3)...p(n), where p(i) is the i-th entry in the linear order given by p. Then (p(i),p(j)) is an inversion of p if i < j but p(i) > p(j). Denote the number of inversions of p by inv(p) and call a 2n-permutation p = p(1)p(2)...p(2n) 2-ordered if p(1) < p(3) < ... < p(2n-1) and p(2) < p(4) < ... < p(2n). Then Sum(inv(p)) = n*4^(n-1), where the sum is taken over all 2-ordered 2n-permutations of p. See Bona reference below. - Ross La Haye, Jan 21 2014
Sum over all peaks of Dyck paths of semilength n of the product of the x and y coordinates. - Alois P. Heinz, May 29 2015
Sum of the number of all edges over all j-dimensional subcubes of the boolean hypercube graph of dimension n, Q_n, for all j, so a(n) = Sum_{j=1..n} binomial(n,j)*2^(n-j) * j*2^(j-1). - Constantinos Kourouzides, Mar 24 2024

Examples

			From _Bernard Schott_, Jan 04 2013: (Start)
See the comment about intersection of X and Y.
If A={b,c}, then in P(A) we have:
{b}Inter{b}={b},
{b}Inter{b,c}={b},
{c}Inter{c}={c},
{c}Inter{b,c}={c},
{b,c}Inter{b}={b},
{b,c}Inter{c}={c},
{b,c}Inter{b,c}={b,c}
and : #{b}+ #{b}+ #{c}+ #{c}+ #{b}+ #{c}+ #{b,c} = 8 = 2*4^(2-1) = a(2).
The other intersections are empty.
(End)
		

References

  • Miklos Bona, Combinatorics of Permutations, Chapman and Hall/CRC, 2004, pp. 1, 43, 64.
  • C. Lanczos, Applied Analysis. Prentice-Hall, Englewood Cliffs, NJ, 1956, p. 516.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Programs

Formula

a(n) = n*4^(n-1).
G.f.: x/(1-4x)^2. a(n+1) is the convolution of powers of 4 (A000302). - Wolfdieter Lang, May 16 2003
Third binomial transform of n. E.g.f.: x*exp(4x). - Paul Barry, Jul 22 2003
a(n) = Sum_{k=0..n} k*binomial(2*n, 2*k). - Benoit Cloitre, Jul 30 2003
For n>=0, a(n+1) = Sum_{i+j+k+l=n} binomial(2i, i)*binomial(2j, j)*binomial(2k, k)*binomial(2l, l). - Philippe Deléham, Jan 22 2004
a(n) = Sum_{k=0..n} 4^(n-k)*binomial(n-k+1, k)*binomial(1, (k+1)/2)*(1-(-1)^k)/2. - Paul Barry, Oct 15 2004
Sum_{n>0} 1/a(n) = 8*log(2) - 4*log(3). - Jaume Oliver Lafont, Sep 11 2009
a(0) = 0, a(n) = 4*a(n-1) + 4^(n-1). - Vincenzo Librandi, Dec 31 2010
a(n+1) is the convolution of A000984 with A002457. - Rui Duarte, Oct 08 2011
a(0) = 0, a(1) = 1, a(n) = 8*a(n-1) - 16*a(n-2). - Harvey P. Dale, Jan 18 2012
a(n) = A002699(n)/2 = A212698(n)/3. - Bernard Schott, Jan 04 2013
G.f.: W(0)*x/2 , where W(k) = 1 + 1/( 1 - 4*x*(k+2)/( 4*x*(k+2) + (k+1)/W(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Aug 19 2013
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*log(5/4). - Amiram Eldar, Oct 28 2020
a(n) = (1/2)*Sum_{k=0..n} k*binomial(2*n, k). Compare this with the formula of Benoit Cloitre above. - Wolfdieter Lang, Nov 12 2021
a(n) = (-1)^(n-1)*det(M(n)) for n > 0, where M(n) is the n X n symmetric Toeplitz matrix whose first row consists of 1, 3, ..., 2*n-1. - Stefano Spezia, Aug 04 2022

A002802 a(n) = (2*n+3)!/(6*n!*(n+1)!).

Original entry on oeis.org

1, 10, 70, 420, 2310, 12012, 60060, 291720, 1385670, 6466460, 29745716, 135207800, 608435100, 2714556600, 12021607800, 52895074320, 231415950150, 1007340018300, 4365140079300, 18839025605400, 81007810103220, 347176329013800, 1483389769422600
Offset: 0

Views

Author

Keywords

Comments

For n >= 1 a(n) is also the number of rooted bicolored unicellular maps of genus 1 on n+2 edges. - Ahmed Fares (ahmedfares(AT)my-deja.com), Aug 20 2001
a(n) is half the number of (n+2) X 2 Young tableaux with a three horizontal walls between the first and second column. If there is a wall between two cells, the entries may be decreasing; see [Banderier, Wallner 2021], A000984 for one horizontal wall, and A002457 for two. - Michael Wallner, Jan 31 2022
From Robert Coquereaux, Feb 12 2024: (Start)
Call B(p,g) the number of genus g partitions of a set with p elements (genus-dependent Bell number). Up to an appropriate shift the given sequence counts the genus 1 partitions of a set: we have a(n) = B(n+4,1), with a(0)= B(4,1)=1.
When shifted with an offset 4 (i.e., defining b(p)=a(p-4), which starts with 0,0,0,1,10,70, etc., and b(4)=1), the given sequence reads b(p) = (1/( 2^4 3 )) * (1/( (2 p - 1) (2 p - 3))) * (1/(p - 4)!) * (2p)!/p!. In this form it appears as a generalization of Catalan numbers (that indeed count the genus 0 partitions).
Call C[p, [alpha], g] the number of partitions of a set with p elements, of cyclic type [alpha], and of genus g (genus g Faa di Bruno coefficients of type [alpha]). Up to an appropriate shift the given sequence also counts the genus 1 partitions of p=2k into k parts of length 2, which is then called C[2k, [2^k], 1], and we have a(n) = C[2k, [2^k], 1] for k=n+2.
The two previous interpretations of this sequence, leading to a(n) = B(n+4, 1) and to a(n) = C[2(n+2), [2^(n+2)], 1] are not related in any obvious way. (End)

Examples

			G.f. = 1 + 10*x + 70*x^2 + 420*x^3 + 2310*x^4 + 12012*x^5 + 60060*x^6 + ...
		

References

  • C. Jordan, Calculus of Finite Differences. Röttig and Romwalter, Budapest, 1939; Chelsea, NY, 1965, p. 449.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A035309, A000108 (for genus 0 maps), A046521 (third column).
Column g=1 of A370235.

Programs

  • GAP
    F:=Factorial;; List([0..25], n-> F(2*n+3)/(6*F(n)*F(n+1)) ); # G. C. Greubel, Jul 20 2019
  • Magma
    F:=Factorial; [F(2*n+3)/(6*F(n)*F(n+1)): n in [0..25]]; // G. C. Greubel, Jul 20 2019
    
  • Maple
    seq(simplify(4^n*hypergeom([-n,-3/2], [1], 1)),n=0..25); # Peter Luschny, Apr 26 2016
  • Mathematica
    Table[(2*n+3)!/(6*n!*(n+1)!), {n, 0, 25}] (* Vladimir Joseph Stephan Orlovsky, Dec 13 2008 *)
  • PARI
    {a(n) = if( n<0, 0, (2*n + 3)! / (6 * n! * (n+1)!))}; /* Michael Somos, Sep 16 2013 */
    
  • PARI
    {a(n) = 2^(n+3) * polcoeff( pollegendre(n+4), n) / 3}; /* Michael Somos, Sep 16 2013 */
    
  • Sage
    f=factorial; [f(2*n+3)/(6*f(n)*f(n+1)) for n in (0..25)] # G. C. Greubel, Jul 20 2019
    

Formula

G.f.: (1 - 4*x)^(-5/2) = 1F0(5/2;;4x).
Asymptotic expression for a(n) is a(n) ~ (n+2)^(3/2) * 4^(n+2) / (sqrt(Pi) * 48).
a(n) = Sum_{a+b+c+d+e=n} f(a)*f(b)*f(c)*f(d)*f(e) with f(n) = binomial(2n, n) = A000984(n). - Philippe Deléham, Jan 22 2004
a(n-1) = (1/4)*Sum_{k=1..n} k*(k+1)*binomial(2*k, k). - Benoit Cloitre, Mar 20 2004
a(n) = A051133(n+1)/3 = A000911(n)/6. - Zerinvary Lajos, Jun 02 2007
From Rui Duarte, Oct 08 2011: (Start)
Also convolution of A000984 with A002697, also convolution of A000302 with A002457.
a(n) = ((2n+3)(2n+1)/(3*1)) * binomial(2n, n).
a(n) = binomial(2n+4, 4) * binomial(2n, n) / binomial(n+2, 2).
a(n) = binomial(n+2, 2) * binomial(2n+4, n+2) / binomial(4, 2).
a(n) = binomial(2n+4, n+2) * (n+2)*(n+1) / 12. (End)
D-finite with recurrence: n*a(n) - 2*(2*n+3)*a(n-1) = 0. - R. J. Mathar, Jan 31 2014
a(n) = 4^n*hypergeom([-n,-3/2], [1], 1). - Peter Luschny, Apr 26 2016
Boas-Buck recurrence: a(n) = (10/n)*Sum_{k=0..n-1} 4^(n-k-1)*a(k), n >= 1, a(0) = 1. Proof from a(n) = A046521(n+2, 2). See a comment there. - Wolfdieter Lang, Aug 10 2017
a(n) = (-4)^n*binomial(-5/2, n). - Peter Luschny, Oct 23 2018
Sum_{n>=0} 1/a(n) = 12 - 2*sqrt(3)*Pi. - Amiram Eldar, Oct 13 2020
E.g.f.: (1/12) exp(2 x) x^2 BesselI[2, 2 x]. - Robert Coquereaux, Feb 12 2024

A002694 Binomial coefficients C(2n, n-2).

Original entry on oeis.org

1, 6, 28, 120, 495, 2002, 8008, 31824, 125970, 497420, 1961256, 7726160, 30421755, 119759850, 471435600, 1855967520, 7307872110, 28781143380, 113380261800, 446775310800, 1761039350070, 6943526580276, 27385657281648, 108043253365600
Offset: 2

Views

Author

Keywords

Comments

Number of lattice paths from (0,0) to (n,n) with steps E=(1,0) and N=(0,1) which touch or cross the line x-y=2. Example: For n=3 there are 6 paths EEENNN, EENENN, EENNEN, EENNNE, ENEENN and NEEENN. - Herbert Kociemba, May 23 2004
Number of dissections of a convex (n+3)-gon by noncrossing diagonals into several regions, exactly n-2 of which are triangular. Example: a(3)=6 because the convex hexagon ABCDEF is dissected by any of the diagonals AC, BD, CE, DF, EA, FB into regions containing exactly 1 triangle. - Emeric Deutsch, May 31 2004
Number of UUU's (triple rises), where U=(1,1), in all Dyck paths of semilength n+1. Example: a(3)=6 because we have UD(UUU)DDD, (UUU)DDDUD, (UUU)DUDDD, (UUU)DDUDD and (U[UU)U]DDDD, the triple rises being shown between parentheses. - Emeric Deutsch, Jun 03 2004
Inverse binomial transform of A026389. - Ross La Haye, Mar 05 2005
Sum of the jump-lengths of all full binary trees with n internal nodes. In the preorder traversal of a full binary tree, any transition from a node at a deeper level to a node on a strictly higher level is called a jump; the positive difference of the levels is called the jump distance; the sum of the jump distances in a given full binary tree is called the jump-length. - Emeric Deutsch, Jan 18 2007
a(n) = number of convex polyominoes (A005436) of perimeter 2n+4 that are directed but not parallelogram polyominoes, because the directed convex polyominoes are counted by the central binomial coefficient binomial(2n,n) and the subset of parallelogram polyominoes is counted by the Catalan number C(n+1) = binomial(2n+2,n+1)/(n+2) and a(n) = binomial(2n,n) - C(n+1). - David Callan, Nov 29 2007
a(n) = number of DUU's in all Dyck paths of semilength n+1. Example: a(3)=6 because we have UU(DUU)DDD, U(DUU)UDDD, U(DUU)DUDD, UDU(DUU)DD, U(DUU)DDUD, UUD(DUU)DD, the DUU's being shown between parentheses and no other Dyck path of semilength 4 contains a DUU. - David Callan, Jul 25 2008
C(2n,n-m) is the number of Dyck-type walks such that their graphs have one marked edge passed 2m times and the other edges are passed 2 times counting "there and back" directions. - Oleksiy Khorunzhiy, Jan 09 2015
Number of paths in the half-plane x >= 0, from (0,0) to (2n,4), and consisting of steps U=(1,1) and D=(1,-1). For example, for n=3, we have the 6 paths: UUUUUD, UUUUDU, UUUDUU, UUDUUU, UDUUUU, DUUUUU, DUUUUU. - José Luis Ramírez Ramírez, Apr 19 2015

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 828.
  • C. Lanczos, Applied Analysis. Prentice-Hall, Englewood Cliffs, NJ, 1956, p. 517.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A006659.
Diagonal 5 of triangle A100257.
Cf. binomial(k*n, n-k): A000027 (k=1), this sequence (k=2), A004321 (k=3), A004334 (k=4), A004347 (k=5), A004361 (k=6), A004375 (k=7), A004389 (k=8), A281580 (k=9).
Cf. binomial(2*n+m, n): A000984 (m = 0), A001700 (m = 1), A001791 (m = 2), A002054 (m = 3), A003516 (m = 5), A002696 (m = 6), A030053 - A030056, A004310 - A004318.

Programs

  • GAP
    List([2..30], n-> Binomial(2*n,n-2)); # G. C. Greubel, Mar 21 2019
  • Haskell
    a002694 n = a007318' (2 * n) (n - 2)  -- Reinhard Zumkeller, Jun 18 2012
    
  • Magma
    [Binomial(2*n, n-2): n in [2..30]]; // Vincenzo Librandi, Apr 20 2015
    
  • Maple
    a:=n->sum(binomial(n,j-1)*binomial(n,j+1),j=1..n): seq(a(n), n=2..25); # Zerinvary Lajos, Nov 26 2006
  • Mathematica
    CoefficientList[ Series[ 16/(((Sqrt[1 - 4 x] + 1)^4)*Sqrt[1 - 4 x]), {x, 0, 23}], x] (* Robert G. Wilson v, Aug 08 2011 *)
    Table[Binomial[2n,n-2],{n,2,30}] (* Harvey P. Dale, Jun 12 2014 *)
  • PARI
    {a(n) = binomial(2*n,n-2)}; \\ G. C. Greubel, Mar 21 2019
    
  • Sage
    [binomial(2*n,n-2) for n in (2..30)] # G. C. Greubel, Mar 21 2019
    

Formula

a(n) = A067310(n, 1) as this is number of ways of arranging n chords on a circle (handshakes between 2n people across a table) with exactly 1 simple intersection. - Henry Bottomley, Oct 07 2002
E.g.f.: exp(2*x) * BesselI(2, 2*x). - Vladeta Jovovic, Aug 21 2003
G.f.: (1-sqrt(1-4*z))^4/(16*z^2*sqrt(1-4*z)). - Emeric Deutsch, Jan 28 2004
a(n) = Sum_{k=0..n} C(n, k)*C(n, k+2). - Paul Barry, Sep 20 2004
D-finite with recurrence: -(n-2)*(n+2)*a(n) + 2*n*(2*n-1)*a(n-1) = 0. - R. J. Mathar, Dec 04 2012
G.f.: z^2*C(z)^4/(1-2*z*C(z)), where C(z) is the g.f. of Catalan numbers. - José Luis Ramírez Ramírez, Apr 19 2015
a(n) = Sum_{k=1..n} binomial(2*n-k,n-k-1). - Vladimir Kruchinin, Oct 22 2016
G.f.: x^2* 2F1(5/2,3;5;4*x). - R. J. Mathar, Jan 27 2020
From Amiram Eldar, May 16 2022: (Start)
Sum_{n>=2} 1/a(n) = 23/6 - 13*Pi/(9*sqrt(3)).
Sum_{n>=2} (-1)^n/a(n) = 106*log(phi)/(5*sqrt(5)) - 37/10, where phi is the golden ratio (A001622). (End)
From Peter Bala, Oct 13 2024: (Start)
a(n) = Integral_{x = 0..4} x^n * w(x) dx, where the weight function w(x) = 1/(2*Pi) * (x^2 - 4*x + 2)/sqrt(x*(4 - x)).
G.f. x^2 * B(x) * C(x)^4, where B(x) = 1/sqrt(1 - 4*x) is the g.f. of the central binomial coefficients A000984 and C(x) = (1 - sqrt(1 - 4*x))/(2*x) is the g.f. of the Catalan numbers A000108. (End)
Previous Showing 11-20 of 146 results. Next