Next Article in Journal
Strong Chromatic Index of Outerplanar Graphs
Previous Article in Journal
On an Integral Equation with the Riemann Function Kernel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Note on the Higher-Order Derivatives of the Hyperharmonic Polynomials and the r-Stirling Polynomials of the First Kind

by
José L. Cereceda
Independent Researcher, Collado Villalba, 28400 Madrid, Spain
Axioms 2022, 11(4), 167; https://doi.org/10.3390/axioms11040167
Submission received: 14 February 2022 / Revised: 31 March 2022 / Accepted: 6 April 2022 / Published: 8 April 2022
(This article belongs to the Section Mathematical Analysis)

Abstract

:
In this paper, we focus on the higher-order derivatives of the hyperharmonic polynomials, which are a generalization of the ordinary harmonic numbers. We determine the hyperharmonic polynomials and their successive derivatives in terms of the r-Stirling polynomials of the first kind and show the relationship between the (exponential) complete Bell polynomials and the r-Stirling numbers of the first kind. Furthermore, we provide a new formula for obtaining the generalized Bernoulli polynomials by exploiting their link with the higher-order derivatives of the hyperharmonic polynomials. In addition, we obtain various identities involving the r-Stirling numbers of the first kind, the Bernoulli numbers and polynomials, the Stirling numbers of the first and second kind, and the harmonic numbers.

1. Introduction

The n-th harmonic number is the sum of the reciprocals of the first n positive integers
H n = 1 + 1 2 + 1 3 + + 1 n ,
with H 0 = 0 . Harmonic numbers play an important role in various branches of mathematics and applications such as number theory, combinatorics, analysis, special functions, and computer science and have been profusely generalized by many authors (see, for instance, [1,2,3,4,5,6,7,8,9,10,11,12,13,14] and references therein). In particular, the importance of the harmonic numbers and their generalizations in the evaluation of the special values of the Riemann zeta function, Hurwitz zeta function, and more generally zeta functions of arithmetical nature (see, for instance, [15,16,17,18,19]) should be stressed.
In this paper, we deal with a generalization of the harmonic numbers known as hyperharmonic numbers, which were introduced by Conway and Guy in 1996 [20] (p. 258). Following the notation in [17,21,22], we denote the hyperharmonic numbers by H n ( r ) . These can be defined recursively as follows.
Definition 1.
For integers n , r 0 , the n-th hyperharmonic number of order r, H n ( r ) , is given by
H n ( r ) = 0 if n = 0 , 1 n if n > 0 and r = 0 , i = 1 n H i ( r 1 ) if n , r 1 .
Clearly, H n ( 1 ) = H n . Equivalently, the hyperharmonic numbers can be defined by the generating function [23]
n = 0 H n ( r ) t n = ln ( 1 t ) ( 1 t ) r .
Moreover, Conway and Guy [20] (p. 258) provided the following formula
H n ( r ) = n + r 1 r 1 H n + r 1 H r 1 , r 1 ,
connecting the hyperharmonic numbers with the ordinary harmonic numbers.
Starting with (2), Mező [21] constructed the analytic extension of the hyperharmonic numbers. Specifically, he defined the hyperharmonic function H z ( w ) as follows:
H z ( w ) = ( z ) w z Γ ( w ) Ψ ( z + w ) Ψ ( w ) ,
involving the Pochhammer symbol ( z ) w , the classical gamma Γ ( w ) , and digamma Ψ ( w ) functions, and where w , z + w C \ Z { 0 } (see, e.g., [24] (Chapter 1) for a detailed description of Γ ( w ) and Ψ ( w ) ). Furthermore, Mező computed the first derivative of H z ( w ) with respect to the variables w and z. Subsequently, Dil [22] presented formulas to calculate special values of H z ( w ) and showed a way to obtain higher derivatives of the hyperharmonic function with respect to w and z. In this regard, he proved that [22] (Equation (15))
d n + 1 d w n + 1 ( w ) z = Γ ( z + 1 ) d n d w n H z ( w ) ,
showing that the higher-order derivatives (with respect to w) of H z ( w ) can be expressed in terms of the higher-order derivatives of ( w ) z . For the special case in which z is a non-negative integer, we have Γ ( z + 1 ) = z ! , and thus
d n d w n H z ( w ) = 1 z ! d n + 1 d w n + 1 ( w ) z = d n + 1 d w n + 1 w + z 1 z ,
where the general binomial coefficient [25] (Equation (1.2))
w + z 1 z = w ( w + 1 ) ( w + 2 ) ( w + z 1 ) z ! ,
applies to any complex number w and non-negative integer z. Evaluating the above derivative (3) at w = 1 for the simplest case n = 0 gives the well-known relation (see, e.g., [2] (Equation (8)))
H z = d d w w + z 1 z w = 1 = d d w w + z z w = 0 .
However, the problem of finding closed-form expressions for the higher-order derivatives of H z ( w ) considered in [21,22] was left open. In this paper, we completely determine the higher-order derivatives (with respect to x) of the class of hyperharmonic polynomials H j ( x ) , where the subscript j denotes a non-negative integer variable, while the superscript x stands for any arbitrary (real or complex) value. Next, we define the hyperharmonic polynomials H j ( x ) through their generating function. Without loss of generality, in what follows, we may restrict x to belong to R .
Definition 2.
For integer j 0 and x R , the sequence { H 0 ( x ) , H 1 ( x ) , H 2 ( x ) , } of the hyperharmonic polynomials H j ( x ) Q [ x ] is determined by the generating function
j = 0 H j ( x ) t j = ln ( 1 t ) ( 1 t ) x , | t | < 1 ,
where, by definition, H 0 ( x ) = 0 .
Note that (4) is just the generating function (1) with r replaced by the continuous variable x, so that H j ( x ) reduces to the hyperharmonic number H j ( r ) when x is the non-negative integer r. In particular, H j ( 1 ) = H j . (See [5,26,27,28] for other types of polynomials associated with the harmonic numbers.) The hyperharmonic polynomials admit, among others, the following representations [29]:
H j ( x ) = d d x x + j 1 j ,
H j + 1 ( x + 1 ) = t = 0 j 1 j + 1 t x + t t ,
H j + 1 ( x ) = t = 0 j ( 1 ) t j + 1 t x t ,
and
H j + 1 ( x ) = 1 j ! t = 0 j j + 1 t + 1 B t ( x ) ,
where j t represents the (unsigned) Stirling numbers of the first kind and B t ( x ) is the t-th Bernoulli polynomial. For j 1 , H j ( x ) is a polynomial in x of degree j 1 with leading coefficient 1 ( j 1 ) ! and constant term 1 j . The first few hyperharmonic polynomials are given explicitly by
H 1 ( x ) = 1 , H 2 ( x ) = x + 1 2 , H 3 ( x ) = 1 2 x 2 + x + 1 3 , H 4 ( x ) = 1 6 x 3 + 3 4 x 2 + 11 12 x + 1 4 , H 5 ( x ) = 1 24 x 4 + 1 3 x 3 + 7 8 x 2 + 5 6 x + 1 5 , H 6 ( x ) = 1 120 x 5 + 5 48 x 4 + 17 36 x 3 + 15 16 x 2 + 137 180 x + 1 6 , H 7 ( x ) = 1 720 x 6 + 1 40 x 5 + 25 144 x 4 + 7 12 x 3 + 29 30 x 2 + 7 10 x + 1 7 , H 8 ( x ) = 1 5040 x 7 + 7 1440 x 6 + 23 480 x 5 + 35 144 x 4 + 967 1440 x 3 + 469 480 x 2 + 363 560 x + 1 8 .
Remark 1.
Notice that (5) immediately implies that 0 1 H j ( x ) d x = 1 for all j 1 . Furthermore, (5) entails that
d i d x i H j ( x ) = d i + 1 d x i + 1 x + j 1 j ,
in accordance with (3).
At this point, we should remark two previous results that are of particular interest for the present work. For that, it is convenient to introduce the symbol P ( i , j + r , r ) (with i, j, and r being non-negative integers), which will be defined later in Section 4 (see Definition 4). In essence, P ( i , j + r , r ) denotes the value of the (exponential) complete Bell polynomial for certain arguments related to the harmonic numbers. The first result that we are interested in was given by Wang [30] (Equation (4.2)) and tells us that (in our notation)
i ! k = 0 j j k k i ( r + 1 ) k i = j ! j + r r P ( i , j + r , r ) = d i d x i ( x ) j x = r + 1 .
Let us recall that the r-Stirling numbers of the first kind, k j r , count the number of permutations of the set { 1 , 2 , , k } having j disjoint (non-empty) cycles, such that the first r elements belong to distinct cycles [31]. In particular, k j 0 = k j and k + 1 j + 1 1 = k + 1 j + 1 , where k j is the ordinary Stirling number of the first kind. Actually, as will be shown in Section 2, the leftmost part in the above double identity is equal to i ! times the r-Stirling number of the first kind j + r + 1 i + r + 1 r + 1 . Wang’s result can then equivalently be stated as
j + r r P ( i , j + r , r ) = i ! j ! j + r + 1 i + r + 1 r + 1 = d i d x i x + j + r j x = 0 .
In particular, for r = 0 , we have
P ( i , j , 0 ) = i ! j ! j + 1 i + 1 = d i d x i x + j j x = 0 ,
and P ( 1 , j , 0 ) = H j .
Remark 2.
The relation d i d x i x + j + r j x = 0 = j + r r P ( i , j + r , r ) was rederived (in a somewhat different form) by Wang and Jia in [32] (Theorem 2).
The second result that we alluded to before appears in Equation (19) of the paper by Kargın et al. [12] and tells us that (in our notation)
d i 1 d x i 1 H j ( x + 1 ) x = r = j + r r P ( i , j + r , r ) ,
relating the higher-order derivatives of the hyperharmonic polynomials and the complete Bell polynomials.
Combining the mentioned result of Wang in Equation (10) and that of Kargın et al. in Equation (11), we immediately obtain
d i 1 d x i 1 H j ( x + 1 ) x = r = i ! j ! j + r + 1 i + r + 1 r + 1 ,
or, equivalently,
d i 1 d x i 1 H j ( x ) x = r = i ! j ! j + r i + r r ,
which applies to any non-negative integer r.
Remark 3.
When i = 1 , (12) reduces to
H j ( r ) = 1 j ! j + r r + 1 r ,
expressing the hyperharmonic numbers in terms of the r-Stirling numbers of the first kind. A combinatorial proof of (13) was given in [23] (Theorem 2).
However, although Equation (12) enables us to evaluate the i-th derivative of H j ( x ) for the particular case in which x is a non-negative integer, we still lack a general expression for the higher-order derivatives of the hyperharmonic polynomials. One of the principal objectives of the present paper is to provide a thorough account of the hyperharmonic polynomials and their successive derivatives in terms of the r-Stirling polynomials of the first kind. The rest of the paper is organized as follows.
In Section 2, we introduce the r-Stirling polynomials of the first kind R m , i ( x ) and R ¯ m , i ( x ) . Based on the properties of these polynomials, in Section 3, we express the higher-order derivatives of the hyperharmonic polynomials, d i d x i H j + 1 ( x ) and d i d x i H j + 1 ( x + 1 ) , as an explicit polynomial in x of degree j i . As anticipated by Equation (12), such derivatives can in turn be expressed in terms of the r-Stirling numbers of the first kind when they are evaluated at the non-negative integer r. In Section 4, we exhibit the relationship between the (exponential) complete Bell polynomials and the r-Stirling numbers of the first kind. Specifically, we show that Wang’s relation
j + r r P ( i , j + r , r ) = i ! j ! j + r + 1 i + r + 1 r + 1
arises as a particular case of a theorem due to Kölbig [33] (Theorem). In Section 5, we study the connection (already established in [22] (Proposition 3.13)) between the generalized Bernoulli polynomials and the higher-order derivatives of the hyperharmonic polynomials and provide a new formula for obtaining the generalized Bernoulli polynomials (see Equation (51) below). In Section 6, we consider a series of identities obtained by Spieß [34], Wang [30], and Wuyungaowa [35] involving the numbers j + r r P ( i , j + r , r ) and recast them, by means of (14), into a form involving the r-Stirling numbers of the first kind.
In addition, throughout this paper, we obtain various identities involving the r-Stirling numbers of the first kind, the Bernoulli numbers and polynomials, the Stirling numbers of the first and second kind, and the harmonic numbers. As a preliminary example of such identities, let us observe that, from (8) and (13), we quickly obtain
j = 0 k k + 1 j + 1 B j ( r ) = 1 k + 1 k + r + 1 r + 1 r ,
for any non-negative integers k and r. Furthermore, the Stirling transform [36] (Appendix A) of (15) yields
B k ( r ) = ( 1 ) k j = 0 k ( 1 ) j j + 1 k + 1 j + 1 j + r + 1 r + 1 r ,
and, in particular,
B k = ( 1 ) k j = 0 k ( 1 ) j j ! j + 1 k + 1 j + 1 ,
where B k = B k ( 0 ) are the Bernoulli numbers and k j are the Stirling numbers of the second kind.

