--- title: Hsu–Shiue Family & Generalized Stirling Numbers url: https://www.emergentmind.com/topics/hsu-shiue-family type: topic --- # Hsu–Shiue Family & Generalized Stirling Numbers Searching arXiv for recent and foundational papers on the Hsu–Shiue family and related generalized Stirling-number frameworks. The expression **Hsu–Shiue family** most commonly denotes the three-parameter family of generalized Stirling numbers introduced by Hsu and Shiue, together with the associated polynomial, functional, and Riordan-array structures built from them. In this usage, the family consists of connection coefficients between generalized factorial bases and includes classical Stirling numbers, Lah-type arrays, and multiple noncentral or weighted variants as specializations. A separate special-function usage applies the Hsu–Shiue name to an elliptic-function family obtained from the inverse of a hypergeometric integral involving ${}_2F_1\!\left(\frac13,\frac23;\frac12;\bullet\right)$; that construction was developed by Li-Chien Shen and later recast with streamlined notation and new proofs. Across these settings, the common theme is a parametric coefficient system that mediates between natural bases or uniformizations [1106.5251; 2207.10224; 1907.09938]. ## 1. Foundational definition as generalized Stirling numbers In the combinatorial literature, the Hsu–Shiue family is introduced through a generalized factorial-basis expansion. One standard form is \[ (z)_{n,-a}=\sum_{k=0}^n S(n,k,a,\beta,r)\,(z-r)_{k,-\beta}, \] with \[ (z)_{n,-a}=z(z-a)(z-2a)\cdots(z-(n-1)a). \] A closely related notation used elsewhere is \[ (x)_{n,a}=\sum_{k=0}^n S_{n,k}(a,b;r)\,(x-r)_{k,b}. \] In both forms, the Hsu–Shiue numbers are connection coefficients between two factorial polynomial bases [1106.5251; 2207.10224]. The family is normalized by apex and boundary data such as \[ S(0,0)=1,\qquad S(n,n)=1, \] and, in the formulation of the unified Stirling framework, \[ S(1,0)=r,\qquad S(n,0)=(r)_{n,-a}. \] These coefficients interpolate classical arrays. Two basic specializations are \[ S(n,k,1,0,0)=s(n,k),\qquad S(n,k,0,1,0)=S(n,k), \] recovering the classical Stirling numbers of the first and second kinds, respectively. Other specializations listed in the literature include Lah numbers, Riordan’s non-central Stirling numbers, Carlitz’s weighted and degenerate variants, Charalambides–Koutras numbers, Gould–Hopper non-central Lah numbers, Tsylova and Todorov families, Ahuja–Enneking associated Lah numbers, and Broder’s \(r\)-Stirling numbers [1106.5251]. This definition places the Hsu–Shiue family in the connection-coefficient tradition: the array is determined not by a single recurrence alone, but by how one factorial basis expands in another. That viewpoint is central in later work on GKP triangles, generalized Eulerian numbers, and operator ordering [2207.10224; 2308.10332]. ## 2. Core identities, finite-difference formulas, and recurrence structure Several papers emphasize that the Hsu–Shiue family admits exact coefficient extraction by generalized divided differences. One form is \[ S(n,k,a,\beta,r)=\left.\Delta_\beta^k (z)_{n,-a}\right|_{z=r} = \begin{cases} \beta^k k!\,(z)_{n,-a}[r,r+\beta,\dots,r+k\beta], & \beta\neq 0,\\[4pt] D^k (z)_{n,-a}\big|_{z=r}, & \beta=0. \end{cases} \] For \(\beta\neq 0\), this yields the explicit alternating sum \[ S(n,k,a,\beta,r)=\frac{1}{\beta^k k!}\sum_{j=0}^{k}(-1)^j\binom{k}{j}\bigl(r+(k-j)\beta\bigr)_{n,-a}. \] Equivalent finite-difference forms also appear as \[ S_{n,k}(a,b;r)=\frac{1}{b^k k!}\sum_{j=0}^k (-1)^{k-j}\binom{k}{j}(bj+r)_{n,a}, \] or \[ S_{n,k}(a,b;r)=\frac{1}{b^k k!}\,\Delta_b^k\!\bigl[(bx+r)_{n,a}\bigr]\Big|_{x=0}, \] depending on the chosen parametrization [1106.5251; 2207.10224]. The defining recurrence is of Graham–Knuth–Patashnik type. In one notation, \[ S(n,k;a,\beta,r)=S(n-1,k-1;a,\beta,r)+\bigl(a(n-1)+\beta k+r\bigr)\,S(n-1,k;a,\beta,r), \] while in another, \[ S_{n+1,k+1}=\bigl[-an+b(k+1)+r\bigr]\,S_{n,k+1}+S_{n,k}. \] The same recurrence governs the Hsu–Shiue coefficients when they are used in generalized geometric-polynomial theory, barred preferential arrangements, and boson-ordering expansions [1412.1110; 2207.10224; 2508.13094]. A further structural property is homogeneity: \[ S_{n,k}(Aa,Ab;Ar)=A^{\,n-k}S_{n,k}(a,b;r). \] The family also admits matrix composition and inversion laws. One explicit identity is \[ S(a,c;r_1+r_2)=S(a,b;r_1)\,S(b,c;r_2), \] with inverse \[ S(a,b;r)^{-1}=S(b,a;-r). \] These formulas place the family within a Riordan-type group structure and explain why binomial transforms and basis changes recur so persistently in applications [2207.10224; 2308.10332]. ## 3. Analytic extensions: generating functions, Riordan arrays, and Stirling functions The Hsu–Shiue numbers generate a natural row-polynomial sequence, and the resulting array is an exponential Riordan array. In the notation of the Sheffer-operator literature, the row polynomials are \[ {_n(t;\,A,B,r)}=\sum_{k=0}^n {_{n,k}(A,B,r)}\,t^k, \] with Riordan pair \[ d_{HS}(z;A,r)=(1+Az)^{r/A},\qquad h_{HS}(z;A,B)=\frac{(1+Az)^{B/A}-1}{B}. \] Their exponential generating function is \[ \sum_{n=0}^\infty {_n(t;\,A,B,r)}\frac{z^n}{n!} = (1+Az)^{r/A}\exp\!\left\{\frac{t}{B}\left[(1+Az)^{B/A}-1\right]\right\}. \] The standard exponential-Riordan coefficient formula \[ s_{n,k}=\frac{n!}{k!}\,[z^n]\left[d(z)h(z)^k\right] \] applies directly after specialization to \(d_{HS}\) and \(h_{HS}\) [2508.13094]. Equivalent vertical generating functions appear in the GKP analysis: \[ \sum_{n=0}^{\infty} S_{n,k}(a,b;r)\,\frac{z^n}{n!} = \frac{(1+az)^{r/a}}{k!} \left(\frac{(1+az)^{b/a}-1}{b}\right)^k. \] This identifies \(S(a,b;r)=(S_{n,k}(a,b;r))\) as an exponential Riordan array and situates the family as the canonical triangle for one of the principal GKP cases [2207.10224]. The analytic generalization extends beyond integer row index. A unified Stirling-function theory defines \[ S(y,n,a,\beta,r;\varepsilon) :=\beta^n\Gamma(n+1)\lim_{z\to r}\Delta_{\beta,\varepsilon}^n\bigl((z)_{y,-a}\bigr), \] with absolutely convergent series representation \[ S(y,n,a,\beta,r;\varepsilon) = \beta^n\Gamma(n+1)\sum_{j\ge 0} (-1)^j\binom{n}{j}e^{(n-j)\varepsilon}\,(r+(n-j)\beta)_{y,-a}. \] For integer \(y=n\) and \(\varepsilon=0\), these recover the generalized Stirling numbers. The same framework relates the family to generalized factorials \((z)_{y,\pm k}\), the \(k\)-gamma function \(\Gamma_k(z)\), and recurrences such as \[ S(y,n,a,\beta,r;\varepsilon) = \bigl(r+n\beta-(y-1)a\bigr)\,S(y-1,n,a,\beta,r;\varepsilon) + S(y-1,n-1,a,\beta,r;\varepsilon). \] This analytic continuation is one reason the Hsu–Shiue family functions as more than a discrete triangle [1106.5251]. ## 4. Combinatorial realizations and polynomial hierarchies A major development is the identification of the Hsu–Shiue numbers as weighted counts of **extended Lah distributions**. In that model, an extended Lah distribution is obtained from a Lah distribution by circling a subset of special elements, with the restriction that \(1\) may be circled only if it starts a block. For \(\lambda\in L_n^*\), the weight is \[ w(\lambda)=a^{\operatorname{nrec}(\lambda)}\,\beta^{\operatorname{rec}^*(\lambda)}\,r^{\operatorname{circ}(\lambda)}, \] and the Hsu–Shiue numbers satisfy \[ S(n,k;a,\beta,r)=\sum_{\lambda\in L_{n,k}^*} w(\lambda). \] This interpretation makes the recurrence transparent: the largest label \(n\) either forms a new true block, enters as a non-record-low element, begins a true block as a non-minimal record low, or appears as a circled element, producing exactly the recurrence coefficients \(1\), \(a(n-1)\), \(\beta k\), and \(r\) [1412.1110]. The same combinatorial infrastructure supports a refinement of Xu’s extension of Spivey’s Bell-number formula. Writing \[ B_{n;a,\beta,r}(x)=\sum_{k=0}^n S(n,k;a,\beta,r)x^k, \] one obtains \[ S(m+n,k;a,\beta,r) = \sum_{i=0}^n \sum_{j=0}^m \binom{n}{i}\, S(m,j;a,\beta,r)\, S(i,k-j;a,\beta,r)\, \prod_{l=0}^{n-i-1}\bigl((m+l)a+j\beta\bigr), \] and hence \[ B_{m+n;a,\beta,r}(x) = \sum_{i=0}^n \sum_{j=0}^m \binom{n}{i}\, S(m,j;a,\beta,r)\, B_{i;a,\beta,r}(x)\, \prod_{l=0}^{n-i-1}\bigl((m+l)a+j\beta\bigr). \] These formulas treat the Hsu–Shiue family as a weighted Stirling system underlying Bell-type convolution identities [1412.1110]. Another line of work uses the Hsu–Shiue numbers as the coefficient system for generalized geometric polynomials and generalized barred preferential arrangements. The generalized geometric polynomials are \[ G_n(a,\beta,\gamma;x)=\sum_{k=0}^n \beta^k k! \, S(n,k,a,\beta,\gamma)\,x^k, \] with generating function \[ \sum_{n=0}^{\infty} G_n(a,\beta,\gamma;x)\frac{t^n}{n!} = \frac{(1+at)^{\gamma/a}}{1-x\left((1+at)^{\beta/a}-1\right)}. \] This generating series is interpreted via barred preferential arrangements whose sections have property 1 or property 2, defined by cyclically ordered labelled compartments of sizes \(\beta\) or \(\gamma\) [1908.05014]. The same negative-parameter mechanism produces a **second type of higher order generalised geometric polynomials** \[ \sum_{n=0}^\infty A_n^{\,x}(a,B,y)\,\frac{t^n}{n!} = \frac{(1-at)^{-y/a}}{\bigl(1-x[(1-at)^{-B/a}-1]\bigr)}, \] with coefficient formula \[ A_n^{\,x}(a,B,y) = \sum_{k=0}^n (-1)^{n+k}\,B^k\,k!\,S(n,k;a,-B,-y)\,x^k. \] These polynomials are counted by barred preferential arrangements in which one section has a single cell with \(y\) compartments and the remaining \(x\) sections have cells with \(B\) compartments, all subject to the same splitting rule. Setting \(x=-\tfrac12\) yields the generalized Euler polynomials discussed in the same paper [2004.01490]. ## 5. Transform relations, generalized Eulerian numbers, and boson-operator ordering The Hsu–Shiue triangles are closely linked to generalized Eulerian triangles by row-wise binomial transforms. In the GKP framework, the generalized Eulerian numbers \(E_{n,k}(a,b;c_0,c_\infty)\) are related to Hsu–Shiue numbers by the upper binomial transform and reflection: \[ E_{n,k}(a,b;c_0,c_\infty) = \sum_{j=k}^n (-1)^{j-k}\binom{j}{k}\, (c_0+c_\infty)^{\,n-j}_{\,b}\, S_{n,n-j}(-a,b;c_\infty), \] with inverse \[ (c_0+c_\infty)^{\,n-k}_{\,b}\,S_{n,n-k}(-a,b;c_\infty) = \sum_{j=k}^n \binom{j}{k}\,E_{n,j}(a,b;c_0,c_\infty). \] The corresponding generalized Worpitzky identity is \[ (c_0+c_\infty)_{n,b}\,(x)_{n,a} = \sum_{k=0}^n E_{n,k}(a,b;c_0,c_\infty)\,(x-c_0)_{k,b}\,(x+c_\infty)_{n-k,b}. \] In this formulation, the Hsu–Shiue family is the one-sided factorial-basis connection array, while the generalized Eulerian family is the two-sided or bifactorial connection array [2207.10224]. A parallel treatment in the Weyl–Heisenberg algebra shows that powers of single-annihilator boson strings are governed by Hsu–Shiue coefficients. For \[ \Omega=a^{\dagger L} a\, a^{\dagger R},\qquad e=L+R-1, \] one normal-ordering identity is \[ \left(a^{\dagger L} a a^{\dagger R}\right)^n = \left(a^{\dagger e}\right)^n \left[\sum_{k=0}^n {_{n,k}(-e,1,R)\,a^{\dagger k} a^k}\right], \] while the anti-normal companion is \[ \left(a^{\dagger L} a a^{\dagger R}\right)^n = \left[\sum_{k=0}^n {_{n,k}(e,-1,-L)\,a^k a^{\dagger k}}\right] \left(a^{\dagger e}\right)^n. \] In the \(x,D\) realization, the key substitution is \(x\mapsto xD\) together with \[ (xD)^n=x^n D^n, \] which converts factorial-basis expansions directly into operator identities [2308.10332; 2508.13094]. The exponential generating functions of the Hsu–Shiue polynomials then yield explicit ordered exponentials. One formula is \[ {\rm e}^{\lambda(a^{\dagger L} a a^{\dagger R})} = :{\mathcal{HS}(a^\dagger a, \lambda a^{\dagger e};\, -e,1,R)}_N = :{(1-e\lambda a^{\dagger e})^{-R/e} \exp\left\{ \left[\bigl(1-e\lambda a^{\dagger e}\bigr)^{-1/e}-1\right] a^\dagger a \right\}}_N, \] with an analogous anti-normal formula. A later extension introduces a **two-point Hsu–Shiue family** \(\mathcal T\) that interpolates between the \(s=-1\) and \(s=+1\) endpoints and gives a unified \(s\)-ordered expression for \[ {\rm e}^{\lambda(a^{\dagger L} a a^{\dagger R})}. \] This places the Hsu–Shiue family simultaneously in combinatorics, Riordan-array theory, and quantum-operator ordering [2508.13094]. ## 6. The distinct Hsu–Shiue/Shen elliptic-function family A separate usage of the name occurs in elliptic-function theory. For \(00, \] so the period lattice is rectangular with fundamental periods \(2K\) and \(2iK'\). The poles of \(d\) occur where \(\wp(u)=-2\), and the commentary proves that \(u=\pm 2iK'\) are poles. Since \(d\) is elliptic of order two, these account for the poles modulo periods [1907.09938]. This theory also clarifies which associated functions are elliptic. From \[ 4k^2s^2=(1-d)(2+d)^2,\qquad c^2=1-s^2, \] it follows that \(s^2\) and \(c^2\) are elliptic; because \(d\) has simple poles, \(s^2\) and \(c^2\) have triple poles. By contrast, \(s\) and \(c\) themselves are not elliptic, although the product \(sc\) is elliptic and satisfies \[ sc=-\frac{1}{2}(2+d)d'=-\frac12\bigl((d+2)^2\bigr)'. \] In modern language, this family gives an explicit hypergeometric uniformization of a genus-one algebraic curve, with the inverse hypergeometric construction birationally equivalent to the Weierstrass model [1907.09938]. Source: https://www.emergentmind.com/topics/hsu-shiue-family