Next Article in Journal
Forbidden Pairs of Disconnected Graphs for Traceability of Block-Chains
Next Article in Special Issue
Optimal Plan of Multi-Stress–Strength Reliability Bayesian and Non-Bayesian Methods for the Alpha Power Exponential Model Using Progressive First Failure
Previous Article in Journal
Recent Advances in Surrogate Modeling Methods for Uncertainty Quantification and Propagation
Previous Article in Special Issue
An Alternate Generalized Odd Generalized Exponential Family with Applications to Premium Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Convolutions for Bernoulli and Euler–Genocchi Polynomials of Order (r,m) and Their Probabilistic Interpretation †

by
Robert Frontczak
1,‡ and
Živorad Tomovski
2,*
1
Landesbank Baden-Württemberg (LBBW), 70173 Stuttgart, Germany
2
Department of Mathematics, Faculty of Science, University of Ostrava, 70103 Ostrava, Czech Republic
*
Author to whom correspondence should be addressed.
Dedicated to Prof. Hari M. Srivastava on the Occasion of his 80th Birthday.
Statements and conclusions made in this article by R. F. are entirely those of the author. They do not necessarily reflect the views of LBBW.
Symmetry 2022, 14(6), 1220; https://doi.org/10.3390/sym14061220
Submission received: 12 May 2022 / Revised: 1 June 2022 / Accepted: 6 June 2022 / Published: 13 June 2022
(This article belongs to the Special Issue Symmetric Distributions, Moments and Applications)

Abstract

:
The main purpose of this article is to derive several convolutions for generalized Bernoulli and Euler–Genocchi polynomials of order ( r , m ) , B n ( r , m ) ( x ) and A n ( r , m ) ( x ) , respectively. These polynomials have been introduced recently and contain the generalized Bernoulli, Euler and Genocchi polynomials as special members. Some of our results extend the results of M. Merca and others concerning Bernoulli numbers and polynomials. Probabilistic interpretations of the presented results are also given.

1. Introduction

Bernoulli numbers ( B n ) n 0 are defined by
B ( z ) = n = 0 B n z n n ! = z e z 1 ( | z | < 2 π ) .
A generalization of Bernoulli numbers are Bernoulli polynomials B n ( x ) , x C , defined by
B ( x , z ) = n = 0 B n ( x ) z n n ! = z e z 1 e x z ( | z | < 2 π ) .
These numbers (polynomials) are fascinating objects, appearing in many mathematical branches such as number theory, combinatorics and analysis. The basic properties of Bernoulli numbers and polynomials are discussed in [1,2].
Closely related to Bernoulli polynomials are the Euler and Genocchi polynomials. These polynomials are defined for | z | < π by
n = 0 E n ( x ) z n n ! = 2 e z + 1 e x z a n d n = 0 G n ( x ) z n n ! = 2 z e z + 1 e x z .
Finding recurrences and convolutions for these polynomials is still an active field of research. Many interesting identities for Bernoulli, Euler and Genocchi polynomials can be found in the articles [3,4,5,6,7,8,9,10] for instance. See [11,12,13,14] for some properties of generalizations of these polynomials.
The popularity and importance of Bernoulli numbers and polynomials in number theory comes also from their connection to the Riemann zeta function
ζ ( 2 n ) = ( 1 ) n + 1 ( 2 π ) 2 n 2 · ( 2 n ) ! B 2 n ,
where
ζ ( s ) = n = 1 1 n s , ( s ) > 1 ,
is the Riemann zeta function [15]. A great deal of proof for this relation has been provided over the years. See [16] for references. Recently, Merca [16] proved the following relation between Bernoulli numbers:
k = 0 n n k 2 k B k = ( 2 2 n ) B n .
This relation can be used to derive a recurrence relation for ζ ( 2 n ) . Moreover, in his next article on the topic, Merca [17] used recurrence relations for Bernoulli polynomials B n ( x ) to derive two new infinite families of linear recurrence relations for the Riemann zeta function at positive even integer arguments. Merca’s elegant results are based on the following relations (Theorems 2.1 and 3.1 in [17]): Let n be a positive integer and x , α C . Then
k = 0 n n k B k ( x ) ( α n k ( 1 ) k ( 2 x 1 + α ) n k ) = 0 ,
and
k = 0 n 2 n 2 k α 2 k B n 2 k ( x ) = k = 0 n n k ( 1 ) n k ( 2 x 1 + α ) k + ( 2 x 1 α ) k 2 B n k ( x ) .
In fact, identities Equations (1) and (2) have been widely known and used already in the 19th century as immediate consequences of the simple functional relations
B ( 2 z ) ( e z + 1 ) = 2 B ( z ) a n d B ( x , z ) e y z = B ( x , z ) e ( 2 x 1 + y ) z .
It may be very difficult to identify the first pioneers who discovered them after all this time. We refer to the notable books by Saalschütz [18], Nielsen [19] and Hansen [20] as a resource in which a large number of classical identities, mainly developed in the 18th and the 19th centuries, can be found. In addition, Equation (1) can be extended to [21]
k = 0 n n k m k B k i = 0 m 1 i n k = m B n ,
valid for all integers m 1 . Thus, Equation (1) is the special case for m = 2 . Next, we can define the function Φ n ( m ) ( x , y ) , m 1 , x , y C , by
Φ n ( m ) ( x , y ) = k = 0 n n k m k B k ( x ) y n k .
Then, the following identity has been known for a long time:
Φ n ( m ) ( x , y ) = k = 0 n n k ( m ) k B k ( x ) ( 2 m x m + y ) n k , m 1 .
The identity Equation (4) reduces to Equation (2) when m = 1 and y = α . It is easy to see that Φ n ( m ) ( x , y ) satisfies the following functional equation.
Proposition 1.
For all m , n N the following functional equation holds:
Φ n ( m ) ( x , ( 2 m x m + y ) ) = ( 1 ) n Φ n ( m ) ( x , y ) .
In addition, we can examine a calculation of the sum Φ n ( m ) ( x , y ) + Φ n ( m ) ( x , y ) . We have
Φ n ( m ) ( x , y ) + Φ n ( m ) ( x , y ) = 2 l = 0 n / 2 n 2 l m n 2 l B n 2 l ( x ) y 2 l , i f   n k = 2 l ; 0 , o t h e r w i s e .
In view of such an identity, we recognize that Equation (3) is an obvious consequence of Equation (2). Finally, we remark that if n k is odd, then Φ n ( m ) ( x , y ) = Φ n ( m ) ( x , y ) , i.e., the function Φ n ( m ) ( x , y ) is an odd or asymmetric function with respect to y.
We conclude this section by recalling the definition of B n ( r , m ) ( x ) , which comes from [22].
Definition 1.
For integers r , m 1 the generalized Bernoulli polynomials of order ( r , m ) are defined by the generating function
B ( r , m ; x , z ) = n = 0 B n ( r , m ) ( x ) z n n ! = z r e z 1 m e x z ( | z | < 2 π ) .
The numbers B n ( r , m ) ( 0 ) = B n ( r , m ) are called generalized Bernoulli numbers of order r and m.
The polynomials B n ( r , m ) ( x ) belong to the family of Appell polynomials. We mention that they are defined for n m ( r 1 ) , i.e., B j ( r , m ) ( x ) = 0 for j < m ( r 1 ) . From the definition, it is obvious that B n ( 1 , 1 ) ( x ) = B n ( x ) and B n ( 1 , m ) ( x ) = B n ( m ) ( x ) are the generalized Bernoulli polynomials of order m. Moreover, B n ( 1 , 1 ) = B n are the Bernoulli numbers.
The goal of the present article is to derive several convolutions for generalized Bernoulli polynomials of order ( r , m ) , B n ( r , m ) ( x ) . First, we will generalize Merca’s results for Bernoulli polynomials to the more general class of polynomials. This will be performed in Section 2. In Section 3, we will prove the analogue identities for generalized Euler–Genocchi polynomials A n ( r , m ) ( x ) . A range of additional convolutions for B n ( r , m ) ( x ) and A n ( r , m ) ( x ) will be given in Section 4. Among other things, we will rediscover identity Equation (1) as a special case of our findings. In Section 5 and Section 6 we will state some additional remarks concerning applications and future work.

