arXiv:2101.11039v1 [math.CO] 26 Jan 2021 yl edr set leaders cycle eodkind second fbt iditoue yBoe [ Broder By introduced kind both of blocks. ohkns h uhrcniee that considered author the kinds, both k e. b h eodkind second the hi lcs .e. i. blocks, their numbers. Example. h ubro l emttoso [ of permutations all of number the 1 b , yls(resp. cycles THE n ftems neetn eeaiaino trignmeswsth was numbers Stirling of generalization interesting most the of One e od n phrases. and words Key 2010 Date Let stesm a,let way, same the As ti elkonta the that known well is It 2 b , . . . , • • aur 8 2021. 28, January : σ ahmtc ujc Classification. Subject Mathematics Abstract. trigo h rtkind. first the of Stirling ucin.As,w iealmtrpeetto ftemlil eafun zeta multiple the of representation limit interpre a combinatorial give we using Also, numbers these functions. of properties various fpruain n attoscle ( called partitions and permutations of { For { For 1 1 eapruaino h e [ set the of permutation a be ( , , ,r l, 2 3 k n n , , A . ) 6 6 n SILN UBR:ACMIAOILAPPROACH. COMBINATORIAL A NUMBERS: -STIRLING ,tepermutation the 6, = k ,tepartition the 7, = } }
k . . r lc edr set leaders block onstenme fpermutations of number the counts ) of lcs uhta the that such blocks) denoted , hswr el ihanwgnrlzto of generalization new a with deals work This σ denoted , EBCI HAC BELBACHIR emttos e attos trignmes ymti functio Symmetric numbers, Stirling partitions, Set Permutations, π cl bl eapriino h e [ set the of partition a be ( ( n k trignmeso h rtkind first the of numbers Stirling σ π
= ) = ) cl on h ubro l attoso [ of partitions all of number the count , of ( π σ { { 1 = ,i h e ftesals lmnso hi yls i. cycles, their on elements smallest the of set the is ), π 1. n σ min min denoted , aigexactly having ] 6 1)25()hv h e fcceleaders cycle of set the have (13)(245)(6) = rmr 17,1B3 eodr 50,0A8 05E05. 05A18, 05A05, Secondary 11B83; 11B73, Primary Introduction .Aaoosyt h lsia trignmesof numbers Stirling classical the to Analogously ]. n N N JMAAYAHIA DJEMMADA AND ENE ` , r r = ] 2 c b Siln ubr ftefis kind first the of numbers -Stirling rteeet 1 elements first , ,r l, 1 1 4 , , min min | -trignmeso ohkns estudy We kinds. both of numbers )-Stirling 3 { 1 , 1 5 , bl , c b 2 ( 7 2 2 n , . . . , π , . . . , , . . . , | aetesto lc leaders block of set the have 6 ,i h e ftesals lmnson elements smallest the of set the is ), σ k min min } rs.partitions (resp. yls and cycles, r , n Siln ubr using numbers -Stirling 2 having = ] r , . . . , b c k k } } . . { ain n symmetric and tations 1 k , lead. 2 cycles n , . . . , trignmesof numbers Stirling to sn ( using ction denoted , e n r aigexactly having ] π Siln numbers -Stirling } aigexactly having ) c 1 c , into n k l -tuple 2 r c , . . . , ,r l, ns, rs.the (resp. n k k )- r bl cl count , -Stirling blocks ( ( π k σ A . = ) = ) k 2 BELBACHIR H. AND DJEMMADA Y.
Dumont, in [8], gives the first interpretation for the ”central factorial” numbers of the second kind U(n, k) given by the recurrence (1.1) U(n, k)= U(n 1,k 1) + k2U(n 1,k), for 0 2. The (l,r)-Stirling numbers of both kinds Let us consider the following generalization, Definition 2.1. n (l) The (l,r)-Stirling number of the first kind k r counts the number of l-tuple of permutations (σ1, σ2,...,σl) of [n] having exactly k cycles such that 1, 2,...,r first elements lead, and cl(σ )= cl(σ )= = cl(σ ). 1 2 ··· l Definition 2.2. n (l) The (l,r)-Stirling number of the second kind k r counts the number of l-tuple of partitions (π1, π2,...,πl) of [n] having exactly k blocks such that 1, 2,...,r first elements lead, and bl(π )= bl(π )= = bl(π ). 1 2 ··· l Theorem 2.3. The (l,r)-Stirling numbers of the first satisfy the following recurrences n (l) n 1 (l) n 1 (l) (2.1) = − +(n 1)l − , for n>r k k 1 − k r − r r and n (l) 1 n (l) n (l) (2.2) = , for n r > 1. k (r 1)l k 1 − k 1 ≥ r − − r−1 − r ! with boundary conditions n (l) for ; and n (l) , for . k r =0, n Proof. The (σ1, σ2,...,σl)-permutations of the set [n] having k cycles such that 1, 2,...,r first elements are in distinct cycles and cl(σ1) = cl(σ2) = = cl(σl) is either obtained from: ··· THE (l,r)-STIRLING NUMBERS: A COMBINATORIAL APPROACH. 3 Inserting the nth elements after any element in each permutation of (σ1, σ2,...,σl)- • permutations of the set [n 1] having k cycles such that 1, 2,...,r first elements are − l n−1 (l) in distinct cycles and cl(σ1)= cl(σ2)= = cl(σl), hence there are (n 1) k r choices. ··· − The nth element forms a cycle in each permutation of (σ , σ ,...,σ )-permutations, • 1 2 l the remaining [n 1] have to be (σ1, σ2,...,σl)-permuted in (k 1) cycles under − (l) − the preceding conditions, hence there are n−1 . k−1 r This correspondence yields the first recurrence. For the second recurrence, we use the double counting principle. Let us count the numbers of (σ , σ ,...,σ )-permutations of the set [n] having (k 1) cycles such that 1 2 l − 1,...,r 1 are cycle leaders but r is not, with cl(σ1)= cl(σ2)= = cl(σl), this is either obtained− from: ··· We count the (σ , σ ,...,σ )-permutations of the set [n] having (k 1) cycles such • 1 2 l − that 1,...,r 1 are cycle leaders then we exclude from them the (σ1, σ2,...,σl)- permutations− having r as cycle leader. That gives n (l) n (l) , k 1 − k 1 − r−1 − r Or we count the (σ1, σ2,...,σl)-permutations of the set [n] having k cycles such • that 1,...,r are cycle leaders then we appending the cycle having r as leader at the end of a cycle having a smaller leader. We have (r 1) choices to do in each permutation. That gives − n (l) (r 1)l , − k r from the two ways of counting we get the result. Theorem 2.4. The (l,r)-Stirling numbers of the second satisfy the following recurrences n (l) n 1 (l) n 1 (l) (2.3) = − + kl − , for n>r k k 1 k r − r r and n (l) n (l) n 1 (l) (2.4) = (r 1)l − , for n r> 1. k k − − k ≥ r r−1 r−1 with boundary conditions n (l) , for ; and n (l) , for . k r =0 n Proof. As in Theorem 2.3, the (π1, π2,...,πl)-partitions of the set [n] into k blocks such that 1, 2,...,r first elements are in distinct blocks and bl(π1) = bl(π2) = = bl(πl) is either obtained from: ··· Inserting the nth elements in a block of each partition of (π1, π2,...,πl)-partitions • of the set [n 1] into k blocks such that 1, 2,...,r first elements are in distinct − 4 BELBACHIR H. AND DJEMMADA Y. (l) blocks and bl(π ) = bl(π ) = = bl(π ), hence there are kl n−1 choices (the 1 2 ··· l k r position of the nth element in a block doesn’t matter). The nth element forms a block in each partition of (π , π ,...,π )-partitions, the • 1 2 l remaining [n 1] have to be (π1, π2,...,πl)-partitioned into (k 1) blocks under − (l) − the preceding conditions, hence there are n−1 . k−1 r For the Identity (2.4), we use the double counting principle to count the numbers of (π1, π2,...,πl)-partitions of [n] into k blocks such that 1, 2,..., (r 1) are block leaders but r is not, with bl(π )= bl(π )= = bl(π ), this is either obtained− from: 1 2 ··· l We count the (π , π ,...,π )-partitions of the set [n] into k blocks such that 1,...,r • 1 2 l − 1 are block leaders then we exclude from them the (π1, π2,...,πl)-partitions having r as block leader, with bl(π )= bl(π )= = bl(π ). That gives 1 2 ··· l n (l) n (l) , k − k r−1 r Or we count the (π1, π2,...,πl)-partitions of the set [n] r into k blocks such • that 1,...,r 1 are block leaders then we include the element\{ } r in any block having a smaller− leader then r. We have (r 1) choices to do in{ each} partition of − (π1, π2,...,πl)-partitions, that gives n 1 (l) (r 1)l − , − k r−1 from the two ways of counting we get the result. Remark 2.5. Using the previous recurrences it is easy to get the following special cases n (l) (2.5) = rl(r + 1)l (n 2)l(n 1)l =(rn−r)l, for n r r ··· − − ≥ r and n (l) (2.6) = rl(n−r), for n r. r ≥ r 3. Orthogonality of (l,r)-Stirling numbers pair Theorem 3.1. For n k 0, for all positive integer l, we have the two orthogonality relations bellow ≥ ≥ n (l) (l) ( 1) δ , for n r; n j n,k (3.1) ( 1)j = − ≥ j k − j r r 0, for n and n (l) (l) ( 1) δ , for n r; j k n,k (3.2) ( 1)j = − ≥ n j − j r r 0, for n < r. X Proof. Let us start by Identity (3.1). The proof goes by induction on n For n 4. Properties via symmetric functions Let x1, x2,...,xn be n random variables. We denote, respectively, by ek(x1, x2,...,xn) and hk(x1, x2,...,xn) the elementary symmetric function and the com- plete homogeneous symmetric function of degree k in n-variables given for n k 1, by ≥ ≥ (4.1) e (x , x ,...,x )= x x k 1 2 n i1 ··· ik 1≤i and n (4.4) H(t)= h (x , x , , x )tk = (1 x t)−1. k 1 2 ··· n − i i=1 Xk≥0 Y For more details about symmetric functions we refer readers to [2, 13, 15] and the references therein. Let us now give some results linked to the symmetric functions and their generating functions. Theorem 4.1. The (l,r)-Stirling of the first kind and the elementary symmetric function are linked as n +1 (l) (4.5) = e (rl,...,nl), n +1 k k − r equivalently n (l) (4.6) = e (rl,..., (n 1)l). k n−k − r Proof. It is clear that in each (σ , σ ,...,σ )-permutation having (n k) cycles with 1 2 l − 1,...,r lead, we have 1, 2,...,r,yr+1,...,yn−k lead a cycle and x1, x2,...,xk el- {ements don’t} lead where r Construct (n k) cycles having only one element from 1, 2,...,r, yr+1,...,yn−k , • i. e. − { } σ = (1)(2) ... (r)(yr+1) ... (yn−k), Insert x after an element of cycles smaller than x , we have (x 1) ways of • 1 1 1 − inserting x1. Then Insert x2 after an element of cycles smaller than x2, we have (x2 1) choices, and so on. We have (x1 1)(x2 1) (xk 1) ways to construct a permutation.− − − ··· − l Repeat the process with each permutation σ σ1,...,σl , so we have (x1 1) (x2 • l l ∈{ } − − 1) (xk 1) ways of construction. Summing··· − over all possible set of numbers x , x ,...,x , hence the total number of • { 1 2 k} ways to construct (σ1, σ2,...,σl)-permutations having (n k) cycles with 1,...,r lead is − { } n (l) = (x 1)l(x 1)l (x 1)l n k 1 − 2 − ··· k − r r Theorem 4.2. The (l,r)-Stirling of the first kind and the complete homogeneous symmetric function are linked as n + k (l) (4.7) = h (rl,...,nl), n k r Proof. Let us count the number of (π1, π2,...,πl)-partitions of [n + k] into n blocks with 1, 2,...,r are leaders. First, we denote, y1,y2,...,yk the elements that are not leaders {where y Construct a partition of n blocks with [n + k] y1,y2,...,yk where 1, 2,...,r are • leaders, i. e. \{ } 1 2 ... r z ... z . { }{ } { }{ r+1} { n} Insert the y ,y ,...,y elements to the n blocks. It is clear that y can belong • { 1 2 k} i only to a block having a leader smaller than yi, we have x1 x2 xk ways to do. · ··· l l l Repeat the process with each partition π π1,...,πl , so we have (x1) (x2) ... (xk) • ways of construction. ∈{ } Summing over all possible set of numbers x , x ,...,x , hence the total number of • { 1 2 k} ways to construct (π1, π2,...,πl)-partitions of [n+k] having n blocks with 1,...,r lead is { } n + k (l) = xl xl xl n 1 2 ··· k r r≤x ≤x ≤···≤n 1 X2 l l = hk(r ,...,n ). 5. Generating functions Now, we can use the symmetric functions to construct the generating functions for the (l,r)-Stirling of both kinds. Theorem 5.1. The generating function for the (l,r)-Stirling numbers of the first kind is n (l) n−1 (5.1) zk = zr z + il = zr z + rl z +(r + 1)l z +(n 1)l , k ··· − k r i=r X Y 8 BELBACHIR H. AND DJEMMADA Y. Proof. From Theorem 4.1 and the generating function (4.3) we obtain (l) n k n l l −1 k z = z ek(r ,..., (n 1) )(z ) k r − Xk Xk n−1 il (5.2) = zn 1+ z i=r nY−1 = zr (z + il). i=r Y Theorem 5.2. The generating function for the (l,r)-Stirling numbers of the second kind is −1 n (l) k zk (5.3) zn = zk (1 zil) = . k − (1 zrl)(1 z(r + 1)l)(1 zkl) r i=r ! Xn=k Y − − − Proof. From Theorem 4.2 and the generating function of homogeneous symmetric function (4.4), we obtain −1 n (l) k + j k zn = zk+j = zk h (rl,...,kl)zj = zk (1 zil) . k k j − r j≥0 j≥0 i=r ! Xn≥k X X Y In the following theorem we investigate the symmetric functions to obtain a convolution formula for the (l,r)-Stirling numbers of both kinds. Theorem 5.3. For all positive integers l, n, k and r with (n k r), we have ≥ ≥ l s n + i (2 ) l−1 n +1 (2 ) n + k (5.4) l = . n n +1 is n l r s=0 r r i0 +2i1 X+2 il = k Y − i , ···, i 0 0 ··· l ≥ THE (l,r)-STIRLING NUMBERS: A COMBINATORIAL APPROACH. 9 Proof. Let us consider the generating function of the complete homogeneous symmetric function (4.4). From that we have n 1 h (x ,...,x )zk = k 1 n (1 x z) i=1 i Xk≥0 Y − n l−1 2s 2s 1 1+ xi z = s s (1 x z) 1+ x2 z2 i=1 i s=0 i Y − Y n l−1 1 2s 2s = l 1+ x z 2 2l i i=1 (1 xi z ) s=0 Y − Y l−1 2l 2l 2lk 2s 2s 2sk = hk(x1 ,...,xn )z ek(x1 ,...,xn )z s=0 Xk≥0 Y Xk≥0 l l l l−1 l−1 l−1 = h (x2 ,...,x2 )z2 k e (x ,...,x )zk e (x2,...,x2 )z2k e (x2 ,...,x2 )z2 k k 1 n k 1 n k 1 n ··· k 1 n Xk≥0 Xk≥0 Xk≥0 Xk≥0 l−1 2l 2l 2s 2s k = hil (x1 ,...,xn ) eis (x1 ,...,xn ) z . k≥0 i +2i + +2li = k; s=0 X 0 1 ···X l Y i ,...,i 0. 0 l ≥ From Theorem 4.1 and Theorem 4.2 and by comparing the coefficients of zk of the two sides the result holds true. The simplest case of the previous theorem is the corollary bellow which generalize the result of Broder [6]. Corollary 5.4. For l =1, we have ⌊k/2⌋ n + i (2) n +1 n + k (5.5) = . n n +1+2i k n i=0 r r r X − 6. The (l,r)-Stirling numbers, the sum powers and Bernoulli polynomials Recall, for every integer n 0, the Bernoulli polynomials, denoted Bn(x), are defined by ≥ ∞ tn text (6.1) B (x) = . n n et 1 n=0 X − The sum of the powers of natural numbers is closely related to the Bernoulli polynomials Bn(x). Jacobi [12, 16] gives the following identity using the sum of powers and Bernoulli polynomials n B (n + 1) B (0) (6.2) jm = m+1 − m+1 . m +1 j=1 X The following theorem gives the relation between (l,r)-Stirling of both kinds and Bernoulli polynomials. 10 BELBACHIR H. AND DJEMMADA Y. Theorem 6.1. For all positive integers n, k and l, we have k n +1 (l) n + k j (l) B (n + 1) B (0) (6.3) ( 1)j(j + 1) − = lk+l+1 − lk+l+1 , − n j n lk + l +1 j=0 X − Proof. In the first hand we have Jacobi’s Identity (6.2) n B (n + 1) B (0) (6.4) (jl)k = lk+1 − lk+1 , lk +1 j=1 X in the second hand, we have n 1 H(t)= h (1l, 2l,...,nl)tk = k (1 jst) j=1 Xk≥0 Y − and n l l l k s E(t)= ek(1 , 2 ,...,n )t = (1 + j t), j=1 Xk≥0 Y from the obvious observation that H(t)=1/E( t), we obtain − d H′(t) (6.5) ln H(t)= = H(t)E′( t) dt H(t) − but d n jl n (6.6) ln H(t)= = js(k+1)tk. dt (1 jlt) j=1 j=1 X − Xk≥0 X Then from equations (6.5) and (6.6), we get n js(k+1)tk = H(t)E′( t) − k≥0 j=1 (6.7) X X l l k k−1 l l k−1 = hk(1 ,...,n )t k( 1) ek(1 ,...,n )t . ! − ! Xk≥0 Xk≥1 Cauchy product and equating coefficient of tk gives n n s(k+1) j l l l l (6.8) j = (j + 1)( 1) ej+1(1 ,...,n )hk−j(1 ,...,n ), − j=1 j≥1 X X replacing symmetric functions by stirling numbers from Theorem 4.1 and Theorem 4.2, and comparing with Equation (6.4) we get the result. THE (l,r)-STIRLING NUMBERS: A COMBINATORIAL APPROACH. 11 7. Multiple zeta function and (l,r)-Stirling numbers of the first kind For any ordered sequence of positive integers i1, i2,...,ik, the multiple zeta function is introduced by Hoffman [11] and independently Zagier [17] by the following infinite sums 1 (7.1) ζ(i1, i2,...,ik)= . ji1 ji2 jik 0 Proposition 7.2. For r =1, we have 1 n +1 (l) (7.7) lim = ζ( l k), n→∞ l k +1 { } (n!) where l =(l,l,...,l). { }n n times Proof. The proposition| {z } follows immediately from the definition of multiple zeta function (7.1) as an infinity sums and Theorem 7.1 for r = 1. Corollary 7.3. For k 1, we have ≥ For l =2 • 1 n +1 (2) π2k (7.8) lim = . n→∞ n 2 k +1 (2k + 1)! ( !) For l =4 • 1 n +1 (4) 4(2π)4k 1 2k+1 (7.9) lim = . n→∞ 4 k +1 (4k + 2)! 2 (n!) For l =6 • 1 n +1 (6) 6(2π)6k (7.10) lim = . n→∞ 6 k +1 (6k + 3)! (n!) For l =8 • 1 n +1 (8) π8k 1 4k+2 1 4k+2 (7.11) lim = 28k+3 1+ + 1 . n→∞ 8 k +1 (8k + 4)! √ − √ (n!) 2 2 ! Proof. Authors in [7] give the following special values of multiple zeta function π2n ζ( 2 ) = , { }n (2n + 1)! 4(2π)4n 1 2n+1 ζ( 4 ) = , { }n (4n + 2)! 2 6(2π)6n ζ( 6 ) = , { }n (6n + 3)! π8n 1 4n+2 1 4n+2 ζ( 8 ) = 28n+3 1+ + 1 , { }n (8n + 4)! √ − √ 2 2 ! the corollary is a consequence of the previous special cases and Proposition 7.2. THE (l,r)-STIRLING NUMBERS: A COMBINATORIAL APPROACH. 13 8. Remarks The (l,r)-Stirling gives another graphical view of Rooks polynomials of higher di- • mensions in triangle boards [3, 18] using set partitions. In this work we gives a limit representation of multiple zeta function using (l,r)- • Stirling numbers. π2 We can obtain the well-known Euler identity ζ(2) = 6 from Equation (7.8) for • k = 1. Acknowledgment. We would like to thank the anonymous reviewers for their suggestions and comments which improved the quality of the present paper. The paper was partially supported by the DGRSDT grant №:C0656701. References [1] Ablinger, J., and Bl¨umlein, J. (2013). Harmonic sums, polylogarithms, special numbers, and their generalizations. In Computer Algebra in Quantum Field Theory (pp. 1-32). Springer, Vienna. [2] Abramowitz, M., and Stegun, I. A. (1964). Handbook of mathematical functions with formulas, graphs, and mathematical tables (Vol. 55). US Government printing office. [3] Alayont, F., and Krzywonos, N. (2013). Rook polynomials in three and higher dimensions. Involve, a Journal of Mathematics, 6(1), 35-52. [4] Arakawa, T., and Kaneko, M. (1999). Multiple zeta values, poly-Bernoulli numbers, and related zeta functions. Nagoya Mathematical Journal, 153, 189-209. [5] Bradley, D. M. (2005). Partition identities for the multiple zeta function. In Zeta functions, topology and quantum physics (pp. 19-29). Springer, Boston, MA. [6] Broder, A. Z. (1984). The r-Stirling numbers. Discrete Mathematics, 49(3), 241-259. [7] Borwein, J. M., and Bradley, D. M. (1997). Evaluations of k-fold Euler/Zagier sums: a compendium of results for arbitrary k. the electronic journal of combinatorics, 4(2), R5. [8] Dumont, D. (1974). Interpr´etations combinatoires des nombres de Genocchi. Duke Mathematical Jour- nal, 41(2), 305-318. [9] Foata, D., and Han, G. N. (2000). Principes de combinatoire classique. Lecture notes, Strasbourg. [10] Gelineau, Y., and Zeng, J. (2010). Combinatorial interpretations of the Jacobi-Stirling numbers. the electronic journal of combinatorics, 17(R70), 1. [11] Hoffman, M. (1992). Multiple harmonic series. Pacific Journal of Mathematics, 152(2), 275-290. [12] Jacobi, C. G. J. (1834). De usu legitimo formulae summatoriae Maclaurinianae. Journal f¨ur die reine und angewandte Mathematik, 1834(12), 263-272. [13] Macdonald, I. G. (1998). Symmetric functions and Hall polynomials. Oxford university press. [14] Merca, M., and Cuza, A. I. (2012). A special case of the generalized Girard-Waring formula. Journal of Integer Sequences, 15(2), 3. [15] Merca, M. (2016). New convolutions for complete and elementary symmetric functions. Integral Trans- forms and Special Functions, 27(12), 965-973. [16] Srivastava, H. M. (2000, July). Some formulas for the Bernoulli and Euler polynomials at rational arguments. In Mathematical Proceedings of the Cambridge Philosophical Society (Vol. 129, No. 1, pp. 77-84). [17] Zagier, D. (1994). Values of zeta functions and their applications. In First European Congress of Mathematics Paris, July 6–10, 1992 (pp. 497-512). Birkh¨auser Basel. [18] Zindle, B. (2007). Rook polynomials for chessboards of two and three dimensions. Master thesis. [19] Zudilin, W. W. (2003). Algebraic relations for multiple zeta values. Russian Mathematical Surveys, 58(1), 1-29. 14 BELBACHIR H. AND DJEMMADA Y. USTHB, Faculty of Mathematics, RECITS Laboratory, P.O. Box 32, El Alia, 16111, Bab Ezzouar, Algiers, Algeria Email address: [email protected], [email protected] USTHB, Faculty of Mathematics, RECITS Laboratory, P.O. Box 32, El Alia, 16111, Bab Ezzouar, Algiers, Algeria Email address: [email protected], [email protected]