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 28 results. Next

A113647 Triangle of numbers related to the generalized Catalan sequence C(2;n+1)=A064062(n+1), n>=0.

Original entry on oeis.org

1, 1, 3, 1, 7, 13, 1, 15, 41, 67, 1, 31, 113, 247, 381, 1, 63, 289, 783, 1545, 2307, 1, 127, 705, 2271, 5361, 9975, 14589, 1, 255, 1665, 6207, 16929, 36879, 66057, 95235, 1, 511, 3841, 16255, 50113, 123871, 255985, 446455, 636925, 1, 1023, 8705, 41215, 141441
Offset: 0

Views

Author

Wolfdieter Lang, Jan 13 2006

Keywords

Comments

This triangle, called Y(2,1), appears in the totally asymmetric exclusion process for the (unphysical) values alpha=2, beta=1. See the Derrida et al. refs. given under A064094, where the triangle entries are called Y_{N,K} for given alpha and beta.
The main diagonal (M=1) gives the generalized Catalan sequence C(2,n):=A064062(n).
The diagonal sequences give A064062(n+1), A115137, A115150-A115153, for n+1>= M=1,..,6.

Examples

			Triangle begins:
  1;
  1,3;
  1,7,13;
  1,15,41,67;
  1,31,113,247,381;
  ...
113=a(4,3)= a(4,2) + 2*a(3,3)= 31 + 2*41.
		

Crossrefs

Row sums give A115136.

Formula

a(n, n+1)=A064062(n+1) (main diagonal with M=1); a(n, n-M+2)= a(n, n-M+1) + 2*a(n-1, n-M+2), M>=2; a(n, 1)=1; n>=0.
G.f. for diagonal sequence M=1: GY(1, x):=(2*c(2*x)-1)/(1+x) with c(x) g.f. of A000108 (Catalan); for M=2: GY(2, x)=(1-2*x)*GY(1, x)-1; for M>=3: GY(M, x)= GY(M-1, x) -2*x*GY(M-2, x) + x^(M-2).
G.f. for diagonal sequence M (solution to the above given recurrence): GY(M, x)= (x^(M-1)/(1+x))*( 2^(M+1)*x*(p(M, 2*x)-(2*x)*p(M+1, 2*x)*c(2*x))+1), with c(x) g.f. of A000108 (Catalan) and p(n, x):= -((1/sqrt(x))^(n+1))*S(n-1, 1/sqrt(x)) with Chebyshev's S(n, x) polynomials given in A049310.

A115195 Triangle of numbers, called Y(1,2), related to generalized Catalan numbers A062992(n) = C(2;n+1) = A064062(n+1).

Original entry on oeis.org

1, 2, 3, 4, 10, 13, 8, 28, 54, 67, 16, 72, 180, 314, 381, 32, 176, 536, 1164, 1926, 2307, 64, 416, 1488, 3816, 7668, 12282, 14589, 128, 960, 3936, 11568, 26904, 51468, 80646, 95235, 256, 2176, 10048, 33184, 86992, 189928, 351220, 541690, 636925, 512, 4864
Offset: 0

Views

Author

Wolfdieter Lang, Feb 23 2006

Keywords

Comments

This triangle Y(1,2) appears in the totally asymmetric exclusion process for the (unphysical) values alpha=1, beta=2. See the Derrida et al. refs. given under A064094, where the triangle entries are called Y_{N,K} for given alpha and beta.
The main diagonal (M=1) gives the generalized Catalan sequence C(2,n+1):=A064062(n+1).
The diagonal sequences give A064062(n+1), 2*A084076, 4*A115194, 8*A115202, 16*A115203, 32*A115204 for n+1>= M=1,..,6.

Examples

			Triangle begins:
   1;
   2,  3;
   4, 10,  13;
   8, 28,  54,  67;
  16, 72, 180, 314, 381;
  ...
		

Crossrefs

Row sums give A084076.

Formula

G.f. m-th diagonal, m>=1: ((1 + 2*x*c(2*x))*(2*x*c(2*x))^m)/(2*x*(1+x)) with c(x) the o.g.f. of A000108 (Catalan).

A115197 Convolution of generalized Catalan numbers A064062 (called C(n;2)).

Original entry on oeis.org

1, 2, 7, 32, 169, 974, 5947, 37820, 247885, 1662890, 11362399, 78806936, 553386097, 3926523782, 28108587139, 202764451700, 1472446595221, 10755543924578, 78973277044903, 582558618222416, 4315238786662585
Offset: 0