2. Notes on Merca’s Identities

Our first result is an extension of Theorem 2.1 of [17].
Theorem 1.
Let r and m be positive integers and x , y C . Then,
k = 0 n n k B k ( r , m ) ( x ) ( y n k ( 1 ) m ( r 1 ) + k ( 2 x m + y ) n k ) = 0 .
Especially, with y = m , we have
k = 0 n n k B k ( r , m ) ( x ) ( m n k ( 1 ) m ( r 1 ) + k ( 2 x ) n k ) = 0 .
Proof. 
For | z | < 2 π we have from Equation (6)
n = 0 B n ( r , m ) ( x ) z n n ! n = 0 y n z n n ! = z r e z 1 m e x z e y z = z r m ( 1 ) m ( e z 1 ) m e x z e ( y + 2 x m ) z = ( 1 ) m ( r 1 ) ( z ) r e z 1 m e x z e ( y + 2 x m ) z = ( 1 ) m ( r 1 ) B ( r , m ; x , z ) e ( y + 2 x m ) z = ( 1 ) m ( r 1 ) n = 0 B n ( r , m ) ( x ) ( 1 ) n z n n ! n = 0 ( y + 2 x m ) n z n n ! = n = 0 k = 0 n n k ( 1 ) m ( r 1 ) + k B k ( r , m ) ( x ) ( y + 2 x m ) n k z n n ! .
On the other hand, using Cauchy’s rule, it is obvious that
n = 0 B n ( r , m ) ( x ) z n n ! n = 0 y n z n n ! = n = 0 k = 0 n n k B k ( r , m ) ( x ) y n k z n n ! .
Comparing the coefficients for z n in the two power series proves the formula. □
For ( r , m ) = ( 1 , 1 ) Theorem 1 reduces to Merca’s Theorem 2.1. For ( r , m ) = ( 1 , m ) our theorem gives a convolutional relation for generalized Bernoulli polynomials:
Corollary 1.
For x , y C the following relation holds for generalized Bernoulli polynomials:
k = 0 n n k B k ( m ) ( x ) ( y n k + ( 1 ) k + 1 ( 2 x m + y ) n k ) = 0 .
Corollary 2.
For x C the following relation holds:
B n ( r , m ) ( x ) = k = 0 n n k ( 1 ) m ( r 1 ) + k B k ( r , m ) ( x ) ( 2 x m ) n k .
Especially for the generalized Bernoulli polynomials we have the identity
B n ( m ) ( x ) = k = 0 n n k ( 1 ) k B k ( m ) ( x ) ( 2 x m ) n k .
Proof. 
Set y = 0 in Theorem 1. □
Next, for x , y C , we define the function
F n ( r , m ) ( x , y ) = k = 0 n n k B k ( r , m ) ( x ) y n k , r , m 1 .
Then, we note the following functional equation:
Proposition 2.
For all r , m , n N the following functional equation holds:
F n ( r , m ) ( x , ( 2 m x m + y ) ) = ( 1 ) n + m ( r 1 ) F n ( r , m ) ( x , y ) .
Proof. 
Replacing y with ( 2 x m + y ) in Equation (12), using Equation (7), we get the result. □
We can also calculate
F n ( r , m ) ( x , y ) + F n ( r , m ) ( x , y ) = 2 l = 0 n / 2 n 2 l B n 2 l ( r , m ) ( x ) y 2 l , i f   n k = 2 l ; 0 , o t h e r w i s e ,
and
F n ( r , m ) ( x , y ) F n ( r , m ) ( x , y ) = 2 l = 0 ( n 1 ) / 2 n 2 l + 1 B n 2 l 1 ( r , m ) ( x ) y 2 l + 1 , i f   n k = 2 l + 1 ; 0 , o t h e r w i s e .
These calculations confirm the following facts.
(i)
If n k is odd, then F n ( r , m ) ( x , y ) = F n ( r , m ) ( x , y ) , i.e., the function F n ( r , m ) ( x , y ) is an odd or asymmetric function with respect to y.
(ii)
If n k is even, then F n ( r , m ) ( x , y ) = F n ( r , m ) ( x , y ) , i.e., the function F n ( r , m ) ( x , y ) is an even or symmetric function with respect to y.
The above observations lead to the next corollary, which provides an extension of Theorem 3.1 of [17] to the class B n ( r , m ) ( x ) .
Corollary 3.
Let r and m be positive integers and x , y C . Then,
k = 0 n 2 n 2 k y 2 k B n 2 k ( r , m ) ( x ) = ( 1 ) m ( r 1 ) k = 0 n n k ( 1 ) n k ( 2 x m + y ) k + ( 2 x m y ) k 2 B n k ( r , m ) ( x )
and
k = 0 n 1 2 n 2 k + 1 y 2 k + 1 B n 2 k 1 ( r , m ) ( x ) = ( 1 ) m ( r 1 ) k = 0 n n k ( 1 ) n k ( 2 x m + y ) k ( 2 x m y ) k 2 B n k ( r , m ) ( x ) .
We conclude this section with the following results.
Corollary 4.
Let r and m be positive integers and x , y C . Then,
k = 0 n 2 n 2 k ( 1 ) k y 2 k B n 2 k ( r , m ) ( x ) = ( 1 ) m ( r 1 ) k = 0 n n k ( 1 ) n k B n k ( r , m ) ( x ) p = 0 k k p y p cos π p 2 ( 2 x m ) k p
and
k = 0 n 1 2 n 2 k + 1 ( 1 ) k y 2 k + 1 B n 2 k 1 ( r , m ) ( x ) = ( 1 ) m ( r 1 ) k = 0 n n k ( 1 ) n k B n k ( r , m ) ( x ) p = 0 k k p y p sin π p 2 ( 2 x m ) k p .
Proof. 
Replace y by i y with i = 1 in Corollary 3 and simplify. □
Corollary 5.
Let r and m be positive integers and x C . Then,
k = 0 n 2 n 2 k ( 1 ) k m 2 k B n 2 k ( r , m ) ( x ) = ( 1 ) m ( r 1 ) k = 0 n n k ( 1 ) n k 2 k 1 ( x k + ( x m ) k ) B n k ( r , m ) ( x )
and
k = 0 n 1 2 n 2 k + 1 ( 1 ) k m 2 k + 1 B n 2 k 1 ( r , m ) ( x ) = ( 1 ) m ( r 1 ) k = 0 n n k ( 1 ) n k 2 k 1 ( x k ( x m ) k ) B n k ( r , m ) ( x ) .
Proof. 
Set y = m in Corollary 3 and simplify. □

