---
title: Generalized Riordan Arrays
url: https://www.emergentmind.com/topics/generalized-riordan-arrays
type: topic
---

# Generalized Riordan Arrays

Generalized Riordan arrays denote several extensions of the classical Riordan-array formalism in which lower-triangular matrices are encoded by formal series and act on generating functions by multiplication and substitution. In the classical setting, a Riordan array \((g(x),f(x))\) has entries \(t_{n,k}=[x^n]\,g(x)f(x)^k\), group law \((g,f)(u,v)=\bigl(g(x)u(f(x)),\,v(f(x))\bigr)\), inverse \(\left(\frac{1}{g(\bar f(x))},\bar f(x)\right)\), and bivariate generating function \(\frac{g(x)}{1-yf(x)}\) [2108.13537]. In current literature, the adjective “generalized” covers several non-equivalent enlargements of this template: weighted and \(q\)-dependent coefficient systems, generalized Pascal-matrix conjugations, higher-periodic groups using several substitution series, stretched and almost Riordan constructions, multivariate analogues, and semi-Laurent versions with bi-infinite matrices [1612.07657][1505.03100][2412.05461][2509.04160]. This suggests that the subject is best viewed as a family of compatible operator theories rather than a single universally fixed definition.

## 1. Classical core and operator-theoretic template

The common point of departure is the ordinary Riordan group. A Riordan array is specified by power series
\[
g(x)=g_0+g_1x+g_2x^2+\cdots,\qquad g_0\neq 0,
\]
and
\[
f(x)=f_1x+f_2x^2+\cdots,\qquad f_1\neq 0,
\]
with entries \(t_{n,k}=[x^n]\,g(x)f(x)^k\). Its action on a series \(h(x)\) is the Fundamental Theorem of Riordan arrays,
\[
(g(x),f(x))\cdot h(x)=g(x)\,h(f(x)),
\]
and the associated production matrix is lower Hessenberg and determined by the \(Z\)- and \(A\)-sequences. In the standard characterization recalled in the production-matrix literature, if \(Z(x)=\sum z_nx^n\) and \(A(x)=\sum a_nx^n\), then
\[
M=\left(\frac{1}{1-xZ(x)},\,\frac{x}{A(x)}\right).
\]
This classical skeleton remains the reference point for most generalized constructions [2011.13985].

A second foundational thread is the operator and umbral viewpoint. In the symbolic treatment of Riordan arrays, exponential Riordan arrays are represented by umbral pairs \((\gamma,\alpha)\), and the generalized \(\omega\)-Riordan array attached to a sequence \((c_n)\) is given by
\[
r_{n,k}=c_n[t^n]\,d(t)\,\frac{h(t)^k}{c_k}.
\]
This recovers exponential Riordan arrays when \(c_n=n!\) and ordinary Riordan arrays when \(c_n=1\). The same framework derives the group law, inverse, Lagrange inversion, nonrecursive entry formulas, and generalized recurrences from Abel-type identities, rather than taking generating-function manipulations as primary [1505.07253].

The persistence of the same two operations—multiplication by an invertible prefactor and substitution by a compositional unit—is the organizing principle across the later generalizations. Even when the ambient algebra, indexing set, or coefficient basis changes, the generalized object is typically designed so that matrix multiplication still represents an operator of “multiply, then substitute” type.

## 2. Weighted, \(q\)-dependent, and generalized Pascal-matrix frameworks

One major route to generalization changes the coefficient algebra underlying the array. A generalized binomial system is built from a sequence \(c_n\) with \(c_0=1\), generalized factorials
\[
c_n! = \prod_{m=1}^n c_m,
\]
and generalized binomial coefficients
\[
\binom{n}{m}_c = \frac{c_n!}{c_m!\,c_{n-m}!}.
\]
The associated generalized Pascal matrix \(P_c\) is central, and generalized Riordan arrays are obtained by diagonal conjugation,
\[
|c(x)|^{-1}(b(x),a(x))|c(x)|=(b(x),a(x))_c.
\]
In this setting, the parameter \(q\) organizes a distinguished family: \(q=0\) gives ordinary Riordan arrays, \(q=1\) gives exponential Riordan arrays, \(q\neq -1\) gives arrays associated with \(q\)-binomial coefficients, and \(q=-1\) is exceptional and leads to a separate zero-generalized-Pascal theory [1612.07657].

