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 10 results.

A180966 Hankel transform of A123164.

Original entry on oeis.org

1, 4, 28, 384, 10496, 573440, 62652416, 13690208256, 5982889443328, 5229277301702656, 9141181343655264256, 31958984107701798174720, 223467104335874481157308416, 3125102257923487167715657908224
Offset: 0

Views

Author

Paul Barry, Sep 29 2010

Keywords

Crossrefs

Programs

  • Magma
    [ 2^Binomial(n,2)*(&+[ (-1)^k*Binomial(n-k,k)*2^(2*n-3*k): k in [0..Floor(n/2)]]): n in [0..20]]; // G. C. Greubel, Apr 06 2021
    
  • Mathematica
    a[n_] := 2^Binomial[n, 2] Sum[Binomial[n-k, k] (-2)^k 4^(n-2k), {k, 0, n/2} ]; Table[a[n], {n, 0, 13}] (* Jean-François Alcover, Jun 17 2019 *)
    Table[2^(n^2/2)*ChebyshevU[n, Sqrt[2]], {n,0,20}] (* G. C. Greubel, Apr 06 2021 *)
  • Sage
    [2^(n^2/2)*chebyshev_U(n, sqrt(2)) for n in (0..20)] # G. C. Greubel, Apr 06 2021

Formula

a(n) = 2^C(n,2)*Sum_{k=0..floor(n/2)} C(n-k,k)*(-2)^k*4^(n-2*k).
a(n) = 2^C(n,2)*[x^n] (1/(1 - 4*x + 2*x^2)).
a(n) = 2^(2*n + ((n-1)*n)/2)*Hyper2F1([(1-n)/2, -n/2], [-n], 1/2) for n > 0. - Peter Luschny, Aug 02 2014
a(n) ~ 2^(n^2/2 - 1) * (1 + sqrt(2))^(n+1). - Vaclav Kotesovec, Feb 14 2021
a(n) = 2^(n^2/2)*ChebyshevU(n, sqrt(2)) = 2^(n*(n-1)/2)*A007070(n). - G. C. Greubel, Apr 06 2021

A006318 Large Schröder numbers (or large Schroeder numbers, or big Schroeder numbers).

Original entry on oeis.org

1, 2, 6, 22, 90, 394, 1806, 8558, 41586, 206098, 1037718, 5293446, 27297738, 142078746, 745387038, 3937603038, 20927156706, 111818026018, 600318853926, 3236724317174, 17518619320890, 95149655201962, 518431875418926, 2832923350929742, 15521467648875090
Offset: 0

Views

Author

Keywords

Comments

For the little Schröder numbers (or little Schroeder numbers, or small Schroeder numbers) see A001003.
The number of perfect matchings in a triangular grid of n squares (n = 1, 4, 9, 16, 25, ...). - Roberto E. Martinez II, Nov 05 2001
a(n) is the number of subdiagonal paths from (0, 0) to (n, n) consisting of steps East (1, 0), North (0, 1) and Northeast (1, 1) (sometimes called royal paths). - David Callan, Mar 14 2004
Twice A001003 (except for the first term).
a(n) is the number of dissections of a regular (n+4)-gon by diagonals that do not touch the base. (A diagonal is a straight line joining two nonconsecutive vertices and dissection means the diagonals are noncrossing though they may share an endpoint. One side of the (n+4)-gon is designated the base.) Example: a(1)=2 because a pentagon has only 2 such dissections: the empty one and the one with a diagonal parallel to the base. - David Callan, Aug 02 2004
a(n) is the number of separable permutations, i.e., permutations avoiding 2413 and 3142 (see Shapiro and Stephens). - Vincent Vatter, Aug 16 2006
Eric W. Weisstein comments that the Schröder numbers bear the same relationship to the Delannoy numbers (A001850) as the Catalan numbers (A000108) do to the binomial coefficients. - Jonathan Vos Post, Dec 23 2004
a(n) is the number of lattice paths from (0, 0) to (n+1, n+1) consisting of unit steps north N = (0, 1) and variable-length steps east E = (k, 0), with k a positive integer, that stay strictly below the line y = x except at the endpoints. For example, a(2) = 6 counts 111NNN, 21NNN, 3NNN, 12NNN, 11N1NN, 2N1NN (east steps indicated by their length). If the word "strictly" is replaced by "weakly", the counting sequence becomes the little Schröder numbers, A001003 (offset). - David Callan, Jun 07 2006
a(n) is the number of dissections of a regular (n+3)-gon with base AB that do not contain a triangle of the form ABP with BP a diagonal. Example: a(1) = 2 because the square D-C | | A-B has only 2 such dissections: the empty one and the one with the single diagonal AC (although this dissection contains the triangle ABC, BC is not a diagonal). - David Callan, Jul 14 2006
a(n) is the number of (colored) Motzkin n-paths with each upstep and each flatstep at ground level getting one of 2 colors and each flatstep not at ground level getting one of 3 colors. Example: With their colors immediately following upsteps/flatsteps, a(2) = 6 counts U1D, U2D, F1F1, F1F2, F2F1, F2F2. - David Callan, Aug 16 2006
The Hankel transform of this sequence is A006125(n+1) = [1, 2, 8, 64, 1024, 32768, ...]; example: Det([1, 2, 6, 22; 2, 6, 22, 90; 6, 22, 90, 394; 22, 90, 394, 1806]) = 64. - Philippe Deléham, Sep 03 2006
Triangle A144156 has row sums equal to A006318 with left border A001003. - Gary W. Adamson, Sep 12 2008
a(n) is also the number of order-preserving and order-decreasing partial transformations (of an n-chain). Equivalently, it is the order of the Schröder monoid, PC sub n. - Abdullahi Umar, Oct 02 2008
Sum_{n >= 0} a(n)/10^n - 1 = (9 - sqrt(41))/2. - Mark Dols, Jun 22 2010
1/sqrt(41) = Sum_{n >= 0} Delannoy number(n)/10^n. - Mark Dols, Jun 22 2010
a(n) is also the dimension of the space Hoch(n) related to Hochschild two-cocycles. - Ph. Leroux (ph_ler_math(AT)yahoo.com), Aug 24 2010
Let W = (w(n, k)) denote the augmentation triangle (as at A193091) of A154325; then w(n, n) = A006318(n). - Clark Kimberling, Jul 30 2011
Conjecture: For each n > 2, the polynomial sum_{k = 0}^n a(k)*x^{n-k} is irreducible modulo some prime p < n*(n+1). - Zhi-Wei Sun, Apr 07 2013
From Jon Perry, May 24 2013: (Start)
Consider a Pascal triangle variant where T(n, k) = T(n, k-1) + T(n-1, k-1) + T(n-1, k), i.e., the order of performing the calculation must go from left to right (A033877). This sequence is the rightmost diagonal.
Triangle begins:
1;
1, 2;
1, 4, 6;
1, 6, 16, 22;
1, 8, 30, 68, 90;
... (End)
a(n) is the number of permutations avoiding 2143, 3142 and one of the patterns among 246135, 254613, 263514, 524361, 546132. - Alexander Burstein, Oct 05 2014
a(n) is the number of semi-standard Young tableaux of shape n x 2 with consecutive entries. That is, j in P and 1 <= i<= j imply i in P. - Graham H. Hawkes, Feb 15 2015
a(n) is the number of unary-rooted size n unary-binary trees (each node has either 1 or 2 degree out). - John Bodeen, May 29 2017
Conjecturally, a(n) is the number of permutations pi of length n such that s(pi) avoids the patterns 231 and 321, where s denotes West's stack-sorting map. - Colin Defant, Sep 17 2018
a(n) is the number of n X n permutation matrices which percolate under the 2-neighbor bootstrap percolation rule (see Shapiro and Stephens). The number of general n X n matrices of weight n which percolate is given in A146971. - Jonathan Noel, Oct 05 2018
a(n) is the number of permutations of length n+1 which avoid 3142 and 3241. The permutations are precisely the permutations that are sortable by a decreasing stack followed by an increasing stack in series. - Rebecca Smith, Jun 06 2019
a(n) is the number of permutations of length n+1 avoiding the partially ordered pattern (POP) {3>1, 4>1, 1>2} of length 4. That is, the number of length n+1 permutations having no subsequences of length 4 in which the second element is the smallest, and the first element is smaller than the third and fourth elements. - Sergey Kitaev, Dec 10 2020
Named after the German mathematician Ernst Schröder (1841-1902). - Amiram Eldar, Apr 15 2021
a(n) is the number of sequences of nonnegative integers (u_1, u_2, ..., u_n) such that (i) u_i <= i for all i, and (ii) the nonzero u_i are weakly increasing. For example, a(2) = 6 counts 00, 01, 02, 10, 11, 12. See link "Some bijections for lattice paths" at A001003. - David Callan, Dec 18 2021
a(n) is the number of separable elements of the Weyl group of type B_n/C_n (see Gaetz and Gao). - Fern Gossow, Jul 31 2023
The number of domino tilings of an Aztec triangle of order n. Dually, the number perfect matchings of the edges in the cellular graph formed by a triangular grid of n squares (n = 1, 4, 9, 16, 25, ...) as in Ciucu (1996). - Michael Somos, Sep 16 2024
a(n) is the number of dissections of a convex (n+3)-sided polygon by non-intersecting diagonals such that none of the dividing diagonals passes through a chosen vertex. - Muhammed Sefa Saydam, Mar 01 2025
a(n) is the number of dissections of a convex (n+m+1)-sided polygon by non-intersecting diagonals such that the selected m consecutive sides of the polygon will be in the same subpolygon. - Muhammed Sefa Saydam, Jul 02 2025