3. Analogue Relations for Generalized Euler–Genocchi Polynomials

The definition of generalized Euler–Genocchi polynomials of order ( r , m ) also comes from the paper [22], where many basic properties of the polynomials are discussed.
Definition 2.
Let r and m be integers with r 0 and m 1 . The generalized Euler–Genocchi polynomials of order ( r , m ) , A n ( r , m ) ( x ) , x C , are defined by the generating function
A ( r , m ; x , z ) = n = 0 A n ( r , m ) ( x ) z n n ! = 2 z r e z + 1 m e x z ( | z | < π )
with A j ( r , m ) ( x ) = 0 for j < r m . The numbers A n ( r , m ) ( 0 ) = A n ( r , m ) are called the generalized Euler–Genocchi numbers of order r and m.
We see that A n ( 0 , m ) ( x ) = E n ( m ) ( x ) and A n ( 1 , m ) ( x ) = G n ( m ) ( x ) are the generalized Euler and Genocchi polynomials, respectively, where
n = 0 E n ( m ) ( x ) z n n ! = 2 e z + 1 m e x z ( | z | < π )
and
n = 0 G n ( m ) ( x ) z n n ! = 2 z e z + 1 m e x z ( | z | < π ) .
Finally, we mention that the degenerated case m = 0 gives A n ( r , 0 ) ( x ) = x n for all r 0 .
The first analogue result of Merca’s identities is stated in the next theorem.
Theorem 2.
Let r and m be integers with r 0 and m 1 , and x , y C . Then,
k = 0 n n k A k ( r , m ) ( x ) ( y n k ( 1 ) m r + k ( 2 x m + y ) n k ) = 0 .
Especially, with y = m , we have
k = 0 n n k A k ( r , m ) ( x ) ( m n k ( 1 ) m r + k ( 2 x ) n k ) = 0 .
Proof. 
Due to the high degree of similarity in the proofs, we only sketch the proofs. The identity basically follows from
A ( r , m ; x , z ) e y z = 2 z r e z + 1 m e x z e y z = ( 1 ) m r 2 ( z ) r e z + 1 m e x z e ( 2 x + y m ) z .
Corollary 6.
For x , y C the following relations hold:
k = 0 n n k E k ( m ) ( x ) ( y n k + ( 1 ) k + 1 ( 2 x m + y ) n k ) = 0
and
k = 0 n n k G k ( m ) ( x ) ( y n k + ( 1 ) m + k + 1 ( 2 x m + y ) n k ) = 0 .
Proof. 
Set r = 0 and r = 1 in Theorem 2. □
Corollary 7.
For x C the following relation holds:
A n ( r , m ) ( x ) = k = 0 n n k ( 1 ) m r + k A k ( r , m ) ( x ) ( 2 x m ) n k .
Proof. 
Set y = 0 in Theorem 2. □
Theorem 3.
Let r and m be integers with r 0 and m 1 , and x , y C . Then,
k = 0 n 2 n 2 k y 2 k A n 2 k ( r , m ) ( x ) = ( 1 ) m r k = 0 n n k ( 1 ) n k ( 2 x m + y ) k + ( 2 x m y ) k 2 A n k ( r , m ) ( x )
and
k = 0 n 1 2 n 2 k + 1 y 2 k + 1 A n 2 k 1 ( r , m ) ( x ) = ( 1 ) m r k = 0 n n k ( 1 ) n k ( 2 x m + y ) k ( 2 x m y ) k 2 A n k ( r , m ) ( x ) .
Proof. 
The proof follows the same arguments as the proof of Theorem 3. □
The special cases where y is replaced by i y and y = m are obvious. We continue skipping the presentation of these explicit results.

