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

A001263 Triangle of Narayana numbers T(n,k) = C(n-1,k-1)*C(n,k-1)/k with 1 <= k <= n, read by rows. Also called the Catalan triangle.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 6, 6, 1, 1, 10, 20, 10, 1, 1, 15, 50, 50, 15, 1, 1, 21, 105, 175, 105, 21, 1, 1, 28, 196, 490, 490, 196, 28, 1, 1, 36, 336, 1176, 1764, 1176, 336, 36, 1, 1, 45, 540, 2520, 5292, 5292, 2520, 540, 45, 1, 1, 55, 825, 4950, 13860, 19404, 13860, 4950, 825
Offset: 1

Views

Author

N. J. A. Sloane

Keywords

Comments

Number of antichains (or order ideals) in the poset 2*(k-1)*(n-k) or plane partitions with rows <= k-1, columns <= n-k and entries <= 2. - Mitch Harris, Jul 15 2000
T(n,k) is the number of Dyck n-paths with exactly k peaks. a(n,k) = number of pairs (P,Q) of lattice paths from (0,0) to (k,n+1-k), each consisting of unit steps East or North, such that P lies strictly above Q except at the endpoints. - David Callan, Mar 23 2004
Number of permutations of [n] which avoid-132 and have k-1 descents. - Mike Zabrocki, Aug 26 2004
T(n,k) is the number of paths through n panes of glass, entering and leaving from one side, of length 2n with k reflections (where traversing one pane of glass is the unit length). - Mitch Harris, Jul 06 2006
Antidiagonal sums given by A004148 (without first term).
T(n,k) is the number of full binary trees with n internal nodes and k-1 jumps. 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. - Emeric Deutsch, Jan 18 2007
From Gary W. Adamson, Oct 22 2007: (Start)
The n-th row can be generated by the following operation using an ascending row of (n-1) triangular terms, (A) and a descending row, (B); e.g., row 6:
A: 1....3....6....10....15
B: 15...10....6.....3.....1
C: 1...15...50....50....15....1 = row 6.
Leftmost column of A,B -> first two terms of C; then followed by the operation B*C/A of current column = next term of row C, (e.g., 10*15/3 = 50). Continuing with the operation, we get row 6: (1, 15, 50, 50, 15, 1). (End)
The previous comment can be upgraded to: The ConvOffsStoT transform of the triangular series; and by rows, row 6 is the ConvOffs transform of (1, 3, 6, 10, 15). Refer to triangle A117401 as another example of the ConvOffsStoT transform, and OEIS under Maple Transforms. - Gary W. Adamson, Jul 09 2012
For a connection to Lagrange inversion, see A134264. - Tom Copeland, Aug 15 2008
T(n,k) is also the number of order-decreasing and order-preserving mappings (of an n-element set) of height k (height of a mapping is the cardinal of its image set). - Abdullahi Umar, Aug 21 2008
Row n of this triangle is the h-vector of the simplicial complex dual to an associahedron of type A_n [Fomin & Reading, p.60]. See A033282 for the corresponding array of f-vectors for associahedra of type A_n. See A008459 and A145903 for the h-vectors for associahedra of type B and type D respectively. The Hilbert transform of this triangle (see A145905 for the definition of this transform) is A145904. - Peter Bala, Oct 27 2008
T(n,k) is also the number of noncrossing set partitions of [n] into k blocks. Given a partition P of the set {1,2,...,n}, a crossing in P are four integers [a, b, c, d] with 1 <= a < b < c < d <= n for which a, c are together in a block, and b, d are together in a different block. A noncrossing partition is a partition with no crossings. - Peter Luschny, Apr 29 2011
Noncrossing set partitions are also called genus 0 partitions. In terms of genus-dependent Stirling numbers of the second kind S2(n,k,g) that count partitions of genus g of an n-set into k nonempty subsets, one has T(n,k) = S2(n,k,0). - Robert Coquereaux, Feb 15 2024
Diagonals of A089732 are rows of A001263. - Tom Copeland, May 14 2012
From Peter Bala, Aug 07 2013: (Start)
Let E(y) = Sum_{n >= 0} y^n/(n!*(n+1)!) = 1/sqrt(y)*BesselI(1,2*sqrt(y)). Then this triangle is the generalized Riordan array (E(y), y) with respect to the sequence n!*(n+1)! as defined in Wang and Wang.
Generating function E(y)*E(x*y) = 1 + (1 + x)*y/(1!*2!) + (1 + 3*x + x^2)*y^2/(2!*3!) + (1 + 6*x + 6*x^2 + x^3)*y^3/(3!*4!) + .... Cf. A105278 with a generating function exp(y)*E(x*y).
The n-th power of this array has a generating function E(y)^n*E(x*y). In particular, the matrix inverse A103364 has a generating function E(x*y)/E(y). (End)
T(n,k) is the number of nonintersecting n arches above the x axis, starting and ending on vertices 1 to 2n, with k being the number of arches starting on an odd vertice and ending on a higher even vertice. Example: T(3,2)=3 [16,25,34] [14,23,56] [12,36,45]. - Roger Ford, Jun 14 2014
Fomin and Reading on p. 31 state that the rows of the Narayana matrix are the h-vectors of the associahedra as well as its dual. - Tom Copeland, Jun 27 2017
The row polynomials P(n, x) = Sum_{k=1..n} T(n, k)*x^(k-1), together with P(0, x) = 1, multiplied by (n+1) are the numerator polynomials of the o.g.f.s of the diagonal sequences of the triangle A103371: G(n, x) = (n+1)*P(n, x)/(1 - x)^{2*n+1}, for n >= 0. This is proved with Lagrange's theorem applied to the Riordan triangle A135278 = (1/(1 - x)^2, x/(1 - x)). See an example below. - Wolfdieter Lang, Jul 31 2017
T(n,k) is the number of Dyck paths of semilength n with k-1 uu-blocks (pairs of consecutive up-steps). - Alexander Burstein, Jun 22 2020
In case you were searching for Narayama numbers, the correct spelling is Narayana. - N. J. A. Sloane, Nov 11 2020
Named after the Canadian mathematician Tadepalli Venkata Narayana (1930-1987). They were also called "Runyon numbers" after John P. Runyon (1922-2013) of Bell Telephone Laboratories, who used them in a study of a telephone traffic system. - Amiram Eldar, Apr 15 2021 The Narayana numbers were first studied by Percy Alexander MacMahon (see reference, Article 495) as pointed out by Bóna and Sagan (see link). - Peter Luschny, Apr 28 2022
From Andrea Arlette España, Nov 14 2022: (Start)
T(n,k) is the degree distribution of the paths towards synchronization in the transition diagram associated with the Laplacian system over the complete graph K_n, corresponding to ordered initial conditions x_1 < x_2 < ... < x_n.
T(n,k) for n=2N+1 and k=N+1 is the number of states in the transition diagram associated with the Laplacian system over the complete bipartite graph K_{N,N}, corresponding to ordered (x_1 < x_2 < ... < x_N and x_{N+1} < x_{N+2} < ... < x_{2N}) and balanced (Sum_{i=1..N} x_i/N = Sum_{i=N+1..2N} x_i/N) initial conditions. (End)
From Gus Wiseman, Jan 23 2023: (Start)
Also the number of unlabeled ordered rooted trees with n nodes and k leaves. See the link by Marko Riedel. For example, row n = 5 counts the following trees:
((((o)))) (((o))o) ((o)oo) (oooo)
(((o)o)) ((oo)o)
(((oo))) ((ooo))
((o)(o)) (o(o)o)
((o(o))) (o(oo))
(o((o))) (oo(o))
The unordered version is A055277. Leaves in standard ordered trees are counted by A358371. (End)

