IJCOM International Journal of Combinatorics 1687-9171 1687-9163 Hindawi Publishing Corporation 301394 10.1155/2014/301394 301394 Research Article Embedding Structures Associated with Riordan Arrays and Moment Matrices http://orcid.org/0000-0003-1909-8725 Barry Paul Mansour Toufik School of Science, Waterford Institute of Technology, Waterford Ireland wit.ie 2014 1732014 2014 17 12 2013 14 02 2014 17 3 2014 2014 Copyright © 2014 Paul Barry. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Every ordinary Riordan array contains two naturally embedded Riordan arrays. We explore this phenomenon, and we compare it to the situation for certain moment matrices of families of orthogonal polynomials.

1. Introduction

Riordan arrays  have been used mainly to prove combinatorial identities [2, 3]. Recently, their links to orthogonal polynomials have been investigated [4, 5], while there is a growing literature surrounding their structural properties . In this paper we investigate an embedding structure, common to all ordinary Riordan arrays. We also look at this embedding structure in the context of moment matrices of families of orthogonal polynomials. In addition to some knowledge of Riordan arrays, we assume that the reader has a basic familiarity with the theory of orthogonal polynomials on the real line , production matrices [14, 15], and continued fractions . We will meet a number of integer sequences and integer triangles in this paper. The On-Line Encyclopedia may be consulted for many of them [17, 18]. In this paper we will understand by an ordinary Riordan array an integer number triangle whose (n,k)-th element Tn,k is defined by a pair of power series g(x) and f(x) over the integers with g(x)=1+g1x+g2x2+, f(x)=x+f2x2+f3x3+, in the following manner: (1)Tn,k=[xn]g(x)f(x)k. The group law for Riordan arrays is given by (2)(g,f)·(h,l)=(g(hf),lf). The identity for this law is I=(1,x) and the inverse of (g,f) is (g,f)-1=(1/(gf¯),f¯), where f¯ is the compositional inverse of f. For a power series f(x)=n=0anxn with f(0)=0, we define the reversion or compositional inverse of f to be the power series f¯(x) such that f(f¯(x))=x. We sometimes write this as f¯=Revf.

If a matrix A is the inverse of the coefficient array of a family of orthogonal polynomials, then we will call it a moment matrix, and we will single out the first column as the moment sequence.

2. The Canonical Embedding

Let (g,f) be an ordinary Riordan array R, with general term (3)Tn,k=[xn]gfk.

Then we observe that there are two naturally associated Riordan arrays “embedded” in the array R as follows.

Beginning at the first column of R, we take every second column, “raising” the columns appropriately to obtain a lower-triangular matrix A. The matrix A is then the Riordan array (4)A=(g,f2x) with general term An,k given by (5)An,k=[xn]g(f2x)k=[xn]gx-kf2k=[xn+k]gf2k=Tn+k,2k. Similarly, starting at the second column of R, taking every second column and “raising” all columns appropriately to obtain a lower-triangular matrix, we obtain a matrix B. This matrix B is then a Riordan array, given by (6)B=(gfx,f2x). We have (7)B=(fx,x)·A. The general term Bn,k of B is given by (8)Bn,k=[xn]gfx(f2x)k=[xn]gx-k-1f2k+1=[xn+k+1]gf2k+1=Tn+k+1,2k+1.

Example 1.

We take the example of the binomial matrix (9)R=(11-x,x1-x). We then have (10)A=(11-x,x(1-x)2)withgeneralterm    (n+k2k),B=(1(1-x)2,x(1-x)2)withgeneralterm    (n+k+12k+1). The following decomposition makes this clear: (11)(10000011000012100013310014641015101051). The matrices A and B are the coefficient arrays of the Morgan-Voyce polynomials bn(x) and Bn(x), respectively.

Example 2.