A parallel formulation replaces the coefficient sequence \(c\) by a general weight
\[
W(t)=\sum_{n\ge 0} w_nt^n,\qquad w_0=1,\quad w_n\neq 0,
\]
and defines \(W\)-Riordan arrays through weighted column series. In the standard form, a \(W\)-Riordan array is represented by
\[
a(y)\,W\!\bigl(x\,\beta(y)\bigr),\qquad v(a)=0,\quad v(\beta)=1,
\]
and the resulting group \(R_W\) is isomorphic to
\[
F[[y]]^\times \rtimes yF[[y]]^\times.
\]
The operator interpretation is equally important: \(W\)-Riordan arrays are exactly the matrices of \(W\)-Sheffer sequences, and all groups arising from different weights are conjugate inside the full lower-triangular group of graded polynomial sequences [1505.03100].

The zero-generalized-Pascal extensions enlarge the subject further by admitting limiting or degenerate central matrices such as \(0,qP\) and the fractal products
\[
[0,q]P = 0,qP \times 0,q^2P \times 0,q^3P \times \cdots.
\]
These constructions include the \(q=-1\) case and Pascal triangle modulo \(2\), and the associated matrix groups preserve a Riordan-like multiplication law with a modified algebra product. An analogue of the Lagrange inversion theorem is developed in this setting, and zero generalized Pascal matrices are interpreted as “points at infinity” relative to the space of generalized Pascal matrices [2111.04049].

| Framework | Defining data | Distinguished specializations |
|---|---|---|
| \((c)\)-Riordan arrays | Sequence \(c_n\), conjugation by \(|c(x)|\) | ordinary, exponential, \(q\)-binomial |
| \(W\)-Riordan arrays | Weight \(W(t)=\sum w_nt^n\) | weighted Sheffer/Riordan theories |
| Zero generalized Pascal setting | \(0,qP\), \(P[o,q]\) | \(q=-1\), Pascal triangle modulo \(2\) |

These formulations are structurally close but not identical. A plausible implication is that “generalized Riordan array” functions as an umbrella term for extensions that preserve some Riordan mechanisms while varying the underlying coefficient geometry.

## 3. Higher-periodic, bisection-based, and multivariate generalizations

Another route to generalization replaces the single substitution series \(f\) by several series interleaved in a periodic pattern. The triple Riordan group \(TR\) consists of \(4\)-tuples
\[
(g,f_1,f_2,f_3),
\]
with
\[
g\in \mathbf{R}[[x^3]],\qquad f_1,f_2,f_3\in x\mathbf{R}[[x^3]],
\]
and product law
\[
(g,f_1,f_2,f_3)\cdot(G,F_1,F_2,F_3)
=
\bigl(g\,G(h),\,f_1F_1(h),\,f_2F_2(h),\,f_3F_3(h)\bigr),
\]
where
\[
h(x)=\bigl(f_1(x)f_2(x)f_3(x)\bigr)^{1/3}.
\]
Its bivariate generating function is
\[
\frac{g\bigl(1+yf_1+y^2f_1f_2\bigr)}{1-y^3f_1f_2f_3},
\]
and the specialization \((g,f,f,f)\) embeds an ordinary Riordan array into the triple theory. The paper presents this as the cubic continuation of the double Riordan construction and as a template for further \(m\)-fold periodic groups [2412.05461].

A different higher-order extension is the Sprugnoli group, whose elements are triples \((g,f_1,f_2)\) with \(g\in F_0\), \(f_1\in F_1\), and \(f_2\in xR[x^2]\). Here the columns alternate according to even and odd bisections, and the fundamental action on a series \(h(x)=\sum a_nx^n\) is
\[
(g,f_1,f_2)\cdot h(x) = g(x)\,h_e(xf_2(x)) + g(x)f_1(x)\,h_o(xf_2(x)),
\]
where
\[
h_e(x)=\sum_{n\ge 0} a_{2n}x^n,\qquad h_o(x)=\sum_{n\ge 0} a_{2n+1}x^n.
\]
The bivariate generating function is
\[
g(x)\,\frac{1+y f_1(x)}{1-y^2 x f_2(x)},
\]
the identity is \((1,x,x)\), and the production matrix involves alternating \(A\)- and \(B\)-stripe data rather than a single \(A\)-sequence [2605.16633].

A genuinely different notion of higher generalization is multivariate. In \(d\) variables, one works in
\[
F_d=K[[x_1,\dots,x_d]],
\]
with units \(F_d^\times\) and the group \(G_d\) of formal maps \(g=(g_1,\dots,g_d)\) fixing the origin and compositional invertible. The higher-dimensional Riordan group is
\[
R_d(K)=F_d^\times \rtimes G_d^{\mathrm{op}},
\]
with law
\[
(f,g)(f',g')=\bigl(f\cdot (f'\circ g),\,g'\circ g\bigr),
\]
and matrices are indexed by monic monomials rather than by integers. The corresponding multivariate Fundamental Theorem states that matrix multiplication exactly matches multiplication in \(F_d\) together with composition by formal automorphisms [2009.14130].

These constructions represent different meanings of “higher-order.” The triple and Sprugnoli groups keep one formal variable but impose periodic or parity structure; the multivariate groups replace the one-variable indexing entirely.

## 4. Structural transformations, characterizations, and near-Riordan objects

Generalized Riordan theory is also driven by transformations that preserve, deform, or reconstruct Riordan structure. The production-matrix approach is central here. For an ordinary Riordan array \(M=(g,f)\), the ordinary production matrix is lower Hessenberg and determined by the \(Z\)- and \(A\)-data; the second production matrix is obtained by deleting the first column of the ordinary production matrix, and the matrix generated by it is again a Riordan array with components expressed explicitly in terms of \(g(x)\), \(f(x)\), \(g(f(x))\), and \(f(x)^2/x\). The same pattern is extended to third production matrices and conjecturally to all \(n\)-th production matrices, yielding a canonical hierarchy attached to a single initial Riordan array [2011.13985].

Shifted central extraction gives another systematic construction. If
\[
A=(g(x),xf(x)),\qquad g(0)=1,\ f(0)=1,
\]
and
\[
c(A;r)=\bigl(a_{2n+r,n+k+r}\bigr)_{n,k\ge 0},
\]
then every \(r\)-shifted central triangle is itself a Riordan array, with factorization
\[
c(A;r)=\bigl(\phi'(x)\,f(\phi(x))^{\,r-1},\,\phi(x)\bigr)\cdot A,
\qquad
\phi(x)=\operatorname{Rev}\!\left(\frac{x}{f(x)}\right).
\]
The family satisfies transition relations independent of \(r\), and the generating function of the \(A\)-sequence of \(c(A;r)\) is \(A(x)^2\) when \(A(x)\) is the \(A\)-sequence generating function of the original array [1906.01328].

Partial summation shows how Riordan structure can persist or fail under natural matrix operations. Column partial sums of \((g,f)\) produce the Riordan array
\[
\left(\frac{g(x)}{1-x},\,f(x)\right),
\]
whereas row partial sums generally leave the Riordan group and yield a rectangular matrix whose inverse is lower Hessenberg. A further variant leads to an almost Riordan array linked in the paper to Fibonacci numbers and phyllotaxis, showing that generalized Riordan behavior can survive outside the strict lower-triangular group [2108.13537].

The \(A\)-matrix viewpoint sharpens these structural statements. Riordan arrays can be characterized by local recurrences encoded by an \(A\)-matrix and a \(\rho\)-sequence, and sparse \(A\)-matrices are used to construct Bell-type matrices, quasi-involutions, and moment matrices of perturbed orthogonal polynomials. In particular, certain two-row \(A\)-matrices are conjectured to produce first-column Hankel transforms that satisfy Somos \(4\) recurrences [1912.01126].

Recent work on lattice paths indicates another structural synthesis. A 2025 note states that it describes Riordan arrays and some of their generalizations, including rectifications and triangulations, links Riordan arrays to lattice paths with steps of type \((a,b)\), treats almost Riordan arrays in lattice-path models, and uses the \(A\)-matrix characterization to analyze path families with downward steps [2504.09719].

## 5. Polynomial triangles, Eulerian transforms, and combinatorial matrix families

Generalized Riordan arrays frequently appear as coefficient arrays of polynomial systems and as transformation devices between classical combinatorial triangles. One important example is the theory of generalized Euler polynomials. For a polynomial sequence \(u_n(x)\) of degree \(n\),
\[
\alpha_n(x)=(1-x)^{n+1}\sum_{m=0}^{\infty}u_n(m)x^m
\]
is realized as the normalized row polynomial of a square or generalized Riordan array \((1,a(x))\). The matrices \(U_n\) and \(V_n\) mediate transformations involving \(\log a(x)\), \(a(x)-1\), Stirling numbers, multinomial coefficients, Bell polynomials, and generalized Lagrange series. In this formulation, classical Euler polynomials arise as the special case \(a(x)=e^x\) [1709.02229].

Pascal-like triangles form another major branch. For the ordinary Riordan family
\[
\left(\frac{1}{1-x},\,\frac{x(1+rx)}{1-x}\right),
\]
the rows are palindromic and admit \(\gamma\)-vector decompositions. The associated \(\gamma\)-matrix is the stretched Riordan array
\[
\left(\frac{1}{1-x},\,rx^2\right),
\]
with entries
\[
Y_{n,k}=\binom{n-k}{k}r^k.
\]
Reversion of these Pascal-like triangles produces generalized Narayana families, and their generating functions are expressed by Jacobi continued fractions [1804.05027].

The related \(f\)-matrix theory studies the binomial multiplication
\[
F=M B,
\]
where \(B\) is the Pascal matrix. For the same ordinary Riordan family, the \(f\)-matrix is
\[
F_r=\left(\frac{1}{1-2x-rx^2},\,\frac{x}{1-2x-rx^2}\right),
\]
while in the exponential case
\[
[e^x,\,x(1+rx/2)]\cdot [e^x,x]=[e^{2x+rx^2/2},\,x(1+rx/2)].
\]
The paper emphasizes that the \(y\)-, \(h\)-, and \(f\)-matrices are linked by systematic transformations of the coefficients in corresponding Jacobi continued fractions [1805.02274].

A further parametric deformation is obtained by multiplying any Riordan array \(R=(d(z),h(z))\) on the right by the generalized Pascal array
\[
P(r)=\left(\frac{1}{1-rz},\,\frac{z}{1-rz}\right),
\]
producing
\[
R(r)=R\,P(r)=\left(\frac{d(z)}{1-rh(z)},\,\frac{h(z)}{1-rh(z)}\right).
\]
For the Catalan array
\[
C=(c(z),zc(z)),\qquad c(z)=\frac{1-\sqrt{1-4z}}{2z},
\]
this yields the \(r\)-Catalan triangle
\[
C(r)=C\,P(r),
\]
which interpolates between the Aigner array at \(r=0\) and the Shapiro array at \(r=1\). The associated Sheffer polynomials satisfy the shift relation
\[
p_n^{R(r)}(x)=p_n^R(x+r),
\]
and specializations lead to identities involving Chebyshev polynomials of the second kind, periodic sequences, and Fibonacci numbers [1505.05568].

## 6. Orthogonal-polynomial, Hankel-transform, and semi-Laurent directions

A substantial part of generalized Riordan theory lies at the interface with orthogonal and biorthogonal polynomials. Laurent biorthogonal polynomials with constant-coefficient three-term recurrence
\[
P_n(x)=(x-a)P_{n-1}(x)-Bx\,P_{n-2}(x)
\]
have coefficient arrays that are Riordan arrays, for example
\[
\left(\frac{1}{1+ax},\,\frac{x(1-Bx)}{1+ax}\right),
\]
or, in the preferred normalization,
\[
\left(\frac{1-Bx}{1+ax},\,\frac{x(1-Bx)}{1+ax}\right).
\]
The inverse Riordan array gives a natural associated family of ordinary orthogonal polynomials, and the moment Hankel transform is
\[
h_n=\bigl(B(a+B)\bigr)^{\binom{n+1}{2}}.
\]
The theory extends to a three-parameter generalized recurrence
\[
P_n(x)=(x-a)P_{n-1}(x)-B(x-y)P_{n-2}(x),
\]
which still admits a Riordan-array description [1311.2292].

Catalan-based recurrence solving is another prominent application. Generalized Catalan-Schroeder convolution recurrences are solved by writing their generating functions as Riordan transforms of the Catalan series
\[
c(x)=\frac{1-\sqrt{1-4x}}{2x}.
\]
The resulting sequences are often moment sequences for orthogonal polynomials, and many of their Hankel transforms are conjectured or shown to satisfy Somos \(4\) recurrences. The same framework is used to connect Riordan arrays with elliptic curves, continued fractions, and orthogonal-polynomial moment matrices [1910.00875].

Stretched Riordan arrays amplify this Hankel-transform direction. In conjectural families built from generalized Jacobi continued fractions and Catalan numbers, coefficient sequences are described by stretched Riordan arrays of the form \((g(x),x^rf(x))\), and their Hankel transforms are conjectured to be Somos \(4\), Somos \(6\), or Somos \(8\) sequences. In this setting, the Riordan array provides the coefficient-extraction mechanism, Catalan numbers provide the seed sequence, and the Hankel transform supplies the bridge to bilinear integrable recurrences [2211.12637].

A more radical generalization replaces formal power series by formal semi-Laurent series. If \(g\) is a formal semi-Laurent series and \(f\) is a formal power series with \(\ord f=1\), the generalized Riordan array, or Riordan Laurent matrix, is the bi-infinite matrix
\[
(g,f)_{m,n}=[z^m](g f^n),\qquad m,n\in\mathbb Z.
\]
The multiplication law is
\[
(g,f)*(h,k)=\bigl(g\cdot(h\circ f),\,k\circ f\bigr),
\]
with inverse
\[
(g,f)^{-1}=\bigl(g^{-1}\circ f^{[-1]},\,f^{[-1]}\bigr).
\]
The order-zero subgroup \(RL_1^0(\mathbb K)\) is identified with the classical Riordan group, but the generalized group is proved not to be isomorphic to the classical one. Over \(\mathbb C\), the same paper constructs an infinite-dimensional Lie group modeled on \(s_L\times s\) and a Lie algebra
\[
rl_1(\mathbb C)=\mathcal L_s(\mathbb C)\times z\mathbb X(\mathbb C)
\]
with bracket
\[
[(\gamma_1,\kappa_1),(\gamma_2,\kappa_2)]
=
(\gamma_2'\kappa_1-\gamma_1'\kappa_2,\ \kappa_2'\kappa_1-\kappa_1'\kappa_2).
\]
This places generalized Riordan arrays in a setting that is simultaneously combinatorial, algebraic, and Lie-theoretic [2509.04160].

One persistent misconception is that all generalized Riordan theories are equivalent reformulations of the same object. The literature does not support that view. Weighted \(W\)-Riordan groups are conjugate for different weights [1505.03100], but semi-Laurent generalized Riordan groups are proved not to be isomorphic to the classical Riordan group [2509.04160], and periodic constructions such as the triple and Sprugnoli groups introduce genuinely different column-generation mechanisms [2412.05461][2605.16633]. The modern subject is therefore unified by method and operator form, not by a single definition.

Source: https://www.emergentmind.com/topics/generalized-riordan-arrays