2. r-Stirling Polynomials of the First Kind

In order to properly define the r-Stirling polynomials of the first kind R m , i ( x ) and R ¯ m , i ( x ) , we next formulate the following proposition.
Proposition 1.
For integers 0 i m and r 0 , we have
m + r i + r r = j = 0 m i i + j i m i + j r j ,
and
1 1 m + r + 1 i + r + 1 r + 1 = j = 0 m i i + j i m + 1 i + j + 1 r j .
Proof. 
We prove only the relation (18), since the proof of relation (17) proceeds in an analogous way. According to [31] (Equation (27)) (see also [37] (p. 224)), m + r + 1 i + r + 1 r + 1 can be expressed as
m + r + 1 i + r + 1 r + 1 = j = i m m j j i ( r + 1 ) m j ¯ ,
where r m ¯ = r ( r + 1 ) ( r + m 1 ) , and r 0 ¯ = 1 . Changing the summation variable from j to t, where t = m j , results in
m + r + 1 i + r + 1 r + 1 = t = 0 m i m t m t i ( r + 1 ) t ¯ .
Now, we have
( r + 1 ) t ¯ = k = 0 t t k ( r + 1 ) k = k = 0 t s = 0 k k s t k r s ,
and then
m + r + 1 i + r + 1 r + 1 = t = 0 m i k = 0 t s = 0 k m t m t i k s t k r s = t = 0 m i s = 0 t k = s t k s t k m t m t i r s .
Using the well-known identity (see, e.g., [38] (Equation (6.16))) k = s t k s t k = t + 1 s + 1 , the preceding equation becomes
m + r + 1 i + r + 1 r + 1 = t = 0 m i s = 0 t m t m t i t + 1 s + 1 r s = j = 0 m i t = j m i m t t + 1 j + 1 m t i r j ,
where we have renamed the variable s as j. Invoking the identity (see, e.g., [31] (Equation (52)) and [37] (p. 224))
j + i i m + r + s j + i + r + s r + s = t = j m i m t t + r j + r r m t + s i + s s ,
and specializing to the case in which r = 1 and s = 0 , it follows that
i + j i m + 1 i + j + 1 = t = j m i m t t + 1 j + 1 m t i .
Subsequently, combining (19) and (20), we obtain (18). □
Remark 4.
From (17), it follows that
j + r + 1 i + r + 1 r + 1 = k = 0 j i i + k i j i + k ( r + 1 ) k = k = 0 j k i j k ( r + 1 ) k i ,
as was noted in the introduction in relation to Wang’s identity (10). On the other hand, from (18), we can see that m + r + 1 i + r + 1 r + 1 can be expressed as a polynomial in r with constant term m + 1 i + 1 . Hence, setting r = 0 and j = 0 in (19) yields the identity (cf. [31] (Equation (30)))
m + 1 i + 1 = t = 0 m i t ! m t m t i .
Remark 5.
For integers k , n 0 and m 1 , the hyper-sums of powers of integers S k ( m ) ( n ) are defined recursively by
S k ( m ) ( n ) = j = 1 n S k ( m 1 ) ( j ) ,
with initial conditions S k ( 0 ) ( n ) S k ( n ) = 1 k + 2 k + + n k , and S k ( m ) ( 0 ) = 0 . In [12] (Equation (27)), Kargın et al. obtained the following formula for S k ( m ) ( n ) (see also the related paper [39])
S k ( m ) ( n ) = 1 m ! i = 0 m ( 1 ) i m + n + 1 i + n + 1 n + 1 S k + i ( n ) .
When k = 0 , the hyper-sum S k ( m ) ( n ) is equal to S 0 ( m ) ( n ) = n + m m + 1 . Hence, putting k = 0 in the last equation gives
i = 0 m ( 1 ) i m + n + 1 i + n + 1 n + 1 S i ( n ) = m ! n + m m + 1 .
In particular, for n = 1 , we find
i = 0 m ( 1 ) i m + 2 i + 2 2 = m ! .
In view of Proposition 1, and following Broder [31] (Equation (56)) and Carlitz [40] (Equation (5.2)), we define the r-Stirling polynomials of the first kind as follows.
Definition 3.
For integers 0 i m , the r-Stirling polynomials of the first kind R m , i ( x ) and R ¯ m , i ( x ) are defined as
R m , i ( x ) = j = 0 m i i + j i m i + j x j ,
and
R ¯ m , i ( x ) = j = 0 m i i + j i m + 1 i + j + 1 x j ,
respectively.
By construction, R m , i ( x ) and R ¯ m , i ( x ) reduce to R m , i ( r ) = m + r i + r r and R ¯ m , i ( r ) = m + r + 1 i + r + 1 r + 1 , respectively, when x is a non-negative integer r, so that R ¯ m , i ( r ) = R m , i ( r + 1 ) for any integer r 0 . Since R ¯ m , i ( x ) R m , i ( x + 1 ) vanishes for infinitely many r N , R ¯ m , i ( x ) = R m , i ( x + 1 ) for an arbitrary x. Hence, we have
j = 0 m i i + j i m + 1 i + j + 1 x j = j = 0 m i i + j i m i + j ( x + 1 ) j .
In particular, when x = 1 , it follows that
j = 0 m i ( 1 ) j i + j i m + 1 i + j + 1 = j = 0 m ( 1 ) j i j i m + 1 j + 1 = m i ,
which corresponds to the identity [38] (Equation (6.18)).