4. More Convolutions for B n ( r , m ) ( x ) and A n ( r , m ) ( x )

In this section, several other convolutions for B n ( r , m ) ( x ) and A n ( r , m ) ( x ) are derived. The first two theorems contain convolutions involving B n ( r , m ) ( x ) and powers of 2.
Theorem 4.
Let r and m be positive integers and x , y C . Then
k = 0 n n k 2 k B k ( r , m ) ( x ) y n k = 2 m ( r 1 ) k = 0 n n k B k ( r , m ) ( 2 x + y ) E n k ( m ) ( 0 ) ,
where E n ( m ) ( x ) is the generalized Euler polynomial of order m.
Proof. 
For | z | < 2 π we have from Equation (6)
B ( r , m ; x , 2 z ) e y z = n = 0 2 n B n ( r , m ) ( x ) z n n ! n = 0 y n z n n ! = n = 0 k = 0 n n k 2 k B k ( r , m ) ( x ) y n k z n n ! .
On the other hand, we observe that
B ( r , m ; x , 2 z ) e y z = 2 r m z r e 2 z 1 m e ( 2 x + y ) z = 2 m ( r 1 ) z r e z 1 m e ( 2 x + y ) z 2 e z + 1 m = 2 m ( r 1 ) n = 0 B n ( r , m ) ( 2 x + y ) z n n ! n = 0 E n ( m ) ( 0 ) z n n !
Comparing the coefficients for z n in the two power series proves the formula. □
Remark 1.
From the above proof, it is clear that we can also write
k = 0 n n k 2 k B k ( r , m ) ( x ) y n k = 2 m ( r 1 ) k = 0 n n k B k ( r , m ) ( x + y ) E n k ( m ) ( x ) = 2 m ( r 1 ) k = 0 n n k B k ( r , m ) ( x ) E n k ( m ) ( x + y ) = 2 m ( r 1 ) k = 0 n n k B k ( r , m ) ( y ) E n k ( m ) ( 2 x ) .
Corollary 8.
Let r and m be positive integers and x C . Then,
2 n B n ( r , m ) ( x ) = 2 m ( r 1 ) k = 0 n n k B k ( r , m ) ( 2 x ) E n k ( m ) ( 0 ) ,
k = 0 n n k ( 1 ) n k 2 k B k ( r , m ) ( x ) x n k = 2 m ( r 1 ) k = 0 n n k B k ( r , m ) ( x ) E n k ( m ) ( 0 )
and
k = 0 n n k ( 1 ) n k B k ( r , m ) ( x ) x n k = 2 m ( r 1 ) n k = 0 n n k B k ( r , m ) ( 0 ) E n k ( m ) ( 0 ) .
Proof. 
Set y = 0 , y = x and y = 2 x in (30), respectively. □
Evidently, the sums in the Corollary contain some interesting special cases. The evaluations with ( r , m ) = ( 1 , 1 ) , and x = 0 and x = 1 / 2 , respectively, yield the following identities for Bernoulli numbers:
k = 0 n 1 n k 2 2 n + 2 k n + 1 k B k B n + 1 k = ( 2 n 1 ) B n ,
k = 0 n 1 n k ( 1 ) k 2 2 n + 2 k n + 1 k B k B n + 1 k = ( 2 2 n ( 1 ) n ) B n ,
k = 0 n 1 n k ( 1 ) n k ( 1 2 k 1 ) B k = ( 2 n 1 ) B n
and
k = 0 n n k ( 1 ) n k 2 2 k ( 2 1 k 1 ) B k = 2 n + 1 k = 0 n n k ( 2 1 k 1 ) ( 2 n + 1 k 1 ) n + 1 k B k B n + 1 k ,
where we have employed the following relations
E n 1 ( x ) = 2 n ( B n ( x ) 2 n B n ( x / 2 ) ) ,
and
B n ( 1 / 2 ) = ( 2 1 n 1 ) B n .
It is difficult to say whether the Bernoulli identities Equations (34)–(37) are original. We could not find them in the book [20]. Hence, they are maybe not classical. However, they may have appeared elsewhere before. Furthermore, setting r = m = 1 , x = 0 and y = 1 in Equation (30) gives
k = 0 n n k 2 k B k = k = 0 n n k B k ( 1 ) E n k ( 0 ) = 2 n B n ( 1 / 2 ) ,
where we have used Equation (31). Hence, we rediscover Merca’s identity Equation (1).
Theorem 5.
Let r and m be positive integers and x , y C . Then,
k = 0 n n k 2 k B k ( r , m ) ( x ) y n k = 2 m ( r 1 ) j = 0 m m j ( 1 ) m j 2 j k = 0 n n k B k ( r , m ) ( 0 ) E n k ( m j ) ( 2 x m + y ) ,
where E n ( m ) ( x ) is the generalized Euler polynomial of order m.
Proof. 
The identity follows from
B ( r , m ; x , 2 z ) e y z = 2 r m z r e 2 z 1 m e ( 2 x + y ) z = 2 r m e ( 2 x m + y ) z z r e z 1 z r e 2 z 1 m = 2 r m e ( 2 x m + y ) z j = 0 m m j ( 1 ) m j z r e z 1 j z r e 2 z 1 m j = 2 r m z r e z 1 m e ( 2 x m + y ) z j = 0 m m j ( 1 ) m j 1 e z + 1 m j = 2 r m z r e z 1 m j = 0 m m j ( 1 ) m j 2 j m 2 e z + 1 m j e ( 2 x m + y ) z = 2 r ( m 1 ) j = 0 m m j ( 1 ) m j 2 j z r e z 1 m 2 e z + 1 m j e ( 2 x m + y ) z .
Corollary 9.
For x , y C the generalized Bernoulli polynomials of order ( r , m ) satisfy the following relation:
k = 0 n n k ( 1 ) k m ( r 1 ) 2 k B k ( r , m ) ( x ) y n k = 2 m ( r 1 ) j = 0 m m j ( 1 ) m j 2 j k = 0 n n k B k ( r , m ) ( 0 ) E n k ( m j ) ( m + y 2 x ) .
Proof. 
Replace x with m x and use the reciprocal relation (see [22])
B n ( r , m ) ( m x ) = ( 1 ) n m ( r 1 ) B n ( r , m ) ( x ) .
Setting r = m = 1 and using the fact that E n ( 0 ) ( x ) = x n , we obtain:
Corollary 10.
For x , y C the Bernoulli polynomials satisfy the following relation:
k = 0 n n k 2 k B k ( x ) y n k = k = 0 n n k B k ( 2 ( 2 x 1 + y ) n k E n k ( 2 x 1 + y ) )
and
k = 0 n n k ( 1 ) k 2 k B k ( x ) y n k = k = 0 n n k B k ( 2 ( 1 + y 2 x ) n k E n k ( 1 + y 2 x ) ) .
Inserting x = ( 1 y ) / 2 in Equation (42) or x = ( 1 + y ) / 2 in Equation (43) yields
k = 0 n n k 2 k B k 1 y 2 y n k = ( 2 2 n ) B n
and
k = 0 n n k ( 1 ) k 2 k B k 1 + y 2 y n k = ( 2 2 n ) B n
for each y C . We can see that y = 1 or y = 1 produce Merca’s identity Equation (1). It is also worth mentioning the special cases of the above results for y being a power of 2.
We now focus on presenting other types of convolutions. Some types follow straightforwardly from the definitions Equations (6) and (20). For instance, it is fairly easy to deduce that for each integer p 0 and x , y C
k = 0 n n k B k ( r , m ) ( x ) B n k ( r , p ) ( y ) = B n ( r , m + p ) ( x + y )
and
k = 0 n n k A k ( r , m ) ( x ) A n k ( r , p ) ( y ) = A n ( r , m + p ) ( x + y ) .
Setting p = x = 0 corresponds to the representation
k = 0 n n k B k ( r , m ) x n k = B n ( r , m ) ( x ) a n d k = 0 n n k A k ( r , m ) x n k = A n ( r , m ) ( x ) .
Theorem 6.
For m , p 0 and x , y C we have
k = 0 n n k B k ( r , m ) ( x ) B n k ( r , p ) ( y ) = k = 0 n n k ( 1 ) ( r 1 ) ( m + p ) + k B k ( r , m + p ) ( x ) ( 2 x + y ( m + p ) ) n k
and
k = 0 n n k A k ( r , m ) ( x ) A n k ( r , p ) ( y ) = k = 0 n n k ( 1 ) r ( m + p ) + k A k ( r , m + p ) ( x ) ( 2 x + y ( m + p ) ) n k .
Proof. 
The first identity follows from
z r e z 1 m e x z · z r e z 1 p e y z = ( 1 ) ( r 1 ) ( m + p ) ( z ) r e z 1 m + p e x z e ( 2 x + y ( m + p ) ) z .
The other one can be proved similarly. □
Theorem 7.
For r , m 1 and x , y C , we have the following convolution
m k = 0 n n k B k ( r , m ) ( x ) B n k ( r , 1 ) ( y + 1 m ) = ( x + y m ) ( n ) r B n r ( r , m ) ( x + y m ) ( n + 1 r ( m + 1 ) ) ( n ) r 1 B n + 1 r ( r , m ) ( x + y m ) ,
where ( n ) r denotes the falling factorial defined by
( n ) r = 1 , r = 0 ; n ( n 1 ) ( n r + 1 ) , r 1 .
Proof. 
Since
d d z e ( x + y m ) z ( e z 1 ) m = ( x + y m ) e ( x + y m ) z ( e z 1 ) m m e ( x + y + 1 m ) z ( e z 1 ) m + 1 ,
we have the relation
m z r e z 1 m + 1 e ( x + y + 1 m ) z = ( x + y m ) z r z r e z 1 m e ( x + y m ) z z r ( m + 1 ) d d z e ( x + y m ) z ( e z 1 ) m .
Replacing x with m x and using (41), we immediately get the alternating version of Theorem 7.
Corollary 11.
For r , m 1 and x , y C we have the following convolution
m k = 0 n n k ( 1 ) k m ( r 1 ) B k ( r , m ) ( x ) B n k ( r , 1 ) ( y + 1 m ) = ( y x ) ( n ) r B n r ( r , m ) ( y x ) ( n + 1 r ( m + 1 ) ) ( n ) r 1 B n + 1 r ( r , m ) ( y x ) .
For ( r , m ) = ( 1 , m ) the above results reduce to convolutions for generalized Bernoulli polynomials:
m k = 0 n n k B k ( m ) ( x ) B n k ( y + 1 m ) = ( x + y m ) n B n 1 ( m ) ( x + y m ) ( n m ) B n ( m ) ( x + y m ) ,
and
m k = 0 n n k ( 1 ) k B k ( m ) ( x ) B n k ( y + 1 m ) = ( y x ) n B n 1 ( m ) ( y x ) ( n m ) B n ( m ) ( y x ) .
These results generalize the known convolutions for Bernoulli polynomials, which are obtained for ( r , m ) = ( 1 , 1 ) . The analogue convolution for A n ( r , m ) ( x ) is given next.
Theorem 8.
For r , m 1 and x , y C , the following convolution result holds:
m k = 0 n n k A k ( r , m ) ( x ) A n k ( r , 1 ) ( y + 1 m ) = 2 ( x + y m ) ( n ) r A n r ( r , m ) ( x + y m ) 2 ( n + 1 r ( m + 1 ) ) ( n ) r 1 A n + 1 r ( r , m ) ( x + y m ) .
For r = 0 , the result becomes
m k = 0 n n k E k ( m ) ( x ) E n k ( y + 1 m ) = 2 ( x + y m ) E n ( m ) ( x + y m ) 2 E n + 1 ( m ) ( x + y m ) .
Proof. 
Use
m 2 z r e z + 1 m + 1 e ( x + y + 1 m ) z = 2 ( x + y m ) z r 2 z r e z + 1 m e ( x + y m ) z 2 m + 1 z r ( m + 1 ) d d z e ( x + y m ) z ( e z + 1 ) m .
Corollary 12.
For r , m 1 and x , y C the following convolution result holds:
m k = 0 n n k ( 1 ) k r m A k ( r , m ) ( x ) A n k ( r , 1 ) ( y + 1 m ) = 2 ( y x ) ( n ) r A n r ( r , m ) ( y x ) 2 ( n + 1 r ( m + 1 ) ) ( n ) r 1 A n + 1 r ( r , m ) ( y x ) .
For r = 0 the result becomes
m k = 0 n n k ( 1 ) k E k ( m ) ( x ) E n k ( y + 1 m ) = 2 ( y x ) E n ( m ) ( y x ) 2 E n + 1 ( m ) ( y x ) .
Proof. 
This result follows from the reciprocal relation for A n ( r , m ) ( x ) (see [22])
A n ( r , m ) ( m x ) = ( 1 ) n m r A n ( r , m ) ( x ) .
Theorem 9.
For r , m 1 and x , y C the following convolution result holds:
2 n B n ( r , m ) x + y 2 = 2 m ( r 1 ) ( n ) r m k = 0 n r m n r m k B k ( r , m ) ( x ) A n r m k ( r , m ) ( y ) .
Proof. 
Use
( 2 z ) r e 2 z 1 m e ( x + y ) z = 2 m ( r 1 ) z r m z r e z 1 m e x z 2 z r e z + 1 m e y z .
Upon replacing x with m x and y with m y and using Equations (41) and (57), we also get the alternating version of the previous result.
Corollary 13.
For r , m 1 and x , y C the following convolution result holds:
2 n B n ( r , m ) m x + y 2 = ( 1 ) n m ( r 1 ) 2 m ( r 1 ) ( n ) r m k = 0 n r m n r m k B k ( r , m ) ( x ) A n r m k ( r , m ) ( y ) .
The case ( r , m ) = ( 1 , 1 ) produces relations involving Bernoulli and Genocchi polynomials:
2 n B n x + y 2 = n k = 0 n 1 n 1 k B k ( x ) G n 1 k ( y )
and
2 n B n 1 x + y 2 = ( 1 ) n n k = 0 n 1 n 1 k B k ( x ) G n 1 k ( y ) .