Examples

			a(3) = 22 since the top row of Q^n = (6, 6, 6, 4, 0, 0, 0, ...); where 22 = (6 + 6 + 6 + 4).
G.f. = 1 + 2*x + 6*x^2 + 22*x^3 + 90*x^4 + 394*x^5 + 1806*x^6 + 8858*x^7 + 41586*x^8 + ...
		

References

  • D. Andrica and E. J. Ionascu, On the number of polynomials with coefficients in [n], An. St. Univ. Ovidius Constanta, 2013, to appear.
  • Paul Barry, On Integer-Sequence-Based Constructions of Generalized Pascal Triangles, Journal of Integer Sequences, Vol. 9 (2006), Article 06.2.4.
  • Paul Barry, Riordan-Bernstein Polynomials, Hankel Transforms and Somos Sequences, Journal of Integer Sequences, Vol. 15 2012, #12.8.2.
  • Paul Barry, Riordan arrays, generalized Narayana triangles, and series reversion, Linear Algebra and its Applications, 491 (2016) 343-385.
  • Paul Barry, Jacobsthal Decompositions of Pascal's Triangle, Ternary Trees, and Alternating Sign Matrices, Journal of Integer Sequences, 19, 2016, #16.3.5.
  • O. Bodini, A. Genitrini, F. Peschanski, and N.Rolin, Associativity for binary parallel processes, CALDAM 2015.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, pages 24, 618.
  • S. Brlek, E. Duchi, E. Pergola, and S. Rinaldi, On the equivalence problem for succession rules, Discr. Math., 298 (2005), 142-154.
  • Xiang-Ke Chang, XB Hu, H Lei, and YN Yeh, Combinatorial proofs of addition formulas, The Electronic Journal of Combinatorics, 23(1) (2016), #P1.8.
  • William Y. C. Chen and Carol J. Wang, Noncrossing Linked Partitions and Large (3, 2)-Motzkin Paths, Discrete Math., 312 (2012), 1918-1922.
  • L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 81, #21, (4), q_n.
  • D. E. Davenport, L. W. Shapiro, and L. C. Woodson, The Double Riordan Group, The Electronic Journal of Combinatorics, 18(2) (2012), #P33.
  • Deng, Eva Y. P.; Dukes, Mark; Mansour, Toufik; and Wu, Susan Y. J.; Symmetric Schröder paths and restricted involutions. Discrete Math. 309 (2009), no. 12, 4108-4115. See p. 4109.
  • E. Deutsch, A bijective proof of an equation linking the Schroeder numbers, large and small, Discrete Math., 241 (2001), 235-240.
  • C. Domb and A. J. Barrett, Enumeration of ladder graphs, Discrete Math. 9 (1974), 341-358.
  • Doslic, Tomislav and Veljan, Darko. Logarithmic behavior of some combinatorial sequences. Discrete Math. 308 (2008), no. 11, 2182--2212. MR2404544 (2009j:05019) - From N. J. A. Sloane, May 01 2012
  • M. Dziemianczuk, Generalizing Delannoy numbers via counting weighted lattice paths, INTEGERS, 13 (2013), #A54.
  • Egge, Eric S., Restricted signed permutations counted by the Schröder numbers. Discrete Math. 306 (2006), 552-563. [Many applications of these numbers.]
  • S. Getu et al., How to guess a generating function, SIAM J. Discrete Math., 5 (1992), 497-499.
  • S. Gire, Arbres, permutations a motifs exclus et cartes planaire: quelques problemes algorithmiques et combinatoires, Ph.D. Thesis, Universite Bordeaux I, 1993.
  • N. S. S. Gu, N. Y. Li, and T. Mansour, 2-Binary trees: bijections and related issues, Discr. Math., 308 (2008), 1209-1221.
  • Guruswami, Venkatesan, Enumerative aspects of certain subclasses of perfect graphs. Discrete Math. 205 (1999), 97-117.
  • Silvia Heubach and Toufik Mansour, Combinatorics of Compositions and Words, CRC Press, 2010.
  • D. E. Knuth, The Art of Computer Programming, Vol. 1, Section 2.2.1, Problem 11.
  • D. Kremer, Permutations with forbidden subsequences and a generalized Schröder number, Discrete Math. 218 (2000) 121-130.
  • Kremer, Darla and Shiu, Wai Chee; Finite transition matrices for permutations avoiding pairs of length four patterns. Discrete Math. 268 (2003), 171-183. MR1983276 (2004b:05006). See Table 1.
  • Laradji, A. and Umar, A. Asymptotic results for semigroups of order-preserving partial transformations. Comm. Algebra 34 (2006), 1071-1075. - Abdullahi Umar, Oct 11 2008
  • L. Moser and W. Zayachkowski, Lattice paths with diagonal steps, Scripta Math., 26 (1961), 223-229.
  • L. Shapiro and A. B. Stephens, Bootstrap percolation, the Schröder numbers and the N-kings problem, SIAM J. Discrete Math., Vol. 4 (1991), pp. 275-280.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see page 178 and also Problems 6.39 and 6.40.
  • Lin Yang and S.-L. Yang, The parametric Pascal rhombus. Fib. Q., 57:4 (2019), 337-346.
  • Sheng-Liang Yang and Mei-yang Jiang, The m-Schröder paths and m-Schröder numbers, Disc. Math. (2021) Vol. 344, Issue 2, 112209. doi:10.1016/j.disc.2020.112209. See Table 1.

