Arxiv:1702.06519V1 [Math.CO] 21 Feb 2017 Ho N Ug[]Soe Httenumbers the That Showed [8] Jung and Cheon Where Ie Sflos Let Follows

arXiv:1702.06519v1 [math.CO] 21 Feb 2017 ho n ug[]soe httenumbers the that showed [8] Jung and Cheon where ie sflos Let follows. as tices feeet fcorank of elements of order lmnsof elements hte ubr ftescn kind, second the of numbers Whitney hte ubr ftescn idb sn ooe e partitions. set colored using by kind second the of numbers Whitney iob n amn 2]fuda neetn obntra interp combinatorial interesting an found [27] Rahmani and Mihoubi The oncigcecet ewe oepriua oyoil.Frno For polynomials. particular k some n, between coefficients connecting e od n phrases. and words Key 2010 Date E PRAHT THE TO APPROACH NEW A and r m Winynmeso h eodkind second the of numbers -Whitney eray2,2017. 22, February : x ahmtc ujc Classification. Subject Mathematics Abstract. fthe of ieeta prtr soitdt h grammar the to associated operators differential eodkn.Aan yseilzn otecase the to specializing by Again, kind. second m ieeta acls ial epoeta the that prove we Finally calculus. differential fteSiln ubro h eodkn n hog hmabinomial a them through and sever show kind the we second involving Moreover, Touchard’s. the generalizes of that number polynomials Stirling the of eaiet h eeaie ocadbnma aiy td hi umb their study family, binomial Touchard [ generalized introduce the to relative hte ubro ohkns enul n ue polynomials. Euler and Bernoulli kinds, both of number Whitney n h ira arcsw iesvrlnwcmiaoilidentitie combinatorial new several give we matrices Riordan the and h offiin of coefficient The . n r eoti loanwcmiaoilitrrtto fthe of interpretation combinatorial new a also obtain we 1 = = with Q x r ( n Winynmeso h eodkn bandfo h combinator the from obtained kind second the of numbers -Whitney x OBNTRA IFRNILCALCULUS DIFFERENTIAL COMBINATORIAL ( G n m − containing ) -trignmeso h rtkn.Fo h eainbtenumbra between relation the From kind. first the of numbers ]-Stirling ≥ ntepeetatcew nrdc w e obntra interpre combinatorial new two introduce we article present the In 1) r Dwigplnmasadthe and polynomials -Dowling k Q · · · k ≥ n IULA M A. MIGUEL ieeta operators, Differential ( ( in G n o n integer any for and 0 x eteDwigltieo rank of lattice Dowling the be ) − Q r ( s n n mx ( ftepolynomial the of n G )for 1) + − [10]. ) + s r 1. ) itntui ucin.Ti eunegnrlzsthe generalizes sequence This functions. unit distinct NE N JOS AND ENDEZ ´ n W rmr 18;Scnay1B3 51,05A19. 05A15, 11B73, Secondary 11B83; Primary Introduction = n m X ≥ ( k r ,k n, =0 n r WINYNMESB USING BY NUMBERS -WHITNEY Winynumber, -Whitney and 1 1 W m r = ) Winynmesuigtecombinatorial the using numbers -Whitney > m r k m,r W DwigplnmasaeaSee family Sheffer a are polynomials -Dowling P r W G m,r eitoueanwgeneralization new a introduce we 0 = ( W := m,r x k n ,k n, =0 0 0 m, ( .RAM L. E ´ { ( ,k n, 1. = W y ,k n, 1 eedfie yMzo[4 sthe as Mez˝o [24] by defined were ) ( → ,k n, m,r r eae oteDwiglat- Dowling the to related are ) ) yx x ( n ,ta on h oa number total the count that ), ,k n, k r m where , , Dwigpolynomial. -Dowling r ´ x , IREZ Siln ubr fthe of numbers -Stirling lwl-nw identities well-known al seult h ubrof number the to equal is ) → x novn the involving s } ealta partition a that Recall a ness and inverses, ral G yspecializing By . yefml of family type -eaieintegers n-negative eainfrthe for retation safiiegopof group finite a is c fthe of ics calculus l tations r - r - 2 MIGUEL A. MENDEZ´ AND JOSE´ L. RAMÍREZ of a set A is a class of disjoint subsets of A such that the union of them covers A. The subsets are called blocks. Let r, n ≥ 0 be integers, and let An,r be the set defined by

An,r := {1, 2,...,r,r +1,...,n + r}. The elements 1, 2,...,r will be called special elements by us. A block of a partition of the above set is called special if it contains special elements. Then Wm,r(n, k) counts the number of the set partitions of An,r in k + r blocks such that the elements 1, 2,...,r are in distinct blocks (i.e., any special block contains exactly one special element). All the elements but the last one in non-special blocks are coloured with one of m colours independently and neither the elements in the special blocks nor the special blocks are coloured. The r-Whitney numbers satisfy the following recurrence relation [24]

Wm,r(n, k)= Wm,r(n − 1,k − 1)+(km + r)Wm,r(n − 1,k). Moreover, these numbers have the following exponential generating function [24]: ∞ zn erz emz − 1 k (1) W (n, k) = . m,r n! k! m Xn=k Note that if (m, r)=(1, 0) we obtain the Stirling numbers of the second kind, if (m, r)= (1,r) we have the r-Stirling (or noncentral Stirling) numbers [4], and if (m, r)=(m, 1) we have the Whitney numbers [1, 2]. Many properties of the r-Whitney numbers and their connections to elementary symmetric functions, matrix theory and special polynomials, combinatorial identities and generalizations can be found in [9, 18, 20, 21, 24, 25, 26, 28, 36].

The combinatorial differential calculus was introduced by Joyal in the framework of combinatorial species (see for example [13, 3, 14]). However, by using directly exponential formal power series and the combinatorics of their coefficients (see [19]), the fundamentals of the approach can be explained without the use of the categorical framework involved in the theory of species. The combinatorial differential calculus is closely related to Chen context-free grammar method [6]. However, they are not equivalent. Differential operators that are not derivations can be combinatorially interpreted, and algebraic formulas obtained from that interpretation. Only derivations have a counterpart in the context of Chen grammars. The substitution rules in a context-free grammar can be translated into a differential operator. By the iterated application of the combinatorial version of it, the associated combinatorial objects emerge, and not infrequently, in a simple and natural way.

Using the combinatorial differential calculus we present a new combinatorial interpretation to the r-Whitney numbers of the second kind. As a special case we get a new combinatorial interpretation for the r-Stirling numbers of the second kind and introduce the [m]-Stirling numbers of the second kind, S[m](n, k). Through them we define the [m]-Touchard polynomials, that generalize the classical family and are also of binomial type. Their umbral inverse gives us the [m]-Stirling numbers of the first kind. In a forthcoming paper we will A NEW APPROACH TO THE r-WHITNEY NUMBERS 3 show that the combinatorial interpretation of the r-Whitney numbers (and of the generalization of the Stirling numbers) of the second kind we present here are very natural. They directly count the generalized Whitney numbers of ranked posets, and by Möbius inversion we get the Whitney numbers of the first kind.

By using classical results of umbral calculus we get that the r-Dowling polynomials are of Sheﬀer type relative to the [m]-Touchard binomial family. Finally, from the relation between umbral calculus and the Riordan matrices we give several new combinatorial identities involving the r-Whitney number of both kinds, Bernoulli polynomials and Euler polynomials.

2. New combinatorial models for r-Whitney, r-Stirling, and new [m]-Stirling numbers Chen [6] introduced a combinatorial method by means of context-free grammar to study exponential structures. A context-free grammar G over an alphabet X, whose symbols are commutative indeterminates, is a set of productions or substitutions rules that replace a symbol of X by a formal function (formal power series) in the set of indeterminates X,

x → φx(X), x ∈ A. Here A is a subset of the alphabet X. The formal derivative D is a linear operator deﬁned with respect to a context-free grammar G such that for a letter x ∈ A acts by substitution by φ(X), and it is extended recursively as a derivation. For any formal functions u and v we have: ∂f(u) D(u + v)= D(u)+ D(v),D(uv)= D(u)v + uD(v),D(f(u)) = D(u), ∂u where f(x) is a formal power series. For more applications of the context-free grammar method see for example [5, 7, 11, 15, 16, 17]. The grammar formal derivative D can be equivalently written as a derivation in the algebra of formal power series over a ring R containing Q, R[[X]].

D = φx(X)∂x. x∈A X For A infinite, we have to make a summability assumption over the family of formal power series {φx(X)}x∈A. The variables in X are thought of as colors. Each operator φx(X)∂x is combinatorially interpreted as a corolla having as root a ghost vertex of color x and weighted with the coefficients of the series φ(X) according with the distribution of colors of the leaves in the corolla. A formal power series F (X) are also represented by corollas in a similar way but without the ghost vertex. The combinatorial representation of the operator acting over a series F (X) is obtained by dropping over the combinatorial representation of F (X) corollas whose ghost vertices replace vertices of the same color on corollas representing F (see Figure 1). Details and proofs of why this combinatorial approach works can be seen in [19]. 4 MIGUEL A. MENDEZ´ AND JOSE´ L. RAMÍREZ

2 2 1 φ * i 2 2 2 2 * i 2 1 3 3 2 φ 0 0 F F

Figure 1. Corolla operator φ(X)∂i applied to F.

As an example let us consider the Stirling grammar ([6]) y → xy G := (x → x. Applying the associated formal operator n times to y it is known that we get n (2) Dny = S(n, k)xky, Xk=1 S(n, k) being the Stirling numbers of the second kind. The diﬀerential operator associated to that grammar is D = xy∂y + x∂x. The variables x and y are represented by vertices in two colors, yellow and white. The operator D has as combinatorial representation the sum of two corollas. The operator xy∂y acts over a white vertex replacing it by a ‘ghost’ vertex ∗ and connecting it to a yellow vertex x and to white one y. The operator x∂x places a ghost vertex in a yellow vertex and connects it to another yellow vertex (see Figure 2).

x y x

∗1 ∗1 xy∂y x∂x

Figure 2. Combinatorial operator D = xy∂y + x∂x.

When applying the operator more than one time, for simplicity, we replace the ghost vertices by numbers to keep track of the order in which the operator was applied. If we apply it n times to y (combinatorially represented as a singleton vertex of color white) we obtain an increasing tree with a path (spine) from the root to the white vertex on the A NEW APPROACH TO THE r-WHITNEY NUMBERS 5 top. Along each node of the spine there are linear branches having in the top of each of them one yellow vertex with weight y. The elements of the branches along the spine form a partition π of [n]. The weight of each such tree is equal to x to the number of branches (blocks of π) times y, the weight of the white vertex on top of the spine (see Figure 3). In this way we get a visual proof of Equation (2).

x y x 6 7 4 x 5 2 3 1

7 7 Figure 3. A tree enumerate by D y =(xy∂y + x∂x) y that corresponds to the partition π = {{1, 3}{2, 5, 7}{4, 6}}.

Let us now consider the context-free grammar G deﬁned in [11]. y → yxm, G := (x → x. Hao et al. [11] proved as a particular case that n n mk+1 D yx = Wm(n, k)yx . Xk=0 This grammar G corresponds to the diﬀerential operator: m D = yx ∂y + x∂x. In Figure 4 we show the combinatorial interpretation of this operator for m = 2. Moreover, by the main theorem of [11] we have n n r mk+r (3) D yx = Wm,r(n, k)yx . Xk=0

∗ ∗

2 yx ∂y x∂x

m Figure 4. Combinatorial interpretation for the operator D = yx ∂y + x∂x, for m = 2. 6 MIGUEL A. MENDEZ´ AND JOSE´ L. RAM´IREZ

In Figures 5 and 6, we show the ﬁrst and second derivative for m = 2 and r = 3, respectively.

Dyx3 = yx5 + 3yx3

1 1 1 1

Figure 5. The ﬁrst derivative.

D2y3x = yx7 + 8yx5 + 9yx3

2 1 2 1

1 2 1

2 2

1 1

1 2 1 2 1 2

1 1 2 1 2

1 2 1 1 2

1 2 1 2 1

Figure 6. The second derivative. A NEW APPROACH TO THE r-WHITNEY NUMBERS 7

By the iterated application of the combinatorial version of the operator D, dropping the corollas over one white vertex and r yellow vertices (yxr) we obtain the following combinatorial structure. A forest consisting of one increasing tree, grown from the initial white vertex y, followed by r linearly ordered branchless increasing trees (grown from the r yellow initial vertices, xr). The trees of this forest are characterized as follows (see Figure 7). (1) The ﬁrst tree has a spine of white vertices with an unlabeled white vertex at the top having weight y. There are m (totally ordered) branches of yellow vertices that sprout from each vertex on the spine. Each branch has an unlabeled yellow vertex at the top, with weight x. (2) Each of the r branchless increasing trees consists of a set of internal vertices (that may be empty) with an unlabeled yellow vertex at the top, having weight x. The total weight of each of this forests is equal to yxmk+r, where k is the number of vertices in the spine. Then, from Equation (3) we have that the r-Whitney numbers Wm,r(n, k) count the number of forests as above, having k vertices in the spine, and n internal vertices in total. Reinterpreting the forests of increasing trees we obtain the following combinatorial interpretation of Wm,r(n, k).

Theorem 1. The r-Whitney numbers Wm,r(n, k) count pairs of the form (({(B, fB}B∈π), B) where (1) The ﬁrst component is an m-colored partial partition of {1, 2,...,n}, having exactly k blocks, ⊎B∈πB = A ⊆{1, 2,...,n}. The function fB : B → [m] has the following restriction. It colors the minimum element of the block with color 1. (2) The second component B =(B1, B2,...,Br) is a weak r-composition of [n] − A r

Br = [n] − A. j=1 ] Proof. We are going to prove that the ﬁrst tree is in bijection with the colored partitions described in Item (1). We assign to each vertex of the spine of a given tree the set of vertices of the m branches attached to it plus the vertex itself. In this way we obtain the partial partition π with k blocks in total. Then color the vertices on the ith branch of each block with color i, i = 1, 2,...,m, and assign color 1 to the vertex on the spine. Observe that, since the tree is increasing, the vertex on the spine has the minimum label of its block. This construction is clearly reversible (see Figure 7). The r branchless trees are naturally associated to the composition B by assigning the ith tree its set of internal vertices Bi.

We represent a colored set (B, fB) by placing the color of each element as its exponent f(b) (B, fB) ≡{b |b ∈ B}. In this way, the colored partition associated to the tree in Figure 7 is {{11, 81, 132}, {21, 31, 52, 72}, {61, 91}}. The composition is equal to: B =({4, 10}, ∅, {11, 12}). 8 MIGUEL A. MENDEZ´ AND JOSE´ L. RAM´IREZ

π x x B x y 7 9 x 2 1 5 6 x 1 x x 3 2 1 2 13 1 1013 12 x 2 x 8 1 4 11 1 1

Figure 7. Whitney forest for m = 2 and r = 3.

Example 2. The r-Whitney number W2,2(2, 2) = 1, the pair being

({{11}, {21}}, (∅, ∅)).

For W2,3(2, 1) = 8, it enumerates pairs, the ﬁrst component is a 2-colored partition on a subset A of {1, 2} having one block. The second component a weak 3-composition of [2] − A. The pairs being (1) ({{11, 21}}, (∅, ∅, ∅)), (2) ({{11, 22}}, (∅, ∅, ∅)), (3) ({{11}}, ({2}, ∅, ∅)), (4) ({{11}}, (∅, {2}, ∅)), (5) ({{11}}, (∅, ∅, {2})), (6) ({{21}}, ({1}, ∅, ∅)), (7) ({{21}}, (∅, {1}, ∅)), (8) ({{21}}, (∅, ∅, {1})).

Remark 3. We get the following combinatorial interpretations for generalized Stirling numbers of the second kind obtained by specializing r and m. (1) For m = 1 we get a combinatorial interpretation for the r-Stirling numbers of the second kind Sr(n, k). They count the pairs (π, B), π a partition of some subset A of [n], B a weak r-composition of [n] − A. (2) For r = 0 we get a new generalization S[m](n, k) of the Stirling numbers of the second kind. It counts the number of colored partitions, as in Theorem 1, Item (1), but over the whole set [n]. A NEW APPROACH TO THE r-WHITNEY NUMBERS 9

[m] Deﬁnition 4. We deﬁne the [m]-Touchard polynomials, Tn (x), by n [m] [m] k Tn (x) := S (n, k)x . Xk=1 [m] ∞ The polynomial family {Tn (x)}n=0 is of binomial type, i.e., [m] T0 (x)=1, n n T [m](x + y)= T [m](x)T [m] (y). n k k n−k Xk=0 We shall prove it in Section 4.

3. The r-Dowling Polynomials Cheon and Jung [8] defined the r-Dowling polynomials of degree n by n k Dm,r(n, u) := Wm,r(n, k)u . Xk=0 They found some combinatorial identities by means of Riordan arrays. Let us define the following generating function ∞ ∞ tn tn n H (t)= H(t; x, y)= etD yxr = D nyxr = W (n, k)yxkm+r. m,r n! n! m,r n=0 n=0 X X Xk=0 It is easy to show the following Lemma 5. We have the identities n r r m D yx = yx Dm,r(n, x ), ∞ tn H (t) = etD yxr = yxr D (n, xm). m,r n! m,r n=0 X Theorem 6. The exponential generating function for the r-Dowling polynomials is ∞ tn emt − 1 D (n, u) = exp rt + u . m,r n! m n=0 X Proof. Since D is a derivation, it is easy to show that the operator etD is multiplicative. Hence ′ tD D r tD m+r r tD r tD m tD r Hm,r(t)= e yx = e (yx + ryx )=(e yx ) · (e x) + re yx m mt = Hm,r(t)(x e + r). From that, ′ Hm,r(t) d m mt = ln(Hm,r(t)) = x e + r. Hm,r(t) dt 10 MIGUEL A. MENDEZ´ AND JOSE´ L. RAMÍREZ

r Integrating and using the initial condition Hm,r(0) = yx , we get: xm H (t)= yxr exp{rt + (emt − 1)}. m,r m By Lemma 5, making u = xm we get the result. Theorem 7. The r-Dowling polynomials satisfy the following relation for any integers r, l ≥ 0 n n D (n, u)= ln−kD (k,u). m,r+l k m,r Xk=0 n km+r Proof. Let Dm,r(n, x, y) := k=0 Wm,r(n, k)yx . Then tD r+l PtD r tD l Hm,r+l(t)=b e yx = e yx e x ∞ n tn n = xleltH (t)= xl ln−jD (j, x, y) . m,r n! j m,r n=0 j=0 ! X X Therefore b n n D (n, x, y)= xl ln−jD (j, x, y). m,r j m,r j=0 X Modifying a bit this equalityb we obtain the desired result.b In particular if l = 1 then n n (4) D (n, u)= D (k,u). m,r+1 k m,r Xk=0 In the following theorem we generalize the beautiful relation given por Spivey [33] for the Bell numbers Bn. The Spivey’s formula says that: n m n B = jn−kS(m, j) B , n+m k k j=0 Xk=0 X where S(n, j) is the Stirling number of the second kind. We generalize this identity for the r-Whitney numbers by using diﬀerential operators. Theorem 8. The r-Dowling polynomials satisfy the following formula n h n D (n + h, u)= D (k,u)W (h, j)ujjn−kmn−k. m,r k m,r m,r j=0 Xk=0 X Proof. We compute the derivative in two ways. First, in a direct way ∞ dh tn H (t)= yxr D (n + h, xm) , dth m,r m,r n! n=0 X A NEW APPROACH TO THE r-WHITNEY NUMBERS 11

dh tD tD D h and secondly, by using the identity dth e = e and Lemma 5 dh dh H (t)= etD yxr = etD D hyxr = etD yxrD (h, xm) dth m,r dth m,r ∞ tn (5) = yxr D (n, xm) D (h, xmemt) . m,r n! m,r n=0 ! X Making the change u = xm and from the generating functions in Equation (3) and in Equation (5) we get ∞ ∞ tn tn (6) D (n + h, u) = D (n, u) D (h, uemt) . m,r n! m,r n! m,r n=0 n=0 ! X X mt Expanding Dm,r(h, ue ),

h mt j mjt Dm,r(h, ue )= Wm,r(h, j)u e j=0 X ∞ ∞ h tk h tk = W (h, j)uj mkjk = W (h, j)mkujjk . m,r k! m,r k! j=0 j=0 ! X Xk=0 Xk=0 X By plugging it into Equation (6), computing the Cauchy product and equating coefficients, we obtain the result. The above identity was proved by using a different approach in [36]. Moreover, this identity is a particular case of the main result of [35]. From Theorem 8 we obtain the following convolution formula. Corollary 9. For 0 ≤ s ≤ n + h, we have n h n W (n + h, s)= W (h, j)W (k,s − j)(jm)n−k. m,r k m,r m,r j=0 Xk=0 X In particular if r = 1 we obtain the Theorem 4.3 of [11]. By setting h = 1 in Theorem 8 and Corollary 9 we obtain the following corollary. Corollary 10. The r-Dowling polynomials satisfy the following recursive formula n n D (n +1,u)= rD (n, u)+ u mn−jD (j, u). m,r m,r j m,r j=0 X The r-Whitney numbers of the second kind satisfy the following recursive formula n n W (n +1,k)= rW (n, k)+ mn−jW (j, k − 1). m,r m,r j m,r − j=Xk 1 The above corollary was proved in [8] by using Riordan arrays. 12 MIGUEL A. MENDEZ´ AND JOSE´ L. RAMÍREZ

Theorem 11. The r-Dowling polynomials satisfy the following formula n n D (n, u)= (r − s)n−jD (j, u). m,r j m,s j=0 X Proof. From Lemma 5 we have

tD r tD s tD r−s r−s (r−s)t (7) Hm,r(t)= e yx = e yx e x = x e Hm,s(t). By the Cauchy product we obtain n n D (n, x, y)= xr−s (r − s)n−jD (j, x, y). m,r j m,s j=0 X Therefore, we get theb desired result. b

In particular if s = 1 we have (see Theorem 5.2 of [8])

n n D (n, x)= (r − 1)n−j D (j, x), m,r j m j=0 X where Dm(n, u) are the Dowling polynomials, i.e.,

n k Dm(n, u)= Wm(n, k)u . Xk=0 Corollary 12. The r-Whitney numbers satisfy the following formula n n (8) W (n, k)= (r − s)n−jW (j, k). m,r j m,s j=0 X In particular if s =1 we have n n (9) W (n, k)= (r − 1)n−jW (j, k). m,r j m j=0 X Proof. From Theorem 11

j n r n W (n, k)yxmk+r = xr−s (r − s)n−j W (j, k)yxkm+s m,r j m,s j=0 Xk=0 X Xk=0 n n n = (r − s)n−jW (j, k)yxkm+r. j m,s j=0 Xk=0 X By equating coeﬃcients, we obtain the result. A NEW APPROACH TO THE r-WHITNEY NUMBERS 13

4. The Sheffer family of [m]-Touchard polynomials and r-Dowling polynomials From the generating function of the r-Dowling polynomials (Theorem 6), making r = 0, we get the generating function of the [m]-Touchard polynomials ∞ tn emt − 1 (10) T [m](x) = exp x . n n! m n=0 X From this is immediate that they are of binomial type ([31]). From Equation (4) we have the following identity relating them with the classical Dowling polynomials n n D (n, x)= D (n, x)= T [m](x). m,1 m k k Xk=0 and hence, n n T [m](x)= (−1)n−kD (k, x). n k m Xk=0 Our goal is to ﬁnd the umbral inverse of [m]-Touchard sequence, i.e, the polynomial family [m] n k Tn (x)= k=1 cn,kx satisfying P n b [m] [m] [m] [m] n Tn (T (x)) := S (n, k)Tk (x)= x , ∀n =0, 1, 2 .... Xk=1 b n b [m] [m] [m] n Tn (T (x)) := cn,kTk (x)= x , ∀n =0, 1, 2 .... Xk=1 Let b emt − 1 O(t) := , m and consider its compositional inverse 1 O(t) = ln(1 + mx) m .

Let ∂x the derivative operator acting on the polynomial ring C[x]. Let O(∂x) and O(∂x) be the shift-invariant operators acting also on the polynomial ring C[x], ∞ em∂x − I 1 ∂k Em − I O(∂ )= = mk x = , x m m k! m k=1 X ∞ ∞ 1 k k 1 k k ∂x k k−1 ∂x O(∂ ) = ln(1 + m∂ ) m = (−1) m (k − 1)! = (−1) m (k − 1)! . x x m k! k! Xk=1 Xk=1 Here we denote by Ea is the shift operator Eap(x) = p(x + a). By the classical theory [m] [m] of umbral calculus, Tn (x) and Tn (x) are respectively the sequences associated to the

b 14 MIGUEL A. MENDEZ´ AND JOSE´ L. RAM´IREZ

m− 1 E I m operators m , and ln(1 + m∂x) , Em − I T [m](x + m) − T [m](x) T [m](x)= n n = nT [m] (x) m n m n−1 1 [m] b m b [m] b ln(1 + m∂x) Tn (x) = nTbn−1(x). The derivatives of the formal power series O(x) and O(x) are respectively O′(x) = emx O ′ 1 and ( (x)) = 1+mx . The recurrence formula for families of binomial type ([31] Corollary 1 of Theorem 8), gives us

[m] −m [m] [m] (11) Tn (x) = xE Tn−1(x)= xTn−1(x − m), [m] [m] [m] [m] (12) T (x) = x(1 + m∂x)T (x)= xT (x)+ mx∂xT (x). bn b n−1 b n−1 n−1 From Equation (11) we obtain the polynomial family

[m] (13) Tn (x)= x(x − m)(x − 2m) ··· (x − (n − 1)m). From Equation (11) we get the recurrence for [m]-Stirling numbers of the second kind b (14) S[m](n, k)= S[m](n − 1,k − 1) + kmS[m](n − 1,k).

[m] We define the [m]-Stirling numbers of the first kind as the coefficients connecting Tn (x) with the powers. b Definition 13. We define the [m]-Stirling numbers of the first kind s[m](n, k) as the coef- [m] ficients in the expansion of Tn (x)= x(x − m)(x − 2m) ··· (x − (n − 1)m) in terms of the power sequence, b n x(x − m)(x − 2m) ··· (x − (n − 1)m)= s[m](n, k)xk. Xk=1 Remark 14. From the generating function in Theorem 6, ∞ tn emt − 1 D (n, x) = erx exp x , m,r n! m n=0 X we get that the r-Dowling polynomials,

[m,r] Dn (x) := Dm,r(n, x) are a Sheﬀer family relative to the [m]-Touchard, associated to the pair of generating 1 1/m functions (1+mx)r/m , ln(1 + mx) , 1 emx − 1 −1 , ln(1 + mx)1/m = erx, , (1 + mx)r/m m A NEW APPROACH TO THE r-WHITNEY NUMBERS 15

rx emx−1 being the inverse of e , m as an exponential Riordan array, Equation (18) bellow. For that, see [31], and [30] Theorem 2.3.4. Hence, we have the binomial identity (see [30], Theorem 2.3.9) n n (15) D[m,r](x + y)= D[m,r](x)T [m] (y). n k k n−k Xk=0 [m,r] [m] Its umbral inverse Dn (x) is Sheffer relative to Tn (x). It is associated to the pair emx − 1 erx, . b m b From that we get next theorem Theorem 15. The umbral inverse of the r-Dowling Sheffer sequence is equal to [m,r] −r [m] (16) Dn (x)= E Tn (x)=(x − r)(x − r − m) ··· (x − r − (k − 1)m). The r-Dowling sequence satisfies the identity b b n r(r − m)(r − 2m) ··· (r − (k − 1)m) (17) D[m,r](x)= ∂kT [m](x). n k! x n Xk=0 [m,r] rx emx−1 Proof. Since Dn (x) is associated to (e , m ), it is equal to the inverse of the operator r∂x r [m] e = E applied to Tn (x) (see [30], Theorem 2.3.6). For the same reason we have that b [m,r] r/m [m] D (x)=(1+ m∂x) T (x). b n n r/m Expanding the operator (1 + m∂x) by the binomial formula we obtain the result.

5. Some Applications from Riordan Arrays The r-Whitney numbers can be defined using exponential Riordan arrays. We recall that an infinite lower triangular matrix is called a Riordan array [32] if its kth column satisfies the generating function g(z)(f(z))k for k ≥ 0, where g(z) and f(z) are formal power series with g(0) =6 0, f(0) = 0 and f ′(0) =6 0. The matrix corresponding to the pair f(z),g(z) is denoted by (g(z), f(z)). The product of two Riordan arrays (g(z), f(z)) and (h(z),l(z)) is defined by (g(z), f(z)) ∗ (h(z),l(z))=(g(z)h (f(z)) ,l (f(z))) . We recall that the set of all Riordan matrices is a group under the operator “ ∗ ” [32]. The identity element is I = (1, z), and the inverse of (g(z), f(z)) is (18) (g(z), f(z))−1 = 1/ g ◦ f (z), f(z) , Sometimes, it is useful to use exponential generating functions instead of ordinary generating functions when we apply a Riordan array method. We call the resulting array an exponential Riordan array and we denote it by hg(z), f(z)i. Its column k has generating function g(z)(f(z))k /k!,k =0, 1, 2,... (cf. [34]). 16 MIGUEL A. MENDEZ´ AND JOSE´ L. RAMÍREZ

The r-Whitney numbers of the second kind are given by the exponential Riordan array:

emx − 1 W := [W (n, k)] = erx, . 2 m,r n,k≥0 m

[m,r] Therefore, as we have already seen in Remark 14, the r-Dowling polynomial Dn (x) is a Sheﬀer sequence for

(1 + mx)−r/m, ln(1 + mx)1/m . In this section we give some explicit relations between the r-Dowling polynomials and the Bernoulli and Euler polynomials by using the connection constants an,k in the expression n rn(x) = k=0 an,ksk(x), where the polynomials rn(x) and sn(x) are Sheﬀer sequences. These constants can be determined by the umbral method [30, pp. 131] or equivalent P by Riordan arrays [34, Theorem 6.4]. In particular, let sn(x) and rn(x) be Sheﬀer for n (g(t), f(t)) and (h(t),l(t)), respectively. If rn(x) = k=0 an,ksk(x), then an,k is the entry (n, k)-th of the Riordan array P

g(l(t)) (19) , f(l(t)) . h(l(t))

The Bernoulli polynomials, Bn(x), are deﬁned by the exponential generating function

∞ tn text B (x) = . n n! et − 1 n=0 X

The Bernoulli numbers, Bn, are define by Bn := Bn(0). Moreover, the polynomials Bn(x) et−1 are Sheffer for t , t (cf. [30, 12]). The r-Whitney numbers of the first kind are defined by exponential generating function [24]:

∞ n k z − r ln (1 + mz) (20) w (n, k) =(1+ mz) m . m,r n! mkk! Xn=k Theorem 16. For n ≥ 0 we have

n n n [m,r] (21) B (x)= B − w (l,k)D (x). n l n l m,r k Xk=0 Xl=k A NEW APPROACH TO THE r-WHITNEY NUMBERS 17

n [m,r] Proof. If Bn(x)= k=0 an,kDk (x), then from (19) and (20) we get P 1 t k a = [tn] (1 + mt)−r/m · ln(1 + mt)1/m n,k k! et − 1 lnk(1 + mt) t = [tn] (1 + mt)−r/m mkk! et − 1 ∞ ∞ tn tn = [tn] w (n, k) B (x) m,r n! n n! ! n=0 !! Xn=k X ∞ n n n n z = [t ] B − w (l,k) . l n l m,r n! n=0 ! X Xl=0 Therefore, it is clear (21).

The Euler polynomials, En(x), are deﬁned by the exponential generating function

∞ tn 2ext E (x) = . n n! et +1 n=0 X

The numbers En, are deﬁne by En := En(0). Moreover, the polynomials En(x) are Sheﬀer et+1 for 2 , t . From a similar argument as in above theorem we get the following theorem.

Theorem 17. For n ≥ 0 we have

n n n [m,r] E (x)= E − w (l,k)D (x). n l n l m,r k Xk=0 Xl=k In the following theorem we analyze the connecting coeﬃcients

n [m,r] Dn (x)= an,kBk(x). Xk=0 Theorem 18. For n ≥ 0 we have

n n−k l [m,r] 1 n +1 l +1 l−s [m] (22) D (x)= W (n − l,k)m T (1)B − B (x). n n +1 l +1 s +1 m,r s+1 l s k s=0 Xk=0 Xl=0 X 18 MIGUEL A. MENDEZ´ AND JOSE´ L. RAM´IREZ

Proof. From Equations (19), (1) and (10) we get

m−1 e −1 mt k 1 e m e − 1 a = [tn] ert · n,k k! emt−1 m m ! m−1 rt mt k e −1 e e − 1 e m mt = [tn] · · k! m t emt−1 ! ∞ ∞ ∞ tn tn tn = [tn] W (n, k) · T [m](1) · mnB m,r n! n (n + 1)! n n! ! n=0 ! n=0 !! Xn=k X X ∞ n ∞ l [m] n n t l l−s Ts+1(1) t = [t ] W (n, k) · m B − m,r n! s l s s +1 n! ! n=0 s=0 !! Xn=k X X ∞ n−k l [m] n n n l l−s Ts+1(1) t = [t ] W (n − l,k)m B − l s m,r l s s +1 n! n=0 s=0 ! ! X Xl=0 X ∞ n−k l n n 1 n +1 l +1 l−s [m] t = [t ] W (n − l,k)m B − T (1) n +1 l +1 s +1 m,r l s s+1 n! n=0 s=0 ! ! X Xl=0 X

Therefore, it is clear (21).

From a similar argument we get the following theorem. Theorem 19. For n ≥ 0 we have

n n−k [m,r] 1 n [m] 1 Dn (x)= Wm,r(n − l,k)Tl (1) + Wm,r(n, k) Ek(x). 2 l 2 ! Xk=0 Xl=0 5.1. Some Recurrence Relations. Rogers [29] introduced the concept of the A-sequence. Speciﬁcally, Rogers observed that every element dn+1,k+1 of a Riordan matrix (not belong- ing to row or column 0) could be expressed as a linear combination of the elements in the preceding row. Merlini et al. [22] introduced the Z-sequence, which characterizes column 0, except for the element d0,0. Therefore, the A-sequence, Z-sequence and the element d0,0 completely characterize a proper Riordan array. Summarizing, we have the following theorem.

Theorem 20 ([22]). An inﬁnite lower triangular array D = [dn,k]n,k∈N is a Riordan array if and only if d0,0 =06 and there are sequences A =(a0 =06 , a1, a2,... ) and Z =(z0, z1, z2,... ) such that if n, k ≥ 0, then

dn+1,k+1 = a0dn,k + a1dn,k+1 + a2dn,k+2 + ··· ,

dn+1,0 = z0dn,0 + z1dn,1 + z2dn,2 + ··· , A NEW APPROACH TO THE r-WHITNEY NUMBERS 19 or equivalently g(0) g(z)= and f(z)= z(A(f(z))), 1 − zZ(g(z)) where A and Z are the generating functions of the A-sequence and Z-sequence, respectivally. The generating function for the A-sequence of the exponential Riordan array of the r- Whitney numbers is given by

∞ k t t mt k t = = = ckm , f(t) ln(1 + mt)1/m ln(1 + mt) k! Xk=0 1 n where ck are the Cauchy numbers of first kind. They are defined by cn = 0 x dx. See [23] for general information about Cauchy numbers. Taking in count that any exponential R Riordan array hg(x), f(x)i =(dn,k)n,k≥0 satisfies the recurrence relations (see [34, Corollary 5.7])

∞ n +1 k + j d = j!a d , n+1,k+1 k +1 j j n,k+j j=0 X n n − 1 d − d˜ − = f − d − − , n,k n 1,k l − 1 n l+1 l 1,k 1 Xl=k n n kd = f − d − − , n,k l − 1 n l+1 l 1,k 1 Xl=k ′ where (aj) is the A-sequence and (d˜n,k)n,k≥0 = hg (x), f(x)i. Then we obtain the following corollary. Corollary 21. The r-Whitney numbers of the second kind satisfy the following recurrence relations ∞ n +1 k + j W (n +1,k +1) = c mjW (n, k + j), m,r k +1 j j m,r j=0 X n n − 1 W (n, k) − rW (n − 1,k)= mn−lW (l − 1,k − 1), n ≥ 1, m,r m,r l − 1 m,r Xl=k n n kW (n, k)= mn−lW (l − 1,k − 1). m,r l − 1 m,r Xl=k From the generating function (20) we obtain that the r-Whitney numbers of the ﬁrst kind are given by the exponential Riordan array

−r/m 1/m W1 := [wm,r(n, k)] = (1 + mz) , ln(1 + mz) .

20 MIGUEL A. MENDEZ´ AND JOSE´ L. RAM´IREZ

In this case, the generating function for the A-sequence of the exponential Riordan array W1 is given by ∞ t t mt tk = = = B mk , emt−1 emt − 1 k k! f(t) m Xk=0 where Bn are the Bernoulli numbers. Therefore we obtain the following corollary. Corollary 22. The r-Whitney numbers of theﬁrst kind satisfy the following recurrence relations ∞ n +1 k + j w (n +1,k +1)= B mjw (n, k + j), m,r k +1 j j m,r j=0 X − n 1 n − 1 w (n, k)+ r (n − l − 1)!w (l,k)(−m)n−l−1 m,r l m,r Xl=0 n n − 1 = (−m)n−l(n − l)!w (l − 1,k − 1), n ≥ 1, l − 1 m,r Xl=k n n kw (n, k)= (−m)n−l(n − l)!w (l − 1,k − 1). m,r l − 1 m,r Xl=k 5.2. Determinantal Identity. The r-Whitney numbers of both kinds satisfy the following orthogonality relation (cf. [24]):

n n

Wm,r(n, i)wm,r(i, s)= wm,r(n, i)Wm,r(i, s)= δs,n, i=s i=s X X where δs,n = 1 if s = n and 0 otherwise. From above relations we obtain the inverse relation: n n

fn = wm,r(n, s)gs ⇐⇒ gn = Wm,r(n, s)fs. s=0 s=0 X X −1 Moreover, we have the identity W1 = W2 , where W2 is the exponential Riordan array for the r-Whitney numbers of the second kind. From deﬁnition of the r-Dowling polynomials we obtain the following equality

W2 · X = Dm,r,

2 T [m,r] [m,r] [m,r] T where X = [1, x, x ,... ] and Dm,r = [D0 (x), D1 (x), D2 (x),... ] . Then X = W1Dm,r and n n [m,r] x = wm,r(n, k)Dk (x). Xk=0 A NEW APPROACH TO THE r-WHITNEY NUMBERS 21

Therefore n−1 [m,r] n [m,r] (23) Dn (x)= x − wm,r(n, k)Dk (x), n ≥ 0. Xk=0 From the above equation we obtain the following determinantal identity. Theorem 23. The r-Dowling polynomials polynomials satisfy 1 x ··· xn−1 xn 1 wm,r(1, 0) ··· wm,r(n − 1, 0) wm,r(n, 0)

D[m,r](x)=(−1)n 0 1 ··· wm,r(n − 1, 1) wm,r(n, 1) n . . . . ··· 0 0 ··· 1 w (n, n − 1) m,r

Proof. This identity follows from Equation (23) and by expanding the determinant by the last column.

References [1] M. Benoumhani. On some numbers related to Whitney numbers of Dowling lattices. Adv. Appl. Math. 19(1997), 106–116. [2] M. Benoumhani. On Whitney numbers of Dowling lattices. Discrete Math. 159(1996), 13–33. [3] F. Bergeron, G. Labelle, P. Leroux. Combinatorial Species and Tree-Like Structures. Encyclopedia of Mathematic and its Applications. Volume 67. Cambridge University Press, Cambridge, 1998. [4] A. Z. Broder. The r-Stirling numbers. Discrete Math. 49(1984), 241–259. [5] D. Callan, S.-M. Ma, T. Mansour. Some combinatorial arrays related to the Lotka-Volterra system. Electron. J. Combin. 22(2) (2015), # P2.22. [6] W. Y. C. Chen. Context-Free grammars, differential operators and formal power series. Theoret. Comput. Sci. 117(1993), 113–129. [7] W. Y. C. Chen, A. M. Fu. Context-free grammars for permutations and increasing trees. Adv. Appl. Math. 82(2017), 58–82. [8] G.-S. Cheon, J.-H. Jung. r-Whitney numbers of Dowling lattices. Discrete Math. 312(15)(2012), 2337–2348. [9] C. B. Corcino, R. B. Corcino, R. J. Gasparin. Equivalent asymptotic formulas of second kind r- Whitney numbers. Integral Transforms Spec. Funct. 26(3)(2015), 192–202. [10] T. A. Dowling. A class of geometric lattices based on finite groups. J. Combin. Theory Ser. B 14(1973) 61–86. erratum; J. Combin. Theory Ser. B 15(1973), 211. [11] R. Hao, L. Wang, H. Yang. Context-free grammars for triangular arrays. Acta Math. Sin. (Engl. Ser.) 31(3)(2015), 445–455. [12] O. Herscovici, T. Mansour. Identities involving Touchard polynomials derived from umbral calculus. Adv. Stud. Contemp. Math. (Kyungshang) 25(1)(2015), 39–46. [13] A. Joyal. Une théorie combinatoire des séries formelles. Adv. Math. 42(1981), 1–82. [14] P. Leroux, G. X. Viennot. Combinatorial resolution of systems of differential equations, I. Ordinary differential equations. In Combinatoire Enumérativé(1986), G. Labelle, P. Leroux, Eds., no. 1234 in Lecture Notes in Mathematics, Springer-Verlag, 210–245. [15] S.-M. Ma. Some combinatorial arrays generated by context-free grammars. European J. Combin. 34(2013), 1081–1091. [16] S.-M. Ma. Derivative polynomials and enumeration of permutations by number of interior and left peaks. Discrete Math. 312(2012), 405–412. 22 MIGUEL A. MENDEZ´ AND JOSE´ L. RAMÍREZ

[17] S.-M. Ma, T. Mansour, M. Schork. Normal ordering problem and the extensions of the Stirling grammar. Russ. J. Math. Phys. 21(2)(2014), 242–255. [18] T. Mansour, J. L. Ram´ırez, M. Shattuck. A generalization of the r-Whitney numbers of the second kind. J. Comb. 8(1)(2017) 29–55. [19] M. Méndez. Combinatorial differential operators in: Faàdi Bruno formula, enumeration of bal- lot paths, enriched rooted trees and increasing rooted trees. ArXiv preprint, available online at http://arxiv.org/pdf/1610.03602v1.pdf [20] M. Merca. A convolution for the complete and elementary symmetric functions. Aequationes Math. 86(3) (2013), 217–229. [21] M. Merca. A note on the r-Whitney numbers of Dowling lattices. C. R. Math. Acad. Sci. Paris. 351(16-17)(2013), 649–655. [22] D. Merlini, D. G. Rogers, R. Sprugnoli, M. Cecilia Verri. On some alternative characterizations of Riordan arrays. Canadian J. Math. 49(1997), 301–320. [23] D. Merlini, R. Sprugnoli, M. C. Verri. The Cauchy numbers. Discrete Math. 306(2006), 1906–1920. [24] I. Mez˝o. A new formula for the Bernoulli polynomials. Result. Math. 58(3)(2010), 329–335. [25] I. Mez˝o, J. L. Ram´ırez. The linear algebra of the r-Whitney matrices. Integral Transforms Spec. Funct. 26(3)(2015), 213–225. [26] I. Mez˝o, J. L. Ram´ırez. Some identities of the r-Whitney numbers. Aequationes Math. 90(2) (2016), 393–406. [27] M. Mihoubi, M. Rahmani. The partial r-Bell polynomials. ArXiv preprint, available online at http://arxiv.org/pdf/1308.0863v1.pdf. [28] M. Mihoubi, M. Tiachachat. Some applications of the r-Whitney numbers. C. R. Math. Acad. Sci. Paris. 352(12)(2014), 965–969. [29] D. G. Rogers. Pascal triangles, Catalan numbers and renewal arrays. Discrete Math. 22 (1978), 301– 310. [30] S. Roman. The Umbral Calculus. Pure and Applied Mathematics, 111. Academic Press Inc. (1984). [31] S. Roman, G.-C. Rota. The Umbral Calculus. Adv. Math. 27(1978), 95–188. [32] L. W. Shapiro, S. Getu, W. Woan, L. Woodson. The Riordan group. Discrete Appl. Math. 34(1991), 229–239. [33] M. Z. Spivey. A generalized recurrence for Bell numbers. J. Integer Seq. 11(2008), Article 08.2.5. [34] W. Wang, T. Wang. Generalized Riordan arrays. Discrete Math. 308(2008), 6466–6500. [35] A. Xu. Extensions of Spivey’s Bell number formula. Electron. J. Combin. 19(2), #P6(2012). [36] A. Xu, T. Zhou. Some identities related to the r-Whitney numbers. Integral Transforms Spec. Funct. 27(11)(2016), 920–929.

Departamento de Matematicas,´ Facultad de Ciencias, Universidad Antono Narino,˜ Bo- gota,´ Colombia E-mail address: [email protected] Departamento de Matematicas,´ UNAL, Bogota,´ COLOMBIA E-mail address: [email protected]