Convolution identities of stirling numbers

Nadia N. Li1, Wenchang Chu2
1School of Mathematics and Statistics Zhoukou Normal University, Henan, China
2Via Dalmazio Birago 9/E, Lecce 73100, Italy

Abstract

By means of the generating function method, a linear recurrence relation is explicitly resolved. The solution is expressed in terms of the Stirling numbers of both the first and the second kind. Two remarkable pairs of combinatorial identities (Theorems 3.1 and 3.3) are established as applications, that contain some well–known convolution formulae on Stirling numbers as special cases.

Keywords: Recurrence relation, Convolution formula, Generating function, Stirling number of the first kind, Stirling number of the second kind

1. Introduction

Denote by N the set of natural numbers with N0=N{0}. For an indeterminate x, define the rising and falling factorials by (x)0=x01 and (x)n=x(x+1)(x+n1)fornN,xn=x(x1)(xn+1)fornN.

Then the unsigned Stirling numbers of the first kind [nk] are determined as the connection coefficients: xn=k=0n(1)nk[nk]xk. Analogously, the Stirling numbers of the second kind {nk} are given by xn=k=0n{nk}xk, which admits the explicit expression {nk}=1k!j=0k(1)kj(kj)jn.

These numbers appear frequently in mathematical literatures and have wide applications in combinatorics and number theory (see [2] Chapter 5, [4] §6.2 and [5], just for example).

Recently, Stenlund [6] introduced an interesting bivariate polynomial sequence {Pn(x,z)}nN through the following recurrence relation: (1)Pn+1(x,z)=x(n+zn)xm=1n(nm+znm+1)Pm(x,z)withP1(x,z)=x.

The same author not only found out an explicit double sum expression (2)Pn(x,z)=k=1nj=1k(1)j1(j1)!(n1)![nk]{kj}xjzk1, but also explored applications to combinatorial identities and probability theory.

Inspired by the above work of Stenlund [6], we shall examine the extended polynomial sequence {Qn}nN0 (with three variables λ, x and y) defined by the recurrence relation (3)Qn=(λn)(1)nxk=1n(1)k(yk)Qnk,withQ0=1. The first three terms are recorded as follows: Q1=xyλ,Q2=xy2(1y+2xy)λ2(1λ+2xy),Q3=xy6(y1)(y2)λ6(λ1)(λ2)xy(1y+xy)(λxy)+λxy2(λy).

These Qn-polynomials in (3) generalize the Pn-polynomials in (1) as we can verify below that Qn becomes Pn+1/x when λ=1z and yz. Denoting by Q~n the resulting Qn-polynomial under the above replacements. Then we have Q~0=1 and the recurrence relation Q~n=(1zn)(1)nxk=1n(1)k(zk)Q~nk=(z+nn)xk=1n(k1+zk)Q~nk=(z+nn)xm=1n(nm+znm+1)Q~m1k1+nm.

Reformulating the above equality by (xQ~n)=x(z+nn)xm=1n(nm+znm+1)(xQ~m1), we can see that “xQ~n” satisfies the same relation (1) as Pn+1(x,z) with the same initial condition. Hence, Qn indeed extends Pn as claimed.

The rest of the paper will be organized as follows. In the next section, we shall derive the generating function and the explicit formulae for the Qn-polynomials. As applications, two pairs of convolution sums containing both kinds of Stirling numbers will be evaluated in closed forms in Section 3.

2. Generating function and explicit expressions

For the sequence {Qn}nN0 defined by (3), we shall derive its generating function and explicit formulae. The informed reader will notice that this treatment is more transparent than the induction approach adopted by Stenlund [6].

Multiplying both sides of (3) by τn and then summing over n for n0, we can manipulate the generating function Ω(τ)=n=0Qnτn=1+n=1Qnτn=1+n=1(λn)(τ)nxn=1k=1n(1)k(ykQnk)τn=(1τ)λxk=1(1)k(yk)τkn=kQnkτnk.

This leads us to the functional equation Ω(τ)=(1τ)λxΩ(τ){(1τ)y1}, and the rational function expression.

Lemma 2.1 (Generating function). Ω(τ)=(1τ)λ1x+x(1τ)y.