5. Probabilistic Interpretation of Bernoulli and Euler–Genocchi Polynomials

In this section, we indicate a possible application of the polynomials discussed in this work.
Consider a probability measure space ( Ω , F , P ) , where Ω is a non-empty space, F is a σ -algebra of events and P is a probability measure on F . Let X be a random variable which is a measurable real function on the probability space ( Ω , F , P ) , and let p ( x ) be a probability density function of X. Then mathematical expectation and moments of order n of X are defined by
E [ X ] = x p ( x ) d x a n d E [ X n ] = x n p ( x ) d x ,
with E [ X 1 ] = E [ X ] . The moment-generating function of X is defined by
μ ( t ) = E [ e t X ] = e t x p ( x ) d x , t R .
The relation between the moment generating function μ ( t ) and and moments E [ X n ] is given by
μ ( t ) = k = 0 t k k ! E [ X k ] .
Suppose that (6), (20), (21) and (22) are moment-generating functions of random variables X 1 , X 2 , , X 5 . Then,
E [ X 1 n ] = d n d z n z r e z 1 m e x z | z = 0 = B n ( r , m ) ( x ) , E [ X 2 n ] = d n d z n z e z 1 m e x z | z = 0 = B n ( m ) ( x ) , E [ X 3 n ] = d n d z n 2 z r e z + 1 m e x z | z = 0 = A n ( r , m ) ( x ) , E [ X 4 n ] = d n d z n 2 e z + 1 m e x z | z = 0 = E n ( m ) ( x ) , E [ X 5 n ] = d n d z n 2 z e z + 1 m e x z | z = 0 = G n ( m ) ( x ) .
Another interpretation is the next example. Let X be a discrete random variable on ( Ω , F , P ) . Then, we can define a probability distribution for X by
P ( X = k ) = n k B k ( r , m ) B n ( r , m ) ( 1 ) , k = 0 , 1 , , n .
Note that
k = 0 n P ( X = k ) = k = 0 n n k B k ( r , m ) B n ( r , m ) ( 1 ) = B n ( r , m ) ( 1 ) B n ( r , m ) ( 1 ) = 1 .
Next, replacing k by n k in the identity
k = 0 n n k B k ( r , m ) x n k = B n ( r , m ) ( x )
we get
k = 0 n n k B n k ( r , m ) x k = B n ( r , m ) ( x ) ,
which, upon differentiating with respect to x, yields
k = 0 n k n k B n k ( r , m ) x k 1 = n B n 1 ( r , m ) ( x ) .
Inserting x = 1 gives
k = 0 n k n k B n k ( r , m ) = n B n 1 ( r , m ) ( 1 ) ,
which shows that
E [ X ] = n B n 1 ( r , m ) ( 1 ) .
In an analogous fashion, the quantity A n ( r , m ) ( x ) can be interpreted probabilistically.