Examples

			The initial rows of the triangle are:
  [1] 1
  [2] 1,  1
  [3] 1,  3,   1
  [4] 1,  6,   6,    1
  [5] 1, 10,  20,   10,    1
  [6] 1, 15,  50,   50,   15,    1
  [7] 1, 21, 105,  175,  105,   21,   1
  [8] 1, 28, 196,  490,  490,  196,  28,  1
  [9] 1, 36, 336, 1176, 1764, 1176, 336, 36, 1;
  ...
For all n, 12...n (1 block) and 1|2|3|...|n (n blocks) are noncrossing set partitions.
Example of umbral representation:
  A007318(5,k)=[1,5/1,5*4/(2*1),...,1]=(1,5,10,10,5,1),
  so A001263(5,k)={1,b(5)/b(1),b(5)*b(4)/[b(2)*b(1)],...,1}
  = [1,30/2,30*20/(6*2),...,1]=(1,15,50,50,15,1).
  First = last term = b.(5!)/[b.(0!)*b.(5!)]= 1. - _Tom Copeland_, Sep 21 2011
Row polynomials and diagonal sequences of A103371: n = 4,  P(4, x) = 1 + 6*x + 6*x^2 + x^3, and the o.g.f. of fifth diagonal is G(4, x) = 5* P(4, x)/(1 - x)^9, namely [5, 75, 525, ...]. See a comment above. - _Wolfdieter Lang_, Jul 31 2017

References

  • Berman and Koehler, Cardinalities of finite distributive lattices, Mitteilungen aus dem Mathematischen Seminar Giessen, 121 (1976), pp. 103-124.
  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 196.
  • P. A. MacMahon, Combinatory Analysis, Vols. 1 and 2, Cambridge University Press, 1915, 1916; reprinted by Chelsea, 1960, Sect. 495.
  • T. V. Narayana, Lattice Path Combinatorics with Statistical Applications. Univ. Toronto Press, 1979, pp. 100-101.
  • A. Nkwanta, Lattice paths and RNA secondary structures, in African Americans in Mathematics, ed. N. Dean, Amer. Math. Soc., 1997, pp. 137-147.
  • T. K. Petersen, Eulerian Numbers, Birkhäuser, 2015, Chapter 2.
  • J. Riordan, Combinatorial Identities, Wiley, 1968, p. 17.
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 6.36(a) and (b).

Links

Crossrefs

Other versions are in A090181 and A131198. - Philippe Deléham, Nov 18 2007
Cf. variants: A181143, A181144. - Paul D. Hanna, Oct 13 2010
Row sums give A000108 (Catalan numbers), n>0.
Columns give A000217, A002415, A006542, A006857, A108679. A084938.
Cf. A000372, A002083, A056932, A056939, A056940, A056941, A065329, A073345.
A145596, A145597, A145598, A145599. - Peter Bala, Oct 22 2008
A008459 (h-vectors type B associahedra), A033282 (f-vectors type A associahedra), A145903 (h-vectors type D associahedra), A145904 (Hilbert transform). - Peter Bala, Oct 27 2008
Cf. A016098 and A189232 for numbers of crossing set partitions.
Cf. A243752.
Cf. A089231, A103371, A135278.
Cf. A002378, A007318, A010790, A089732, A105278, A119900, A132710, A134264.
Triangles of generalized binomial coefficients (n,k)_m (or generalized Pascal triangles) for m = 1,...,12: A007318 (Pascal), A001263, A056939, A056940, A056941, A142465, A142467, A142468, A174109, A342889, A342890, A342891.
Cf. A000081, A005043, A032027, A055277, A358371.

Programs

  • GAP
    Flat(List([1..11],n->List([1..n],k->Binomial(n-1,k-1)*Binomial(n,k-1)/k))); # Muniru A Asiru, Jul 12 2018
  • Haskell
    a001263 n k = a001263_tabl !! (n-1) !! (k-1)
a001263_row n = a001263_tabl !! (n-1)
a001263_tabl = zipWith dt a007318_tabl (tail a007318_tabl) where
   dt us vs = zipWith (-) (zipWith (*) us (tail vs))
                          (zipWith (*) (tail us ++ [0]) (init vs))
-- Reinhard Zumkeller, Oct 10 2013
  • Magma
    /* triangle */ [[Binomial(n-1,k-1)*Binomial(n,k-1)/k : k in [1..n]]: n in [1.. 15]]; // Vincenzo Librandi, Oct 19 2014
  • Maple
    A001263 := (n,k)->binomial(n-1,k-1)*binomial(n,k-1)/k;
a:=proc(n,k) option remember; local i; if k=1 or k=n then 1 else add(binomial(n+i-1, 2*k-2)*a(k-1,i),i=1..k-1); fi; end:
# Alternatively, as a (0,0)-based triangle:
R := n -> simplify(hypergeom([-n, -n-1], [2], x)): Trow := n -> seq(coeff(R(n,x),x,j), j=0..n): seq(Trow(n), n=0..9); # Peter Luschny, Mar 19 2018
  • Mathematica
    T[n_, k_] := If[k==0, 0, Binomial[n-1, k-1] Binomial[n, k-1] / k];