3. Hyperharmonic Polynomials and Their Derivatives

From [31] (Theorem 28), we know that R m , i ( x ) = 1 i ! d i d x i x m ¯ . This follows immediately from the horizontal generating function of the Stirling numbers of the first kind, namely,
x m ¯ = j = 0 m m j x j .
Indeed, differentiating i times (with respect to x) both sides of (25), we obtain
1 i ! d i d x i x m ¯ = j = 0 m j i m j x j i = j = 0 m i i + j i m i + j x j .
Thus,
d i d x i x m ¯ = i ! R m , i ( x ) .
Now, we present the following theorem, which shows the explicit expression of the higher-order derivatives of the hyperharmonic polynomials H j + 1 ( x ) and H j + 1 ( x + 1 ) .
Theorem 1.
For integers i , j 0 , the i-th derivative of H j + 1 ( x ) and H j + 1 ( x + 1 ) with respect to x is given, respectively, by
d i d x i H j + 1 ( x ) = ( i + 1 ) ! ( j + 1 ) ! t = 0 j i i + t + 1 i + 1 j + 1 i + t + 1 x t ,
and
d i d x i H j + 1 ( x + 1 ) = ( i + 1 ) ! ( j + 1 ) ! t = 0 j i i + t + 1 i + 1 j + 2 i + t + 2 x t .
Proof. 
Relation (27) readily follows from (26) and the representation (9) for the i-th derivative of H j ( x ) . Thus, combining (9) and (26), we obtain
d i d x i H j + 1 ( x ) = d i + 1 d x i + 1 x + j j + 1 = 1 ( j + 1 ) ! d i + 1 d x i + 1 x j + 1 ¯ = ( i + 1 ) ! ( j + 1 ) ! R j + 1 , i + 1 ( x ) .
Now, using the expression for R j + 1 , i + 1 ( x ) that is obtained from (22), we obtain (27). On the other hand, we have
d i d x i H j + 1 ( x + 1 ) = ( i + 1 ) ! ( j + 1 ) ! R j + 1 , i + 1 ( x + 1 ) = ( i + 1 ) ! ( j + 1 ) ! R ¯ j + 1 , i + 1 ( x ) ,
and, using the expression for R ¯ j + 1 , i + 1 ( x ) that is obtained from (23), we then obtain (28). □
Letting i = 0 in (27) and (28) gives us, respectively, H j + 1 ( x ) and H j + 1 ( x + 1 ) as polynomials in x of degree j, namely,
H j + 1 ( x ) = R j + 1 , 1 ( x ) ( j + 1 ) ! = 1 ( j + 1 ) ! t = 0 j ( t + 1 ) j + 1 t + 1 x t ,
and
H j + 1 ( x + 1 ) = R ¯ j + 1 , 1 ( x ) ( j + 1 ) ! = 1 ( j + 1 ) ! t = 0 j ( t + 1 ) j + 2 t + 2 x t .
Additionally, setting x = 1 and j j 1 in the first of the two equations above produces the well-known identity H j = 1 j ! t = 0 j t j t .
The following corollary is a direct consequence of Theorem 1 and Definition 3.
Corollary 1.
For non-negative integer r, we have
d i d x i H j + 1 ( x ) x = r = ( i + 1 ) ! ( j + 1 ) ! j + r + 1 i + r + 1 r ,
and
d i d x i H j + 1 ( x + 1 ) x = r = ( i + 1 ) ! ( j + 1 ) ! j + r + 2 i + r + 2 r + 1 .
Note that, obviously, any of (29) or (30) is equivalent to (12).
Example 1.
As a simple application of Theorem 1, we may use the representation (6) in combination with (28) to obtain
t = 0 j 1 j + 1 t d i d x i x + t t = ( i + 1 ) ! ( j + 1 ) ! t = 0 j i i + t + 1 i + 1 j + 2 i + t + 2 x t .
Therefore, it follows from (30) that
t = 0 j 1 j + 1 t d i d x i x + t t x = r = ( i + 1 ) ! ( j + 1 ) ! j + r + 2 i + r + 2 r + 1 .
Since d i d x i x + t t x = r = i ! t ! t + r + 1 i + r + 1 r + 1 , this can be expressed as (after replacing r by r 1 )
t = 0 j 1 t ! ( j + 1 t ) t + r i + r r = i + 1 ( j + 1 ) ! j + r + 1 i + r + 1 r ,
which holds for any integers 0 i j and r 0 . In particular, for r = 0 , we have
t = 0 j 1 t ! ( j + 1 t ) t i = i + 1 ( j + 1 ) ! j + 1 i + 1 .
Furthermore, setting i = r = 1 in (31) yields the identity
t = 0 j H t j + 1 t = 2 ( j + 1 ) ! j + 2 3 = H j + 1 2 H j + 1 [ 2 ] ,
which can also be found in [41] (p. 544), and where the notation H j [ 2 ] means t = 1 j 1 / t 2 .
Remark 6.
Using (25) and the recurrence relation j + 1 i + 1 = j j i + 1 + j i , it is easy to show that
d i + 1 d x i + 1 x j + 1 ¯ = ( x + j ) d i + 1 d x i + 1 x j ¯ + ( i + 1 ) d i d x i x j ¯ ,
or, in view of (26),
R j + 1 , i + 1 ( x ) = ( x + j ) R j , i + 1 ( x ) + R j , i ( x ) .
Hence, setting x = r (where r is a non-negative integer), we obtain the corresponding recurrence for the r-Stirling numbers of the first kind
j + r + 1 i + r + 1 r = ( r + j ) j + r i + r + 1 r + j + r i + r r .
Furthermore, letting i = 0 in (33) and recalling (13) leads to the following recurrence for the hyperharmonic numbers:
H j + 1 ( r ) = r + j j + 1 H j ( r ) + 1 j + 1 j + r 1 j ,
which holds for any integers r , j 0 . In particular, when r = 1 , we recover the recurrence relation defining the harmonic numbers, namely, H j + 1 = H j + 1 j + 1 , with H 0 = 0 .

4. Complete Bell Polynomials and r -Stirling Numbers of the First Kind