Crossrefs

Apart from leading term, twice A001003 (the small Schroeder numbers). Cf. A025240.
Sequences A085403, A086456, A103137, A112478 are essentially the same sequence.
Main diagonal of A033877.
Row sums of A104219. Bisections give A138462, A138463.
Row sums of A175124.
The sequences listed in Yang-Jiang's Table 1 appear to be A006318, A001003, A027307, A034015, A144097, A243675, A260332, A243676. - N. J. A. Sloane, Mar 28 2021

Programs

  • GAP
    Concatenation([1],List([1..25],n->(1/n)*Sum([0..n],k->2^k*Binomial(n,k)*Binomial(n,k-1)))); # Muniru A Asiru, Nov 29 2018
  • Haskell
    a006318 n = a004148_list !! n
    a006318_list = 1 : f [1] where
       f xs = y : f (y : xs) where
         y = head xs + sum (zipWith (*) xs $ reverse xs)
    -- Reinhard Zumkeller, Nov 13 2012
    
  • Maple
    Order := 24: solve(series((y-y^2)/(1+y),y)=x,y); # then A(x)=y(x)/x
    BB:=(-1-z-sqrt(1-6*z+z^2))/2: BBser:=series(BB, z=0, 24): seq(coeff(BBser, z, n), n=1..23); # Zerinvary Lajos, Apr 10 2007
    A006318_list := proc(n) local j, a, w; a := array(0..n); a[0] := 1;
    for w from 1 to n do a[w] := 2*a[w-1]+add(a[j]*a[w-j-1], j=1..w-1) od; convert(a,list)end: A006318_list(22); # Peter Luschny, May 19 2011
    A006318 := n-> add(binomial(n+k, n-k) * binomial(2*k, k)/(k+1), k=0..n): seq(A006318(n), n=0..22); # Johannes W. Meijer, Jul 14 2013
    seq(simplify(hypergeom([-n,n+1],[2],-1)), n=0..100); # Robert Israel, Mar 23 2015
  • Mathematica
    a[0] = 1; a[n_Integer] := a[n] = a[n - 1] + Sum[a[k]*a[n - 1 - k], {k, 0, n - 1}]; Array[a[#] &, 30]
    InverseSeries[Series[(y - y^2)/(1 + y), {y, 0, 24}], x] (* then A(x) = y(x)/x *) (* Len Smiley, Apr 11 2000 *)
    CoefficientList[Series[(1 - x - (1 - 6x + x^2)^(1/2))/(2x), {x, 0, 30}], x] (* Harvey P. Dale, May 01 2011 *)
    a[ n_] := 2 Hypergeometric2F1[ -n + 1, n + 2, 2, -1]; (* Michael Somos, Apr 03 2013 *)
    a[ n_] := With[{m = If[ n < 0, -1 - n, n]}, SeriesCoefficient[(1 - x - Sqrt[ 1 - 6 x + x^2])/(2 x), {x, 0, m}]]; (* Michael Somos, Jun 10 2015 *)
    Table[-(GegenbauerC[n+1, -1/2, 3] + KroneckerDelta[n])/2, {n, 0, 30}] (* Vladimir Reshetnikov, Nov 12 2016 *)
    CoefficientList[Nest[1+x(#+#^2)&, 1+O[x], 20], x] (* Oliver Seipel, Dec 21 2024 *)
  • PARI
    {a(n) = if( n<0, n = -1-n); polcoeff( (1 - x - sqrt( 1 - 6*x + x^2 + x^2 * O(x^n))) / 2, n+1)}; /* Michael Somos, Apr 03 2013 */
    
  • PARI
    {a(n) = if( n<1, 1, sum( k=0, n, 2^k * binomial( n, k) * binomial( n, k-1)) / n)};
    
  • Python
    from gmpy2 import divexact
    A006318 = [1, 2]
    for n in range(3,10**3):
        A006318.append(int(divexact(A006318[-1]*(6*n-9)-(n-3)*A006318[-2],n)))
    # Chai Wah Wu, Sep 01 2014
    
  • Sage
    # Generalized algorithm of L. Seidel
    def A006318_list(n) :
        D = [0]*(n+1); D[1] = 1
        b = True; h = 1; R = []
        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]
                R.append(D[h-1]);
            b = not b
        return R
    A006318_list(23) # Peter Luschny, Jun 02 2012
    

Formula

G.f.: (1 - x - (1 - 6*x + x^2)^(1/2))/(2*x).
a(n) = 2*hypergeom([-n+1, n+2], [2], -1). - Vladeta Jovovic, Apr 24 2003
For n > 0, a(n) = (1/n)*Sum_{k = 0..n} 2^k*C(n, k)*C(n, k-1). - Benoit Cloitre, May 10 2003
The g.f. satisfies (1 - x)*A(x) - x*A(x)^2 = 1. - Ralf Stephan, Jun 30 2003
For the asymptotic behavior, see A001003 (remembering that A006318 = 2*A001003). - N. J. A. Sloane, Apr 10 2011
From Philippe Deléham, Nov 28 2003: (Start)
Row sums of A088617 and A060693.
a(n) = Sum_{k = 0..n} C(n+k, n)*C(n, k)/(k+1). (End)
With offset 1: a(1) = 1, a(n) = a(n-1) + Sum_{i = 1..n-1} a(i)*a(n-i). - Benoit Cloitre, Mar 16 2004
a(n) = Sum_{k = 0..n} A000108(k)*binomial(n+k, n-k). - Benoit Cloitre, May 09 2004
a(n) = Sum_{k = 0..n} A011117(n, k). - Philippe Deléham, Jul 10 2004
a(n) = (CentralDelannoy(n+1) - 3 * CentralDelannoy(n))/(2*n) = (-CentralDelannoy(n+1) + 6 * CentralDelannoy(n) - CentralDelannoy(n-1))/2 for n >= 1, where CentralDelannoy is A001850. - David Callan, Aug 16 2006
From Abdullahi Umar, Oct 11 2008: (Start)
A123164(n+1) - A123164(n) = (2*n+1)*a(n) (n >= 0).
and 2*A123164(n) = (n+1)*a(n) - (n-1)*a(n-1) (n > 0). (End)
Define the general Delannoy numbers d(i, j) as in A001850. Then a(k) = d(2*k, k) - d(2*k, k-1) and a(0) = 1, Sum_{j=0..n} ((-1)^j * (d(n, j) + d(n-1, j-1)) * a(n-j)) = 0. - Peter E John, Oct 19 2006
Given an integer t >= 1 and initial values u = [a_0, a_1, ..., a_{t-1}], we may define an infinite sequence Phi(u) by setting a_n = a_{n-1} + a_0*a_{n-1} + a_1*a_{n-2} + ... + a_{n-2}*a_1 for n >= t. For example, Phi([1]) is the Catalan numbers A000108. The present sequence is (essentially) Phi([2]). - Gary W. Adamson, Oct 27 2008
G.f.: 1/(1-2x/(1-x/(1-2x/(1-x/(1-2x/(1-x/(1-2x/(1-x/(1-2x/(1-x.... (continued fraction). - Paul Barry, Dec 08 2008
G.f.: 1/(1 - x - x/(1 - x - x/(1 - x - x/(1 - x - x/(1 - x - x/(1 - ... (continued fraction). - Paul Barry, Jan 29 2009
a(n) ~ ((3 + 2*sqrt(2))^n)/(n*sqrt(2*Pi*n)*sqrt(3*sqrt(2) - 4))*(1-(9*sqrt(2) + 24)/(32*n) + ...). - G. Nemes (nemesgery(AT)gmail.com), Jan 25 2009
Logarithmic derivative yields A002003. - Paul D. Hanna, Oct 25 2010
a(n) = the upper left term in M^(n+1), M = the production matrix:
1, 1, 0, 0, 0, 0, ...
1, 1, 1, 0, 0, 0, ...
2, 2, 1, 1, 0, 0, ...
4, 4, 2, 1, 1, 0, ...
8, 8, 8, 2, 1, 1, ...
... - Gary W. Adamson, Jul 08 2011
a(n) is the sum of top row terms in Q^n, Q = an infinite square production matrix as follows:
1, 1, 0, 0, 0, 0, ...
1, 1, 2, 0, 0, 0, ...
1, 1, 1, 2, 0, 0, ...
1, 1, 1, 1, 2, 0, ...
1, 1, 1, 1, 1, 2, ...
... - Gary W. Adamson, Aug 23 2011
From Tom Copeland, Sep 21 2011: (Start)
With F(x) = (1 - 3*x - sqrt(1 - 6*x + x^2))/(2*x) an o.g.f. (nulling the n = 0 term) for A006318, G(x) = x/(2 + 3*x + x^2) is the compositional inverse.
Consequently, with H(x) = 1/ (dG(x)/dx) = (2 + 3*x + x^2)^2 / (2 - x^2),
a(n) = (1/n!)*[(H(x)*d/dx)^n] x evaluated at x = 0, i.e.,
F(x) = exp[x*H(u)*d/du] u, evaluated at u = 0. Also, dF(x)/dx = H(F(x)). (End)
a(n-1) = number of ordered complete binary trees with n leaves having k internal vertices colored black, the remaining n - 1 - k internal vertices colored white, and such that each vertex and its rightmost child have different colors ([Drake, Example 1.6.7]). For a refinement of this sequence see A175124. - Peter Bala, Sep 29 2011
D-finite with recurrence: (n-2)*a(n-2) - 3*(2*n-1)*a(n-1) + (n+1)*a(n) = 0. - Vaclav Kotesovec, Oct 05 2012
G.f.: A(x) = (1 - x - sqrt(1 - 6*x + x^2))/(2*x) = (1 - G(0))/x; G(k) = 1 + x - 2*x/G(k+1); (continued fraction, 1-step). - Sergei N. Gladkovskii, Jan 04 2012
G.f.: A(x) = (1 - x - sqrt(1 - 6*x + x^2))/(2*x) = (G(0) - 1)/x; G(k) = 1 - x/(1 - 2/G(k+1)); (continued fraction, 2-step). - Sergei N. Gladkovskii, Jan 04 2012
a(n+1) = a(n) + Sum_{k=0..n} a(k)*(n-k). - Reinhard Zumkeller, Nov 13 2012
G.f.: 1/Q(0) where Q(k) = 1 + k*(1 - x) - x - x*(k+1)*(k+2)/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Mar 14 2013
a(-1-n) = a(n). - Michael Somos, Apr 03 2013
G.f.: 1/x - 1 - U(0)/x, where U(k) = 1 - x - x/U(k+1); (continued fraction). - Sergei N. Gladkovskii, Jul 16 2013
G.f.: (2 - 2*x - G(0))/(4*x), where G(k) = 1 + 1/( 1 - x*(6 - x)*(2*k - 1)/(x*(6 - x)*(2*k - 1) + 2*(k + 1)/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jul 16 2013
a(n) = 1/(n + 1) * (Sum_{j=0..n} C(n+j, j)*C(n+j+1, j+1)*(Sum_{k=0..n-j} (-1)^k*C(n+j+k, k))). - Graham H. Hawkes, Feb 15 2015
a(n) = hypergeom([-n, n+1], [2], -1). - Peter Luschny, Mar 23 2015
a(n) = sqrt(2) * LegendreP(n, -1, 3) where LegendreP is the associated Legendre function of the first kind (in Maple's notation). - Robert Israel, Mar 23 2015
G.f. A(x) satisfies: A(x) = Sum_{j>=0} x^j * Sum_{k=0..j} binomial(j,k)*A(x)^k. - Ilya Gutkovskiy, Apr 11 2019
From Peter Bala, May 13 2024: (Start)
a(n) = 2 * Sum_{k = 0..floor(n/2)} binomial(n, 2*k)*binomial(2*n-2*k, n)/(n-2*k+1) for n >= 1.
a(n) = Integral_{x = 0..1} Legendre_P(n, 2*x+1) dx. (End)
G.f. A(x) = 1/(1 - x) * c(x/(1-x)^2), where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the g.f. of the Catalan numbers A000108. - Peter Bala, Aug 29 2024

Extensions

Edited by Charles R Greathouse IV, Apr 20 2010

A122542 Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, 2, -1, 0, 0, 0, 0, 0, ...] 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, 2, 1, 0, 2, 4, 1, 0, 2, 8, 6, 1, 0, 2, 12, 18, 8, 1, 0, 2, 16, 38, 32, 10, 1, 0, 2, 20, 66, 88, 50, 12, 1, 0, 2, 24, 102, 192, 170, 72, 14, 1, 0, 2, 28, 146, 360, 450, 292, 98, 16, 1, 0, 2, 32, 198, 608, 1002, 912, 462, 128, 18, 1
Offset: 0

Views

Author

Philippe Deléham, Sep 19 2006, May 28 2007

Keywords

Comments

Riordan array (1, x*(1+x)/(1-x)). Rising and falling diagonals are the tribonacci numbers A000213, A001590.

Examples

			Triangle begins:
  1;
  0, 1;
  0, 2,  1;
  0, 2,  4,   1;
  0, 2,  8,   6,   1;
  0, 2, 12,  18,   8,    1;
  0, 2, 16,  38,  32,   10,   1;
  0, 2, 20,  66,  88,   50,  12,   1;
  0, 2, 24, 102, 192,  170,  72,  14,   1;
  0, 2, 28, 146, 360,  450, 292,  98,  16,  1;
  0, 2, 32, 198, 608, 1002, 912, 462, 128, 18, 1;
		

Crossrefs

Other versions: A035607, A113413, A119800, A266213.
Sums include: A000007, A001333 (row), A001590 (diagonal), A007483, A057077 (signed row), A078016 (signed diagonal), A086901, A091928, A104934, A122558, A122690.

Programs

  • Haskell
    a122542 n k = a122542_tabl !! n !! k
    a122542_row n = a122542_tabl !! n
    a122542_tabl = map fst $ iterate
       (\(us, vs) -> (vs, zipWith (+) ([0] ++ us ++ [0]) $
                          zipWith (+) ([0] ++ vs) (vs ++ [0]))) ([1], [0, 1])
    -- Reinhard Zumkeller, Jul 20 2013, Apr 17 2013
    
  • Magma
    function T(n, k) // T = A122542
      if k eq 0 then return 0^n;
      elif k eq n then return 1;
      else return T(n-1,k) + T(n-1,k-1) + T(n-2,k-1);
      end if;
    end function;
    [T(n, k): k in [0..n], n in [0..12]]; // G. C. Greubel, Oct 27 2024
  • Mathematica
    CoefficientList[#, y]& /@ CoefficientList[(1-x)/(1 - (1+y)x - y x^2) + O[x]^11, x] // Flatten (* Jean-François Alcover, Sep 09 2018 *)
    (* Second program *)
    T[n_, k_]:= T[n, k]= If[k==n, 1, If[k==0, 0, T[n-1,k-1] +T[n-1,k] +T[n-2,k- 1] ]]; (* T = A122542 *)
    Table[T[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Oct 27 2024 *)
  • Sage
    def A122542_row(n):
        @cached_function
        def prec(n, k):
            if k==n: return 1
            if k==0: return 0
            return prec(n-1,k-1)+2*sum(prec(n-i,k-1) for i in (2..n-k+1))
        return [prec(n, k) for k in (0..n)]
    for n in (0..10): print(A122542_row(n)) # Peter Luschny, Mar 16 2016
    

Formula

Sum_{k=0..n} x^k*T(n,k) = A000007(n), A001333(n), A104934(n), A122558(n), A122690(n), A091928(n) for x = 0, 1, 2, 3, 4, 5. - Philippe Deléham, Jan 25 2012
Sum_{k=0..n} 3^(n-k)*T(n,k) = A086901(n).
Sum_{k=0..n} 2^(n-k)*T(n,k) = A007483(n-1), n >= 1. - Philippe Deléham, Oct 08 2006
T(2*n, n) = A123164(n).
T(n, k) = T(n-1,k) + T(n-1,k-1) + T(n-2,k-1), n > 1. - Philippe Deléham, Jan 25 2012
G.f.: (1-x)/(1-(1+y)*x-y*x^2). - Philippe Deléham, Mar 02 2012
From G. C. Greubel, Oct 27 2024: (Start)
Sum_{k=0..n} (-1)^k*T(n, k) = A057077(n).
Sum_{k=0..floor(n/2)} T(n-k, k) = A001590(n+1).
Sum_{k=0..floor(n/2)} (-1)^k*T(n-k, k) = A078016(n). (End)

A103885 a(n) = [x^(2*n)] ((1 + x)/(1 - x))^n.

Original entry on oeis.org

1, 2, 16, 146, 1408, 14002, 142000, 1459810, 15158272, 158611106, 1669752016, 17664712562, 187641279616, 2000029880786, 21380213588848, 229129634462146, 2460955893981184, 26482855453375042, 285475524009208720, 3082024598888203090, 33319523640218177408
Offset: 0

Views

Author

Ralf Stephan, Feb 20 2005

Keywords

Comments

From Peter Bala, Mar 01 2020: (Start)
The recurrence given below can be rewritten in the form
(2*n+1)*(2*n+2)*P(2,n)*a(n+1) - (2*n-1)*(2*n-2)*P(2,-n)*a(n-1) = Q(2,n^2)*a(n), where the polynomial Q(2,n) = 4*(55*n^2 - 34*n + 3) and the polynomial P(2,n) = 5*n^2 - 5*n + 1 satisfies the symmetry condition P(2,n) = P(2,1-n) and has real zeros.
More generally, for fixed m = 1,2,3,..., we conjecture that the sequence b(n) := a(m*n) satisfies a recurrence of the form ( Product_{k = 1..2*m} (2*m*n + k) ) * P(2*m,n)*b(n+1) + (-1)^m*( Product_{k = 1..2*m} (2*m*n - k) ) * P(2*m,-n)*b(n-1) = Q(2*m,n^2)*b(n), where the polynomials P(2*m,n) and Q(2*m,n) have degree 2*m. Conjecturally, the polynomial P(2*m,n) = P(2*m,1-n) and has real zeros in the interval [0, 1]. The 4*m zeros of the polynomial Q(2*m,n^2) seem to belong to the interval [-1, 1] and 4*m - 2 of these zeros appear to be approximated by the rational numbers +- k/(3*m), where 1 <= k <= 3*m - 2, k not a multiple of 3. (End)

Crossrefs

Programs

  • Magma
    A103885:= func< n | n eq 0 select 1 else (&+[ Binomial(n, k)*Binomial(2*n+k-1, n-1): k in [0..n]]) >;
    [A103885(n): n in [0..40]]; // G. C. Greubel, Oct 27 2024
    
  • Maple
    a := n -> `if`(n=0, 1, 2*n*hypergeom([1 - 2*n, 1 - n], [2], 2)):
    seq(simplify(a(n)), n=0..17); # Peter Luschny, Dec 30 2019
    # Alternative (after Peter Bala ):
    gf := n -> ( (1 + x)/(1 - x) )^n: ser := n -> series(gf(n), x, 40):
    seq(coeff(ser(n), x, 2*n), n=0..17); # Peter Luschny, Mar 20 2020
  • Mathematica
    Prepend[Table[Sum[2^i Binomial[n, i] Binomial[2n-1, i-1], {i, 1, 2n}], {n,1,20}], 1] (* Vaclav Kotesovec, Jul 01 2015 *)
  • PARI
    a(n) = if (n==0, 1, sum(i=0, n, 2^i * binomial(n, i) * binomial(2*n-1, i-1))); \\ Michel Marcus, Mar 21 2020
    
  • SageMath
    def A103885(n): return 1 if n==0 else sum(binomial(n, k)*binomial(2*n+k-1, n-1) for k in range(n+1))
    [A103885(n) for n in range(41)] # G. C. Greubel, Oct 27 2024

Formula

a(n) = Sum_{i=0..n} 2^i * binomial(n,i) * binomial(2*n-1,i-1). [Original definition, with summation range {i=1..n}.]
a(n) = A103884(n, n).
G.f.: A(x) = x*B(x)'/B(x), where B(x) is g.f. of A027307. - Vladimir Kruchinin, Jun 30 2015
From Vaclav Kotesovec, Jul 01 2015: (Start)
Recurrence: n*(2*n-1)*(5*n^2 - 15*n + 11)*a(n) = 2*(55*n^4 - 220*n^3 + 296*n^2 - 152*n + 24)*a(n-1) + (n-2)*(2*n-3)*(5*n^2 - 5*n + 1)*a(n-2).
a(n) ~ ((11 + 5*sqrt(5))/2)^n / (2 * 5^(1/4) * sqrt(Pi*n)). (End)
a(n) = [x^n] (1/(1 - x - x/(1 - x - x/(1 - x - x/(1 - x - x/(1 - ...))))))^n, a continued fraction. - Ilya Gutkovskiy, Sep 29 2017
a(n) = 2*n*hypergeom([1 - 2*n, 1 - n], [2], 2) for n >= 1. - Peter Luschny, Dec 30 2019
From Peter Bala, Mar 01 2020: (Start)
a(n) = Sum_{k = 0..n} C(n, k)*C(2*n+k-1, n-1), with a(0) = 1.
a(n) = Sum_{k = 0..n} C(2*n, 2*k)*C(2*n-k-1, n-1), with a(0) = 1.
a(n) = (1/2)*Sum_{k = 0..n} C(2*n, n-k)*C(2*n+k-1, k). Cf. A156894.
a(n) = [x^n] S(x)^n, where S(x) = (1 - x - sqrt(1 - 6*x + x^2))/(2*x) is the o.g.f. of the sequence of large Schröder numbers A006318.
a(n) = (1/2) * [x^(n)] ( (1 + x)/(1 - x) )^(2*n). Cf. A002003(n) = [x^n] ( (1 + x)/(1 - x) )^n.
Conjecture: a(n) = - [x^n] G(x)^(-n), where G(x) = 1 + 2*x + 14*x^2 + 134*x^3 + 1482*x^4 + ... is the o.g.f. of A144097.
a(p) == 2 ( mod p^3 ) for prime p >= 5. (End)
From Peter Bala, Sep 22 2021: (Start)
a(n) = Sum_{k = 0..n} 4^k*binomial(n+k-1,n)*binomial(n,k)^2 / binomial(2*k,k).
Equivalently, a(n) = [x^n] T(n,(1+x)/(1-x)), where T(n,x) is the n-th Chebyshev polynomial of the first kind. Cf. A103882. (End)
For n>0, a(n) = (1/3) * [x^n] (1/S(-x))^(3*n), where S(x) = (1 - x - sqrt(1 - 6*x + x^2))/(2*x) is the o.g.f. of the sequence of large Schröder numbers A006318. Cf. A370102. - Peter Bala, Jul 29 2024

Extensions

a(0) = 1 added and new name by Peter Bala, Mar 01 2020

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

Original entry on oeis.org

1, 3, 19, 138, 1059, 8378, 67582, 552576, 4563235, 37972290, 317894394, 2674398268, 22590697614, 191475925332, 1627653567916, 13870754053388, 118464647799075, 1013709715774130, 8689197042438274, 74594573994750972, 641252293546113434, 5519339268476249676, 47558930664216470628
Offset: 0

Views

Author

Paul Barry, Feb 17 2009

Keywords

Crossrefs

Programs

  • Magma
    A156894:= func< n | (&+[ Binomial(n,k)*Binomial(2*n+k-1,k): k in [0..n]]) >;
    [A156894(n): n in [0..30]]; // G. C. Greubel, Jan 06 2022
    
  • Maple
    a := n -> hypergeom([-n, 2*n], [1], -1);
    seq(round(evalf(a(n),32)), n=0..19); # Peter Luschny, Aug 02 2014
  • Mathematica
    Table[Sum[Binomial[n,k]Binomial[2n+k-1,k],{k,0,n}],{n,0,20}] (* Harvey P. Dale, Nov 12 2014 *)
  • PARI
    a(n) = if (n < 1, 1, sum(k=0, n, binomial(n,k)*binomial(2*n+k-1,k)));
    vector(50, n, a(n-1)) \\ Altug Alkan, Oct 05 2015
    
  • Sage
    [round( hypergeometric([-n, 2*n], [1], -1) ) for n in (0..30)] # G. C. Greubel, Jan 06 2022

Formula

a(n) = [x^n] ((1+x)/(1-x)^2)^n.
a(n) = (4*(n+1)*(2*n+1)*A003169(n+1) - (5*n+1)*(2*n-1)*A003169(n))/(17*n + 5) for n>0. - Mark van Hoeij, Jul 14 2010
a(n) = Hypergeometric2F1([-n, 2*n], [1], -1). - Peter Luschny, Aug 02 2014
Conjecture: 64*n*(2*n-1)*a(n) -16*(89*n^2 -134*n +63)*a(n-1) +4*(661*n^2 -2619*n +2576)*a(n-2) -3*(119*n^2 -713*n +1092)*a(n-3) +6*(2*n-7)*(n-4)*a(n-4) = 0. - R. J. Mathar, Feb 05 2015
Conjecture: 16*n*(782*n +5365)*(2*n-1)*a(n) +8*(3128*n^3 -362053*n^2 +593930*n -290328)*a(n-1) -3*(726869*n^3 -5105981*n^2 +11667946*n -8715544)*a(n-2) +158*(2*n-5)*(n-3)*(391*n -764)*a(n-3) = 0. - R. J. Mathar, Feb 05 2015
Conjecture: 4*n*(2*n-1)*(17*n^2 -52*n +39)*a(n) -(1207*n^4 -4899*n^3 +6692*n^2 -3504*n +576)*a(n-1) +2*(n-2)*(2*n-3)*(17*n^2 -18*n +4)*a(n-2) = 0. - R. J. Mathar, Feb 05 2015 [the Maple command sumrecursion (binomial(n,k) * binomial(2*n+k-1,k), k, a(n)) verifies this recurrence. - Peter Bala, Oct 05 2015 ]
a(n) ~ sqrt(578 + 306*sqrt(17)) * (71 + 17*sqrt(17))^n / (17 * sqrt(Pi*n) * 2^(4*n+2)). - Vaclav Kotesovec, Feb 05 2015
exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + 3*x + 14*x^2 + 79*x^3 + ... is the o.g.f. of A003169 (taken with offset 0). - Peter Bala, Oct 05 2015
From Peter Bala, Mar 20 2020: (Start)
a(p) == 3 ( mod p^3 ) for prime p >= 5. Cf. A002003, A103885 and A119259.
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)

A123160 Triangle read by rows: T(n,k) = n!*(n+k-1)!/((n-k)!*(n-1)!*(k!)^2) for 0 <= k <= n, with T(0,0) = 1.

Original entry on oeis.org

1, 1, 1, 1, 4, 3, 1, 9, 18, 10, 1, 16, 60, 80, 35, 1, 25, 150, 350, 350, 126, 1, 36, 315, 1120, 1890, 1512, 462, 1, 49, 588, 2940, 7350, 9702, 6468, 1716, 1, 64, 1008, 6720, 23100, 44352, 48048, 27456, 6435, 1, 81, 1620, 13860, 62370, 162162, 252252, 231660, 115830, 24310
Offset: 0

Views

Author

Roger L. Bagula, Oct 02 2006

Keywords

Comments

T(n,k) is also the number of order-preserving partial transformations (of an n-element chain) of width k (width(alpha) = |Dom(alpha)|). - Abdullahi Umar, Aug 25 2008

Examples

			Triangle begins:
  1;
  1,  1;
  1,  4,   3;
  1,  9,  18,  10;
  1, 16,  60,  80,  35;
  1, 25, 150, 350, 350, 126;
  ...
		

References

  • Frederick T. Wall, Chemical Thermodynamics, W. H. Freeman, San Francisco, 1965 pages 296 and 305

Crossrefs

Programs

  • Magma
    [Binomial(n,k)*Binomial(n+k-1,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Jun 19 2022
    
  • Maple
    T:=proc(n,k) if k=0 and n=0 then 1 elif k<=n then n!*(n+k-1)!/(n-k)!/(n-1)!/(k!)^2 else 0 fi end: for n from 0 to 10 do seq(T(n,k),k=0..n) od; # yields sequence in triangular form
  • Mathematica
    T[n_, m_]= If [n==m==0, 1, n!*(n+m-1)!/((n-m)!*(n-1)!(m!)^2)];
    Table[T[n, m], {n,0,10}, {m,0,n}]//Flatten
    max = 9; s = (x+1)/(2*Sqrt[(1-x)^2-4*y])+1/2 + O[x]^(max+2) + O[y]^(max+2); T[n_, k_] := SeriesCoefficient[s, {x, 0, n}, {y, 0, k}]; Table[T[n-k, k], {n, 0, max}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 18 2015, after Vladimir Kruchinin *)
  • SageMath
    def A123160(n,k): return binomial(n, k)*binomial(n+k-1, k)
    flatten([[A123160(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Jun 19 2022

Formula

T(n, m) = n!*(n + m - 1)!/((n - m)!*(n - 1)!(m!)^2), with T(0, 0) = 1.
T(n, k) = binomial(n,k)*binomial(n+k-1,k). The row polynomials (except the first) are (1+x)*P(n,0,1,2x+1), where P(n,a,b,x) denotes the Jacobi polynomial. The columns of this triangle give the diagonals of A122899. - Peter Bala, Jan 24 2008
T(n, k) = binomial(n,k)*(n+k-1)!/((n-1)!*k!).
T(n, k)= binomial(n,k)*binomial(n+k-1,n-1). - Abdullahi Umar, Aug 25 2008
G.f.: (x+1)/(2*sqrt((1-x)^2-4*y)) + 1/2. - Vladimir Kruchinin, Jun 16 2015
From _Peter Bala, Jul 20 2015: (Start)
O.g.f. (1 + x)/( 2*sqrt((1 - x)^2 - 4*x*y) ) + 1/2 = 1 + (1 + y)*x + (1 + 4*y + 3*y^2)*x^2 + ....
For n >= 1, the n-th row polynomial R(n,y) = (1 + y)*r(n-1,y), where r(n,y) is the n-th row polynomial of A178301.
exp( Sum_{n >= 1} R(n,y)*x^n/n ) = 1 + (1 + y)*x + (1 + 3*y + 2*y^2)*x^2 + ... is the o.g.f for A088617. (End)
From G. C. Greubel, Jun 19 2022: (Start)
T(n, n) = A088218(n).
T(n, n-1) = A037965(n).
T(n, n-2) = A085373(n-2).
Sum_{k=0..n} T(n, k) = A123164(n).
Sum_{k=0..floor(n/2)} T(n-k, k) = A005773(n). (End)

Extensions

Edited by N. J. A. Sloane, Oct 26 2006 and Jul 03 2008

A336521 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where T(n,k) is the coefficient of x^(k*n) in expansion of ( (1 + x)/(1 - x) )^n.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 8, 1, 1, 2, 16, 38, 1, 1, 2, 24, 146, 192, 1, 1, 2, 32, 326, 1408, 1002, 1, 1, 2, 40, 578, 4672, 14002, 5336, 1, 1, 2, 48, 902, 11008, 69002, 142000, 28814, 1, 1, 2, 56, 1298, 21440, 216002, 1038984, 1459810, 157184, 1, 1, 2, 64, 1766, 36992, 525002, 4320608, 15856206, 15158272, 864146, 1
Offset: 0

Views

Author

Seiichi Manyama, Jul 24 2020

Keywords

Examples

			Square array begins:
  1,    1,     1,     1,      1,      1, ...
  1,    2,     2,     2,      2,      2, ...
  1,    8,    16,    24,     32,     40, ...
  1,   38,   146,   326,    578,    902, ...
  1,  192,  1408,  4672,  11008,  21440, ...
  1, 1002, 14002, 69002, 216002, 525002, ...
		

Crossrefs

Column k=0-3 give A000012, A123164, A103885, A333715.
Main diagonal gives A336522.

Programs

  • Mathematica
    T[n_, 0] := 1; T[n_, k_] := Sum[Binomial[n, j] * Binomial[k*n + j - 1, n - 1], {j, 0, n}]; Table[T[k, n - k], {n, 0, 10}, {k, 0, n}] // Flatten (* Amiram Eldar, Jul 24 2020 *)

Formula

T(n,k) = (1/k) * [x^n] ( (1 + x)/(1 - x) )^(k*n) for k > 0 and n > 0.
T(n,k) = Sum_{j=0..n} binomial(n,j) * binomial(k*n+j-1,n-1).
T(n,k) = (1/k) * Sum_{j=0..n} binomial(k*n,n-j) * binomial(k*n+j-1,j) for k > 0 and n > 0.
T(n,k) = Sum_{j=1..n} 2^j * binomial(n,j) * binomial(k*n-1,j-1) for n > 0.
T(n,k) = binomial(k*n-1, n-1)*hypergeom([-n, k*n], [1+(k-1)*n], -1) for k > 0. - Stefano Spezia, Aug 09 2025

A111516 Triangle G(n,k) read by rows: number of order-preserving partial transformations (of an n-element totally ordered set) of waist k (waist(alpha) = max(Im(alpha))).

Original entry on oeis.org

1, 1, 1, 1, 3, 4, 1, 7, 12, 18, 1, 15, 32, 56, 88, 1, 31, 80, 160, 280, 450, 1, 63, 192, 432, 832, 1452, 2364, 1, 127, 448, 1120, 2352, 4424, 7700, 12642, 1, 255, 1024, 2816, 6400, 12896, 23872, 41456, 68464, 1, 511, 2304, 6912, 16896, 36288, 71136, 130176, 225648, 374274
Offset: 0

Views

Author

Abdullahi Umar, Aug 25 2008

Keywords

Examples

			Triangle G(n,k) (with rows n >= 0 and columns k = 0..n) begins:
  1;
  1,  1;
  1,  3,   4;
  1,  7,  12,  18;
  1, 15,  32,  56,  88;
  1, 31,  80, 160, 280,  450;
  1, 63, 192, 432, 832, 1452, 2364;
  ...
G(2,2) = 4 because there are exactly 4 order-preserving partial transformations (on a 2-element chain) of waist 2, namely: (1)->(2), (2)->(2), (1,2)->(1,2), and (1,2)->(2,2) - the mappings are coordinate-wise.
		

Crossrefs

Cf. A050146 (main diagonal), A055807, A123164 (row sums).

Formula

G(n,k) = Sum_{j=1..n} C(n,j)*C(k+j-2,j-1) for 1 <= k <= n. [Corrected by Petros Hadjicostas, Feb 13 2021]
G(n,k) = 2*G(n-1,k) - G(n-1,k-1) + G(n,k-1) for n >= 2 and 1 <= k <= n-1 with initial conditions G(n,0) = 1 for n >= 0, G(n,1) = 2^n - 1 for n >= 1, and G(n,n) = A050146(n) for n >= 2.
Sum_{k=0..n} G(n,k) = A123164(n).
From Petros Hadjicostas, Feb 13 2021: (Start)
G(n,k) = A055807(n+k,k) for 0 <= k <= n.
Bivariate o.g.f.: Sum_{n,k>=0} G(n,k)*x^n*y^k = ((2 - 2*y - 2*x*y + x*y^2) + 2*x*(y - 1)/(1 - x) + x*y*(2 - 3*y - 2*x*y + x*y^2)/sqrt(1 - 6*x*y + x^2*y^2))/(2*(1 - 2*x - y + x*y)). (End)

Extensions

G(7,5) corrected by and more terms from Petros Hadjicostas, Feb 13 2021

A378378 a(n) = Sum_{k=0..n} binomial(n,k) * binomial(n+3*k-1,3*k).

Original entry on oeis.org

1, 2, 16, 170, 1920, 22402, 266800, 3222634, 39328768, 483752258, 5987236816, 74474238698, 930212870784, 11659157743170, 146567181170160, 1847198697449770, 23332153206562816, 295286370825453442, 3743540075432798608, 47532529217041519658, 604366048841146280320
Offset: 0

Views

Author

Seiichi Manyama, Nov 24 2024

Keywords

Crossrefs

Main diagonal of A378318.

Programs

  • Mathematica
    a[n_]:=HypergeometricPFQ[{(1+n)/3,(2+n)/3,-n,n/3},{1/3,2/3,1},-1]; Array[a,21,0] (* Stefano Spezia, Nov 24 2024 *)
  • PARI
    a(n) = sum(k=0, n, binomial(n, k)*binomial(n+3*k-1, 3*k));

Formula

a(n) = hypergeom([(1+n)/3, (2+n)/3, -n, n/3], [1/3, 2/3, 1], -1). - Stefano Spezia, Nov 24 2024

A144066 T(n, k) is the number of order-preserving partial transformations (of an n-element chain) of height k (height(alpha) = |Im(alpha)|); triangle T read by rows.

Original entry on oeis.org

1, 1, 1, 1, 6, 1, 1, 21, 15, 1, 1, 60, 102, 28, 1, 1, 155, 490, 310, 45, 1, 1, 378, 1935, 2220, 735, 66, 1, 1, 889, 6741, 12285, 7315, 1491, 91, 1
Offset: 0

Views

Author

Abdullahi Umar, Sep 09 2008

Keywords

Comments

T(n, k) is also the number of elements in the Green's J-classes of the monoid of order-preserving partial transformations (of an n-element chain). Sum of rows of T(n, k) is A123164.

Examples

			Triangle T(n,k) (with rows n >= 0 and columns k = 0..n) begins:
  1;
  1,   1;
  1,   6,    1;
  1,  21,   15,     1;
  1,  60,  102,    28,    1;
  1, 155,  490,   310,   45,    1;
  1, 378, 1935,  2220,  735,   66,  1;
  1, 889, 6741, 12285, 7315, 1491, 91, 1;
  ...
T(2,1) = 6 because there are exactly 6 order-preserving partial transformations (on a 2-element chain) of height 1, namely: (1)->(1), (1)->(2), (2)->(1), (2)->(2), (1,2)->(1,1), and (1,2)->(2,2) -- the mappings are coordinate-wise.
		

Crossrefs

Formula

T(n,k) = C(n,k)*A112857(n,k) for 0 <= k <= n.
C(n-1,k-1)*T(n,k) = 2((n-k+1)/(n-k))*T(n-1,k) + C(n,k)*T(n-1,k-1). [This is wrong.]
From Petros Hadjicostas, Feb 14 2021: (Start)
T(n,k) = 2*n*T(n-1,k)/(n-k) + n*T(n-1,k-1)/k for 1 <= k <= n-1 with T(n,0) = T(n,n) = 1 for n >= 0.
T(n,1) = n*(2^n - 1) = A066524(n) for n >= 1.
T(n,n-1) = n*(2*n - 1) = A000384(n) for n >= 1.
T(n,n-2) = A076454(n-1) for n >= 2. (End)
Showing 1-10 of 10 results.