6. Conclusions and Future Work

In this paper, mainly focusing on convolutions, we established additional properties of the generalized Bernoulli and Euler–Genocchi polynomials B n ( r , m ) ( x ) and A n ( r , m ) ( x ) , respectively. These properties provide generalizations of some known facts about generalized Bernoulli and Euler polynomials, respectively. In the future, we intend to work in two different directions. First, it seem desirable to find some new kinds of closed-form expressions for our polynomials (such as combinatorial, integral, hypergeometric and determinantal ones). Such expressions will provide us with new and significant properties of these polynomials. Second, it is possible to study the Apostol-type generalized polynomials B n ( r , m ) ( x ; λ ) associated with the complex parameter λ 0 , which are defined by the generating function
B ( r , m ; x , z ; λ ) = z r λ e z 1 m e x z = n = 0 B n ( r , m ) ( x ; λ ) z n n ! .
Then, it is fairly easy to identify the relations B ( r , m ; x , z ) = B ( r , m ; x , z ; 1 ) and A ( r , m ; x , z ) = ( 2 ) m B ( r , m ; x , z ; 1 ) . Therefore, by means of B n ( r , m ) ( x ; λ ) it is possible to discuss both B n ( r , m ) ( x ) and A n ( r , m ) ( x ) , at once.