Let Y n ( x 1 , x 2 , , x n ) be the n-th (exponential) complete Bell polynomial defined by Y 0 = 1 and (see, e.g., [42] (p. 134))
exp j = 1 x j t j j ! = 1 + p = 1 Y p ( x 1 , x 2 , , x p ) t p p ! ,
and let H ( j , k ; r ) be the function
H ( j , k ; r ) = t = 1 j 1 ( t + r ) k ,
for integers r 0 and j , k 1 . Following Spieß [34] and Kargın et al. [12], we define the symbol P ( i , j + r , r ) as follows.
Definition 4.
For non-negative integers i, j, and r, we have
P ( i , j + r , r ) = 1 , f o r i = 0 , a n d j , r 0 , P i H ( j , 1 ; r ) , H ( j , 2 ; r ) , , H ( j , i ; r ) , f o r 1 i j , j 1 , a n d r 0 , 0 , f o r 0 j < i , i 1 , a n d r 0 ,
where the polynomial P i ( x 1 , x 2 , , x i ) is defined by
P i ( x 1 , x 2 , , x i ) = ( 1 ) i Y i ( 0 ! x 1 , 1 ! x 2 , , ( i 1 ) ! x i ) ,
or, equivalently,
P i ( x 1 , x 2 , , x i ) = Y i ( 0 ! x 1 , 1 ! x 2 , , ( 1 ) i 1 ( i 1 ) ! x i ) .
The first five polynomials P i ( x 1 , x 2 , , x i ) are given by
P 1 ( x 1 ) = x 1 , P 2 ( x 1 , x 2 ) = x 1 2 x 2 , P 3 ( x 1 , x 2 , x 3 ) = x 1 3 3 x 1 x 2 + 2 x 3 , P 4 ( x 1 , x 2 , x 3 , x 4 ) = x 1 4 6 x 1 2 x 2 + 8 x 1 x 3 + 3 x 2 2 6 x 4 , P 5 ( x 1 , x 2 , x 3 , x 4 , x 5 ) = x 1 5 10 x 1 3 x 2 + 20 x 1 2 x 3 + 15 x 1 x 2 2 30 x 1 x 4 20 x 2 x 3 + 24 x 5 .
Thus, we have, for example, P ( 1 , j + r , r ) = H ( j , 1 ; r ) , P ( 2 , j + r , r ) = H ( j , 1 ; r ) 2 H ( j , 2 ; r ) , P ( 3 , j + r , r ) = H ( j , 1 ; r ) 3 3 H ( j , 1 ; r ) H ( j , 2 ; r ) + 2 H ( j , 3 ; r ) , etc.
Remark 7.
By letting i = 1 and r r 1 in (11), and noting that P ( 1 , j + r 1 , r 1 ) = H ( j , 1 ; r 1 ) , we obtain (assuming that r 1 )
H j ( x + 1 ) x = r 1 = H j ( r ) = j + r 1 r 1 H ( j , 1 ; r 1 ) = j + r 1 r 1 t = 1 j 1 t + r 1 = j + r 1 r 1 H j + r 1 H r 1 ,
thus recovering Conway and Guy’s formula (2).
With these ingredients at hand, next we show that Wang’s relation (14) (which we reproduce here for convenience)
j + r r P ( i , j + r , r ) = i ! j ! j + r + 1 i + r + 1 r + 1 ,
is a direct consequence of the following theorem set forth by Kölbig [33] (Theorem).
Theorem 2
(Kölbig, 1994). Let α R with α 1 , 2 , , j , and
H ( j , k ; α ) = t = 1 j 1 ( t + α ) k ,
for integers j , k 1 . We then have
P q H ( j , 1 ; α ) , H ( j , 2 ; α ) , , H ( j , q ; α ) = q ! ( 1 + α ) j ¯ S ( j , q ; α ) , q j , 0 , q > j ,
where
S ( j , q ; α ) = t = q j t q j t ( 1 + α ) t q .
Indeed, taking α = r in Kölbig’s theorem (with integer r 0 ), we have
j + r r P ( i , j + r , r ) = j + r r P i H ( j , 1 ; r ) , H ( j , 2 ; r ) , , H ( j , i ; r ) = j + r r i ! ( 1 + r ) j ¯ S ( j , i ; r ) ,
where it is assumed that i j , with i = 1 , 2 , . Furthermore, it turns out that j + r r = 1 j ! ( 1 + r ) j ¯ ; thus,
j + r r P ( i , j + r , r ) = i ! j ! S ( j , i ; r ) ,
where
S ( j , i ; r ) = t = i j t i j t ( 1 + r ) t i .
On the other hand, we have
j + r + 1 i + r + 1 r + 1 = R j , i ( r + 1 ) = s = 0 j i i + s i j i + s ( r + 1 ) s = t = i j t i j t ( r + 1 ) t i .
Therefore, S ( j , i ; r ) = j + r + 1 i + r + 1 r + 1 ; thus, we obtain (34).
Remark 8.
By setting r = 0 in (34), we recover the well-known result (see [42] (Equation (7b)))
j + 1 i + 1 = j ! i ! P ( i , j , 0 ) = j ! i ! P i H j [ 1 ] , H j [ 2 ] , , H j [ i ] ,
where the notation H j [ i ] means t = 1 j 1 / t i .
Remark 9.
According to [34] (Theorem 16), it turns out that, for m , r 0 ,
k = 0 m P ( r , k , 0 ) k + 1 = P ( r + 1 , m + 1 , 0 ) r + 1 .
Noting that P ( r , k , 0 ) = r ! k ! k + 1 r + 1 , the above relation is equivalent to the identity (cf. [38] (Equation (6.21)))
m + 1 r + 1 = m ! k = 0 m 1 k ! k r .
Remark 10.
The generating function of the numbers j + r r P ( i , j + r , r ) is given by (see, e.g., [30] (Equation (1.6)))
j = i j + r r P ( i , j + r , r ) t j = ( ln ( 1 t ) ) i ( 1 t ) r + 1 .
Therefore, taking r r 1 in (34) and using (35), we obtain the exponential generating function of the r-Stirling numbers of the first kind, namely,
j = i j + r i + r r t j j ! = 1 i ! ( ln ( 1 t ) ) i ( 1 t ) r .
In particular, for r = 0 and i 1 , it follows that
( ln ( 1 t ) ) i = j = i i ! j ! j i t j = j = i i j P ( i 1 , j 1 , 0 ) t j ,
in agreement with [34] (Theorem 9). Note that setting i = 1 in the last equation yields the Maclaurin series of the natural logarithm
ln ( 1 t ) = j = 1 t j j = t t 2 2 t 3 3 .
Remark 11.
Using (34) into (21) leads to the relation
i = 0 m ( 1 ) i i ! P ( i , m + n , n ) S i ( n ) = n m + 1 ,
and, in particular,
i = 0 m ( 1 ) i i ! P ( i , m + 1 , 1 ) = 1 m + 1 .
We conclude this section by noting that, by virtue of (34), the theorem [12] (Theorem 6) can be reformulated as follows:
k = l n n + r k + r r k l B k l ( q ) = l + 1 n + 1 n + q + r l + q + r q + r 1 ,
where B k l ( q ) is a Bernoulli polynomial. Equation (36) may also be written as
k = l n n + r + 1 k + r + 1 r + 1 k l B k l ( q ) = l + 1 n + 1 n + q + r + 1 l + q + r + 1 q + r ,
or else,
k = l n n + r k + r r k l B k l ( q + 1 ) = l + 1 n + 1 n + q + r + 1 l + q + r + 1 q + r ,
for non-negative integers l , q , r , and n l . As an example, putting l = 2 , q = 0 , and r = 1 in (38), we obtain
k = 2 n ( 1 ) k n + 1 k + 1 k ( k 1 ) B k 2 = 6 n + 1 n + 2 4 = n ! H n + 1 3 3 H n + 1 H n + 1 [ 2 ] + 2 H n + 1 [ 3 ] ,
in accordance with the particular identity found in [12] (p. 8). Moreover, from (37) and (38), we find that
k = l n n + r + 1 k + r + 1 r + 1 k l B k l ( x ) = k = l n n + r k + r r k l B k l ( x + 1 ) ,
which holds for an arbitrary x. In particular, for r = 0 , and renaming the indices k l j , l i , and n m , the preceding identity becomes
j = 0 m i i + j i m + 1 i + j + 1 B j ( x ) = j = 0 m i i + j i m i + j B j ( x + 1 ) ,
which may be compared with (24).
Moreover, (36) can be generalized as follows:
k = l n n + r k + r r k l B k l ( x ) = l + 1 n + 1 k = l n k + 1 l + 1 n + 1 k + 1 ( x + r 1 ) k l ,
which holds for an arbitrary x. For l = 0 , (39) can be compactly written as
k = 0 n n + r k + r r B k ( x ) = n ! H n + 1 ( x + r 1 ) .
Conversely, (39) can be obtained by performing the l-th derivative with respect to x of both sides of (40).

5. Connection with the Generalized Bernoulli Polynomials

For a complex parameter α , the generalized Bernoulli polynomials B n ( α ) ( x ) are generated by the relation (see, e.g., [43,44,45])
z e z 1 α e x z = n = 0 B n ( α ) ( x ) z n n ! , | z | < 2 π ,
where B 0 ( α ) = 1 . Note that the value α = 1 corresponds to the classical Bernoulli polynomials B n ( x ) = B n ( 1 ) ( x ) . In the same way, B n ( α ) = B n ( α ) ( 0 ) are the generalized Bernoulli numbers (also known as Bernoulli numbers of order α or Nörlund polynomials), which are rational polynomials in α of degree n. In particular, B n = B n ( 1 ) are the ordinary Bernoulli numbers.
There is a close relationship between the generalized Bernoulli polynomials and the higher-order derivatives of the hyperharmonic polynomials. This relationship stems from the fact that (see, e.g., [46] (Equation (2.5)) and [22] (Equation (35)))
B ν ( m + 1 ) ( x ) = ( 1 ) m ν ! m ! d m ν d x m ν ( 1 x ) m ¯ ,
where m and ν are non-negative integers. Equivalently to (41), we have
B ν ( m + 1 ) ( x + 1 ) = ν ! d m ν d x m ν x m ,
in accordance with the result found by Gould in [25] (Equation (13.2)). In particular, for ν = m , from (42), we quickly obtain
B m ( m + 1 ) ( x + 1 ) = m ! x m = ( 1 ) m j = 0 m ( 1 ) j m j x j .
As indicated in [25] (see the sentence immediately following Equation (13.2) of [25]), a negative exponent in d m ν d x m ν x m has to be interpreted as integrating x m ν m times. Therefore, for example, when ν = m + 1 , from (42), we have
B m + 1 ( m + 1 ) ( x + 1 ) = ( m + 1 ) ! x m d x = B m + 1 ( m + 1 ) ( 1 ) + ( 1 ) m ( m + 1 ) j = 0 m ( 1 ) j j + 1 m j x j + 1 ,
where B m + 1 ( m + 1 ) ( 1 ) is the integration constant. Likewise, when ν = m + 2 , from (42) and the preceding equation, we have
B m + 2 ( m + 1 ) ( x + 1 ) = ( m + 2 ) B m + 1 ( m + 1 ) ( x + 1 ) d x = B m + 2 ( m + 1 ) ( 1 ) + ( m + 2 ) B m + 1 ( m + 1 ) ( 1 ) x + ( 1 ) m ( m + 1 ) ( m + 2 ) j = 0 m ( 1 ) j ( j + 1 ) ( j + 2 ) m j x j + 2 ,
where B m + 2 ( m + 1 ) ( 1 ) is the integration constant. In general, for an arbitrary m 0 and ν 0 , we have
B m + ν ( m + 1 ) ( x + 1 ) = t = 0 ν 1 m + ν t B m + ν t ( m + 1 ) ( 1 ) x t + ( 1 ) m m + ν ν j = 0 m ( 1 ) j j + ν ν 1 m j x j + ν ,
where it is understood that the first summation on the right hand side of (43) is zero when ν = 0 . For the special case in which m = 0 , (43) yields the Bernoulli polynomials evaluated at x + 1 , namely,
B ν ( x + 1 ) = ( 1 ) ν t = 0 ν ( 1 ) t ν t B ν t x t ,
with B ν t being the Bernoulli numbers.
On the other hand, as can be easily verified, (41) can be expressed in terms of the hyperharmonic polynomials as follows (see [22] (Proposition 3.13))
B ν ( m + 1 ) ( 1 x ) = ( 1 ) ν ν ! d m ν 1 d x m ν 1 H m ( x ) .
(Note the corrected factor ( 1 ) ν in (44) instead of the original factor ( 1 ) m + 1 appearing in [22] (Proposition 3.13).) Furthermore, (44) can in turn be written as
j i B j i ( j + 1 ) ( 1 x ) = ( 1 ) j i j ! i ! d i 1 d x i 1 H j ( x ) .
Remark 12.
Recalling (12), it follows from (45) that
B j i ( j + 1 ) ( 1 r ) = ( 1 ) j i j i 1 j + r i + r r ,
which holds for any integers 0 i j and r 0 . The following special cases can be obtained from (46):
  • B j ( j + 1 ) ( 1 r ) = ( 1 ) j j + r r r = ( 1 ) j r j ¯ . In particular, B j ( j + 1 ) = ( 1 ) j j ! .
  • B j 1 ( j + 1 ) ( 1 r ) = ( 1 ) j 1 ( j 1 ) ! H j ( r ) , j 1 .
  • B j i ( j + 1 ) = ( 1 ) j i j i 1 j + 1 i + 1 . In particular, B j ( j + 2 ) = ( 1 ) j j ! H j + 1 .
  • B j i ( j + 1 ) ( 1 ) = ( 1 ) j i j i 1 j i . In particular, B j ( j + 2 ) ( 1 ) = ( 1 ) j j ! j + 1 .