Views

Author

Wolfdieter Lang, Feb 23 2006

Keywords

Comments

Row sums of triangle A115193, called C(1,2).
The o.g.f. given below follows from the Riordan matrix structure of the triangle A115193. See the o.g.f. for the row polynomials of A115193.

Formula

a(n)= sum(A115193(n,m),m=0..n), n>=0.
G.f.: ((1+2*x*c(2*x))/(1+x))^2 = ((1-2*x) + 6*x*c(2*x))/(1+x)^2, with the o.g.f. c(x) of Catalan numbers A000108.
a(n)= sum(C(2;n-k)*C(2;k),k=0..n), n>=0, with C(2;n):= A064062(n).
a(n)=4*A178792(n)-3*(n+1)*A064062(n+1) [From Joseph Abate, Jun 21 2010]
n*a(n) +(-7*n+13)*a(n-1) +4*(-2*n+1)*a(n-2)=0. - R. J. Mathar, Aug 09 2017

A130649 Generalized ballot transform of generalized Catalan numbers C(2;n)=A064062(n).

Original entry on oeis.org

1, 2, 8, 42, 260, 1796, 13372, 104922, 854964, 7165500
Offset: 0

Views

Author

Philippe Deléham, Jun 20 2007

Keywords

Formula

a(n)=Sum_{k, 0<=k<=n}A039599(n,k)*A064062(k).

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

A002416 a(n) = 2^(n^2).

Original entry on oeis.org

1, 2, 16, 512, 65536, 33554432, 68719476736, 562949953421312, 18446744073709551616, 2417851639229258349412352, 1267650600228229401496703205376, 2658455991569831745807614120560689152, 22300745198530623141535718272648361505980416, 748288838313422294120286634350736906063837462003712
Offset: 0

Views

Author

Keywords

Comments

For n >= 1, a(n) is the number of n X n (0, 1) matrices.
Also number of directed graphs on n labeled nodes allowing self-loops (cf. A053763).
1/2^(n^2) is the Hankel transform of C(n, n/2)*(1 + (-1)^n)/(2*2^n), or C(2n, n)/4^n with interpolated zeros. - Paul Barry, Sep 27 2007
Hankel transform of A064062. - Philippe Deléham, Nov 19 2007
a(n) is also the order of the semigroup (monoid) of all binary relations on an n-set. - Abdullahi Umar, Sep 14 2008
With offset = 1, a(n) is the number of n X n (0, 1) matrices with an even number of 1's in every row and in every column. - Geoffrey Critzer, May 23 2013
a(n) is the number of functions from an n-set to its power set (by definition of function including the empty function only when n = 0). - Rick L. Shepherd, Dec 27 2014

Examples

			G.f. = 1 + 2*x + 16*x^2 + 512*x^3 + 65536*x^4 + 33554432*x^5 + ...
		

References

  • John M. Howie, Fundamentals of semigroup theory. Oxford: Clarendon Press, (1995). - Abdullahi Umar, Sep 14 2008

Crossrefs

Bisection of A060656.

Programs

Formula

G.f. satisfies: A(x) = 1 + 2*x*A(4x). - Paul D. Hanna, Dec 04 2009
a(n) = 2^n * Sum_{i = 0..C(n, 2)} C(C(n, 2), i)*3^i. The summation conditions on the number of symmetric pairs (a,b) with aA027465, A013610. - Geoffrey Critzer, Nov 05 2024
G.f.: 1 / (1 - 2^1*x / (1 - 2^1*(2^2-1)*x / (1 - 2^5 * x / (1 - 2^3*(2^4-1)*x / (1 - 2^9*x / (1 - 2^5*(2^6-1)*x / ...)))))). - Michael Somos, May 12 2012
a(n) = [x^n] 1/(1 - 2^n*x). - Ilya Gutkovskiy, Oct 10 2017
Sum_{n>=0} 1/a(n) = A319015. - Amiram Eldar, Oct 14 2020

A106566 Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 1, 1, 1, 1, 1, 1, 1, ... ] DELTA [1, 0, 0, 0, 0, 0, 0, 0, ... ] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 2, 1, 0, 5, 5, 3, 1, 0, 14, 14, 9, 4, 1, 0, 42, 42, 28, 14, 5, 1, 0, 132, 132, 90, 48, 20, 6, 1, 0, 429, 429, 297, 165, 75, 27, 7, 1, 0, 1430, 1430, 1001, 572, 275, 110, 35, 8, 1, 0, 4862, 4862, 3432, 2002, 1001, 429, 154, 44, 9, 1
Offset: 0

