A066822 -id:A066822 - cp's OEIS frontend

A038622 Triangular array that counts rooted polyominoes.

1, 2, 1, 5, 3, 1, 13, 9, 4, 1, 35, 26, 14, 5, 1, 96, 75, 45, 20, 6, 1, 267, 216, 140, 71, 27, 7, 1, 750, 623, 427, 238, 105, 35, 8, 1, 2123, 1800, 1288, 770, 378, 148, 44, 9, 1, 6046, 5211, 3858, 2436, 1296, 570, 201, 54, 10, 1, 17303, 15115, 11505, 7590, 4302, 2067, 825, 265

Offset: 0

N. J. A. Sloane, torsten.sillke(AT)lhsystems.com

The PARI program gives any row k and any n-th term for this triangular array in square or right triangle array format. - Randall L Rathbun, Jan 20 2002
Triangle T(n,k), 0 <= k <= n, read by rows given by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = 2*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + T(n-1,k) + T(n-1,k+1) for k >= 1. - Philippe Deléham, Mar 27 2007
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = x*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + y*T(n-1,k) + T(n-1,k+1) for k >= 1. Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
Triangle read by rows = partial sums of A064189 terms starting from the right. - Gary W. Adamson, Oct 25 2008
Column k has e.g.f. exp(x)*(Bessel_I(k,2x)+Bessel_I(k+1,2x)). - Paul Barry, Mar 08 2011

			From _Paul Barry_, Mar 08 2011: (Start)
Triangle begins
     1;
     2,    1;
     5,    3,    1;
    13,    9,    4,   1;
    35,   26,   14,   5,   1;
    96,   75,   45,  20,   6,   1;
   267,  216,  140,  71,  27,   7,  1;
   750,  623,  427, 238, 105,  35,  8, 1;
  2123, 1800, 1288, 770, 378, 148, 44, 9, 1;
Production matrix is
  2, 1,
  1, 1, 1,
  0, 1, 1, 1,
  0, 0, 1, 1, 1,
  0, 0, 0, 1, 1, 1,
  0, 0, 0, 0, 1, 1, 1,
  0, 0, 0, 0, 0, 1, 1, 1,
  0, 0, 0, 0, 0, 0, 1, 1, 1,
  0, 0, 0, 0, 0, 0, 0, 1, 1, 1
(End)

Reinhard Zumkeller, Rows n = 0..120 of triangle, flattened
D. Gouyou-Beauchamps and G. Viennot, Equivalence of the two-dimensional directed animal problem to a one-dimensional path problem, Adv. in Appl. Math. 9 (1988), no. 3, 334-357. See Table 1 on page 340.

Cf. A005773 (1st column), A005774 (2nd column), A005775, A066822, A000244 (row sums).

Cf. A064189, A007318, A061554, A092392.

import Data.List (transpose)
a038622 n k = a038622_tabl !! n !! k
a038622_row n = a038622_tabl !! n
a038622_tabl = iterate (\row -> map sum $
   transpose [tail row ++ [0,0], row ++ [0], [head row] ++ row]) [1]
-- Reinhard Zumkeller, Feb 26 2013

T := (n,k) -> simplify(GegenbauerC(n-k,-n+1,-1/2)+GegenbauerC(n-k-1,-n+1,-1/2)):
for n from 1 to 9 do seq(T(n,k),k=1..n) od; # Peter Luschny, May 12 2016

nmax = 10; t[n_ /; n > 0, k_ /; k >= 1] := t[n, k] = t[n-1, k-1] + t[n-1, k] + t[n-1, k+1]; t[0, 0] = 1; t[0, ] = 0; t[?Negative, ?Negative] = 0; t[n, 0] := 2 t[n-1, 0] + t[n-1, 1]; Flatten[ Table[ t[n, k], {n, 0, nmax}, {k, 0, n}]](* Jean-François Alcover, Nov 09 2011 *)

s=[0,1]; {A038622(n,k)=if(n==0,1,t=(2*(n+k)*(n+k-1)*s[2]+3*(n+k-1)*(n+k-2)*s[1])/((n+2*k-1)*n); s[1]=s[2]; s[2]=t; t)}