Let us now invoke the so-called harmonic polynomials  H j ( x ) , which were introduced by Cheon and El-Mikkawy in [5] (Section 5). They are defined by the generating function
j = 0 H j ( x ) t j = ln ( 1 t ) t ( 1 t ) 1 x ,
where H j ( 0 ) = H j + 1 . The first few harmonic polynomials are given explicitly by
H 0 ( x ) = 1 , H 1 ( x ) = x + 3 2 , H 2 ( x ) = 1 2 x 2 2 x + 11 6 , H 3 ( x ) = 1 6 x 3 + 5 4 x 2 35 12 x + 25 12 , H 4 ( x ) = 1 24 x 4 1 2 x 3 + 17 8 x 2 15 4 x + 137 60 , H 5 ( x ) = 1 120 x 5 + 7 48 x 4 35 36 x 3 + 49 16 x 2 203 45 x + 49 20 , H 6 ( x ) = 1 720 x 6 1 30 x 5 + 23 72 x 4 14 9 x 3 + 967 240 x 2 469 90 x + 363 140 , H 7 ( x ) = 1 5040 x 7 + 1 160 x 6 13 160 x 5 + 9 16 x 4 1069 480 x 3 + 801 160 x 2 29531 5040 x + 761 280 .
Upon comparing (4) and (47), it is clear that
H j ( x ) = H j + 1 ( 1 x ) , j 0 .
Therefore, by virtue of (48), we can write (44) as
B ν ( m + 2 ) ( x ) = ( 1 ) m ν ! d m ν d x m ν H m ( x ) ,
where now we assume that ν m . In particular,
B m ( m + 2 ) ( x ) = ( 1 ) m m ! H m ( x ) .
As with the hyperharmonic polynomials, the harmonic polynomials enjoy different representations [5,29]. Among them, we highlight the following one:
H m ( x ) = 1 m ! t = 0 m ( 1 ) t m + 1 t + 1 B t ( x ) ,
expressing the harmonic polynomials in terms of the Stirling numbers of the first kind and the Bernoulli polynomials. Therefore, combining (49) and (50) yields
B ν ( m + 2 ) ( x ) = ( 1 ) m m ν 1 t = m ν m ( 1 ) t t m v m + 1 t + 1 B t + ν m ( x ) .
Thus,
B ν ( m + 2 ) = ( 1 ) m m ν 1 t = m ν m ( 1 ) t t m v m + 1 t + 1 B t + ν m ,
and
B ν ( m + 2 ) ( 1 ) = ( 1 ) ν m ν 1 t = m ν m t m v m + 1 t + 1 B t + ν m .
It should be emphasized that, by using (51), one can obtain (at least in principle) the explicit expression of the generalized Bernoulli polynomial B ν ( α ) ( x ) for any ν 0 and arbitrary parameters α and x. For example, setting ν = 1 , , 6 in (51), we obtain in the first place
B 1 ( m + 2 ) ( x ) = x 1 2 ( m + 2 ) , B 2 ( m + 2 ) ( x ) = x 2 ( m + 2 ) x + 1 12 ( m + 2 ) ( 3 m + 5 ) , B 3 ( m + 2 ) ( x ) = x 3 3 2 ( m + 2 ) x 2 + 1 4 ( m + 2 ) ( 3 m + 5 ) x 1 8 ( m + 1 ) ( m + 2 ) 2 , B 4 ( m + 2 ) ( x ) = x 4 2 ( m + 2 ) x 3 + 1 2 ( m + 2 ) ( 3 m + 5 ) x 2 1 2 ( m + 1 ) ( m + 2 ) 2 x + 1 240 ( m + 2 ) ( 15 m 3 + 60 m 2 + 65 m + 12 ) , B 5 ( m + 2 ) ( x ) = x 5 5 2 ( m + 2 ) x 4 + 5 6 ( m + 2 ) ( 3 m + 5 ) x 3 5 4 ( m + 1 ) ( m + 2 ) 2 x 2 + 1 48 ( m + 2 ) ( 15 m 3 + 60 m 2 + 65 m + 12 ) x 1 96 ( m + 1 ) ( m + 2 ) 2 ( 3 m 2 + 5 m 4 ) , B 6 ( m + 2 ) ( x ) = x 6 3 ( m + 2 ) x 5 + 5 4 ( m + 2 ) ( 3 m + 5 ) x 4 5 2 ( m + 1 ) ( m + 2 ) 2 x 3 + 1 16 ( m + 2 ) ( 15 m 3 + 60 m 2 + 65 m + 12 ) x 2 1 16 ( m + 1 ) ( m + 2 ) 2 ( 3 m 2 + 5 m 4 ) x + 1 4032 ( m + 2 ) ( 63 m 5 + 315 m 4 + 315 m 3 539 m 2 938 m 240 ) .
Subsequently, making the transformation m α 2 in the above expressions for B ν ( m + 2 ) ( x ) , ν = 1 , , 6 , we obtain the corresponding generalized Bernoulli polynomials (cf. [45] (p. 143))
B 1 ( α ) ( x ) = x 1 2 α , B 2 ( α ) ( x ) = x 2 α x + 1 12 α ( 3 α 1 ) , B 3 ( α ) ( x ) = x 3 3 2 α x 2 + 1 4 α ( 3 α 1 ) x 1 8 α 2 ( α 1 ) , B 4 ( α ) ( x ) = x 4 2 α x 3 + 1 2 α ( 3 α 1 ) x 2 1 2 α 2 ( α 1 ) x + 1 240 α ( 15 α 3 30 α 2 + 5 α + 2 ) , B 5 ( α ) ( x ) = x 5 5 2 α x 4 + 5 6 α ( 3 α 1 ) x 3 5 4 α 2 ( α 1 ) x 2 + 1 48 α ( 15 α 3 30 α 2 + 5 α + 2 ) x 1 96 α 2 ( α 1 ) ( 3 α 2 7 α 2 ) , B 6 ( α ) ( x ) = x 6 3 α x 5 + 5 4 α ( 3 α 1 ) x 4 5 2 α 2 ( α 1 ) x 3 + 1 16 α ( 15 α 3 30 α 2 + 5 α + 2 ) x 2 1 16 α 2 ( α 1 ) ( 3 α 2 7 α 2 ) x + 1 4032 α ( 63 α 5 315 α 4 + 315 α 3 + 91 α 2 42 α 16 ) .
Note that, when α = 1 , the above generalized polynomials B ν ( α ) ( x ) , where ν = 1 , , 6 , become the ordinary Bernoulli polynomials
B 1 ( x ) = x 1 2 , B 2 ( x ) = x 2 x + 1 6 , B 3 ( x ) = x 3 3 2 x 2 + 1 2 x , B 4 ( x ) = x 4 2 x 3 + x 2 1 30 , B 5 ( x ) = x 5 5 2 x 4 + 5 3 x 3 1 6 x , B 6 ( x ) = x 6 3 x 5 + 5 2 x 4 1 2 x 2 + 1 42 .
Furthermore, when x = 0 , the above generalized polynomials B ν ( α ) ( x ) , where ν = 1 , , 6 , reduce to the following Bernoulli numbers of order α :
B 1 ( α ) = 1 2 α , B 2 ( α ) = 1 4 α 2 1 12 α , B 3 ( α ) = 1 8 α 3 + 1 8 α 2 , B 4 ( α ) = 1 16 α 4 1 8 α 3 + 1 48 α 2 + 1 120 α , B 5 ( α ) = 1 32 α 5 + 5 48 α 4 5 96 α 3 1 48 α 2 , B 6 ( α ) = 1 64 α 6 5 64 α 5 + 5 64 α 4 + 13 576 α 3 1 96 α 2 1 252 α .
Of course, by setting α = 1 in B ν ( α ) , ν = 1 , , 6 , we recover the ordinary Bernoulli numbers B 1 = 1 2 , B 2 = 1 6 , B 3 = 0 , B 4 = 1 30 , B 5 = 0 , and B 6 = 1 42 .