By writing as a geometric series Ω(τ)=11x×(1τ)λ1xx1(1τ)y=k=0(x)k(1x)k+1(1τ)λ+ky, and then extracting the coefficient of τn, we find the infinite series expression.

Proposition 2.2 (Single sum formula). Qn=k=0(1)n+k(λ+kyn)xk(1x)k+1.

Expanding the rightmost fraction into binomial series, we have further Qn=k=0(1)n+k(λ+kyn)j=0(k+jjxk+j)=m=0(1)nxmk=0m(1)k(mk)(λ+kyn), where the last line is justified by the replacement k+j=m. Observe that the inner sum with respect to k results substantially in the differences of order m about a polynomial of degree n, which is equal to zero when m>n (cf. [3] Equation 3.150). Therefore, we have derived a finite double sum expression.

Proposition 2.3 (Double sum formula). Qn=m=0nk=0m(1)n+k(mk)(λ+kyn)xm.

Expressing the last binomial coefficient in terms of unsigned Stirling numbers and then making use of the binomial theorem, we have (λ+kyn)=λ+kynn!=i=0n(1)nin![ni](λ+yk)i=i=0n(1)nin![ni]j=0i(ij)λij(ky)j.

By substitutions, Qn can be written as a four–fold sum Qn=m=0nk=0m(mk)xmi=0n(1)kin![ni]j=0i(ij)λij(ky)j,=m=0nxmi=0n(1)min![ni]j=0i(ij)λijyjk=0m(1)mk(mk)kj.

Writing the rightmost sum in terms of Stirling number of the second kind, we get a triple sum expression (4)Qn=m=0nxmi=0n(1)mim!n![ni]j=0i(ij){jm}λijyj.

Taking into account that {jm}=0 for 0j<m and (ij)=0 for 0i<j, we can restate the last formula as follows.

Proposition 2.4 (Triple sum formula) Qn=m=0ni=mnj=mi(1)mim!n![ni](ij){jm}xmλijyj.

By comparing the coefficients of xm between the two expressions in Proposition 2.3 and Proposition 2.4, we derive the following summation formula.

Theorem 2.5 (Convolution formula). i=mnj=mi(1)n+i[ni](ij){jm}λijyj=n!m!k=0m(1)m+k(mk)(λ+kyn).

As kindly pointed out by an anonymous referee, letting λ=0 in Theorem 2.5 gives a nicer formula than the corresponding ones found by Stenlund [6] Theorem 4 and Corollary 5).

3. Applications to convolution formulae

By specifying y to particular values, we shall derive from Theorem 2.5 two pairs of combinatorial identities as in Theorems 3.1 and 3.3 that extend, with an extra parameter “λ“, the related known convolution formulae on Stirling numbers.

The first pair of identities are given in the following theorem.

Theorem 3.1 (1 ≤m≤n). (a)i=mnj=mi(1)n+i[ni](ij){jm}λij=n!m!(λnm),(b)i=mnj=mi(1)i+j[ni](ij){jm}λij=n!m!(λmnm)(1)m+n.

This theorem may be considered as significant extensions of the two well–known results (cf. [4]§6.1) recorded in the corollary below, that correspond to the special cases of λ=0.

Corollary 3.2 (1 ≤m≤n). (a)k=mn(1)nk[nk{km={1,m=n;0,mn;(b)k=mn[nk{km=n!m!(n1m1.{Lah number}

Proof. Proof of Theorem 3.1. When y=1 in Theorem 2.5, the corresponding equality becomes i=mnj=mi(1)n+i[ni](ij){jm}λij=n!m!k=0m(1)m+k(mk)(λ+kn).

Then the first formula “(a)” follows by evaluating the above sum on the right (cf. [3] Eq.  3.47) k=0m(1)m+k(mk)(λ+kn)=k=0m(1)m+k(mk)=0n(λn)(k)==0n(λn)(m)k=m(1)m+k(mk)==0n(λn)(m(11)m)=(λnm), where the Chu–Vandemonde convolution formula and the binomial theorem have been utilized.

Alternatively, for y=1, the corresponding equality in Theorem 2.5 reads as i=mnj=mi(1)i+j[ni](ij){jm}λij=(1)mnn!m!k=0m(1)k(mk)(λkn).

Then we can analogously evaluate the sum on the right (cf. [3]Equation 3.49) as follows: k=0m(1)k(mk)(λkn)=k=0m(1)k(mk)=0n(λmn)(mk)==0n(λmn)(m)k=0m(1)k(mk)==0n(λmn)(m)(11)m=(λmnm).

This proves the second formula “(b)” in the theorem. ◻

Furthermore, setting y=2 and y=1/2 in Theorem 2.5 yields another pair of double sum identities of which the case λ=0 is again known.

Theorem 3.3 (1 ≤m≤n). (a)i=mnj=mi(1)n+i[ni](ij){jm)λij2nj=n!m!=mnm(λn)(m)2m+n2,(b)i=mnj=mi(1)m+i[ni](ij){jm)λij2jm=n!m!=mn(1)m+n(λn)(mm2m).

Proof. Proof of Theorem 3.3. For y=12, the equality in Theorem 2.5 reduces to i=mnj=mi(1)n+i[ni](ij){jm)λij2nj=2nn!m!k=0m(1)m+k(mk)(λ+k2n).

In this case, the first formula “(a)” is confirmed by first reformulating the above sum on the right k=0m(1)m+k(mk)(λ+k2n)=k=0m(1)m+k(mk)=0n(λn)(k2)==0n(λn)k=0m(1)m+k(mk)(k2), and then evaluating the inner sum by the binomial formula (cf. [3]Equation 3.164) k=0m(1)mk(mk)(k2)=m(m2m2).

When y=2, the corresponding equality in Theorem 2.5 can be stated as i=mnj=mi(1)m+i[ni](ij){jm}λij2jm=2mn!m!k=0m(1)n+k(mk)(λ+2kn).

The above sum on the right can be rewritten as k=0m(1)n+k(mk)(λ+2kn)=k=0m(1)n+k(mk)=0n(λn)(2k)==0n(1)n(λn)k=0m(1)k(mk(2k).

Denote by [xn]f(x) the coefficient of xn in the formal power series f(x). Observe that the last inner sum (cf. [3]Equation 3.64) can be expressed and then evaluated as k=0m(1)k(mk)(2k)=[x]k=0m(1)k(mk)(1+x)2k=[x]{1(1+x)2}m=(1)m[xm](2+x)m=(1)m(mm22m). After substitutions, the second formula “(b)” in the theorem is done. ◻

Finally letting λ=0 in Theorem 3.3, we recover the following two elegant identities due to Yang and Qiao [7], who discovered them by employing the Riordan array and expressed the results in terms of Bessel numbers [1].

Corollary 3.4 (1 ≤m≤n). (a)k=mn[nk]{km}(2)nk=(n1)!(m1)!(nnm2mn),(b)k=mn[nk]{km}(2)km=n!m!(mnm)(2)mn).

Declarations

The authors declare no conflict of interest.

References:

  1. G. Cheon, J. Jung, and L. Shapiro. Generalized bessel numbers and some combinatorial settings. Discrete Mathematics, 313(20):2127–2138, 2013. https://doi.org/10.1016/j.disc.2013.05.001.
  2. L. Comtet. Advanced Combinatorics. Dordrecht–Holland, The Netherlands, 1974.
  3. H. W. Gould. Combinatorial Identities. W. V. Morgantown, 1972.
  4. R. L. Graham, D. E. Knuth, and O. Patashnik. Concrete Mathematics. Addison–Wesley Publishing Company, Reading, Massachusetts, 1989.
  5. D. E. Knuth. Two notes on notation. American Mathematical Monthly, 99:403–422, 1992.
  6. D. Stenlund. Some observations on the connection between stirling numbers and bessel numbers. arXiv preprint arXiv:2007.11557, 2020. https://doi.org/10.48550/arXiv.2007.11557.
  7. S. L. Yang and Z. K. Qiao. The bessel numbers and bessel matrices. Journal of Mathematical Research and Exposition, 31(4):627–636, 2011. http://doi.org/10.3770/j.issn:1000-341X.2011.04.006.