Flatten[Table[Binomial[n-1,k-1] Binomial[n,k-1]/k,{n,15},{k,n}]] (* Harvey P. Dale, Feb 29 2012 *)
TRow[n_] := CoefficientList[Hypergeometric2F1[1 - n, -n, 2, x], x];
Table[TRow[n], {n, 1, 11}] // Flatten (* Peter Luschny, Mar 19 2018 *)
aot[n_]:=If[n==1,{{}},Join@@Table[Tuples[aot/@c],{c,Join@@Permutations/@IntegerPartitions[n-1]}]];
Table[Length[Select[aot[n],Length[Position[#,{}]]==k&]],{n,2,9},{k,1,n-1}] (* Gus Wiseman, Jan 23 2023 *)
T[1, 1] := 1; T[n_, k_]/;1<=k<=n := T[n, k] = (2n/k-1) T[n-1,k-1] + T[n-1, k]; T[n_, k_] := 0; Flatten@Table[T[n, k], {n, 1, 11}, {k, 1, n}] (* Oliver Seipel, Dec 31 2024 *)
  • PARI
    {a(n, k) = if(k==0, 0, binomial(n-1, k-1) * binomial(n, k-1) / k)};
  • PARI
    {T(n,k)=polcoeff(polcoeff(exp(sum(m=1,n,sum(j=0,m,binomial(m,j)^2*y^j)*x^m/m) +O(x^(n+1))),n,x),k,y)} \\ Paul D. Hanna, Oct 13 2010
  • Sage
    @CachedFunction
def T(n, k):
    if k == n or k == 1: return 1
    if k <= 0 or k > n: return 0
    return binomial(n, 2) * (T(n-1, k)/((n-k)*(n-k+1)) + T(n-1, k-1)/(k*(k-1)))
for n in (1..9): print([T(n, k) for k in (1..n)])  # Peter Luschny, Oct 28 2014

Formula

a(n, k) = C(n-1, k-1)*C(n, k-1)/k for k!=0; a(n, 0)=0.
Triangle equals [0, 1, 0, 1, 0, 1, ...] DELTA [1, 0, 1, 0, 1, 0, 1, 0, ...] where DELTA is Deléham's operator defined in A084938.
0Mike Zabrocki, Aug 26 2004
T(n, k) = C(n, k)*C(n-1, k-1) - C(n, k-1)*C(n-1, k) (determinant of a 2 X 2 subarray of Pascal's triangle A007318). - Gerald McGarvey, Feb 24 2005
T(n, k) = binomial(n-1, k-1)^2 - binomial(n-1, k)*binomial(n-1, k-2). - David Callan, Nov 02 2005
a(n,k) = C(n,2) (a(n-1,k)/((n-k)*(n-k+1)) + a(n-1,k-1)/(k*(k-1))) a(n,k) = C(n,k)*C(n,k-1)/n. - Mitch Harris, Jul 06 2006
Central column = A000891, (2n)!*(2n+1)! / (n!*(n+1)!)^2. - Zerinvary Lajos, Oct 29 2006
G.f.: (1-x*(1+y)-sqrt((1-x*(1+y))^2-4*y*x^2))/(2*x) = Sum_{n>0, k>0} a(n, k)*x^n*y^k.
From Peter Bala, Oct 22 2008: (Start)
Relation with Jacobi polynomials of parameter (1,1):
Row n+1 generating polynomial equals 1/(n+1)*x*(1-x)^n*Jacobi_P(n,1,1,(1+x)/(1-x)). It follows that the zeros of the Narayana polynomials are all real and nonpositive, as noted above. O.g.f for column k+2: 1/(k+1) * y^(k+2)/(1-y)^(k+3) * Jacobi_P(k,1,1,(1+y)/(1-y)). Cf. A008459.
T(n+1,k) is the number of walks of n unit steps on the square lattice (i.e., each step in the direction either up (U), down (D), right (R) or left (L)) starting from the origin and finishing at lattice points on the x axis and which remain in the upper half-plane y >= 0 [Guy]. For example, T(4,3) = 6 counts the six walks RRL, LRR, RLR, UDL, URD and RUD, from the origin to the lattice point (1,0), each of 3 steps. Compare with tables A145596 - A145599.
Define a functional I on formal power series of the form f(x) = 1 + ax + bx^2 + ... by the following iterative process. Define inductively f^(1)(x) = f(x) and f^(n+1)(x) = f(x*f^(n)(x)) for n >= 1. Then set I(f(x)) = lim_{n -> infinity} f^(n)(x) in the x-adic topology on the ring of formal power series; the operator I may also be defined by I(f(x)) := 1/x*series reversion of x/f(x).
The o.g.f. for this array is I(1 + t*x + t*x^2 + t*x^3 + ...) = 1 + t*x + (t + t^2)*x^2 + (t + 3*t^2 + t^3)*x^3 + ... = 1/(1 - x*t/(1 - x/(1 - x*t/(1 - x/(1 - ...))))) (as a continued fraction). Cf. A108767, A132081 and A141618. (End)
G.f.: 1/(1-x-xy-x^2y/(1-x-xy-x^2y/(1-... (continued fraction). - Paul Barry, Sep 28 2010
E.g.f.: exp((1+y)x)*Bessel_I(1,2*sqrt(y)x)/(sqrt(y)*x). - Paul Barry, Sep 28 2010
G.f.: A(x,y) = exp( Sum_{n>=1} [Sum_{k=0..n} C(n,k)^2*y^k] * x^n/n ). - Paul D. Hanna, Oct 13 2010
With F(x,t) = (1-(1+t)*x-sqrt(1-2*(1+t)*x+((t-1)*x)^2))/(2*x) an o.g.f. in x for the Narayana polynomials in t, G(x,t) = x/(t+(1+t)*x+x^2) is the compositional inverse in x. Consequently, with H(x,t) = 1/ (dG(x,t)/dx) = (t+(1+t)*x+x^2)^2 / (t-x^2), the n-th Narayana polynomial in t is given by (1/n!)*((H(x,t)*D_x)^n)x evaluated at x=0, i.e., F(x,t) = exp(x*H(u,t)*D_u)u, evaluated at u = 0. Also, dF(x,t)/dx = H(F(x,t),t). - Tom Copeland, Sep 04 2011
With offset 0, A001263 = Sum_{j>=0} A132710^j / A010790(j), a normalized Bessel fct. May be represented as the Pascal matrix A007318, n!/[(n-k)!*k!], umbralized with b(n)=A002378(n) for n>0 and b(0)=1: A001263(n,k)= b.(n!)/{b.[(n-k)!]*b.(k!)} where b.(n!) = b(n)*b(n-1)...*b(0), a generalized factorial (see example). - Tom Copeland, Sep 21 2011
With F(x,t) = {1-(1-t)*x-sqrt[1-2*(1+t)*x+[(t-1)*x]^2]}/2 a shifted o.g.f. in x for the Narayana polynomials in t, G(x,t)= x/[t-1+1/(1-x)] is the compositional inverse in x. Therefore, with H(x,t)=1/(dG(x,t)/dx)=[t-1+1/(1-x)]^2/{t-[x/(1-x)]^2}, (see A119900), the (n-1)-th Narayana polynomial in t is given by (1/n!)*((H(x,t)*d/dx)^n)x evaluated at x=0, i.e., F(x,t) = exp(x*H(u,t)*d/du) u, evaluated at u = 0. Also, dF(x,t)/dx = H(F(x,t),t). - Tom Copeland, Sep 30 2011
T(n,k) = binomial(n-1,k-1)*binomial(n+1,k)-binomial(n,k-1)*binomial(n,k). - Philippe Deléham, Nov 05 2011
A166360(n-k) = T(n,k) mod 2. - Reinhard Zumkeller, Oct 10 2013
Damped sum of a column, in leading order: lim_{d->0} d^(2k-1) Sum_{N>=k} T(N,k)(1-d)^N=Catalan(n). - Joachim Wuttke, Sep 11 2014
Multiplying the n-th column by n! generates the revert of the unsigned Lah numbers, A089231. - Tom Copeland, Jan 07 2016
Row polynomials: (x - 1)^(n+1)*(P(n+1,(1 + x)/(x - 1)) - P(n-1,(1 + x)/(x - 1)))/((4*n + 2)), n = 1,2,... and where P(n,x) denotes the n-th Legendre polynomial. - Peter Bala, Mar 03 2017
The coefficients of the row polynomials R(n, x) = hypergeom([-n,-n-1], [2], x) generate the triangle based in (0,0). - Peter Luschny, Mar 19 2018
Multiplying the n-th diagonal by n!, with the main diagonal n=1, generates the Lah matrix A105278. With G equal to the infinitesimal generator of A132710, the Narayana triangle equals Sum_{n >= 0} G^n/((n+1)!*n!) = (sqrt(G))^(-1) * I_1(2*sqrt(G)), where G^0 is the identity matrix and I_1(x) is the modified Bessel function of the first kind of order 1. (cf. Sep 21 2011 formula also.) - Tom Copeland, Sep 23 2020
T(n,k) = T(n,k-1)*C(n-k+2,2)/C(k,2). - Yuchun Ji, Dec 21 2020
From Sergii Voloshyn, Nov 25 2024: (Start)
G.f.: F(x,y) = (1-x*(1+y)-sqrt((1-x*(1+y))^2-4*y*x^2))/(2*x) is the solution of the differential equation x^3 * d^2(x*F(x,y))/dx^2 = y * d^2(x*F(x,y))/dy^2.
Let E be the operator x*D*D, where D denotes the derivative operator d/dx. Then (1/(n! (1 + n)!)) * E^n(x/(1 - x)) = (row n generating polynomial)/(1 - x)^(2*n+1) = Sum_{k >= 0} C(n-1, k-1)*C(n, k-1)/k*x^k. For example, when n = 4 we have (1/4!/5!)*E^3(x/(1 - x)) = x (1 + 6 x + 6 x^2 + x^3)/(1 - x)^9. (End)

Extensions

Deleted certain dangerous or potentially dangerous links. - N. J. A. Sloane, Jan 30 2021

A108625 Square array, read by antidiagonals, where row n equals the crystal ball sequence for the A_n lattice.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 7, 5, 1, 1, 13, 19, 7, 1, 1, 21, 55, 37, 9, 1, 1, 31, 131, 147, 61, 11, 1, 1, 43, 271, 471, 309, 91, 13, 1, 1, 57, 505, 1281, 1251, 561, 127, 15, 1, 1, 73, 869, 3067, 4251, 2751, 923, 169, 17, 1, 1, 91, 1405, 6637, 12559, 11253, 5321, 1415, 217, 19, 1
Offset: 0

Views

Author

Paul D. Hanna, Jun 12 2005

Keywords

Comments

Compare to the corresponding array A108553 of crystal ball sequences for D_n lattice.
From Peter Bala, Jul 18 2008: (Start)
Row reverse of A099608.
This array has a remarkable relationship with the constant zeta(2). The row, column and diagonal entries of the array occur in series acceleration formulas for zeta(2).
For the entries in row n we have zeta(2) = 2*(1 - 1/2^2 + 1/3^2 - ... + (-1)^(n+1)/n^2) + (-1)^n*Sum_{k >= 1} 1/(k^2*T(n,k-1)*T(n,k)). For example, n = 4 gives zeta(2) = 2*(1 - 1/4 + 1/9 - 1/16) + 1/(1*21) + 1/(4*21*131) + 1/(9*131*471) + ... . See A142995 for further details.
For the entries in column k we have zeta(2) = (1 + 1/4 + 1/9 + ... + 1/k^2) + 2*Sum_{n >= 1} (-1)^(n+1)/(n^2*T(n-1,k)*T(n,k)). For example, k = 4 gives zeta(2) = (1 + 1/4 + 1/9 + 1/16) + 2*(1/(1*9) - 1/(4*9*61) + 1/(9*61*309) - ... ). See A142999 for further details.
Also, as consequence of Apery's proof of the irrationality of zeta(2), we have a series acceleration formula along the main diagonal of the table: zeta(2) = 5 * Sum_{n >= 1} (-1)^(n+1)/(n^2*T(n,n)*T(n-1,n-1)) = 5*(1/3 - 1/(2^2*3*19) + 1/(3^2*19*147) - ...).
There also appear to be series acceleration results along other diagonals. For example, for the main subdiagonal, calculation supports the result zeta(2) = 2 - Sum_{n >= 1} (-1)^(n+1)*(n^2+(2*n+1)^2)/(n^2*(n+1)^2*T(n,n-1)*T(n+1,n)) = 2 - 10/(2^2*7) + 29/(6^2*7*55) - 58/(12^2*55*471) + ..., while for the main superdiagonal we appear to have zeta(2) = 1 + Sum_{n >= 1} (-1)^(n+1)*((n+1)^2 + (2*n+1)^2)/(n^2*(n+1)^2*T(n-1,n)*T(n,n+1)) = 1 + 13/(2^2*5) - 34/(6^2*5*37) + 65/(12^2*37*309) - ... .
Similar series acceleration results hold for Apery's constant zeta(3) involving the crystal ball sequences for the product lattices A_n x A_n; see A143007 for further details. Similar results also hold between the constant log(2) and the crystal ball sequences of the hypercubic lattices A_1 x...x A_1 and between log(2) and the crystal ball sequences for lattices of type C_n ; see A008288 and A142992 respectively for further details. (End)
This array is the Hilbert transform of triangle A008459 (see A145905 for the definition of the Hilbert transform). - Peter Bala, Oct 28 2008

Examples

			Square array begins:
  1,   1,    1,     1,      1,       1,       1, ... A000012;
  1,   3,    5,     7,      9,      11,      13, ... A005408;
  1,   7,   19,    37,     61,      91,     127, ... A003215;
  1,  13,   55,   147,    309,     561,     923, ... A005902;
  1,  21,  131,   471,   1251,    2751,    5321, ... A008384;
  1,  31,  271,  1281,   4251,   11253,   25493, ... A008386;
  1,  43,  505,  3067,  12559,   39733,  104959, ... A008388;
  1,  57,  869,  6637,  33111,  124223,  380731, ... A008390;
  1,  73, 1405, 13237,  79459,  350683, 1240399, ... A008392;
  1,  91, 2161, 24691, 176251,  907753, 3685123, ... A008394;
  1, 111, 3191, 43561, 365751, 2181257, ...      ... A008396;
  ...
As a triangle:
  [0]  1
  [1]  1,  1
  [2]  1,  3,   1
  [3]  1,  7,   5,    1
  [4]  1, 13,  19,    7,    1
  [5]  1, 21,  55,   37,    9,    1
  [6]  1, 31, 131,  147,   61,   11,   1
  [7]  1, 43, 271,  471,  309,   91,  13,   1
  [8]  1, 57, 505, 1281, 1251,  561, 127,  15,  1
  [9]  1, 73, 869, 3067, 4251, 2751, 923, 169, 17, 1
       ...
Inverse binomial transform of rows yield rows of triangle A063007:
  1;
  1,  2;
  1,  6,   6;
  1, 12,  30,  20;
  1, 20,  90, 140,  70;
  1, 30, 210, 560, 630, 252; ...
Product of the g.f. of row n and (1-x)^(n+1) generates the symmetric triangle A008459:
  1;
  1,  1;
  1,  4,   1;
  1,  9,   9,   1;
  1, 16,  36,  16,  1;
  1, 25, 100, 100, 25, 1;
  ...

Links

Crossrefs

Rows include: A003215 (row 2), A005902 (row 3), A008384 (row 4), A008386 (row 5), A008388 (row 6), A008390 (row 7), A008392 (row 8), A008394 (row 9), A008396 (row 10).
Cf. A063007, A099601 (n-th term of A_{2n} lattice), A108553.
Cf. A008459 (h-vectors type B associahedra), A145904, A145905.
Cf. A005258 (main diagonal), A108626 (antidiagonal sums).
Cf. A108628, A208675, A363867, A363868, A363869, A363870, A363871.

Programs

  • Magma
    T:= func< n,k | (&+[Binomial(n,j)^2*Binomial(n+k-j,k-j): j in [0..k]]) >; // array
A108625:= func< n,k | T(n-k,k) >; // antidiagonals
[A108625(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Oct 05 2023
  • Maple
    T := (n,k) -> binomial(n, k)*hypergeom([-k, k - n, k - n], [1, -n], 1):
seq(seq(simplify(T(n,k)),k=0..n),n=0..10); # Peter Luschny, Feb 10 2018
  • Mathematica
    T[n_, k_]:= HypergeometricPFQ[{-n, -k, n+1}, {1, 1}, 1] (* Michael Somos, Jun 03 2012 *)
  • PARI
    T(n,k)=sum(i=0,k,binomial(n,i)^2*binomial(n+k-i,k-i))
  • SageMath
    def T(n,k): return sum(binomial(n,j)^2*binomial(n+k-j, k-j) for j in range(k+1)) # array
def A108625(n,k): return T(n-k, k) # antidiagonals
flatten([[A108625(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Oct 05 2023

Formula

T(n, k) = Sum_{i=0..k} C(n, i)^2 * C(n+k-i, k-i).
G.f. for row n: (Sum_{i=0..n} C(n, i)^2 * x^i)/(1-x)^(n+1).
Sum_{k=0..n} T(n-k, k) = A108626(n) (antidiagonal sums).
From Peter Bala, Jul 23 2008 (Start):
O.g.f. row n: 1/(1 - x)*Legendre_P(n,(1 + x)/(1 - x)).
G.f. for square array: 1/sqrt((1 - x)*((1 - t)^2 - x*(1 + t)^2)) = (1 + x + x^2 + x^3 + ...) + (1 + 3*x + 5*x^2 + 7*x^3 + ...)*t + (1 + 7*x + 19*x^2 + 37*x^3 + ...)*t^2 + ... . Cf. A142977.
Main diagonal is A005258.
Recurrence relations:
Row n entries: (k+1)^2*T(n,k+1) = (2*k^2+2*k+n^2+n+1)*T(n,k) - k^2*T(n,k-1), k = 1,2,3,... ;
Column k entries: (n+1)^2*T(n+1,k) = (2*k+1)*(2*n+1)*T(n,k) + n^2*T(n-1,k), n = 1,2,3,... ;
Main diagonal entries: (n+1)^2*T(n+1,n+1) = (11*n^2+11*n+3)*T(n,n) + n^2*T(n-1,n-1), n = 1,2,3,... .
Series acceleration formulas for zeta(2):
Row n: zeta(2) = 2*(1 - 1/2^2 + 1/3^2 - ... + (-1)^(n+1)/n^2) + (-1)^n*Sum_{k >= 1} 1/(k^2*T(n,k-1)*T(n,k));
Column k: zeta(2) = 1 + 1/2^2 + 1/3^2 + ... + 1/k^2 + 2*Sum_{n >= 1} (-1)^(n+1)/(n^2*T(n-1,k)*T(n,k));
Main diagonal: zeta(2) = 5 * Sum_{n >= 1} (-1)^(n+1)/(n^2*T(n-1,n-1)*T(n,n)).
Conjectural result for superdiagonals: zeta(2) = 1 + 1/2^2 + ... + 1/k^2 + Sum_{n >= 1} (-1)^(n+1) * (5*n^2 + 6*k*n + 2*k^2)/(n^2*(n+k)^2*T(n-1,n+k-1)*T(n,n+k)), k = 0,1,2... .
Conjectural result for subdiagonals: zeta(2) = 2*(1 - 1/2^2 + ... + (-1)^(k+1)/k^2) + (-1)^k*Sum_{n >= 1} (-1)^(n+1)*(5*n^2 + 4*k*n + k^2)/(n^2*(n+k)^2*T(n+k-1,n-1)*T(n+k,n)), k = 0,1,2... .
Conjectural congruences: the main superdiagonal numbers S(n) := T(n,n+1) appear to satisfy the supercongruences S(m*p^r - 1) = S(m*p^(r-1) - 1) (mod p^(3*r)) for all primes p >= 5 and all positive integers m and r. If p is prime of the form 4*n + 1 we can write p = a^2 + b^2 with a an odd number. Then calculation suggests the congruence S((p-1)/2) == 2*a^2 (mod p). (End)
From Michael Somos, Jun 03 2012: (Start)
T(n, k) = hypergeom([-n, -k, n + 1], [1, 1], 1).
T(n, n-1) = A208675(n).
T(n+1, n) = A108628(n). (End)
T(n, k) = binomial(n, k)*hypergeom([-k, k - n, k - n], [1, -n], 1). - Peter Luschny, Feb 10 2018
From Peter Bala, Jun 23 2023: (Start)
T(n, k) = Sum_{i = 0..k} (-1)^i * binomial(n, i)*binomial(n+k-i, k-i)^2.
T(n, k) = binomial(n+k, k)^2 * hypergeom([-n, -k, -k], [-n - k, -n - k], 1). (End)
From Peter Bala, Jun 28 2023; (Start)
T(n,k) = the coefficient of (x^n)*(y^k)*(z^n) in the expansion of 1/( (1 - x - y)*(1 - z ) - x*y*z ).
T(n,k) = B(n, k, n) in the notation of Straub, equation 24.
The supercongruences T(n*p^r, k*p^r) == T(n*p^(r-1), k*p^(r-1)) (mod p^(3*r)) hold for all primes p >= 5 and positive integers n and k.
The formula T(n,k) = hypergeom([n+1, -n, -k], [1, 1], 1) allows the table indexing to be extended to negative values of n and k; clearly, we find that T(-n,k) = T(n-1,k) for all n and k. It appears that T(n,-k) = (-1)^n*T(n,k-1) for n >= 0, while T(n,-k) = (-1)^(n+1)*T(n,k-1) for n <= -1 [added Sep 10 2023: these follow from the identities immediately below]. (End)
T(n,k) = Sum_{i = 0..n} (-1)^(n+i) * binomial(n, i)*binomial(n+i, i)*binomial(k+i, i) = (-1)^n * hypergeom([n + 1, -n, k + 1], [1, 1], 1). - Peter Bala, Sep 10 2023
From G. C. Greubel, Oct 05 2023: (Start)
Let t(n,k) = T(n-k, k) (antidiagonals).
t(n, k) = Hypergeometric3F2([k-n, -k, n-k+1], [1,1], 1).
T(n, 2*n) = A363867(n).
T(3*n, n) = A363868(n).
T(2*n, 2*n) = A363869(n).
T(n, 3*n) = A363870(n).
T(2*n, 3*n) = A363871(n). (End)
T(n, k) = Sum_{i = 0..n} binomial(n, i)*binomial(n+i, i)*binomial(k, i). - Peter Bala, Feb 26 2024
Sum_{k = 0..n} (-1)^(n+k)*binomial(n, k)*binomial(n+k, k)*T(n, k) = A005259(n), the Apéry numbers associated with zeta(3). - Peter Bala, Jul 18 2024
From Peter Bala, Sep 21 2024: (Start)
Sum_{k = 0..n} (-1)^(n+k)*binomial(n, k)*T(n, k) = binomial(2*n, n) = A000984(n).
Sum_{k = 0..n} (-1)^(n+k)*binomial(n, k)*binomial(n+k, k)*T(n-1, n-k) = A376458(n).
Sum_{k = 0..n} (-1)^(n+k)*binomial(n, k)*binomial(n+k, k)*T(i, k) = A143007(n, i). (End)
From Peter Bala, Oct 12 2024: (Start)
The square array = A063007 * transpose(A007318).
Conjecture: for positive integer m, Sum_{k = 0..n} (-1)^(n+k) * binomial(n, k) * T(m*n, k) = ((m+1)*n)!/( ((m-1)*n)!*n!^2) (verified up to m = 10 using the MulZeil procedure in Doron Zeilberger's MultiZeilberger package). (End)

A047969 Square array of nexus numbers a(n,k) = (n+1)^(k+1) - n^(k+1) (n >= 0, k >= 0) read by upwards antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 7, 1, 1, 7, 19, 15, 1, 1, 9, 37, 65, 31, 1, 1, 11, 61, 175, 211, 63, 1, 1, 13, 91, 369, 781, 665, 127, 1, 1, 15, 127, 671, 2101, 3367, 2059, 255, 1, 1, 17, 169, 1105, 4651, 11529, 14197, 6305, 511, 1, 1, 19, 217, 1695, 9031
Offset: 0

Views

Author

N. J. A. Sloane

Keywords

Comments

If each row started with an initial 0 (i.e., a(n,k) = (n+1)^k - n^k) then each row would be the binomial transform of the preceding row. - Henry Bottomley, May 31 2001
a(n-1, k-1) is the number of ordered k-tuples of positive integers such that the largest of these integers is n. - Alford Arnold, Sep 07 2005
From Alford Arnold, Jul 21 2006: (Start)
The sequences in A047969 can also be calculated using the Eulerian Array (A008292) and Pascal's Triangle (A007318) as illustrated below: (cf. A101095).
1 1 1 1 1 1
1 1 1 1 1 1
-----------------------------------------
1 2 3 4 5 6
1 2 3 4 5
1 3 5 7 9 11
-----------------------------------------
1 3 6 10 15 21
4 12 24 40 60
1 3 6 10
1 7 19 37 61 91
-----------------------------------------
1 4 10 20 35 56
11 44 110 220 385
11 44 110 220
1 4 10
1 15 65 175 369 671
----------------------------------------- (End)
From Peter Bala, Oct 26 2008: (Start)
The above remarks of Alford Arnold may be summarized by saying that (the transpose of) this array is the Hilbert transform of the triangle of Eulerian numbers A008292 (see A145905 for the definition of the Hilbert transform). In this context, A008292 is best viewed as the array of h-vectors of permutohedra of type A. See A108553 for the Hilbert transform of the array of h-vectors of type D permutohedra. Compare this array with A009998.
The polynomials n^k - (n-1)^k, k = 1,2,3,..., which give the nonzero entries in the columns of this array, satisfy a Riemann hypothesis: their zeros lie on the vertical line Re s = 1/2 in the complex plane. See A019538 for the connection between the polynomials n^k - (n-1)^k and the Stirling polynomials of the simplicial complexes dual to the type A permutohedra.
(End)
Empirical: (n+1)^(k+1) - n^(k+1) is the number of first differences of length k+1 arrays of numbers in 0..n, k > 0. - R. H. Hardin, Jun 30 2013
a(n-1, k-1) is the number of bargraphs of width k and height n. Examples: a(1,2) = 7 because we have [1,1,2], [1,2,1], [2,1,1], [1,2,2], [2,1,2], [2,2,1], and [2,2,2]; a(2,1) = 5 because we have [1,3], [2,3], [3,1], [3,2], and [3,3] (bargraphs are given as compositions). This comment is equivalent to A. Arnold's Sep 2005 comment. - Emeric Deutsch, Jan 30 2017

Examples

			Array a begins:
  [n\k][0  1   2    3    4   5  6  ...
  [0]   1  1   1    1    1   1  1  ...
  [1]   1  3   7   15   31  63  ...
  [2]   1  5  19   65  211  ...
  [3]   1  7  37  175  ...
  ...
Triangle T begins:
  n\m   0   1    2     3     4      5      6      7      8     9  10 ...
  0:    1
  1:    1   1
  2:    1   3    1
  3:    1   5    7     1
  4:    1   7   19    15     1
  5:    1   9   37    65    31      1
  6:    1  11   61   175   211     63      1
  7:    1  13   91   369   781    665    127      1
  8:    1  15  127   671  2101   3367   2059    255      1
  9:    1  17  169  1105  4651  11529  14197   6305    511     1
  10:   1  19  217  1695  9031  31031  61741  58975  19171  1023   1
  ...  - _Wolfdieter Lang_, May 07 2021

References

  • J. H. Conway and R. K. Guy, The Book of Numbers, Copernicus Press, NY, 1996, p. 54.

Links

Crossrefs

Cf. A047970.
Cf. A009998, A108553 (Hilbert transform of array of h-vectors of type D permutohedra), A145904, A145905.
Row n sequences of array a: A000012, A000225(k+1), A001047(k+1), A005061(k+1), A005060(k+1), A005062(k+1), A016169(k+1), A016177(k+1), A016185(k+1), A016189(k+1), A016195(k+1), A016197(k+1).
Column k sequences of array a: (nexus numbers): A000012, A005408, A003215, A005917(n+1), A022521, A022522, A022523, A022524, A022525, A022526, A022527, A022528.
Cf. A343237 (row reversed triangle).

Programs

  • Mathematica
    Flatten[Table[n = d - e; k = e; (n + 1)^(k + 1) - n^(k + 1), {d, 0, 100}, {e, 0, d}]] (* T. D. Noe, Feb 22 2012 *)
  • Maxima
    T(n,m):=if m=0 then 1 else sum(k!*(-1)^(m+k)*stirling2(m,k)*binomial(n+k-1,n),k,0,m); /* Vladimir Kruchinin, Jan 28 2018 */

Formula

From Vladimir Kruchinin: (Start)
O.g.f. of e.g.f of rows of array: ((1-x)*exp(y))/(1-x*exp(y))^2.
T(n,m) = Sum_{k=0..m} k!*(-1)^(m+k)*Stirling2(m,k)*C(n+k-1,n), T(n,0)=1.(End)
From Wolfdieter Lang, May 07 2021: (Start)
T(n,m) = a(n-m,m) = (n-m+1)^(m+1) - (n-m)^(m+1), n >= 0, m = 0, 1,..., n.
O.g.f. column k of the array: polylog(-(k+1), x)*(1-x)/x. See the Peter Bala comment above, and the Eulerian triangle A008292 formula by Vladeta Jovovic, Sep 02 2002.
E.g.f. of e.g.f. of row of the array: exp(y)*(1 + x*(exp(y) - 1))*exp(x*exp(y)).
O.g.f. of triangle's exponential row polynomials R(n, y) = Sum_{m=0} T(n, m)*(y^m)/m!: G(x, y) = exp(x*y)*(1 - x)/(1 - x*exp(x*y))^2. (End)

A271942 Triangle read by rows: T(n,k) is the number of bargraphs of semiperimeter n having width k (n>=2, k>=1).

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 6, 1, 1, 7, 16, 10, 1, 1, 9, 31, 40, 15, 1, 1, 11, 51, 105, 85, 21, 1, 1, 13, 76, 219, 295, 161, 28, 1, 1, 15, 106, 396, 771, 721, 280, 36, 1, 1, 17, 141, 650, 1681, 2331, 1582, 456, 45, 1, 1, 19, 181, 995, 3235, 6083, 6244, 3186, 705, 55, 1, 1, 21, 226, 1445, 5685, 13663, 19348, 15156, 5985, 1045, 66, 1, 1, 23, 276, 2014, 9325, 27483, 50464, 55308, 33903, 10615, 1496, 78, 1
Offset: 2

Views

Author

Emeric Deutsch, May 21 2016

Keywords

Comments

Sum of entries in row n = A082582(n).
Sum(k*T(n,k), k>=1) = A271943(n).
Connection with A145904 should be explored.

Examples

			Row 4 is 1,3,1 because the 5 (=A082582(4)) bargraphs of semiperimeter 4 correspond to the compositions [1,1,1], [1,2], [2,1], [2,2], [3] which, clearly, have widths 3,2,2,2,1.
Triangle starts:
                               1
                              1, 1
                            1, 3, 1
                           1, 5, 6, 1
                        1, 7, 16, 10, 1
                      1, 9, 31, 40, 15, 1
                   1, 11, 51, 105, 85, 21, 1
                1, 13, 76, 219, 295, 161, 28, 1
             1, 15, 106, 396, 771, 721, 280, 36, 1
         1, 17, 141, 650, 1681, 2331, 1582, 456, 45, 1

Links

Crossrefs

Cf. A082582, A145904, A271943.

Programs

  • Maple
    eq := x*z*G^2-(1-x*z-z-x*z^2)*G+x*z^2 = 0: G := RootOf(eq, G): Gser := simplify(series(G, z = 0, 23)): for n from 2 to 20 do P[n] := sort(expand(coeff(Gser, z, n))) end do: for n from 2 to 20 do seq(coeff(P[n], x, j), j = 1 .. degree(P[n])) end do; # yields sequence in triangular form
# second Maple program:
b:= proc(n, y, t) option remember; expand(`if`(n=0, (1-t),
      `if`(t<0, 0, b(n-1, y+1, 1))+`if`(t>0 or y<2, 0,
       b(n, y-1, -1))+`if`(y<1, 0, b(n-1, y, 0)*z)))
    end:
T:= n-> (p-> seq(coeff(p, z, i), i=1..n-1))(b(n, 0$2)):
seq(T(n), n=2..20);  # Alois P. Heinz, Jun 06 2016
# Alternative, (assuming offset (0,0)):
T := (n, k) -> simplify(hypergeom([-k, k + 3, k - n], [1, 2], 1)):
seq(seq(T(n, k), k=0..n), n=0..9); # Peter Luschny, Oct 18 2020
  • Mathematica
    b[n_, y_, t_] := b[n, y, t] = Expand[If[n == 0, {1 - t}, If[t < 0, 0, b[n - 1, y + 1, 1]] + If[t > 0 || y < 2, 0, b[n, y - 1, -1]] + If[y < 1, 0, b[n - 1, y, 0]*z]]];
T[n_] := Function[p, Table[Coefficient[p, z, i], {i, 1, n-1}]][b[n, 0, 0] ];
Table[T[n], {n, 2, 20}] // Flatten (* Jean-François Alcover, Jul 21 2016, after Alois P. Heinz *)
T[n_, k_] := Sum[(Binomial[n - k - 1, j]*Binomial[n - k, j]*Binomial[2*n - k - 2*j, 2*n - 2*k])/(j + 1), {j, 0, n - k}]; Flatten[Table[T[n, k], {n, 0, 12}, {k, 0, n}]] (* , provided one bases the offset in (0, 0). Detlef Meya, Jan 07 2023 *)

Formula

G.f.: G(x,z) satisfies xzG^2-(1-xz-z-xz^2)G+xz^2=0 (z marks semiperimeter, x marks width).
T(n, k) = hypergeom([-k, k + 3, k - n], [1, 2], 1), provided one bases the offset in (0, 0). - Peter Luschny, Oct 18 2020
T(n, k) = Sum_{j=0..n - k} (binomial(n - k - 1, j)*binomial(n - k, j)*binomial(2*n - k - 2*j, 2*n - 2*k))/(j + 1), provided one bases the offset in (0, 0). - Detlef Meya, Jan 07 2023

A273350 Triangle read by rows: T(n,k) is the number of bargraphs of semiperimeter n having k up steps (n >= 2, k >= 1).

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 6, 5, 1, 1, 10, 16, 7, 1, 1, 15, 40, 31, 9, 1, 1, 21, 85, 105, 51, 11, 1, 1, 28, 161, 295, 219, 76, 13, 1, 1, 36, 280, 721, 771, 396, 106, 15, 1, 1, 45, 456, 1582, 2331, 1681, 650, 141, 17, 1, 1, 55, 705, 3186, 6244, 6083, 3235, 995, 181, 19, 1
Offset: 2

Views

Author

Emeric Deutsch, Jun 02 2016

Keywords

Comments

Number of entries in row n is n-1.
Sum of entries in row n = A082582(n).
Sum_{k>=1} k*T(n,k) = A273351(n).
Connection with A145904 should be explored.
This seems to be a duplicate of A145904. - Alois P. Heinz, Jun 04 2016. [This could probably be proved by showing that the g.f.s are the same. - N. J. A. Sloane, Jul 02 2016]

Examples

			Row 4 is 1,3,1 because the 5 (=A082582(4)) bargraphs of semiperimeter 4 correspond to the compositions [1,1,1], [1,2], [2,1], [2,2], [3] which, clearly, have 1,2,2,2,3 up steps.
Triangle starts
  1;
  1,  1;
  1,  3,  1;
  1,  6,  5,  1;
  1, 10, 16,  7,  1;
  1, 15, 40, 31,  9,  1.

Links

Crossrefs

Cf. A082582, A145904, A273351.

Programs

  • Maple
    eq := z*G^2-(1-z-y*z-y*z^2)*G+y*z^2 = 0: G := RootOf(eq, G): Gser := simplify(series(G, z = 0, 18)): for n from 2 to 16 do P[n] := sort(expand(coeff(Gser, z, n))) end do: for n from 2 to 16 do seq(coeff(P[n], y, j), j = 1 .. n-1) end do; # yields sequence in triangular form
  • Mathematica
    G = (-1 + z + y z + y z^2 + Sqrt[-4y z^3 + (-1 + z + y z + y z^2)^2])/(2z);
Rest @ CoefficientList[#, y]& /@ Drop[-CoefficientList[G + O[z]^13, z], 2] // Flatten (* Jean-François Alcover, Oct 24 2019 *)
  • Maxima
    T(n,m):=if n<=m then 0 else 1/(n-m)*sum(binomial(n-m,i-1)*binomial(n-m,i)*binomial(n-i-1,m-i),i,0,m); /* Vladimir Kruchinin, Jan 15 2018 */

Formula

G.f.: G=G(t,z) satisfies zG^2-(1-z-yz-yz^2)G+yz^2=0 (z marks semiperimeter, y marks up-steps).
T(n,m) = (1/(n-m))*Sum_{i=0..m} C(n-m,i-1)*C(n-m,i)*C(n-i-1,m-i). - Vladimir Kruchinin, Jan 15 2018
Showing 1-5 of 5 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.