6. Further Identities Involving the r-Stirling Numbers of the First Kind

In this section, we will provide further identities involving the r-Stirling numbers of the first kind, the Bernoulli numbers and polynomials, the ordinary Stirling numbers of the first and second kind, and the harmonic numbers.
The following proposition exceptionally involves the r-Stirling numbers of the second kind k + r j + r r [31]. It provides a representation of the Bernoulli polynomials B k ( x ) evaluated at x = r (with r being a non-negative integer).
Proposition 2.
For any non-negative integers k , r , we have
B k ( r ) = ( 1 ) k k + r r r + j = 1 k ( 1 ) j + 1 ( j 1 ) ! j + 1 k + r j + r r ,
where it is understood that the summation on the right hand side of (52) is zero when k = 0 . In particular,
B k = δ k , 0 + ( 1 ) k + 1 j = 1 k ( 1 ) j ( j 1 ) ! j + 1 k j .
Proof. 
First note that, for x = r , (40) can be written as
j = 0 k k + r j + r r B j ( r ) = k ! H k + 1 ( 1 ) ,
whose Stirling transform is given by
B k ( r ) = j = 0 k ( 1 ) k j j ! k + r j + r r H j + 1 ( 1 ) .
On the other hand, it can be seen from (7) that
H j + 1 ( 1 ) = 1 j ( j + 1 ) , j 1 , 1 , j = 0 .
Therefore, by combining the last two equations, we arrive at (52). □
Remark 13.
The above representation (52) for B k ( r ) complements the formula for B k ( r ) given in (16) as well as the following formula established in [47] (Theorem 1.1):
B k ( r ) = j = 0 k ( 1 ) j j ! j + 1 k + r j + r r .
Moreover, as was noted in the introduction, Spieß [34], Wang [30], and Wuyungaowa [35] derived a number of identities involving the numbers j + r r P ( i , j + r , r ) . By making use of (34), we next restate some of these identities in terms of the r-Stirling numbers of the first kind.
  • From [34] (Theorem 10) and (34), we obtain
    k = 0 m ( 1 ) k k ! r m k k + r i + r r = ( 1 ) m m ! m i ,
    which holds for any integers r , i 0 . When r = 1 , (53) gives the recurrence m + 1 i + 1 = m i + m m i + 1 . On the other hand, when i = 1 , we can write (53) in terms of the hyperharmonic numbers H k ( r ) as
    k = 0 m ( 1 ) k r m k H k ( r ) = ( 1 ) m m , m 1 .
  • From [34] (Theorem 13) and (34), we obtain
    k = r m + 1 s r ! s ! k ! ( m + 1 k ) ! k r m + 1 k s = ( r + s ) ! ( m + 1 ) ! m + 1 r + s ,
    which holds for any integers r , s 0 . In particular, for s = 1 , we retrieve the identity in (32).
  • From [34] (Theorem 15) and (34), and after some rearrangements, we obtain
    q r k = 1 m 1 ( k + q ) ! k r = 1 q ! 1 ( m + q ) ! j = 1 r q j 1 m + 1 j ,
    which is valid for any integers 1 r m and q 1 . Setting r = 1 , 2 , and 3 in (54) gives
    k = 1 m 1 k ( k + 1 ) ( k + q ) = 1 q · q ! 1 q ( m + 1 ) ( m + q ) , k = 1 m H k 1 k ( k + 1 ) ( k + q ) = 1 q 2 · q ! 1 + q H m q 2 ( m + 1 ) ( m + q ) ,
    and
    k = 1 m H k 1 2 H k 1 [ 2 ] k ( k + 1 ) ( k + q ) = 2 q 3 · q ! 2 + 2 q H m + q 2 H m 2 H m [ 2 ] q 3 ( m + 1 ) ( m + q ) ,
    respectively. The above three identities are to be compared with the corresponding Examples 1, 2, and 3 previously obtained in [34] (p. 849). Moreover, it is worth pointing out that, for the case in which r = m , (54) yields the horizontal generating function for the Stirling numbers m + 1 j + 1 , namely,
    j = 0 m m + 1 j + 1 q j = ( q + 1 ) ( q + 2 ) ( q + m ) ,
    which holds for an arbitrary q.
  • From [30] (Equation (3.4)) and (34), we obtain
    k = 0 m ( 1 ) k k ! m k k + r r 1 k + r + 1 i + r + 1 r + 1 = ( 1 ) i m ! m + r r 1 m i .
    Setting here r = i = 1 , we obtain
    k = 0 m ( 1 ) k m k H k + 1 = 1 m ( m + 1 ) , m 1 ,
    which may be compared with the more commonly known identity (see, for example, [32] (Equation (3.2))) k = 0 m ( 1 ) k m k H k = 1 m , where m 1 .
  • From [30] (Equation (3.22)) and (34), we obtain
    m + r + 1 i + r + 1 r + 1 = m ! k = 0 m 1 k ! m k + r r k i = m ! k = 0 m ( 1 ) k i 1 k ! m + r k + r k i .
    In particular, for r = 0 , we have (cf. [30] (Equation (3.25)))
    m + 1 i + 1 = m ! k = 0 m 1 k ! k i = m ! k = 0 m ( 1 ) k i 1 k ! m k k i .
  • From [30] (Equations (3.32) and (3.33)) and (34), we obtain
    k = 0 m 1 k ! r + m k 1 m k k + s i + s s = 1 m ! m + r + s i + r + s r + s ,
    and
    m + r i + r r = m ! k = 0 m 1 k ! r + m k 2 m k k + 1 i + 1 ,
    respectively, where r and s are non-negative integers. Notice that, from (55) and (56), one quickly obtains
    k = 0 m 1 k ! r + m k 1 m k k + s i + s s = k = 0 m 1 k ! r + s + m k 2 m k k + 1 i + 1 .
    Furthermore, regarding (56), for r = 2 , it reads
    m + 2 i + 2 2 = m ! k = 0 m 1 k ! k + 1 i + 1 ,
    or, equivalently,
    P i H ( m , 1 ; 1 ) , H ( m , 2 ; 1 ) , , H ( m , i ; 1 ) = i ! m + 1 k = 0 m 1 k ! k + 1 i + 1 .
    In particular, when i = 1 , we recover the well-known identity k = 0 m H k = ( m + 1 ) H m m .
  • From [30] (Equations (4.1) and (4.3)) and (34), we obtain
    k = 0 m ( 1 ) m k m k k + r i + r r = m i r m i ,
    and
    k = i m ( r ) k i k i m + r k + r r = m i ,
    respectively, where i and r are non-negative integers. When r = 1 , (57) and (58) reduce to
    k = 0 m ( 1 ) m k m k k + 1 i + 1 = m i ,
    and
    k = i m ( 1 ) k i k i m + 1 k + 1 = m i ,
    respectively. Furthermore, when i = 1 and r = 2 , from (57) we obtain
    k = 0 m ( 1 ) m k ( k + 1 ) ! m k H k + 1 1 = m 2 m 1 .
    On the other hand, putting i = 1 in [30] (Equation (4.5)) yields
    k = 0 m ( 1 ) m k ( k + 1 ) ! m k = 2 m .
    Therefore, combining the last two identities, we obtain
    k = 0 m ( 1 ) m k ( k + 1 ) ! m k H k + 1 = ( m + 2 ) 2 m 1 ,
    which may be compared with [30] (Equation (4.10)).
  • From [30] (Equation (4.40)) and (34), we obtain
    k = 0 m ( 1 ) k i i ! k i m + r + 1 k + r + 1 r + 1 = m ! r + m i m i ,
    which holds for any integers 0 i m and r 0 . For r = 0 , the above identity becomes
    k = 0 m ( 1 ) k i i ! k i m + 1 k + 1 = m ! .
  • As we saw at the beginning of Section 5, the generalized Bernoulli numbers B k ( i ) are defined by the generating function
    t e t 1 i = k = 0 B k ( i ) t k k ! ,
    where B k ( 1 ) = B k corresponds to the ordinary Bernoulli numbers. In what follows, we restrict the superscript i to be a positive integer i 1 . Then, from [30] (Equations (5.2)) and (34), we obtain
    k = 0 m ( 1 ) k m + r k + r r B k ( i ) = m + i i 1 m + i + r i + r r .
    In particular, for i = 1 , we obtain
    k = 0 m ( 1 ) k m + r k + r r B k = 1 m + 1 m + r + 1 r + 1 r ,
    which, for r = 1 , reduces to the well-known identity (see, e.g., [30] (Equation (5.6)))
    k = 0 m ( 1 ) k m + 1 k + 1 B k = m ! H m + 1 .
    Moreover, substituting the representation of B k ( i ) given by Kim et al. [48]
    B k ( i ) = j = 0 k ( 1 ) j i + j i 1 i + j i k j ,
    into (59), we obtain
    k = 0 m j = 0 k ( 1 ) k + j i + j i 1 i + j i k j m + r k + r r = m + i i 1 m + i + r i + r r .
    In particular, setting i = 2 and r = 0 in (60) leads to
    k = 0 m j = 0 k ( 1 ) k + j j ! j + 2 k j m k H j + 1 = m ! m + 2 H m + 1 .
    Likewise, employing the representation of B k ( i ) given in [43] (Equation (15)),
    B k ( i ) = j = 0 k ( 1 ) j k + i i + j i + j 1 i 1 k + j j 1 k + j j ,
    we obtain from (59)
    k = 0 m j = 0 k ( 1 ) k + j k + i i + j i + j 1 i 1 k + j j k + j j m + r k + r r = m + i i 1 m + i + r i + r r .
    In particular, when i = 3 and r = 0 , (61) implies the relation
    k = 0 m j = 0 k ( 1 ) k + j ( j + 1 ) ( j + 2 ) k + 3 j + 3 k + j j k + j j m k = 6 · m ! m + 3 H m + 2 2 H m + 2 [ 2 ] .
  • From [30] (Equations (5.3)) and (34), we obtain
    k = 0 m m + r + 1 k + r + 1 r + 1 B k ( i + 1 ) = m ! m + i + r + 1 i + r H m + i + r + 1 H i + r .
    Taking r r 1 and i i 1 in (62) gives
    k = 0 m m + r k + r r B k ( i ) = m ! m + i + r 1 i + r 2 H m + i + r 1 H i + r 2 = m ! H m + 1 ( i + r 1 ) ,
    which is just the identity (40) with x replaced by i.
  • In [30] (p. 1508), we find the identity
    k = 0 m P ( k , m + r + i , r + i ) B k ( i ) k ! = r + i i m + i i 1 P ( i , m + r + i , r ) .
    Using (34), and after some minor manipulations, we can write this as
    k = 0 m m + r + i k + r + i r + i B k ( i ) = m + i i 1 m + i + r i + r r .
    In particular, if r = i = 1 here, then we have
    k = 0 m P k H ( m , 1 ; 1 ) , H ( m , 2 ; 1 ) , , H ( m , k ; 1 ) B k k ! = H m + 1 m + 1 .
    Furthermore, from (59) and (63), we find
    k = 0 m ( 1 ) k m + r k + r r B k ( i ) = k = 0 m m + r + i k + r + i r + i B k ( i ) .
    Lastly, it is to be noted that, by setting i = 1 in (59) and (63) and (66), and comparing the resulting equations with (15), we obtain the double identity
    k = 0 m ( 1 ) k m + r k + r r B k = k = 0 m m + r + 1 k + r + 1 r + 1 B k = k = 0 m m + 1 k + 1 B k ( r ) ,
    which holds for any non-negative integers m and r.
  • Combining the recurrence of the numbers j + r r P ( i , j + r , r ) appearing in [30] (p. 1505)
    j + r r P ( i , j + r , r ) = j + r 1 r 1 P ( i , j + r 1 , r 1 ) + j + r 1 r P ( i , j + r 1 , r )
    as well as (34), we readily obtain the corresponding recurrence for the r-Stirling numbers of the first kind:
    j + r + 1 i + r + 1 r + 1 = j + r i + r r + j j + r i + r + 1 r + 1 ,
    which may be compared with (33).
  • From [35] (Equation (4)) and (34), we obtain the identity
    j = i m 1 m + 1 j 1 ( j 1 ) ! j + r i + r r + 1 = i m ! m + r + 1 i + r + 1 r + 1 ,
    which holds for any integers 1 i m and r 0 . Note the close resemblance of this identity to (31). In particular, for r = 0 , we find
    j = i m 1 m + 1 j 1 ( j 1 ) ! j i = i m ! m + 1 i + 1 ,
    which may be compared with (32).