a(n, k) = a(n-1, k-1) + a(n-1, k) + a(n-1, k+1) for k>0, a(n, k) = 2*a(n-1, k) + a(n-1, k+1) for k=0.
Riordan array ((sqrt(1-2x-3x^2)+3x-1)/(2x(1-3x)),(1-x-sqrt(1-2x-3x^2))/(2x)). Inverse of Riordan array ((1-x)/(1+x+x^2),x/(1+x+x^2)). First column is A005773(n+1). Row sums are 3^n (A000244). If L=A038622, then L*L' is the Hankel matrix for A005773(n+1), where L' is the transpose of L. - Paul Barry, Sep 18 2006
T(n,k) = GegenbauerC(n-k,-n+1,-1/2) + GegenbauerC(n-k-1,-n+1,-1/2). In this form also the missing first column of the triangle 1,1,1,3,7,19,... (cf. A002426) can be computed. - Peter Luschny, May 12 2016
From Peter Bala, Jul 12 2021: (Start)
T(n,k) = Sum_{j = k..n} binomial(n,j)*binomial(j,floor((j-k)/2)).
Matrix product of Riordan arrays ( 1/(1 - x), x/(1 - x) ) * ( (1 - x*c(x^2))/(1 - 2*x), x*c(x^2) ) = A007318 * A061554 (triangle version), where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the g.f. of the Catalan numbers A000108.
Triangle equals A007318^(-1) * A092392 * A007318. (End)
The n-th row polynomial R(n,x) equals the n-th degree Taylor polynomial of the function (1 + x)*(1 + x + x^2)^n expanded about the point x = 0. - Peter Bala, Sep 06 2022

More terms from David W. Wilson

A005774 Number of directed animals of size n (k=1 column of A038622); number of (s(0), s(1), ..., s(n)) such that s(i) is a nonnegative integer and |s(i) - s(i-1)| <= 1 for i = 1,2,...,n, where s(0) = 2; also sum of row n+1 of array T in A026323.

Original entry on oeis.org

0, 1, 3, 9, 26, 75, 216, 623, 1800, 5211, 15115, 43923, 127854, 372749, 1088283, 3181545, 9312312, 27287091, 80038449, 234988827, 690513030, 2030695569, 5976418602, 17601021837, 51869858544, 152951628725, 451271872701, 1332147482253

Offset: 0

Simon Plouffe

Number of ordered trees with n+1 edges, having root degree at least 2 and nonroot outdegrees at most 2. - Emeric Deutsch, Aug 02 2002
From Petkovsek's algorithm, this recurrence does not have any closed form solutions. So there is no hypergeometric closed form for a(n). - Herbert S. Wilf
Sum of two consecutive trinomial coefficients starting two positions before central one. Example: a(4) = 10+16 and (1 + x + x^2)^4 = ... + 10*x^2 + 16*x^3 + 19*x^4 + ... - David Callan, Feb 07 2004
Image of n (A001477) under the Motzkin related matrix A107131. Binomial transform of A037952. - Paul Barry, May 12 2005
a(n) = total number of ascents (maximal runs of consecutive upsteps) in all Motzkin (n+1)-paths. For example, the 9 Motzkin 4-paths are FFFF, FFUD, FUDF, FUFD, UDFF, UDUD, UFDF, UFFD, UUDD and they contain a total of 9 ascents and so a(3)=9 (U=upstep, D=downstep, F=flatstep). - David Callan, Aug 16 2006
Image of the sequence (0,1,2,3,3,3,...) under the array A122896. - Paul Barry, Sep 18 2006
This is some kind of Motzkin transform of A079978 because the substitution x-> x*A001006(x) in the independent variable of the g.f. A079978(x) yields 1,0 followed by this sequence here. - R. J. Mathar, Nov 08 2008

			G.f.: x + 3*x^2 + 9*x^3 + 26*x^4 + 75*x^5 + 216*x^6 + 623*x^7 + ...

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

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
P. Barry, Riordan arrays, generalized Narayana triangles, and series reversion, Linear Algebra and its Applications, 491 (2016) 343-385.
D. Gouyou-Beauchamps, G. Viennot, Equivalence of the two-dimensional directed animal problem to a one-dimensional path problem, Adv. in Appl. Math. 9 (1988), no. 3, 334-357.
Christian Krattenthaler, Daniel Yaqubi, Some determinants of path generating functions, II, arXiv:1802.05990 [math.CO], 2018; Adv. Appl. Math. 101 (2018), 232-265.
Simon Plouffe, Approximations de Séries Génératrices et Quelques Conjectures, Dissertation, Université du Québec à Montréal, 1992.
Simon Plouffe, Une méthode pour obtenir la fonction génératrice d'une série, FPSAC 1993, Florence. Formal Power Series and Algebraic Combinatorics.

Cf. A038622, A098494, A026010, A005773, A005775, A066822.

a005774 0 = 0
a005774 n = a038622 n 1 -- Reinhard Zumkeller, Feb 26 2013

seq( add(binomial(i,k+1)*binomial(i-k,k), k=0..floor(i/2)), i=0..30 ); # Detlef Pauly (dettodet(AT)yahoo.de), Nov 09 2001
seq(simplify(GegenbauerC(n-2,-n,-1/2) + GegenbauerC(n-1,-n,-1/2)), n=0..27); # Peter Luschny, May 12 2016

CoefficientList[Series[(1-x-Sqrt[1-2x-3x^2])/(x(1-3x+Sqrt[1-2x-3x^2])),{x,0,30}],x] (* Harvey P. Dale, Sep 20 2011 *)
RecurrenceTable[{a[0]==0, a[1]==1,a[n]==(2n(n+1)a[n-1]+3n(n-1)a[n-2])/ ((n+2)(n-1))},a,{n,30}] (* Harvey P. Dale, Nov 09 2012 *)

s=[0,1]; {A005774(n)=k=(2*(n+2)*(n+1)*s[2]+3*(n+1)*n*s[1])/((n+3)*n); s[1]=s[2]; s[2]=k; k}

{a(n) = if( n<2, n>0, (2 * (n+1) * n *a(n-1) + 3 * (n-1) * n * a(n-2)) / (n+2) / (n-1))}; /* Michael Somos, May 01 2003 */

Inverse binomial transform of [ 0, 1, 5, 21, 84, ... ] (A002054). - John W. Layman
D-finite with recurrence (n+2)*(n-1)*a(n) = 2*n*(n+1)*a(n-1) + 3*n*(n-1)*a(n-2) for all n in Z. - Michael Somos, May 01 2003
E.g.f.: exp(x)*(BesselI(1, 2*x)+BesselI(2, 2*x)). - Vladeta Jovovic, Jan 01 2004
G.f.: (1-x-sqrt(1-2x-3x^2))/(x(1-3x+sqrt(1-2x-3x^2))); a(n)= Sum_{k=0..n} C(k+1, n-k+1)*C(n, k)*k/(k+1); a(n) = Sum_{k=0..n} C(n, k)*C(k, floor((k-1)/2)). - Paul Barry, May 12 2005
Starting (1, 3, 9, 26, ...) = binomial transform of A026010: (1, 2, 4, 7, 14, 25, 50, 91, ...). - Gary W. Adamson, Oct 22 2007
a(n)*(2+n) = (4+4*n)*a(n-1) - n*a(n-2) + (12-6*n)*a(n-3). - Simon Plouffe, Feb 09 2012
a(n) ~ 3^(n+1/2) / sqrt(Pi*n). - Vaclav Kotesovec, Mar 10 2014
0 = a(n)*(+36*a(n+1) + 18*a(n+2) - 96*a(n+3) + 30*a(n+4)) + a(n+1)*(-6*a(n+1) + 49*a(n+2) - 26*a(n+3) + 3*a(n+4)) + a(n+2)*(+15*a(n+3) - 8*a(n+4)) + a(n+3)*(a(n+4)) if n >= 0. - Michael Somos, Aug 06 2014
a(n) = GegenbauerC(n-2,-n,-1/2) + GegenbauerC(n-1,-n,-1/2). - Peter Luschny, May 12 2016

Further descriptions from Clark Kimberling

A005775 Number of compact-rooted directed animals of size n having 3 source points.

Original entry on oeis.org

1, 4, 14, 45, 140, 427, 1288, 3858, 11505, 34210, 101530, 300950, 891345, 2638650, 7809000, 23107488, 68375547, 202336092, 598817490, 1772479905, 5247421410, 15538054455, 46019183840, 136325212750, 403933918375, 1197131976846, 3548715207534, 10521965227669

Offset: 3

Simon Plouffe

Binomial transform of A037955. - Paul Barry, Dec 28 2006
Apparently, the number of Dyck paths of semilength n that contain at least one UUU but avoid UUU's starting above level 0. - David Scambler, Jul 02 2013
a(n) = number of paths in the half-plane x >= 0 from (0,0) to (n-1,2) or (n-1,-3), and consisting of steps U=(1,1), D=(1,-1) and H=(1,0). For example, for n=5, we have the 14 paths: HHUU, UUHH, UHHU, HUUH, HUHU, UHUH, UDUU, UUDU, UUUD, DUUU, DDDH, HDDD, DHDD, DDHD. - José Luis Ramírez Ramírez, Apr 19 2015

			G.f. = x^3 + 4*x^4 + 14*x^5 + 45*x^6 + 140*x^7 + 427*x^8 + 1288*x^9 + 3858*x^10 + ...

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