Author

Philippe Deléham, May 30 2005

Keywords

Comments

Catalan convolution triangle; g.f. for column k: (x*c(x))^k with c(x) g.f. for A000108 (Catalan numbers).
Riordan array (1, xc(x)), where c(x) the g.f. of A000108; inverse of Riordan array (1, x*(1-x)) (see A109466).
Diagonal sums give A132364. - Philippe Deléham, Nov 11 2007

Examples

			Triangle begins:
  1;
  0,   1;
  0,   1,   1;
  0,   2,   2,  1;
  0,   5,   5,  3,  1;
  0,  14,  14,  9,  4,  1;
  0,  42,  42, 28, 14,  5, 1;
  0, 132, 132, 90, 48, 20, 6, 1;
From _Paul Barry_, Sep 28 2009: (Start)
Production array is
  0, 1,
  0, 1, 1,
  0, 1, 1, 1,
  0, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1, 1, 1,
  0, 1, 1, 1, 1, 1, 1, 1, 1, 1 (End)
		

Crossrefs

The three triangles A059365, A106566 and A099039 are the same except for signs and the leading term.
See also A009766, A033184, A059365 for other versions.
The following are all versions of (essentially) the same Catalan triangle: A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

  • Magma
    A106566:= func< n,k | n eq 0 select 1 else (k/n)*Binomial(2*n-k-1, n-k) >;
    [A106566(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Sep 06 2021
    
  • Maple
    A106566 := proc(n,k)
        if n = 0 then
            1;
        elif k < 0 or k > n then
            0;
        else
            binomial(2*n-k-1,n-k)*k/n ;
        end if;
    end proc: # R. J. Mathar, Mar 01 2015
  • Mathematica
    T[n_, k_] := Binomial[2n-k-1, n-k]*k/n; T[0, 0] = 1; Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Feb 18 2017 *)
    (* The function RiordanArray is defined in A256893. *)
    RiordanArray[1&, #(1-Sqrt[1-4#])/(2#)&, 11] // Flatten (* Jean-François Alcover, Jul 16 2019 *)
  • PARI
    {T(n, k) = if( k<=0 || k>n, n==0 && k==0, binomial(2*n - k, n) * k/(2*n - k))}; /* Michael Somos, Oct 01 2022 */
  • Sage
    def A106566(n, k): return 1 if (n==0) else (k/n)*binomial(2*n-k-1, n-k)
    flatten([[A106566(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Sep 06 2021
    

Formula

T(n, k) = binomial(2n-k-1, n-k)*k/n for 0 <= k <= n with n > 0; T(0, 0) = 1; T(0, k) = 0 if k > 0.
T(0, 0) = 1; T(n, 0) = 0 if n > 0; T(0, k) = 0 if k > 0; for k > 0 and n > 0: T(n, k) = Sum_{j>=0} T(n-1, k-1+j).
Sum_{j>=0} T(n+j, 2j) = binomial(2n-1, n), n > 0.
Sum_{j>=0} T(n+j, 2j+1) = binomial(2n-2, n-1), n > 0.
Sum_{k>=0} (-1)^(n+k)*T(n, k) = A064310(n). T(n, k) = (-1)^(n+k)*A099039(n, k).
Sum_{k=0..n} T(n, k)*x^k = A000007(n), A000108(n), A000984(n), A007854(n), A076035(n), A076036(n), A127628(n), A126694(n), A115970(n) for x = 0,1,2,3,4,5,6,7,8 respectively.
Sum_{k>=0} T(n, k)*x^(n-k) = C(x, n); C(x, n) are the generalized Catalan numbers.
Sum_{j=0..n-k} T(n+k,2*k+j) = A039599(n,k).
Sum_{j>=0} T(n,j)*binomial(j,k) = A039599(n,k).
Sum_{k=0..n} T(n,k)*A000108(k) = A127632(n).
Sum_{k=0..n} T(n,k)*(x+1)^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. - Philippe Deléham, Aug 25 2007
Sum_{k=0..n} T(n,k)*A000108(k-1) = A121988(n), with A000108(-1)=0. - Philippe Deléham, Aug 27 2007
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 = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Oct 27 2007
T(n,k)*2^(n-k) = A110510(n,k); T(n,k)*3^(n-k) = A110518(n,k). - Philippe Deléham, Nov 11 2007
Sum_{k=0..n} T(n,k)*A000045(k) = A109262(n), A000045: Fibonacci numbers. - Philippe Deléham, Oct 28 2008
Sum_{k=0..n} T(n,k)*A000129(k) = A143464(n), A000129: Pell numbers. - Philippe Deléham, Oct 28 2008
Sum_{k=0..n} T(n,k)*A100335(k) = A002450(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A100334(k) = A001906(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A099322(k) = A015565(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A106233(k) = A003462(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A151821(k+1) = A100320(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A082505(k+1) = A144706(n). - Philippe Deléham, Oct 30 2008
Sum_{k=0..n} T(n,k)*A000045(2k+2) = A026671(n). - Philippe Deléham, Feb 11 2009
Sum_{k=0..n} T(n,k)*A122367(k) = A026726(n). - Philippe Deléham, Feb 11 2009
Sum_{k=0..n} T(n,k)*A008619(k) = A000958(n+1). - Philippe Deléham, Nov 15 2009
Sum_{k=0..n} T(n,k)*A027941(k+1) = A026674(n+1). - Philippe Deléham, Feb 01 2014
G.f.: Sum_{n>=0, k>=0} T(n, k)*x^k*z^n = 1/(1 - x*z*c(z)) where c(z) the g.f. of A000108. - Michael Somos, Oct 01 2022

Extensions

Formula corrected by Philippe Deléham, Oct 31 2008
Corrected by Philippe Deléham, Sep 17 2009
Corrected by Alois P. Heinz, Aug 02 2012

A064063 Generalized Catalan numbers C(3; n).

Original entry on oeis.org

1, 1, 4, 25, 190, 1606, 14506, 137089, 1338790, 13403950, 136846144, 1419257434, 14911016596, 158363649640, 1697452010230, 18338919413425, 199496184219910, 2183299541440150, 24021874198331080, 265559590979820910, 2948253066186839140, 32857382497018933060
Offset: 0

Author

Wolfdieter Lang, Sep 13 2001

Keywords

Comments

a(n+1) = Y_{n}(n+1) = Z_{n}, n >= 0, in the Derrida et al. 1992 reference (see A064094) for alpha=3, beta =1 (or alpha=1, beta=3).
Hankel transform is A060722. - Paul Barry, Jan 30 2009

References

  • S. Ramanujan, Modular Equations and Approximations to pi, pp. 23-39 of Collected Papers of Srinivasa Ramanujan, Ed. G. H. Hardy et al., AMS Chelsea 2000. See page 39, equation (50).

Crossrefs

Cf. A064062 (C(2; n)).

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 30); Coefficients(R!( 6/(5+Sqrt(1-12*x)) )); // G. C. Greubel, May 02 2019
  • Mathematica
    CoefficientList[Series[6/(5+Sqrt[1-12 x]),{x,0,50}],x]  (* Harvey P. Dale, Mar 11 2011 *)
  • PARI
    a(n)=if(n<0,0,polcoeff(serreverse((x-2*x^2)/(1+x)^2+O(x^(n+1))),n)) \\ Ralf Stephan, Jun 12 2004
    
  • PARI
    {a(n)= if(n<1, n==0, polcoeff( serreverse( (x-2*x^2)/ (1+x)^2 +x*O(x^n)), n))} /* Michael Somos, Apr 11 2007 */
    
  • Sage
    def a(n):
        if n==0: return 1
        return hypergeometric([1-n, n], [-n], 3).simplify()
    [a(n) for n in range(24)] # Peter Luschny, Nov 30 2014
    

Formula

G.f.: (1+3*x*c(3*x)/2)/(1+x/2) = 1/(1-x*c(3*x)) with c(x) g.f. of Catalan numbers A000108.
a(n) = Sum_{m=0..n-1} (n-m)*binomial(n-1+m, m)*(3^m)/n.
a(n) = (-1/2)^n * (1 - 3*Sum_{k=0..n-1} C(k)*(-6)^k), n >= 1, a(0) = 1, with C(n) = A000108(n) (Catalan).
a(n) = Sum_{k=0..n} A059365(n, k)*3^(n-k). - Philippe Deléham, Jan 19 2004
Given the semi-axes a,b of an ellipse, then Ramanujan gave the highly accurate formula for the perimeter p = Pi((a+b) + (3(a-b)^2)/(10(a+b) + sqrt(a^2 + 14ab + b^2))). If we let h = ((a-b)/(a+b))^2, then (p/(Pi(a+b))-1)/4 = (3/4)* h/(10 + sqrt(4 - 3*h)) = 1*(h/16) + 1*(h/16)^2 + 4*(h/16)^3 + 25*(h/16)^4 + ... . - Michael Somos, Apr 11 2007
G.f.: 1/(1-x/(1-3*x/(1-3*x/(1-3*x/(1-.... = 1/(1-x-3*x^2/(1-6*x-9*x^2/(1-6*x-9*x^2/(1-.... (continued fractions). - Paul Barry, Jan 30 2009
G.f.: 6/(5+sqrt(1-12*x)). - Harvey P. Dale, Mar 11 2011
From Gary W. Adamson, Jul 12 2011: (Start)
a(n) = upper left term in M^n, M = the infinite square production matrix:
1, 1, 0, 0, 0, 0, ...
3, 3, 3, 0, 0, 0, ...
3, 3, 3, 3, 0, 0, ...
3, 3, 3, 3, 3, 0, ...
3, 3, 3, 3, 3, 3, ...
... (End)
D-finite with recurrence: 2*n*a(n) + (-23*n+36)*a(n-1) + 6*(-2*n+3)*a(n-2) = 0. - R. J. Mathar, Dec 03 2012 (Formula verified and used for computations. - Fung Lam, Mar 05 2014)
a(n) ~ 3^(n+1) * 4^n / (25*n^(3/2)*sqrt(Pi)). - Vaclav Kotesovec, Mar 05 2014
a(n) = hypergeometric([1-n, n], [-n], 3) for n>0. - Peter Luschny, Nov 30 2014

A064094 Triangle composed of generalized Catalan numbers.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 5, 3, 1, 1, 1, 14, 13, 4, 1, 1, 1, 42, 67, 25, 5, 1, 1, 1, 132, 381, 190, 41, 6, 1, 1, 1, 429, 2307, 1606, 413, 61, 7, 1, 1, 1, 1430, 14589, 14506, 4641, 766, 85, 8, 1, 1, 1, 4862, 95235, 137089, 55797, 10746, 1279, 113, 9, 1, 1
Offset: 0

Author

Wolfdieter Lang, Sep 13 2001

Keywords

Comments

The column m sequence (without leading zeros and the first 1) appears in the Derrida et al. 1992 reference as Z_{N}=Y_{N}(N+1), N >=0, for alpha = m, beta = 1 (or alpha = 1, beta = m). In the Derrida et al. 1993 reference the formula in eq. (39) gives Z_{N}(alpha,beta)/(alpha*beta)^N for N>=1.

Examples

			Triangle begins:
  1;
  1,    1;
  1,    1,     1;
  1,    2,     1,     1;
  1,    5,     3,     1,    1;
  1,   14,    13,     4,    1,   1;
  1,   42,    67,    25,    5,   1,   1;
  1,  132,   381,   190,   41,   6,   1,   1;
  1,  429,  2307,  1606,  413,  61,   7,   1,   1;
  1, 1430, 14589, 14506, 4641, 766,  85,   8,   1,   1;
		

Crossrefs

Columns (without leading zeros): A000012 (k=0), A000108 (k=1), A064062 (k=2), A064063 (k=3), A064087 (k=4), A064088 (k=5), A064089 (k=6), A064090 (k=7), A064091 (k=8), A064092 (k=9), A064093 (k=10).
Cf. A064095 (row sums).

Programs

  • Magma
    function A064094(n,k)
      if k eq 0 or k eq n then return 1;
      else return (&+[(n-k-j)*Binomial(n-k-1+j, j)*k^j: j in [0..n-k-1]])/(n-k);
      end if;
    end function;
    [A064094(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Sep 27 2024
    
  • Mathematica
    T[n_, 0] = 1; T[n_, 1] := CatalanNumber[n - 1]; T[n_, n_] = 1; T[n_, m_] := (1/(1 - m))^(n - m)*(1 - m*Sum[ CatalanNumber[k]*(m*(1 - m))^k, {k, 0, n - m - 1}]); Table[ T[n, m], {n, 0, 10}, {m, 0, n}] // Flatten (* Jean-François Alcover, Jul 05 2013 *)
  • SageMath
    def A064094(n,k):
        if (k==0 or k==n): return 1
        else: return sum((n-k-j)*binomial(n-k-1+j,j)*k^j for j in range(n-k))//(n-k)
    flatten([[A064094(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Sep 27 2024

Formula

G.f. for column m: (x^m)/(1-x*c(m*x)) = (x^m)*((m-1)+m*x*c(m*x))/(m-1+x) with the g.f. c(x) of Catalan numbers A000108.
T(n, m) = Sum_{j=0..n-m-1} (n-m-j)*binomial(n-m-1+j, j)*(m^j)/(n-m) or T(n, m) = (1/(1-m))^(n-m)*(1 - m*Sum_{j=0..n-m-1} C(j)*(m*(1-m))^j ), for n - m >= 1, T(n, n) = 1, T(n, m) = 0 if nA000108(k) (Catalan).

A119259 Central terms of the triangle in A119258.

Original entry on oeis.org

1, 3, 17, 111, 769, 5503, 40193, 297727, 2228225, 16807935, 127574017, 973168639, 7454392321, 57298911231, 441739706369, 3414246490111, 26447737520129, 205272288591871, 1595964714385409, 12427568655368191, 96905907580960769, 756583504975757311, 5913649000782757889
Offset: 0

Author

Reinhard Zumkeller, May 11 2006

Keywords

Comments

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. - Peter Bala, Jan 06 2022

References

  • R. P. Stanley, Enumerative Combinatorics Volume 2, Cambridge Univ. Press, 1999, Theorem 6.33, p. 197.

Programs

  • Haskell
    a119259 n = a119258 (2 * n) n  -- Reinhard Zumkeller, Aug 06 2014
    
  • Mathematica
    Table[Binomial[2k - 1, k] Hypergeometric2F1[-2k, -k, 1 - 2k, -1], {k, 0, 10}] (* Vladimir Reshetnikov, Feb 16 2011 *)
  • Python
    from itertools import count, islice
    def A119259_gen(): # generator of terms
        yield from (1,3)
        a, c = 2, 1
        for n in count(1):
            yield (a<>1
    A119259_list = list(islice(A119259_gen(),20)) # Chai Wah Wu, Apr 26 2023

Formula

a(n) = A119258(2*n,n).
a(n) = Sum_{k=0..n} C(2*n,k)*C(2*n-k-1,n-k). - Paul Barry, Sep 28 2007
a(n) = Sum_{k=0..n} C(n+k-1,k)*2^k. - Paul Barry, Sep 28 2007
2*a(n) = A064062(n)+A178792(n). - Joseph Abate, Jul 21 2010
G.f.: (4*x^2+3*sqrt(1-8*x)*x-5*x)/(sqrt(1-8*x)*(2*x^2+x-1)-8*x^2-7*x+1). - Vladimir Kruchinin, Aug 19 2013
a(n) = (-1)^n - 2^(n+1)*binomial(2*n,n-1)*hyper2F1([1,2*n+1],[n+2],2). - Peter Luschny, Jul 25 2014
a(n) = (-1)^n + 2^(n+1)*A014300(n). - Peter Luschny, Jul 25 2014
a(n) = [x^n] ( (1 + x)^2/(1 - x) )^n. Exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + 3*x + 13*x^2 + 67*x^3 + ... is essentially the o.g.f. for A064062. - Peter Bala, Oct 01 2015
The o.g.f. is the diagonal of the bivariate rational function 1/(1 - t*(1 + x)^2/(1 - x)) and hence is algebraic by Stanley 1999, Theorem 6.33, p.197. - Peter Bala, Aug 21 2016
n*(3*n-4)*a(n) +(-21*n^2+40*n-12)*a(n-1) -4*(3*n-1)*(2*n-3)*a(n-2)=0. - R. J. Mathar, Aug 09 2017
From Peter Bala, Mar 23 2020: (Start)
a(p) == 3 ( mod p^3 ) for prime p >= 5. Cf. A002003, A103885 and A156894.
More generally, we conjecture that a(n*p^k) == a(n*p^(k-1)) ( mod p^(3*k) ) for prime p >= 5 and positive integers n and k. (End)
G.f.: (8*x)/(sqrt(1-8*x)*(1+4*x)-1+8*x). - Fabian Pereyra, Jul 20 2024
a(n) = 2^(n+1)*binomial(2*n,n) - A178792(n). - Akiva Weinberger, Dec 06 2024
a(n) = Sum_{k=0..n} 2^k * (-1)^(n-k) * binomial(2*n,k). - Seiichi Manyama, Jul 31 2025
Showing 1-10 of 28 results. Next