We take the Riordan array (12)R=(c(x),xc(x)), where (13)c(x)=1-1-4x2x. Then we find that (14)A=(c(x),xc(x)2),B=(c(x)2,xc(x)2). The matrix R begins with (15)(100000110000221000553100141494104242281451). We note that the matrix A=(c(x),xc(x)2) is the moment array for the family of orthogonal polynomials with coefficient array given by (16)A-1=(c(x),xc(x)2)-1=(11+x,x(1+x)2). Denoting this family by Pn(x), we have (17)Pn(x)=(x-2)Pn-1(x)-Pn-2(x), with P0(x)=1 and P1(x)=x-1. Similarly the matrix B=(c(x)2,xc(x)2) is the moment array for the family of orthogonal polynomials with coefficient array given by (18)B-1=(1(1+x)2,x(1+x)2). Denoting this family by Qn(x), we have (19)Qn(x)  =(x-2)Qn-1(x)-Qn-2(x), with Q0(x)=1 and Q1(x)=x-2.

The inverse matrix R-1 is given by (20)(100000-1100000-2100001-3100003-41000-16-51), which is the Riordan array (1-x,x(1-x)). In it we see the elements of A-1 and B-1 in staggered fashion.

3. A Counter Example

It is natural to ask the question: is a matrix that contains two embedded Riordan arrays as above itself a Riordan array? The following example shows that this is not a sufficient condition on an array to be Riordan.

Example 3.

We will construct an invertible integer lower-triangular matrix which has two embedded Riordan arrays in the fashion above, but which is not itself a Riordan array. We start with the essentially two-period sequence (an)n0(21)1,2,3,2,3,2,3,. We form the matrix (22)(100000-2100000-3100000-2100000-3100000-21). The inverse of this matrix begins with (23)(100000210000631000126210036186310723612621), where we note an alternating pattern of constant columns (with generating functions (1+2x)/(1-6x2) and (1+3x)/(1-6x2), resp.). Removing the first row of this matrix provides us with a production matrix, which is not of the form that produces a Riordan array (after the first column, subsequent columns would be shifted versions of the second column [14, 15]). Thus the resulting matrix will not be a Riordan array. This resulting matrix begins with (24)(1000002100001051000623171004302155110103194159738987121). We now observe that for this matrix, we have (25)A=(11+2x,x1+5x+6x2)-1,B=(11+5x+6x2,x1+5x+6x2)-1. We notice that the sequence (26)1,2,10,62,430,3194, has generating function given by the continued fraction (27)11-2x1-3x1-2x1-, and secondly that (28)1+5x+6x2=1+(2+3)x+2.3x2=(1+2x)(1+3x). This construction is easily generalized.

4. Embedding a Riordan Array

Another natural question to ask is: if we are given a Riordan array A, is it possible to embed it as above into a Riordan array R? For this, we let (29)A=(u,v), and seek to determine (30)R=(g,f) such that A embeds into R. For this, we need (31)u=g,v=f2x. Thus we require that (32)f=xv=xvx. Since we are working in the context of integer valued Riordan arrays, we require that v/x generates an integer sequence. We can state our result as follows.

Proposition 4.

The Riordan array (33)A=(u,v) can be embedded in the Riordan array (34)R=(g,xvx) on condition that v/x is the generating function of an integer sequence.

Example 5.

The Riordan array (35)A=(11-4x,x1-4x) can be embedded in the Riordan array (36)R=(11-4x,x1-4x). For this example, the matrix A begins with (37)(10000021000066100020301010070140701410252630420126181), while R begins with (38)(1000002100006410002016610070643081025225614048101).

5. A Cascading Decomposition

We note that we can “cascade" this embedding process, in the sense that, given a Riordan array R, with embedded Riordan arrays A and B, we can consider decomposing A and B in their turns and then continue this process. For instance, we can decompose (39)A=(g,f2x) into the two matrices (40)AA=(g,f4x3),BA=(gf2x2,f4x3). In their turn AA and BA can be decomposed and so on.

6. Embedding and Orthogonal Polynomials

The phenomenon of embedding as described above is not confined to Riordan arrays, as the continued fraction example above shows. To further amplify this point, we give another example involving a continued fraction. Although we take a particular case, the general case can be inferred easily from it. Thus we take the particular case of the continued fraction (41)11-2x1-3x1-5x1-2x1-3x1-5x1-. This continued fraction is equal to (42)11-2x-6x21-8x-10x21-5x-15x21-7x-6x21-8x-10x21-5x-. By the theory of orthogonal polynomials, the power series expressed by both continued fractions is the generating function for the moment sequence of the family of orthogonal polynomials whose moment matrix (the inverse of the coefficient array of the orthogonal polynomials) has production matrix given by (43)(210000681000010510000157100006810000105). This production matrix generates the moment matrix A that begins with (44)(100000210000101010008010015100760103019022107700109002310350301). In order to produce an embedding for this matrix, we proceed as follows. We form the matrix (45)(100000-2100000-3100000-5100000-2100000-31). We invert this matrix, remove the first row of the resulting matrix, and use this new matrix as a production matrix. The generated matrix R then begins with (46)(10000021000010510008040101007603801001210770038501030130151). The moment matrix A is evidently embedded in the matrix R. We can show that the corresponding matrix B is the moment matrix for the family of orthogonal polynomials whose moments have generating function given by (47)11-5x-15x21-7x-6x21-8x-10x21-5x-. The matrix B begins with (48)(100000510000401210003801302010038501410300251040400155204060440321), and it has production matrix (49)(510000157100006810000105100001571000068).

The matrix R-1 can now be characterized as the coefficient array of a family of polynomials Rn(x) defined as follows. We let Pn(x) be the family of orthogonal polynomials with coefficient array A-1, and we let Qn(x) be the family of orthogonal polynomials with coefficient B-1. Then we have (50)Rn(x)={Qn/2(x)xn/2,if  n  is  even;Pn/2(x)xn/2,  otherwise.

In the general case of a moment sequence generated by the continued fraction (51)11-αx1-βx1-γx1-αx1-, the matrix A will be generated by the production matrix (52)(α10000αββ+γ10000αγα+β10000βγα+γ10000αββ+γ10000αγα+β), while the matrix B is generated by the production matrix (53)(α+β10000βγα+γ10000αββ+γ10000αγα+β10000βγα+γ10000αββ+γ).

7. Embedding, Orthogonal Polynomials, and Riordan Arrays

In this section, we consider the case of two related families of orthogonal polynomials Pn(x) and Qn(x) defined by (54)Pn(x)=(x-7)Pn-1(x)-12Pn-2, with P0(x)=1, P1(x)=x-3, and (55)Qn(x)=(x-7)Qn-1(x)-12Qn-2, with Q0(x)=1, Q1(x)=x-7. We note that (56)(1+3x)(1+4x)=1+7x+12x2. The coefficient array of the polynomials Pn(x) is then given by the Riordan array (57)A-1=(11+3x,x1+7x+12x2), while that of Qn(x) is given by (58)B-1=(11+7x+12x2,x1+7x+12x2). The moment matrix A begins with (59)(100000310000211010001831031710017851108234241018651123493034414311), while the moment matrix B starts with (60)(100000710000611410005951712110062172044330281068047244854690538351). Letting (61)Rn(x)={Qn/2(x)  xn/2,  ifniseven;  P(n+1)/2(x)xn/2,otherwise, we find that the inverse R of the coefficient array of the family Rn(x) is given by (62)(10000031000021710001836110100178559510314101865162171108171171). Thus the two Riordan arrays A and B, which are the moment arrays of the two families of orthogonal polynomials Pn(x) and Qn(x), respectively, embed into the generalized moment array R for the family of polynomials Rn(x). Now the production matrix of R begins with (63)(3100001241000361231001444812410432144361231172857614448124), where the columns have generating functions (1+4x)/(1-12x2), (1+3)/(1-12x2), respectively. We now observe that this production matrix is obtained by removing the first row of the inverse of the matrix (64)(100000-3100000-4100000-3100000-4100000-31). We note finally that the sequence (65)1,3,21,183,1785,18651,204141, has generating function given by (66)g(x)=11-3x1-4x1-3x1-. We have in fact that (67)g(x)=1xRevx(1-4x)1-x. On the other hand, if we let (68)Pn(x)=(x-7)Pn-1(x)-12Pn-2(x), but this time take P0(x)=1 and P1(x)=x-4, then the matrix A becomes (69)A=(11+4x,x1+7x+12x2)-1. The matrix A then begins with (70)(100000410000281110002441171810023801279255251024868143933364442321), where the moment sequence (71)1,4,28,244,2380, has generating function (72)11-4x1-3x1-4x1-, or equivalently (73)11-4x-12x21-7x-12x21-7x-. In this case, the matrix R begins with (74)(10000041000028710002446111100238059511714102486862171279171181). This matrix is then associated with the matrix (75)(100000-4100000-3100000-4100000-3100000-41).

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Shapiro L. W. Getu S. Woan W.-J. Woodson L. C. The Riordan group Discrete Applied Mathematics 1991 34 1–3 229 239 2-s2.0-0000483447 10.1016/0166-218X(91)90088-E MR1137996 ZBL0754.05010 Corsani C. Merlini D. Sprugnoli R. Left-inversion of combinatorial sums Discrete Mathematics 1998 180 1–3 107 122 2-s2.0-0041856467 Sprugnoli R. Riordan arrays and combinatorial sums Discrete Mathematics 1994 132 1–3 267 290 2-s2.0-0001731038 10.1016/0012-365X(92)00570-H MR1297386 ZBL0814.05003 Barry P. Hennessy A. Meixner-type results for Riordan arrays and associated integer sequences Journal of Integer Sequences 2010 13 9, article 10.9.4 2-s2.0-78650364353 MR2746252 ZBL1201.42017 Barry P. Riordan arrays, orthogonal polynomials as moments, and Hankel transforms Journal of Integer Sequences 2011 14 2, article 11.2.2 2-s2.0-79952941711 MR2783956 ZBL1213.42086 Cheon G.-S. Kim H. Shapiro L. W. Riordan group involutions Linear Algebra and Its Applications 2008 428 4 941 952 2-s2.0-37249034114 10.1016/j.laa.2007.09.003 MR2382102 ZBL1131.05012 Davenport D. E. Shapiro L. W. Woodson L. C. The double Riordan group The Electronic Journal of Combinatorics 2012 18 2, Paper 33 MR2880683 ZBL1243.05009 He T.-X. Sprugnoli R. Sequence characterization of Riordan arrays Discrete Mathematics 2009 309 12 3962 3974 2-s2.0-67349097924 10.1016/j.disc.2008.11.021 MR2537389 ZBL1228.05014 Jin S.-T. A characterization of the Riordan Bell subgroup by C-sequences The Korean Journal of Mathematics 2009 17 2 147 154 Merlini D. Rogers D. G. Sprugnoli R. Verri M. C. On some alternative characterizations of Riordan arrays Canadian Journal of Mathematics 1997 49 2 301 320 2-s2.0-0031508649 10.4153/CJM-1997-015-x MR1447493 ZBL0886.05013 Chihara T. S. An Introduction to Orthogonal Polynomials 1978 New York, NY, USA Gordon and Breach MR0481884 Gautschi W. Orthogonal Polynomials: Computation and Approximation 2004 Oxford, UK Clarendon Press MR2061539 Szegö G. Orthogonal Polynomials 1975 4th Providence, RI, USA American Mathematical Society Deutsch E. Ferrari L. Rinaldi S. Production matrices Advances in Applied Mathematics 2005 34 1 101 122 2-s2.0-10344247127 10.1016/j.aam.2004.05.002 MR2102277 ZBL1064.05012 Deutsch E. Ferrari L. Rinaldi S. Production matrices and Riordan arrays Annals of Combinatorics 2009 13 1 65 85 2-s2.0-70449528778 10.1007/s00026-009-0013-1 MR2529720 ZBL1229.05015 Wall H. S. Analytic Theory of Continued Fractions 1967 New York, NY, USA AMS Chelsea Sloane N. J. A. The On-Line Encyclopedia of Integer Sequences 2011, http://oeis.org Sloane N. J. A. The on-line Encyclopedia of Integer Sequences Notices of the American Mathematical Society 2003 50 8 912 915 MR1992789 ZBL1044.11108