Reinhard Zumkeller, Table of n, a(n) for n = 3..1000
Jean-Luc Baril and José Luis Ramírez, Descent distribution on Catalan words avoiding ordered pairs of Relations, arXiv:2302.12741 [math.CO], 2023.
D. Gouyou-Beauchamps and G. Viennot, Equivalence of the two-dimensional directed animal problem to a one-dimensional path problem, Adv. in Appl. Math. 9 (1988), no. 3, 334-357.

Cf. A005773.
k=2 column of array in A038622.
Cf. A005774, A066822.

a005775 = flip a038622 2 . (subtract 1)  -- Reinhard Zumkeller, Feb 26 2013

seq(simplify(GegenbauerC(n-4,-n+1,-1/2) + GegenbauerC(n-3,-n+1,-1/2)),n=3..28); # Peter Luschny, May 12 2016

nmax = 28; t[n_ /; n > 0, k_ /; k >= 1] := t[n, k] = t[n-1, k-1] + t[n-1, k] + t[n-1, k+1]; t[0, 0] = 1; t[0, ] = 0; t[?Negative, ?Negative] = 0; t[n, 0] := 2*t[n-1, 0] + t[n-1, 1]; a[n_] := t[n-1, 2]; Table[a[n], {n, 3, nmax} ] (* Jean-François Alcover, Jul 03 2013, from A038622 *)

{a(n) = polcoeff( (x^2 + x - 1 + (x^2 - 3*x + 1) * sqrt((1 + x) / (1 - 3*x) + x^3 * O(x^n))) / (2*x^2), n)};

{a(n) = n--; sum(k=0, n, binomial(n, k) * binomial(k, k\2 -1))}; /* Michael Somos, May 12 2016 */

D-finite with recurrence (n+2)*(n-3)*a(n) = 2*n*(n-1)*a(n-1) + 3*(n-1)*(n-2)*a(n-2), a(2)=0, a(3)=1. - Michael Somos, Feb 02 2002
G.f.: (x^2 + x - 1 +(x^2 - 3*x + 1)*sqrt((1+x)/(1-3*x)))/(2*x^2).
From Paul Barry, Dec 28 2006: (Start)
E.g.f.: exp(x)*(Bessel_I(2,2*x) + Bessel_I(3,2*x));
a(n+1) = Sum_{k=0..n} C(n,k)*C(k,floor(k/2)-1). (End)
a(n) ~ 3^(n-1/2) / sqrt(Pi*n). - Vaclav Kotesovec, Feb 25 2014
G.f.: (z^3*M(z)^2+z^4*M(z)^3)/(1-z-2*z^2*M(z)), where M(z) is the g.f. of Motzkin paths. - José Luis Ramírez Ramírez, Apr 19 2015
a(n) = GegenbauerC(n-4,-n+1,-1/2) + GegenbauerC(n-3,-n+1,-1/2). - Peter Luschny, May 12 2016
0 = a(n)*(+9*a(n+1) - 63*a(n+2) - 54*a(n+3) + 87*a(n+4) - 21*a(n+5))+ a(n+1)*(+21*a(n+1) + 79*a(n+2) + 13*a(n+3) - 118*a(n+4) + 35*a(n+5)) + a(n+2)*(-14*a(n+2) + 79*a(n+3) - 67*a(n+4) + 14*a(n+5)) + a(n+3)*(+6*a(n+3) + 19*a(n+4) - 11*a(n+5)) + a(n+4)*(+a(n+4) + a(n+5)) if n >= 0. - Michael Somos, May 12 2016
a(n) = A005773(n) - A001006(n) for n >= 3. - John Keith, Nov 20 2020

More terms from Randall L Rathbun, Jan 19 2002

Edited by Michael Somos, Feb 02 2002

cp's OEIS Frontend

A038622 Triangular array that counts rooted polyominoes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

Extensions

A005774 Number of directed animals of size n (k=1 column of A038622); number of (s(0), s(1), ..., s(n)) such that s(i) is a nonnegative integer and |s(i) - s(i-1)| <= 1 for i = 1,2,...,n, where s(0) = 2; also sum of row n+1 of array T in A026323.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Formula

Extensions

A005775 Number of compact-rooted directed animals of size n having 3 source points.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Formula

Extensions