Riordan Array
   HOME

TheInfoList



OR:

A Riordan array is an
infinite Infinite may refer to: Mathematics *Infinite set, a set that is not a finite set *Infinity, an abstract concept describing something without any limit Music Performers *Infinite (group), a South Korean boy band *Infinite (rapper), Canadian ra ...
lower
triangular matrix In mathematics, a triangular matrix is a special kind of square matrix. A square matrix is called if all the entries ''above'' the main diagonal are zero. Similarly, a square matrix is called if all the entries ''below'' the main diagonal are z ...
, D, constructed from two
formal power series In mathematics, a formal series is an infinite sum that is considered independently from any notion of convergence, and can be manipulated with the usual algebraic operations on series (addition, subtraction, multiplication, division, partial su ...
, d(t) of order 0 and h(t) of order 1, such that d_ = ^n(t) h(t)^k. A Riordan array is an element of the Riordan group. It was defined by mathematician Louis W. Shapiro and named after John Riordan. The study of Riordan arrays is a field influenced by and contributing to other areas such as
combinatorics Combinatorics is an area of mathematics primarily concerned with counting, both as a means and as an end to obtaining results, and certain properties of finite structures. It is closely related to many other areas of mathematics and has many ...
,
group theory In abstract algebra, group theory studies the algebraic structures known as group (mathematics), groups. The concept of a group is central to abstract algebra: other well-known algebraic structures, such as ring (mathematics), rings, field ( ...
,
matrix theory In mathematics, a matrix (: matrices) is a rectangular array or table of numbers, symbols, or expressions, with elements or entries arranged in rows and columns, which is used to represent a mathematical object or property of such an object. ...
,
number theory Number theory is a branch of pure mathematics devoted primarily to the study of the integers and arithmetic functions. Number theorists study prime numbers as well as the properties of mathematical objects constructed from integers (for example ...
,
probability Probability is a branch of mathematics and statistics concerning events and numerical descriptions of how likely they are to occur. The probability of an event is a number between 0 and 1; the larger the probability, the more likely an e ...
,
sequences In mathematics, a sequence is an enumerated collection of objects in which repetitions are allowed and order matters. Like a set, it contains members (also called ''elements'', or ''terms''). The number of elements (possibly infinite) is call ...
and
series Series may refer to: People with the name * Caroline Series (born 1951), English mathematician, daughter of George Series * George Series (1920–1995), English physicist Arts, entertainment, and media Music * Series, the ordered sets used i ...
,
Lie groups In mathematics, a Lie group (pronounced ) is a group that is also a differentiable manifold, such that group multiplication and taking inverses are both differentiable. A manifold is a space that locally resembles Euclidean space, whereas ...
and
Lie algebras In mathematics, a Lie algebra (pronounced ) is a vector space \mathfrak g together with an operation called the Lie bracket, an alternating bilinear map \mathfrak g \times \mathfrak g \rightarrow \mathfrak g, that satisfies the Jacobi identi ...
,
orthogonal polynomials In mathematics, an orthogonal polynomial sequence is a family of polynomials such that any two different polynomials in the sequence are orthogonal In mathematics, orthogonality (mathematics), orthogonality is the generalization of the geom ...
,
graph theory In mathematics and computer science, graph theory is the study of ''graph (discrete mathematics), graphs'', which are mathematical structures used to model pairwise relations between objects. A graph in this context is made up of ''Vertex (graph ...
,
networks Network, networking and networked may refer to: Science and technology * Network theory, the study of graphs as a representation of relations between discrete objects * Network science, an academic field that studies complex networks Mathematics ...
,
unimodal In mathematics, unimodality means possessing a unique mode. More generally, unimodality means there is only a single highest value, somehow defined, of some mathematical object. Unimodal probability distribution In statistics, a unimodal p ...
sequences, combinatorial identities,
elliptic curves In mathematics, an elliptic curve is a Smoothness, smooth, Projective variety, projective, algebraic curve of Genus of an algebraic curve, genus one, on which there is a specified point . An elliptic curve is defined over a field (mathematics), ...
,
numerical approximation Numerical analysis is the study of algorithms that use numerical approximation (as opposed to symbolic manipulations) for the problems of mathematical analysis (as distinguished from discrete mathematics). It is the study of numerical methods th ...
,
asymptotic analysis In mathematical analysis, asymptotic analysis, also known as asymptotics, is a method of describing Limit (mathematics), limiting behavior. As an illustration, suppose that we are interested in the properties of a function as becomes very larg ...
, and
data analysis Data analysis is the process of inspecting, Data cleansing, cleansing, Data transformation, transforming, and Data modeling, modeling data with the goal of discovering useful information, informing conclusions, and supporting decision-making. Da ...
. Riordan arrays also unify tools such as
generating functions In mathematics, a generating function is a representation of an infinite sequence of numbers as the coefficients of a formal power series. Generating functions are often expressed in closed form (rather than as a series), by some expression invo ...
,
computer algebra systems A computer algebra system (CAS) or symbolic algebra system (SAS) is any mathematical software with the ability to manipulate mathematical expressions in a way similar to the traditional manual computations of mathematicians and scientists. The d ...
,
formal languages In logic, mathematics, computer science, and linguistics, a formal language is a set of string (computer science), strings whose symbols are taken from a set called "#Definition, alphabet". The alphabet of a formal language consists of symbol ...
, and path models. Books on the subject, such as ''The Riordan Array'' (Shapiro et al., 1991), have been published.


Formal definition

A formal power series a(x)=a_0+a_1x+a_2x^2+ \cdots =\sum_a_j x^j \in \mathbb x (where \mathbb x is the
ring (The) Ring(s) may refer to: * Ring (jewellery), a round band, usually made of metal, worn as ornamental jewelry * To make a sound with a bell, and the sound made by a bell Arts, entertainment, and media Film and TV * ''The Ring'' (franchise), a ...
of formal power series with complex coefficients) is said to have
order Order, ORDER or Orders may refer to: * A socio-political or established or existing order, e.g. World order, Ancien Regime, Pax Britannica * Categorization, the process in which ideas and objects are recognized, differentiated, and understood ...
r if a_0=\cdots=a_ =0\neq a_r. Write \mathcal_r for the set of formal power series of order r. A power series a(x) has a ''multiplicative inverse'' (i.e. 1/a(x) is a power series) if and only if it has order 0, i.e. if and only if it lies in \mathcal_0; it has a ''composition inverse'' that is there exists a power series \bar a such that \bar a(a(x))=x if and only if it has order 1, i.e. if and only if it lies in \mathcal_1. As mentioned previously, a Riordan array is usually defined via a pair of power series (d(t),h(t))\in \mathcal_0 \times \mathcal_1. The "array" part in its name stems from the fact that one associates to (d(t),h(t)) the array of complex numbers defined by d_:= ^nd(t)h(t)^k, n,k\in \mathbb (here " ^ncdots" means "coefficient of t^n in \cdots"). Thus column k of the array consists of the sequence of coefficients of the power series d(t)h(t)^k; in particular, column 0 determines and is determined by the power series d(t). Because d(t) is of order 0, it has a multiplicative inverse, and it follows that from the array's column 1 we can recover h(t) as h(t)=d(t)^ d(t)h(t). Since h(t) has order 1, h(t)^k is of order k and so is d(t)h(t)^k. It follows that the array d_ is lower triangular and exhibits a geometric progression (d_)_=(d_0 h_1^k)_ on its main diagonal. It also follows that the map sending a pair of power series (d(t),h(t))\in \mathcal_0 \times \mathcal_1 to its triangular array is
injective In mathematics, an injective function (also known as injection, or one-to-one function ) is a function that maps distinct elements of its domain to distinct elements of its codomain; that is, implies (equivalently by contraposition, impl ...
.


Example

An example of a Riordan array is given by the pair of power series \left(\frac, \frac\right)=\left(\sum_ x^j, \sum_ x^\right)\in \mathcal F_0\times \mathcal F_1. It is not difficult to show that this pair generates the infinite triangular array of binomial coefficients d_=\binom, also called the
Pascal matrix Pascal, Pascal's or PASCAL may refer to: People and fictional characters * Pascal (given name), including a list of people with the name * Pascal (surname), including a list of people and fictional characters with the name ** Blaise Pascal, Fren ...
: P=\left(\begin 1 & & & & & & \\ 1 &1 & & & & & \\ 1 &2 &1 & & & & \cdots \\ 1 &3 &3 &1 & & & \\ 1 &4 &6 &4 &1 & & \\ & & \vdots & & & & \ddots \end\right). Proof: If q(x)=\sum_ q_j x^j is a power series with associated coefficient sequence (q_0,q_1,q_2,\dotsc), then, by Cauchy multiplication of power series, q(x)\frac=\sum_ (0+q_0+q_1+ \cdots+q_)x^j. Thus, the latter series has the coefficient sequence (0,q_0,q_0+q_1,q_0+q_1+q_2,\dotsc), and hence ^nq(x)\frac= q_0+\cdots+q_. Fix any k\in \mathbb_. If q_n=\binom, so that (q_n)_ represents column k of the Pascal array, then \sum_^ q_j= \sum_^ \binom=\binom. This argument shows by induction on k that \frac\left(\frac\right)^k has column k of the Pascal array as coefficient sequence.


Properties

Below are some often-used facts about Riordan arrays. Note that the matrix multiplication rules applied to infinite lower triangular matrices lead to finite sums only and the product of two infinite lower triangular matrices is infinite lower triangular. The next two theorems were first stated and proved by Shapiro et al., which describes them as derived from results in papers by
Gian-Carlo Rota Gian-Carlo Rota (April 27, 1932 – April 18, 1999) was an Italian-American mathematician and philosopher. He spent most of his career at the Massachusetts Institute of Technology, where he worked in combinatorics, functional analysis, proba ...
and the book of Roman. Theorem: a. Let (a(x),b(x)) and (c(x),d(x)) be Riordan arrays, viewed as infinite lower triangular matrices. Then the product of these matrices is the array associated to the pair (a(x) c(b(x)), d(b(x)) ) of formal power series, which is itself a Riordan array. b. This fact justifies the definition of the multiplication '*' of Riordan arrays viewed as pairs of power series by Proof: Since a(x) and c(x) have order 0, it is clear that a(x) c(b(x)) has order 0. Similarly, b(x),d(x)\in \mathcal F_1 implies d(b(x))\in \mathcal F_1. Therefore, (a(x) c(b(x)), d(b(x)) ) is a Riordan array. Define a matrix M as the Riordan array (a(x),b(x)). By definition, its j-th column M_ is the sequence of coefficients of the power series a(x)b(x)^j. If we multiply this matrix from the right with the sequence (r_0,r_1,r_2, . . . )^T we get as a result a linear combination of columns of M which we can read as a linear combination of power series, namely \sum_ r_\nu M_ = \sum_ r_\nu a(x)b(x)^\nu = a(x) \sum_ r_\nu b(x)^\nu. Thus, viewing sequence (r_0,r_1,r_2, . . . )^T as codified by the power series r(x), we showed that(a(x),b(x))*r(x) = a(x) r(b(x)). Here the * is the symbol for indicating correspondence on the power series level with matrix multiplication. We multiplied a Riordan array (a(x),b(x)) with a single power series. Now let (c(x),d(x)) be another Riordan array viewed as a matrix. One can form the product (a(x),b(x))(c(x),d(x)). The j-th column of this product is just (a(x),b(x)) multiplied with the j-th column of (c(x),d(x)). Since the latter corresponds to the power series c(x) d(x)^j, it follows by the above that the j-th column of (a(x),b(x))(c(x),d(x)) corresponds to a(x) c(b(x)) d(b(x))^j. As this holds for all column indices j occurring in (c(x),d(x)) we have shown part a. Part b is now clear. \Box Theorem: The family of Riordan arrays endowed with the product '*' defined above forms a group: the Riordan group. Proof: The associativity of the multiplication '*' follows from associativity of matrix multiplication. Next note (1,x)*(c(x),d(x))=(1\cdot c(x), d(x))=(c(x),d(x)). So (1,x) is a left neutral element. Finally, we claim that (c(\bar d(x))^, \bar d(x)) is the left inverse to the power series (c(x),d(x)). For this check the computation (c(\bar d(x))^, \bar d(x))* (c(x),d(x)) =( (c(\bar d(x) )^ c(d(x)) , d(\bar d(x)) )= (1,x). As is well known, an associative structure which has a left neutral element and where each element has a left inverse is a group. \Box Of course, not all invertible infinite lower triangular arrays are Riordan arrays. Here is a useful characterization for the arrays that are Riordan. The following result is apparently due to Rogers. Theorem: An infinite lower triangular array D=(d_)_ is a Riordan array if and only if there exist a sequence traditionally called the A-sequence, A=(a_0\neq 0,a_1, . . . ) such that Proof. \Rightarrow: Let D be the Riordan array stemming from (d(t),h(t)). Since d(t)\in \mathcal F_0, d_\neq 0. Since h(t) has order 1, it follows that (d(t)h(t)/t, h(t)) is a Riordan array and by the group property there exists a Riordan array (A(t),B(t)) such that (d(t),h(t))*(A(t),B(t))=(d(t)h(t)/t, h(t)). Computing the left-hand side yields (d(t)A(h(t)),B(h(t)), and hence, comparison yields B(h(t))= h(t). Thus, B(t)=t is a solution to this equation; it is unique since B is composition invertible. Thus, we can rewrite the equation as (d(t),h(t))*(A(t), t)=(d(t)h(t)/t, h(t)). From the matrix multiplication law, the (n,k)-entry of the left-hand side of this latter equation is At the other hand the (n,k)-entry of the right-hand side of the equation above is so that i results. From *_1 we also get d_=a_0 d_ for all n\geq 0 and since we know that the diagonal elements are nonzero, we have a_0\neq 0. Note that using equation *_1 one can compute all entries knowing the entries (d_)_. \Leftarrow : Now assume that, for a triangular array, we have the equations *_1 for some sequence (a_j)_. Let A(t) be the generating function of that sequence and define h(t) from the equation tA(h(t))=h(t). Check that it is possible to solve the resulting equations for the coefficients of h; and since a_0\neq 0 one gets that h(t) has order 1. Let d(t) be the generating function of the sequence (d_,d_,d_,. . . ). Then for the pair (d(t), h(t)) we find (d(t),h(t))*(A(t),t)= (d(t)A(h(t)), h(t))= (d(t)h(t)/t, h(t)). This is the same equations we found in the first part of the proof, and going through its reasoning, we find equations as in *_1 . Since d(t) (or the sequence of its coefficients) determines the other entries, we see that the initial array is the array we deduced. Thus, the array in *_1 is a Riordan array. \Box Clearly, the A-sequence alone does not contain all the information about a Riordan array. Indeed, it only determines h(t) and places no restriction on d(t). To determine d(t) "horizontally", a similarly defined Z-sequence is used. Theorem. Let (d_)_ be an infinite lower triangular array whose diagonal sequence (d_)_ does not contain zeroes. Then there exists a unique sequence Z=(z_0,z_1,z_2, . . .) such that Proof: By the triangularity of the array, the equation claimed is equivalent to d_=\sum_^n z_j d_. For n=0, this equation is d_=z_0 d_ and, as d_\neq 0, it allows computing z_0 uniquely. In general, if z_0,z_1, . . ., z_ are known, then d_-\sum_^ z_j d_ = z_n d_ allows computing z_n uniquely. \Box


References

{{reflist Combinatorics