Author Contributions

Conceptualization, R.F. and Ž.T.; Methodology, R.F. and Ž.T.; validation, R.F. and Ž.T.; formal analysis, R.F. and Ž.T.; writing—original draft preparation, R.F.; writing—review and editing, R.F. and Ž.T.; supervision, R.F. and Ž.T. All authors have read and agreed to the published version of the manuscript.

Funding

Ž.T. was supported of the Department of Mathematics, Faculty of Sciences, University of Ostrava.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; Dover Publications: New York, NY, USA, 1972. [Google Scholar]
  2. Apostol, T.M. Introduction to Analytic Number Theory; Springer: Heidelberg/Berlin, Germany, 1976. [Google Scholar]
  3. Agoh, T. Convolution identities for Bernoulli and Genocchi polynomials. Electron. J. Combin. 2014, 21, 65. [Google Scholar] [CrossRef]
  4. Agoh, T.; Dilcher, K. Convolution identities and lacunary recurrences for Bernoulli numbers. J. Number Theory 2007, 124, 105–122. [Google Scholar] [CrossRef] [Green Version]
  5. Agoh, T.; Dilcher, K. Higher-order recurrences for Bernoulli numbers. J. Number Theory 2009, 129, 1837–1847. [Google Scholar] [CrossRef] [Green Version]
  6. Agoh, T.; Dilcher, K. Shortened recurrence relations for Bernoulli numbers. Discret. Math. 2009, 309, 887–898. [Google Scholar] [CrossRef] [Green Version]
  7. Chu, W.; Wang, C.Y. Convolution formulae for Bernoulli numbers. Integral Transforms Spec. Funct. 2010, 21, 437–457. [Google Scholar] [CrossRef]
  8. Chu, W.; Zhou, R. Convolutions of Bernoulli and Euler polynomials. Sarajevo J. Math. 2010, 6, 147–163. [Google Scholar]
  9. Dilcher, K.; Vignat, C. General convolution identities for Bernoulli and Euler polynomials. J. Math. Anal. Appl. 2016, 435, 1478–1498. [Google Scholar] [CrossRef]
  10. Srivastava, H.M.; Vignat, C. Probabilistic proofs of some relationships between the Bernoulli and Euler polynomials. Eur. J. Pure Appl. Math. 2012, 5, 97–107. [Google Scholar]
  11. Dere, R.; Simsek, Y.; Srivastava, H.M. A unified presentation of three families of generalized Apostol-type polynomials based upon the theory of the umbral calculus and the umbral algebra. J. Number Theory 2013, 133, 3245–3263. [Google Scholar] [CrossRef]
  12. Ozden, H.; Simsek, Y.; Srivastava, H.M. A unified presentation of the generating functions of the generalized Bernoulli, Euler and Genocchi polynomials. Comput. Math. Appl. 2010, 60, 2779–2787. [Google Scholar] [CrossRef] [Green Version]
  13. Srivastava, H.M.; He, Y.; Araci, S.; Abdel-Aty, M. Higher-order convolutions for Apostol-Bernoulli, Apostol-Euler and Apostol-Genocchi polynomials. Mathematics 2018, 6, 329. [Google Scholar]
  14. Srivastava, H.M.; Masjed-Jamei, M.; Beyki, M.R. Some new generalizations and applications of the Apostol-Bernoulli, Apostol-Euler and Apostol-Genocchi polynomials. Rocky Mountain J. Math. 2019, 49, 681–690. [Google Scholar] [CrossRef]
  15. Srivastava, H.M. The Zeta and related functions: Recent developments. J. Adv. Engrg. Comput. 2019, 3, 329–354. [Google Scholar] [CrossRef] [Green Version]
  16. Merca, M. On the Song recurrence relation for the Riemann zeta function. Miskolc Math. Notes 2017, 17, 941–945. [Google Scholar] [CrossRef] [Green Version]
  17. Merca, M. On families of linear recurrence relations for the special values of the Riemann zeta function. J. Number Theory 2017, 170, 55–65. [Google Scholar] [CrossRef]
  18. Saalschütz, L. Vorlesungen Über die Bernoullischen Zahlen, Ihren Zusammenhang mit den Secanten-Coefficienten und Ihre Wichtigeren Anwendungen; Springer: Heidelberg/Berlin, Germany, 1893. [Google Scholar]
  19. Nielsen, N. Traité Úlúmentaire des Nombres de Bernoulli; Gauthier-Villars: Paris, France, 1923. [Google Scholar]
  20. Hansen, E.R. Bernoulli, Euler, and Stirling Polynomials and Numbers and Neumann Polynomials. In A Table of Series and Products; Prentice Hall: Englewood Cliffs, NJ, USA, 1975; pp. 331–353. [Google Scholar]
  21. Deeba, E.Y.; Rodriguez, D.M. Stirling’s series and Bernoulli numbers. Amer. Math. Mon. 1991, 98, 423–426. [Google Scholar] [CrossRef]
  22. Frontczak, R.; Tomovski, Z. Generalized Euler-Genocchi polynomials and Lucas numbers. Integers 2020, 20, A52. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Frontczak, R.; Tomovski, Ž. Convolutions for Bernoulli and Euler–Genocchi Polynomials of Order (r,m) and Their Probabilistic Interpretation. Symmetry 2022, 14, 1220. https://doi.org/10.3390/sym14061220

AMA Style

Frontczak R, Tomovski Ž. Convolutions for Bernoulli and Euler–Genocchi Polynomials of Order (r,m) and Their Probabilistic Interpretation. Symmetry. 2022; 14(6):1220. https://doi.org/10.3390/sym14061220

Chicago/Turabian Style

Frontczak, Robert, and Živorad Tomovski. 2022. "Convolutions for Bernoulli and Euler–Genocchi Polynomials of Order (r,m) and Their Probabilistic Interpretation" Symmetry 14, no. 6: 1220. https://doi.org/10.3390/sym14061220

APA Style

Frontczak, R., & Tomovski, Ž. (2022). Convolutions for Bernoulli and Euler–Genocchi Polynomials of Order (r,m) and Their Probabilistic Interpretation. Symmetry, 14(6), 1220. https://doi.org/10.3390/sym14061220

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