7. Conclusions

In this paper, we have unified and generalized some previous, unconnected results obtained in [12,30] (see also [21,22,23,29]) concerning the higher-order derivatives of the hyperharmonic polynomials and their relationship with the complete Bell polynomials. Specifically, we have fully characterized the hyperharmonic polynomials and their successive derivatives in terms of the r-Stirling polynomials of the first kind. In particular, when evaluated at some non-negative integer r, such derivatives can be expressed in terms of the r-Stirling numbers of the first kind or, equivalently, in terms of the numbers P ( i , j + r , r ) , which, in turn, can be expressed in terms of the complete Bell polynomials. Moreover, by exploiting the link between the generalized Bernoulli polynomials B n ( α ) ( x ) and the higher-order derivatives of the hyperharmonic polynomials, we have derived a new formula that enables us to obtain B n ( α ) ( x ) for any non-negative integer n. Finally, we have provided a number of identities involving the r-Stirling numbers of the first kind, the Bernoulli numbers and polynomials, the ordinary Stirling numbers of both kinds, and the harmonic numbers.
It is to be noted that some of the identities obtained in this paper may be new, including the identity considered in the following proposition.
Proposition 3.
For any integers 0 i m and r 1 , we have
j = 0 m 1 j ! j + r i + r r = ( 1 ) i m ! r i + 1 k = 0 i ( r ) k m + r + 1 k + r r .
Proof. 
We start with the following elementary identity
j = 0 m x + j + r j = 1 + m x + r + 1 x + m + r m .
On the other hand, from (10), we have
d i d x i x + j + r j x = 0 = i ! j ! j + r + 1 i + r + 1 r + 1 .
Hence, it follows that
j = 0 m i ! j ! j + r + 1 i + r + 1 r + 1 = d i d x i 1 + m x + r + 1 x + m + r m x = 0 .
Now, applying the Leibniz rule for the i-th derivative of the product of two functions yields
d i d x i 1 + m x + r + 1 x + m + r m x = 0 = k = 0 i i k d i k d x i k 1 + m x + r + 1 x = 0 d k d x k x + m + r m x = 0 = 1 + m r + 1 i ! m ! m + r + 1 i + r + 1 r + 1 + k = 0 i 1 i k d i k d x i k 1 + m x + r + 1 x = 0 d k d x k x + m + r m x = 0 = 1 + m r + 1 i ! m ! m + r + 1 i + r + 1 r + 1 + k = 0 i 1 ( 1 ) i k i ! ( m 1 ) ! ( r + 1 ) i k + 1 m + r + 1 k + r + 1 r + 1 .
Thus, from (65) and (66), we have
j = 0 m 1 j ! j + r + 1 i + r + 1 r + 1 = 1 + m r + 1 1 m ! m + r + 1 i + r + 1 r + 1 + k = 0 i 1 ( 1 ) i k ( m 1 ) ! ( r + 1 ) i k + 1 m + r + 1 k + r + 1 r + 1 ,
which can be expressed more compactly as
j = 0 m 1 1 j ! j + r + 1 i + r + 1 r + 1 = k = 0 i ( 1 ) i k ( m 1 ) ! ( r + 1 ) i k + 1 m + r + 1 k + r + 1 r + 1 .
Finally, making m m + 1 and r r 1 in the last equation, we obtain (64). □
Remark 14.
For i = 0 , 1 , and m, (64) leads to the relations
j = 0 m r + j 1 j = r + m m , j = 0 m H j ( r ) = 1 m ! r m + r + 1 r + 1 r 1 r r + m m ,
and
r m + 1 = ( 1 ) m k = 0 m ( r ) k m + r + 1 k + r r ,
respectively, where r stands for any arbitrary integer r 1 .

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The author would like to express their deep appreciation to Guest Editor Taras Goy. The author also gratefully acknowledges the anonymous referees for useful comments that enhanced the quality of this paper.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Adamchik, V. On Stirling numbers and Euler sums. J. Comput. Appl. Math. 1997, 79, 119–130. [Google Scholar] [CrossRef] [Green Version]
  2. Paule, P.; Schneider, C. Computer proofs of a new family of harmonic number identities. Adv. Appl. Math. 2003, 31, 359–378. [Google Scholar] [CrossRef] [Green Version]
  3. Zheng, D.-Y. Further summation formulae related to generalized harmonic numbers. J. Math. Anal. Appl. 2007, 335, 692–706. [Google Scholar] [CrossRef] [Green Version]
  4. Cheon, G.-S.; El-Mikkawy, M.E.A. Generalized harmonic number identities and a related matrix representation. J. Korean Math. Soc. 2007, 44, 487–498. [Google Scholar] [CrossRef] [Green Version]
  5. Cheon, G.-S.; El-Mikkawy, M.E.A. Generalized harmonic numbers with Riordan arrays. J. Number Theory 2008, 128, 413–425. [Google Scholar] [CrossRef] [Green Version]
  6. Feng, C.-J.; Zhao, F.-Z. Some results for generalized harmonic numbers. Integers 2009, 9, 605–619. [Google Scholar] [CrossRef]
  7. Coffey, M.W.; Lubbers, N. On generalized harmonic number sums. Appl. Math. Comput. 2010, 217, 689–698. [Google Scholar] [CrossRef]
  8. Choi, J.; Srivastava, H.M. Some summation formulas involving harmonic numbers and generalized harmonic numbers. Math. Comput. Model. 2011, 54, 2220–2234. [Google Scholar] [CrossRef]
  9. Liang, H. Identities involving generalized harmonic numbers and other special combinatorial sequences. J. Integer Seq. 2012, 15, 12–19. [Google Scholar]
  10. Choi, J. Finite summation formulas involving binomial coefficients, harmonic numbers and generalized harmonic numbers. J. Inequal. Appl. 2013, 2013, 49. [Google Scholar] [CrossRef] [Green Version]
  11. Wang, J.; Wei, C. Derivative operator and summation formulae involving generalized harmonic numbers. J. Math. Anal. Appl. 2016, 434, 315–341. [Google Scholar] [CrossRef]
  12. Kargın, L.; Cenkci, M.; Dil, A.; Can, M. Generalized harmonic numbers via poly-Bernoulli polynomials. arXiv 2020, arXiv:2008.00284. To appear in Publ. Math. Debrecen (private communication from L. Kargın to the present author). [Google Scholar]
  13. Dil, A.; Muniroǧlu, E. Applications of derivative and difference operators on some sequences. Appl. Anal. Discret. Math. 2020, 14, 406–430. [Google Scholar] [CrossRef]
  14. Guo, D. Some combinatorial identities concerning harmonic numbers and binomial coefficients. Discret. Math. Lett. 2022, 8, 41–48. [Google Scholar] [CrossRef]
  15. Flajolet, P.; Salvy, B. Euler sums and contour integral representations. Exp. Math. 1998, 7, 15–35. [Google Scholar] [CrossRef]
  16. Choi, J.; Srivastava, H.M. Explicit evaluation of Euler and related sums. Ramanujan J. 2005, 10, 51–70. [Google Scholar] [CrossRef]
  17. Mezo, I.; Dil, A. Hyperharmonic series involving Hurwitz zeta function. J. Number Theory 2010, 130, 360–369. [Google Scholar] [CrossRef] [Green Version]
  18. Szabłowski, P.J. On connection between values of Riemann zeta function at rationals and generalized harmonic numbers. Math. Aeterna 2015, 5, 583–592. [Google Scholar]
  19. Hajli, M. On a formula for the regularized determinant of zeta functions with application to some Dirichlet series. Q. J. Math. 2020, 71, 843–865. [Google Scholar] [CrossRef]
  20. Conway, J.H.; Guy, R.K. The Book of Numbers; Copernicus: New York, NY, USA, 1996. [Google Scholar]
  21. Mezo, I. Analytic extension of hyperharmonic numbers. Online J. Anal. Comb. 2009, 4, 1. [Google Scholar]
  22. Dil, A. On the hyperharmonic function. Süleyman Demirel Univ. J. Nat. Appl. Sci. 2019, 23, 187–193. [Google Scholar] [CrossRef]
  23. Benjamin, A.T.; Gaebler, D.; Gaebler, R. A combinatorial approach to hyperharmonic numbers. Integers 2003, 3, #A15. [Google Scholar]
  24. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F.G. Higher Transcendental Functions; McGraw-Hill Book Company: New York, NY, USA, 1953; Volume I. [Google Scholar]
  25. Quaintance, J.; Gould, H.W. Combinatorial Identities for Stirling Numbers. The Unpublished Notes of H.W. Gould; World Scientific Publishing: Singapore, 2016. [Google Scholar]
  26. Cheon, G.-S.; Hwang, S.-G.; Lee, S.-G. Several polynomials associated with the harmonic numbers. Discret. Appl. Math. 2007, 155, 2573–2584. [Google Scholar] [CrossRef] [Green Version]
  27. Dil, A.; Kurt, V. Polynomials related to harmonic numbers and evaluation of harmonic number series I. Integers 2012, 12, #A38. [Google Scholar] [CrossRef] [Green Version]
  28. Dil, A.; Kurt, V. Polynomials related to harmonic numbers and evaluation of harmonic number series II. Appl. Anal. Discret. Math. 2011, 5, 212–229. [Google Scholar] [CrossRef] [Green Version]
  29. Cereceda, J.L. Daehee, hyperharmonic, and power sum polynomials. arXiv 2021, arXiv:2108.03651. [Google Scholar]
  30. Wang, W. Riordan arrays and harmonic numbers identities. Comput. Math. Appl. 2010, 60, 1494–1509. [Google Scholar] [CrossRef] [Green Version]
  31. Broder, A.Z. The r-Stirling numbers. Discret. Math. 1984, 49, 241–259. [Google Scholar] [CrossRef] [Green Version]
  32. Wang, W.; Jia, C. Harmonic number identities via the Newton–Andrews method. Ramanujan J. 2014, 35, 263–285. [Google Scholar] [CrossRef]
  33. Kölbig, K.S. The complete Bell polynomials for certain arguments in terms of Stirling numbers of the first kind. J. Comput. Appl. Math. 1994, 51, 113–116. [Google Scholar] [CrossRef] [Green Version]
  34. Spieß, J. Some identities involving harmonic numbers. Math. Comp. 1990, 55, 839–863. [Google Scholar] [CrossRef] [Green Version]
  35. Wuyungaowa, S.W. Sums of involving the harmonic numbers and the binomial coefficients. Amer. J. Comput. Math. 2015, 5, 96–105. [Google Scholar] [CrossRef] [Green Version]
  36. Boyadzhiev, K.N. Notes on the Binomial Transform: Theory and Table with Appendix on Stirling Transform; World Scientific: Singapore, 2018. [Google Scholar]
  37. Mezo, I. Combinatorics and Number Theory of Counting Sequences; CRC Press Taylor & Francis Group: Boca Raton, FL, USA, 2020. [Google Scholar]
  38. Graham, R.L.; Knuth, D.E.; Patashnik, O. Concrete Mathematics: A Foundation for Computer Science; Addison-Wesley: Reading, MA, USA, 1989. [Google Scholar]
  39. Cereceda, J.L. Iterative procedure for hypersums of powers of integers. J. Integer Seq. 2014, 17, 14–15. [Google Scholar]
  40. Carlitz, L. Weighted Stirling numbers of the first and second kind—I. Fibonacci Q. 1980, 18, 147–162. [Google Scholar]
  41. Kargın, L.; Can, M. Harmonic number identities via polynomials with r-Lah coefficients. C. R. Math. Acad. Sci. Paris 2020, 358, 535–550. [Google Scholar] [CrossRef]
  42. Comtet, L. Advanced Combinatorics: The Art of Finite and Infinite Expansions; Reidel Publishing Company: Dordrecht, The Netherlands, 1974. [Google Scholar]
  43. Srivastava, H.M.; Todorov, P.G. An explicit formula for the generalized Bernoulli polynomials. J. Math. Anal. Appl. 1988, 130, 509–513. [Google Scholar] [CrossRef] [Green Version]
  44. Liu, G.-D.; Srivastava, H.M. Explicit formulas for the Nörlund polynomials B n ( x ) and b n ( x ) . Comput. Math. Appl. 2006, 51, 1377–1384. [Google Scholar] [CrossRef] [Green Version]
  45. Elezović, N. Generalized Bernoulli polynomials and numbers, revisited. Mediterr. J. Math. 2016, 13, 141–151. [Google Scholar] [CrossRef]
  46. Coffey, M.W. Series representations of the Riemann and Hurwitz zeta functions and series and integral representations of the first Stieltjes constant. arXiv 2011, arXiv:1106.5146. [Google Scholar]
  47. Guo, B.-N.; Mezo, I.; Qi, F. An explicit formula for Bernoulli polynomials in terms of r-Stirling numbers of the second kind. Rocky Mt. J. Math. 2016, 46, 1919–1923. [Google Scholar] [CrossRef] [Green Version]
  48. Kim, D.S.; Kim, T.; Lee, S.-H.; Seo, J.-J. Higher-order Daehee numbers and polynomials. Int. J. Math. Anal. Ruse 2014, 8, 273–283. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cereceda, J.L. Note on the Higher-Order Derivatives of the Hyperharmonic Polynomials and the r-Stirling Polynomials of the First Kind. Axioms 2022, 11, 167. https://doi.org/10.3390/axioms11040167

AMA Style

Cereceda JL. Note on the Higher-Order Derivatives of the Hyperharmonic Polynomials and the r-Stirling Polynomials of the First Kind. Axioms. 2022; 11(4):167. https://doi.org/10.3390/axioms11040167

Chicago/Turabian Style

Cereceda, José L. 2022. "Note on the Higher-Order Derivatives of the Hyperharmonic Polynomials and the r-Stirling Polynomials of the First Kind" Axioms 11, no. 4: 167. https://doi.org/10.3390/axioms11040167

APA Style

Cereceda, J. L. (2022). Note on the Higher-Order Derivatives of the Hyperharmonic Polynomials and the r-Stirling Polynomials of the First Kind. Axioms, 11(4), 167. https://doi.org/10.3390/axioms11040167

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop