SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research

New identities involving Bernoulli, Euler and Genocchi numbers

Su Hu1, Daeyeoul Kim2 and Min-Soo Kim3*

Author Affiliations

1 Department of Mathematical Sciences, Korea Advanced Institute of Science and Technology (KAIST), 373-1 Guseong-dong, Yuseong-gu, Daejeon, 305-701, Republic of Korea

2 National Institute for Mathematical Sciences, Yuseong-daero 1689-gil, Yuseong-gu, Daejeon, 305-811, Republic of Korea

3 Division of Cultural Education, Kyungnam University, 7 (Woryeong-dong) kyungnamdaehak-ro, Masanhappo-gu, Changwon-si, Gyeongsangnam-do, 631-701, Republic of Korea

For all author emails, please log on.

Advances in Difference Equations 2013, 2013:74  doi:10.1186/1687-1847-2013-74

The electronic version of this article is the complete one and can be found online at: http://www.advancesindifferenceequations.com/content/2013/1/74


Received:26 October 2012
Accepted:5 March 2013
Published:26 March 2013

© 2013 Hu et al.; licensee Springer

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Using p-adic integral, many new convolution identities involving Bernoulli, Euler and Genocchi numbers are given.

MSC: 11B68, 11S80.

Keywords:
Bernoulli numbers; Euler numbers; sums of products; fermionic p-adic integral; Volkenborn integral

1 Introduction

The study of the identities involving Bernoulli numbers and polynomials, Euler numbers and polynomials has a long history. More than 250 years ago, Euler discovered the following sums of products identity involving Bernoulli numbers:

l = 1 n 1 ( 2 n 2 l ) B 2 l B 2 n 2 l = ( 2 n + 1 ) B 2 n ( n > 1 ) (1.1)

(see [1]).

This identity has been extended by many authors in different directions (see [2-21]). Recently, Kim and Hu [22] obtained an analogy of Euler’s sums of products identity for Apostol-Bernoulli numbers.

In 1978, Miki [23] obtained the following well-known identity, which is now known as Miki’s identity:

k = 2 n 2 B k B n k k ( n k ) k = 2 n 2 ( n k ) B k B n k k ( n k ) = 2 H n B n n

for every n = 4 , 5 ,  , where H n = 1 + 1 / 2 + + 1 / n .

In 1997, Matiyasevich [24] found another identity of this type

( n + 2 ) k = 2 n 2 B k B n k 2 k = 2 n 2 ( n + 2 k ) B k B n k = n ( n + 1 ) B n

for any n = 4 , 5 ,  .

In 2004, Dunne and Schubert [25] obtained the convolution identities for sums of products of Bernoulli numbers motivated by the role of these identities in quantum field theory and string theory. In 2006, Crabb [26] showed that Gessel’s generalization of Miki’s identity [27] is a direct consequence of a functional equation for the generating function. During the same year, Sun and Pan [28] established two general identities involving Bernoulli and Euler polynomials, which imply both Miki and Matiyasevich’s identities.

According to statement by Cohen in the first paragraph of [[29], Chapter 11], the p-adic functions with nice properties are powerful tools for studying many results of classical number theory in a straightforward manner, for instance, strengthening of almost all the arithmetic results on Bernoulli numbers.

For example, in Exercises 2 and 3 of [[29], Chapter 11], Cohen demonstrated a method to prove the following reciprocity formula for Bernoulli numbers using p-adic integral:

m ! j = 0 m B m j ( m j ) ! B n + j + 1 ( j + 1 ) ! + n ! j = 0 n B n j ( n j ) ! B m + j + 1 ( j + 1 ) ! = B m + n

(see Exercise 3(c) of [[29], Chapter 11]).

Recently, using p-adic integral, Kim et al.[8] proved several identities on Bernoulli and Euler numbers. More comprehensive coverage can be found in the monographs by Chio et al.[2], Kim et al.[4], Kim et al.[5], Kim et al.[6], Kim et al.[7], Kim and Kim [9], Kim et al.[13], Lee and Kim [14].

In this paper, following the methods of [8], we shall further provide many new convolution identities involving Bernoulli, Euler and Genocchi numbers.

Let B n ( x ) , E n ( x ) and G n ( x ) be the nth Bernoulli, Euler and Genocchi polynomials, respectively. In what follows, we use E n to denote the special value of E n ( x ) at 0, that is, E n = E n ( 0 ) .

Let δ i , j be the Kronecker symbol defined by δ i , i = 1 and δ i , j = 0 for i j .

In fact, the following identities are shown in this paper:

l = 0 2 n + 1 ( 2 n + 1 l ) B 2 n l + 1 E l = δ 0 , n (Theorem 3.3) , (1.2)

k = 0 2 m + 1 l = 0 2 n + 2 ( 2 m + 1 k ) ( 2 n + 2 l ) B 2 m k + 1 E 2 n l + 2 E k + l = 0 (Theorem 3.4) , (1.3)

k = 0 2 m + 2 l = 0 2 n + 1 ( 2 m + 2 k ) ( 2 n + 1 l ) B 2 m k + 2 E 2 n l + 1 E k + l = 2 B 2 m + 2 E 2 n + 1 (Theorem 3.4) , (1.4)

l = 0 2 n + 2 ( 2 n + 2 l ) E 2 n l + 2 ( E l 2 E l + 1 ) = 0 (Corollary 3.5) , (1.5)

k = 0 2 m + 1 ( 2 m + 1 k ) B 2 m k + 1 ( E k + 2 E k + 1 ) = 0 (Corollary 3.5) , (1.6)

k = 0 2 m ( 2 m k ) B 2 m k ( E k 2 E k + 1 ) = 2 B 2 m (Corollary 3.5) , (1.7)

l = 0 2 n + 1 ( 2 n + 1 l ) E 2 n l + 1 ( E l 6 E l + 1 + 6 E l + 2 ) = 2 E 2 n + 1 (Corollary 3.5) , (1.8)

l = 0 2 n + 1 ( 2 n + 1 l ) E 2 n l + 1 E l = 2 E 2 n + 1 if  n N  (Remark 3.6) , (1.9)

(1.10)

(1.11)

(1.12)

(1.13)

This paper is organized as follows. In the next section, we recall the fundamental results between p-adic integral and Bernoulli and Euler numbers. Then using these results we prove new identities (identities (1.2), (1.3), (1.4), (1.5), (1.7), (1.8), (1.9)) involving Bernoulli and Euler numbers in Section 3, prove new identities (identities (1.10), (1.11), (1.12)) involving higher-order Bernoulli and Euler numbers and polynomials in Section 4, and also prove a new identity (identity (1.13)) involving both Bernoulli, Euler and Genocchi numbers in the final section.

2 p-adic analysis

In this section, we recall the fundamental results between p-adic integral and Bernoulli and Euler numbers, and we see that the properties of p-adic integrals may imply most well-known facts on Bernoulli numbers and polynomials, Euler numbers and polynomials.

We assume p is an odd prime number. The symbols Z p , Q p and C p denote the rings of p-adic integers, the field of p-adic numbers and the field of p-adic completion of the algebraic closure of Q p , respectively. ℕ denotes the set of natural numbers and N 0 denotes N { 0 } .

Let UD ( Z p ) be the space of uniformly (or strictly) differentiable function on Z p . Then the Volkenborn integral of f is defined by

I ( f ( z ) ) = Z p f ( z ) d μ ( z ) = lim N 1 p N a = 0 p N 1 f ( a ) (2.1)

and this limit always exists when f UD ( Z p ) (see [[30], p.264]). For such functions we have

I ( f 1 ) I ( f ) = f ( 0 ) , (2.2)

where f ( 0 ) = ( d f ( z ) / d z ) | z = 0 and f 1 ( z ) = f ( z + 1 ) . By (2.2), we have the following Volkenborn p-adic integral representation of the generating function of Bernoulli polynomials B n ( x ) :

I ( e ( x + z ) t ) = ( t e t 1 ) e x t = n = 0 B n ( x ) t n n ! , (2.3)

which is equivalent to

I ( ( x + z ) n ) = B n ( x ) (2.4)

(see [12]).

Setting x = 0 in (2.4), we obtain the Volkenborn p-adic integral representation of the nth Bernoulli numbers

B n = B n ( 0 ) = I ( z n ) (2.5)

(see [12,30]).

Thus from the binomial theorem and (2.4), we get

B n ( x ) = I ( ( x + z ) n ) = l = 0 n ( n l ) x l I ( z n l ) = l = 0 n ( n l ) x l B n l , (2.6)

where

( n l ) = n ( n 1 ) ( n l + 1 ) l ! . (2.7)

The recurrence formula of B n is given by

B 0 = 1 , ( B + 1 ) n B n = { 1 , n = 1 , 0 , n > 1 , (2.8)

with the usual convention of denoting B n by B n .

By (2.6) and (2.8), we have

B 0 = 1 , B n ( 1 ) B n = δ 1 , n if  n N . (2.9)

The fermionic p-adic integral I 1 ( f ) on Z p is defined by

I 1 ( f ( z ) ) = Z p f ( z ) d μ 1 ( z ) = lim N a = 0 p N 1 f ( a ) ( 1 ) a , (2.10)

where f UD ( Z p ) . This integral was introduced by Kim [12] in order to derive useful formulas involving the Euler numbers and polynomials. It has also been defined independently by Shiratani and Yamamoto [31] in order to interpolate the Euler numbers p-adically. Osipov [32] gave a new proof of the existence of the Kubota-Leopoldt p-adic zeta function by using the integral representation

I ε ( f ( z ) ) = Z p f ( z ) d μ ε ( z ) = lim N a = 0 p l N 1 f ( a ) ε a , (2.11)

where f UD ( Z p ) , ε k = 1 , ε 1 , ( k , p ) = 1 , and k ( p l 1 ) . Note that when k = 2 , I ε ( f ) is the fermionic p-adic integral I 1 ( f ) on Z p . Recently, the fermionic p-adic integral I 1 ( f ) on Z p has been used by the third author to give a brief proof of Stein’s classical result on Euler numbers modulo power of two [10]; it has also been used by Maïga [33] to give some new identities and congruences concerning Euler numbers and polynomials.

From the definition of I 1 ( f ) , we have

I 1 ( f 1 ) + I 1 ( f ) = 2 f ( 0 ) , (2.12)

where f 1 ( z ) = f ( z + 1 ) (see [[12], Lemma 1]). Thus from induction, we have the following fermionic p-adic integral equation:

I 1 ( f α ) + ( 1 ) α 1 I 1 ( f ) = 2 i = 0 α 1 ( 1 ) α i 1 f ( i ) , (2.13)

where α N and

f α ( z ) = f ( z + α ) (2.14)

(see [[12], Theorem 2]).

Using formula (2.12), we have the following fermionic p-adic integral representation of the generating function of Euler polynomials E n ( x ) :

I 1 ( e ( x + z ) t ) = 2 e t + 1 e x t = n = 0 E n ( x ) t n n ! , (2.15)

namely

I 1 ( ( x + z ) n ) = E n ( x ) (2.16)

(see [[34], Proposition 2.1]).

Setting x = 0 in (2.16), we obtain the fermionic p-adic integral representation of the nth Euler numbers

E n = E n ( 0 ) = I 1 ( z n ) (2.17)

(comparing with (2.5)).

Therefore, by the binomial theorem, (2.16), and (2.17), we have

E n ( x ) = l = 0 n ( n l ) x l I 1 ( z n l ) = l = 0 n ( n l ) x l E n l . (2.18)

The recurrence formula of E n is given by

E 0 = 1 , ( E + 1 ) n + E n = 0 if  n N , (2.19)

with the usual convention of denoting E n by E n .

By (2.18) and (2.19), we have

E n ( 1 ) + E n = 2 δ 0 , n if  n N 0 . (2.20)

The p-adic integral also implies the following reflection symmetric relations for the Bernoulli and the Euler polynomials (see formulas (2.12) and (2.13) in [34]):

B n ( x ) = ( 1 ) n B n ( 1 x ) , where  n N 0 , (2.21)

E n ( x ) = ( 1 ) n E n ( 1 x ) , where  n N 0 . (2.22)

3 Identities involving Bernoulli and Euler numbers

In this section, we prove new identities involving Bernoulli and Euler numbers.

Lemma 3.1For l N 0 , we have

I 1 ( ( 1 x ) l ) = 2 + E l 2 δ 0 , l .

Proof Putting f ( x ) = ( 1 x ) l in (2.12), by (2.17), we have

I 1 ( ( 1 x ) l ) = 2 I 1 ( ( x ) l ) = 2 + ( 1 ) l + 1 E l . (3.1)

By (2.22), with x = 0 , and (2.20), we obtain

( 1 ) l E l = ( 1 ) l E l ( 0 ) = E l ( 1 ) = 2 δ 0 , l E l ,

substituting it to (3.1), we have the desired result. □

The following well-known fact on Euler polynomials may also be established using p-adic integral. We refer to Corollary 1.1 in [35] for another proof.

Theorem 3.2 E 2 n = 0 for n 1 .

Proof We prove our result following the method of the proof for Theorem 2.1 in [8].

By (2.17) and (2.18), we have

I 1 ( E n ( x ) ) = l = 0 n ( n l ) E n l I 1 ( x l ) = l = 0 n ( n l ) E n l E l . (3.2)

On the other hand, using (2.18), (2.22) and Lemma 3.1, we obtain

I 1 ( E n ( x ) ) = ( 1 ) n I 1 ( E n ( 1 x ) ) = ( 1 ) n l = 0 n ( n l ) E n l I 1 ( ( 1 x ) l ) = ( 1 ) n l = 0 n ( n l ) E n l ( 2 + E l 2 δ 0 , l ) = ( 1 ) n 2 ( l = 0 n ( n l ) E n l E n ) + ( 1 ) n l = 0 n ( n l ) E n l E l , (3.3)

where δ 0 , l is the Kronecker symbol.

Comparing (3.2) and (3.3), then using (2.18), with x = 1 , and (2.20), we have

( 1 + ( 1 ) n + 1 ) l = 0 n ( n l ) E n l E l = ( 1 ) n 2 ( l = 0 n ( n l ) E n l E n ) = ( 1 ) n 2 ( E n ( 1 ) E n ) = ( 1 ) n 4 ( δ 0 , n E n ) if  n N 0 . (3.4)

Putting n = 0 on the left-hand side of (3.4), we have 0 = 1 E 0 , hence E 0 = 1 . Replacing n by 2n ( n 1 ), the right-hand side becomes E 2 n = 0 , which completes the proof. □

Theorem 3.3For n N 0 , we have

l = 0 2 n + 1 ( 2 n + 1 l ) B 2 n l + 1 E l = δ 0 , n .

Proof The proof goes the same way as in the above theorem.

By (2.5) and (2.6), we have

I 1 ( B n ( x ) ) = l = 0 n ( n l ) B n l I 1 ( x l ) = l = 0 n ( n l ) B n l E l . (3.5)

On the other hand, using (2.6), (2.21) and Lemma 3.1, we obtain

I 1 ( B n ( x ) ) = ( 1 ) n I 1 ( B n ( 1 x ) ) = ( 1 ) n l = 0 n ( n l ) B n l I 1 ( ( 1 x ) l ) = ( 1 ) n l = 0 n ( n l ) B n l ( 2 + E l 2 δ 0 , l ) = ( 1 ) n 2 ( l = 0 n ( n l ) B n l B n ) + ( 1 ) n l = 0 n ( n l ) B n l E l . (3.6)

Comparing (3.5) and (3.6), then using (2.6), with x = 1 , and (2.9), we have

( 1 + ( 1 ) n + 1 ) l = 0 n ( n l ) B n l E l = ( 1 ) n 2 ( l = 0 n ( n l ) B n l B n ) = ( 1 ) n 2 ( B n ( 1 ) B n ) = ( 1 ) n 2 δ 1 , n .

Letting n = 2 m + 1 ( m 0 ) in this identity, we get our result. □

Theorem 3.4For m , n N 0 , we have

k = 0 2 m + 1 l = 0 2 n + 2 ( 2 m + 1 k ) ( 2 n + 2 l ) B 2 m k + 1 E 2 n l + 2 E k + l = 0 , (3.7)

k = 0 2 m + 2 l = 0 2 n + 1 ( 2 m + 2 k ) ( 2 n + 1 l ) B 2 m k + 2 E 2 n l + 1 E k + l = 2 B 2 m + 2 E 2 n + 1 . (3.8)

Taking m = 0 , n = 0 in (3.7) and (3.8), respectively, we obtain the following corollary.

Corollary 3.5For m , n N 0 , we have

l = 0 2 n + 2 ( 2 n + 2 l ) E 2 n l + 2 ( E l 2 E l + 1 ) = 0 , (3.9)

k = 0 2 m + 1 ( 2 m + 1 k ) B 2 m k + 1 ( E k + 2 E k + 1 ) = 0 , (3.10)

k = 0 2 m ( 2 m k ) B 2 m k ( E k 2 E k + 1 ) = 2 B 2 m , (3.11)

l = 0 2 n + 1 ( 2 n + 1 l ) E 2 n l + 1 ( E l 6 E l + 1 + 6 E l + 2 ) = 2 E 2 n + 1 . (3.12)

Proof of Theorem 3.4 We prove our result following the method of the proof for Theorem 2.3 in [8].

Putting x = 1 in (2.6), (2.18), and using (2.9), (2.20), we get

B n ( 1 ) = l = 0 n ( n l ) B n l = B n + δ 1 , n , E n ( 1 ) = l = 0 n ( n l ) E n l = 2 δ 0 , n E n . (3.13)

By the linear property of p-adic integral and (2.17), we have the following fermionic p-adic integral representation of the product of Bernoulli and Euler polynomials:

I 1 ( B m ( x ) E n ( x ) ) = k = 0 m l = 0 n ( m k ) ( n l ) B m k E n l I 1 ( x k + l ) = k = 0 m l = 0 n ( m k ) ( n l ) B m k E n l E k + l . (3.14)

Further, by (2.21), (2.22), Lemma 3.1 and (3.13), we have

I 1 ( B m ( x ) E n ( x ) ) = I 1 ( ( 1 ) m + n B m ( 1 x ) E n ( 1 x ) ) = ( 1 ) m + n k = 0 m l = 0 n ( m k ) ( n l ) B m k E n l I 1 ( ( 1 x ) k + l ) = ( 1 ) m + n k = 0 m l = 0 n ( m k ) ( n l ) B m k E n l ( 2 + E k + l 2 δ 0 , k + l ) = ( 1 ) m + n 2 k = 0 m l = 0 n ( m k ) ( n l ) B m k E n l + ( 1 ) m + n + 1 2 B m E n + ( 1 ) m + n k = 0 m l = 0 n ( m k ) ( n l ) B m k E n l E k + l = ( 1 ) m + n 2 ( k = 0 m ( m k ) B m k ) ( l = 0 n ( n l ) E n l ) + ( 1 ) m + n + 1 2 B m E n + ( 1 ) m + n k = 0 m l = 0 n ( m k ) ( n l ) B m k E n l E k + l = ( 1 ) m + n 2 B m ( 1 ) E n ( 1 ) + ( 1 ) m + n + 1 2 B m E n + ( 1 ) m + n k = 0 m l = 0 n ( m k ) ( n l ) B m k E n l E k + l . (3.15)

Comparing (3.14) and (3.15), then using (3.13), we have

( 1 + ( 1 ) m + n + 1 ) k = 0 m l = 0 n ( m k ) ( n l ) B m k E n l E k + l = ( 1 ) m + n 2 B m ( 1 ) E n ( 1 ) + ( 1 ) m + n + 1 2 B m E n = ( 1 ) m + n 2 ( B m + δ 1 , m ) ( 2 δ 0 , n E n ) + ( 1 ) m + n + 1 2 B m E n = ( 1 ) m + n 2 ( 2 B m δ 0 , n 2 B m E n + 2 δ 1 , m δ 0 , n E n δ 1 , m ) . (3.16)

Replacing m by 2 m + 1 and n by 2 n + 1 in (3.16), respectively, we have

( 1 + ( 1 ) n ) k = 0 2 m + 1 l = 0 n ( 2 m + 1 k ) ( n l ) B 2 m k + 1 E n l E k + l = ( 1 ) n + 1 2 ( 2 B 2 m + 1 δ 0 , n 2 B 2 m + 1 E n + 2 δ 1 , 2 m + 1 δ 0 , n E n δ 1 , 2 m + 1 ) (3.17)

and

( 1 + ( 1 ) m ) k = 0 m l = 0 2 n + 1 ( m k ) ( 2 n + 1 l ) B m k E 2 n l + 1 E k + l = ( 1 ) m 2 ( 2 B m E 2 n + 1 + E 2 n + 1 δ 1 , m ) . (3.18)

Finally, replacing n by 2 n + 2 in (3.17) and m by 2m in (3.18), using Theorem 3.2, we get our result. □

Remark 3.6 Putting n = 0 in (3.17), we recover Theorem 3.3. Setting m = 0 in (3.18), we obtain

l = 0 2 n + 1 ( 2 n + 1 l ) E 2 n l + 1 E l = E 2 n + 1 E 0 + E 2 n E 1 + + E 1 E 2 n ( 2 n 1 ) - terms + E 0 E 2 n + 1 = 2 E 2 n + 1 if  n N .

This identity is trivial, because by Theorem 3.2, the ( 2 n 1 ) th terms are identical to zero.

4 Identities involving higher-order Bernoulli and Euler numbers

In this section, we prove new identities involving higher-order Bernoulli and Euler numbers and polynomials.

The higher-order Bernoulli polynomials B n ( α ) ( x ) and the higher-order Euler polynomials E n ( α ) ( x ) are defined by the following generating functions:

( t e t 1 ) α e x t = n = 0 B n ( α ) ( x ) t n n ! (4.1)

and

( 2 e t + 1 ) α e x t = n = 0 E n ( α ) ( x ) t n n ! , (4.2)

respectively. For n N 0 , we have B n ( α ) = B n ( α ) ( 0 ) and E n ( α ) = E n ( α ) ( 0 ) , where B n ( α ) and E n ( α ) denote the higher-order Bernoulli and Euler numbers, respectively (see [[20], Section 1.6] and [[21], (16), (17)]). Moreover, from (4.1) and (4.2), we have

B n ( α ) ( x ) = l = 0 n ( n l ) B n l ( α ) x l , (4.3)

B n ( α ) ( x ) = ( 1 ) n B n ( α ) ( α x ) , (4.4)

E n ( α ) ( x ) = l = 0 n ( n l ) E n l ( α ) x l , (4.5)

E n ( α ) ( x ) = ( 1 ) n E n ( α ) ( α x ) . (4.6)

Lemma 4.1For α N and l N 0 , we have

I 1 ( ( α x ) l ) = ( 1 ) α ( 2 δ 0 , l E l ) + 2 i = 0 α 1 ( 1 ) i ( α i ) l .

Proof Letting f ( x ) = ( α x ) l ( α N ) in (2.14), we have f α ( x ) = ( x ) l .

By (2.13), we have

I 1 ( ( α x ) l ) = ( 1 ) α I 1 ( ( x ) l ) + 2 i = 0 α 1 ( 1 ) i f ( i ) = ( 1 ) α ( 1 ) l E l + 2 i = 0 α 1 ( 1 ) i ( α i ) l . (4.7)

By (2.20), (2.22) with x = 0 , we have ( 1 ) l E l = E l ( 1 ) = 2 δ 0 , l E l , and we obtain the assertion of the lemma. □

Theorem 4.2For α N and n N 0 , we have

Proof From (2.17), (4.3) and the linear property of p-adic integral, we obtain the following fermionic p-adic integral representation of the higher-order Bernoulli polynomials:

I 1 ( B n ( α ) ( x ) ) = l = 0 n ( n l ) B n l ( α ) I 1 ( x l ) = l = 0 n ( n l ) B n l ( α ) E l . (4.8)

On the other hand, by (4.4) and Lemma 4.1, we have

I 1 ( B n ( α ) ( x ) ) = ( 1 ) n l = 0 n ( n l ) B n l ( α ) I 1 ( ( α x ) l ) = ( 1 ) n l = 0 n ( n l ) B n l ( α ) ( ( 1 ) α ( 2 δ 0 , l E l ) + 2 i = 0 α 1 ( 1 ) i ( α i ) l ) = ( 1 ) n + α 2 B n ( α ) + ( 1 ) n 2 l = 0 n ( n l ) B n l ( α ) i = 0 α 1 ( 1 ) i ( α i ) l ( 1 ) n + α l = 0 n ( n l ) B n l ( α ) E l . (4.9)

Clearly,

i = 0 α 1 ( 1 ) i ( α i ) l = i = 1 α ( 1 ) α i i l . (4.10)

Comparing (4.8) and (4.9), then using (4.10) and (4.3) with x = i , we obtain

( 1 + ( 1 ) n + α ) l = 0 n ( n l ) B n l ( α ) E l = ( 1 ) n + α 2 B n ( α ) + ( 1 ) n 2 l = 0 n ( n l ) B n l ( α ) i = 0 α 1 ( 1 ) i ( α i ) l = ( 1 ) n + α 2 B n ( α ) + ( 1 ) n 2 l = 0 n ( n l ) B n l ( α ) i = 1 α ( 1 ) α i i l = ( 1 ) n + α 2 B n ( α ) + ( 1 ) n 2 i = 1 α ( 1 ) α i ( l = 0 n ( n l ) B n l ( α ) i l ) = ( 1 ) n + α 2 B n ( α ) + ( 1 ) n 2 i = 1 α ( 1 ) α i B n ( α ) ( i ) = ( 1 ) n + α 2 ( B n ( α ) + i = 1 α ( 1 ) i B n ( α ) ( i ) ) . (4.11)

Since B n ( α ) = B n ( α ) ( 0 ) , therefore, by (4.11), we have

( 1 ) n + α 2 i = 0 α ( 1 ) i B n ( α ) ( i ) = ( 1 + ( 1 ) n + α ) l = 0 n ( n l ) B n l ( α ) E l ,

which is the required result. □

Theorem 4.3For α N and n N 0 , we have

Proof This follows from the same process as in the proof of Theorem 4.2 by using (2.17), (4.5), (4.6) and Lemma 4.1. □

Theorem 4.4For m , n N 0 and α N , we have

Proof First we have the following fermionic p-adic integral representation of the product of higher-order Bernoulli and Euler polynomials:

I 1 ( B m ( α ) ( x ) E n ( α ) ( x ) ) = k = 0 m l = 0 n ( m k ) ( n l ) B m k ( α ) E n l ( α ) I 1 ( x k + l ) = k = 0 m l = 0 n ( m k ) ( n l ) B m k ( α ) E n l ( α ) E k + l . (4.12)

On the other hand, by (4.4), (4.6), Lemma 4.1, by (4.3) and (4.5) with x = i , we have

I 1 ( B m ( α ) ( x ) E n ( α ) ( x ) ) = I 1 ( ( 1 ) m + n B m ( α ) ( 1 x ) E n ( α ) ( 1 x ) ) = ( 1 ) m + n k = 0 m l = 0 n ( m k ) ( n l ) B m k ( α ) E n l ( α ) I 1 ( ( α x ) k + l ) = ( 1 ) m + n + α k = 0 m l = 0 n ( m k ) ( n l ) B m k ( α ) E n l ( α ) ( 2 δ 0 , k + l E k + l + 2 i = 1 α ( 1 ) i i k + l ) = ( 1 ) m + n + α 2 ( B m ( α ) E n ( α ) + i = 1 α ( 1 ) i k = 0 m l = 0 n ( m k ) ( n l ) B m k ( α ) i k E n l ( α ) i l ) ( 1 ) m + n + α k = 0 m l = 0 n ( m k ) ( n l ) B m k ( α ) E n l ( α ) E k + l = ( 1 ) m + n + α 2 ( B m ( α ) E n ( α ) + i = 1 α ( 1 ) i ( k = 0 m ( m k ) B m k ( α ) i k ) ( l = 0 n ( n l ) E n l ( α ) i l ) ) ( 1 ) m + n + α k = 0 m l = 0 n ( m k ) ( n l ) B m k ( α ) E n l ( α ) E k + l = ( 1 ) m + n + α 2 ( B m ( α ) E n ( α ) + i = 1 α ( 1 ) i B m ( α ) ( i ) E n ( α ) ( i ) ) ( 1 ) m + n + α k = 0 m l = 0 n ( m k ) ( n l ) B m k ( α ) E n l ( α ) E k + l . (4.13)

Comparing (4.12) and (4.13), and also noticing that B n ( α ) = B n ( α ) ( 0 ) and E n ( α ) = E n ( α ) ( 0 ) , we have

( 1 + ( 1 ) m + n + α ) k = 0 m l = 0 n ( m k ) ( n l ) B m k ( α ) E n l ( α ) E k + l = ( 1 ) m + n + α 2 i = 0 α ( 1 ) i B m ( α ) ( i ) E n ( α ) ( i ) ,

which is the assertion of the theorem. □

5 Further remarks and observations

In this section, we show that the same methods as in Sections 3 and 4 can be used to obtain a new identity involving both Bernoulli, Euler and Genocchi numbers.

The higher-order Genocchi polynomials G n ( α ) ( x ) are defined by the generating function

( 2 t e t + 1 ) α e x t = n = 0 G n ( α ) ( x ) t n n ! . (5.1)

The Genocchi polynomials G n ( x ) are given by G n ( x ) = G n ( 1 ) ( x ) . For x = 0 , we have the Genocchi numbers G n , i.e., G n = G n ( 0 ) . Letting x = 0 , we also have G n ( α ) = G n ( α ) ( 0 ) , where G n ( α ) denotes the higher-order Genocchi numbers.

The generating function of Genocchi polynomials is similar to those of Bernoulli and Euler polynomials, so it may be expected that the Genocchi numbers also satisfy similar identities as those established in Sections 3 and 4.

From the generating function (5.1), it is easy to deduce that

G n ( α ) ( x ) = l = 0 n ( n l ) G n l ( α ) x l , (5.2)

G n ( α ) ( x ) = ( 1 ) n + α G n ( α ) ( α x ) . (5.3)

Therefore, by the fermionic p-adic integral representation of G n ( α ) ( x ) , (5.2) and (5.3), we have

( 1 + ( 1 ) n ) l = 0 n ( n l ) G n l ( α ) E l = ( 1 ) n 2 i = 0 α ( 1 ) i G n ( α ) ( i ) . (5.4)

For , by (5.4), we have i = 0 α ( 1 ) i G n ( α ) ( i ) = 0 . From (5.4), it is easily seen that

l = 0 2 n ( 2 n l ) G 2 n l ( α ) E l = i = 0 α ( 1 ) i G 2 n ( α ) ( i ) , n N . (5.5)

Also, by the fermionic p-adic representations of B n ( α ) ( x ) and G n ( α ) ( x ) , we have

( 1 + ( 1 ) m + n ) k = 0 m l = 0 n ( m k ) ( n l ) B m k ( α ) G n l ( α ) E k + l = ( 1 ) m + n 2 i = 0 α ( 1 ) i B m ( α ) ( i ) G n ( α ) ( i ) , (5.6)

which is equivalent to

where m , n N 0 .

Competing interests

The authors declare that they have no competing interest.

Authors’ contributions

All authors contributed equally in this paper. They read and approved the final manuscript.

Acknowledgements

This work was supported by the Kyungnam University Foundation Grant, 2013.

References

  1. Sitaramachandrarao, R, Davis, B: Some identities involving the Riemann zeta function II. Indian J. Pure Appl. Math.. 17, 1175–1186 (1986)

  2. Choi, J, Kim, DS, Kim, T, Kim, YH: A note on some identities of Frobenius-Euler numbers and polynomials. Int. J. Math. Math. Sci.. 2012, (2012) Article ID 861797

  3. Dilcher, K: Sums of products of Bernoulli numbers. J. Number Theory. 60, 23–41 (1996). Publisher Full Text OpenURL

  4. Kim, DS, Dolgy, DV, Kim, H-M, Lee, S-H, Kim, T: Integral formulae of Bernoulli polynomials. Discrete Dyn. Nat. Soc.. 2012, (2012) Article ID 269847

  5. Kim, DS, Kim, T, Choi, J, Kim, YH: Some identities on Bernoulli and Euler numbers. Discrete Dyn. Nat. Soc.. 2012, (2012) Article ID 486158

  6. Kim, DS, Kim, T, Choi, J, Kim, YH: Identities involving q-Bernoulli and q-Euler numbers. Abstr. Appl. Anal.. 2012, (2012) Article ID 674210

  7. Kim, DS, Kim, T, Dolgy, DV, Lee, SH, Rim, S-H: Some properties and identities of Bernoulli and Euler polynomials associated with p-adic integral on Z p . Abstr. Appl. Anal.. 2012, (2012) Article ID 847901

  8. Kim, DS, Kim, T, Lee, S-H, Dolgy, DV, Rim, S-H: Some new identities on the Bernoulli and Euler numbers. Discrete Dyn. Nat. Soc.. 2011, (2011) Article ID 856132

  9. Kim, H-M, Kim, DS: Arithmetic identities involving Bernoulli and Euler numbers. Int. J. Math. Math. Sci.. 2012, (2012) Article ID 689797

  10. Kim, M-S: On Euler numbers, polynomials and related p-adic integrals. J. Number Theory. 129, 2166–2179 (2009). Publisher Full Text OpenURL

  11. Kim, M-S: A note on sums of products of Bernoulli numbers. Appl. Math. Lett.. 24(1), 55–61 (2011). Publisher Full Text OpenURL

  12. Kim, T: On the analogs of Euler numbers and polynomials associated with p-adic q-integral on Z p at q = 1 . J. Math. Anal. Appl.. 331, 779–792 (2007). Publisher Full Text OpenURL

  13. Kim, T, Kim, DS, Bayad, A, Rim, S-H: Identities on the Bernoulli and the Euler numbers and polynomials. Ars Comb.. CVII, 455–463 (2012)

  14. Lee, I, Kim, DS: Derivation of identities involving Bernoulli and Euler numbers. Int. J. Math. Math. Sci.. 2012, (2012) Article ID 598543

  15. Petojević, A: New sums of products of Bernoulli numbers. Integral Transforms Spec. Funct.. 19, 105–114 (2008). Publisher Full Text OpenURL

  16. Petojević, A, Srivastava, HM: Computation of Euler’s type sums of the products of Bernoulli numbers. Appl. Math. Lett.. 22, 796–801 (2009). Publisher Full Text OpenURL

  17. Raabe, JL: Zurückführung einiger Summen und bestmmtiem Integrale auf die Jacob-Bernoullische Function. J. Reine Angew. Math.. 42, 348–367 (1851)

  18. Simsek, Y: q-analogue of twisted l-series and q-twisted Euler numbers. J. Number Theory. 110, 267–278 (2005). Publisher Full Text OpenURL

  19. Srivastava, HM: Some formulas for the Bernoulli and Euler polynomials at rational arguments. Math. Proc. Camb. Philos. Soc.. 129, 77–84 (2000). Publisher Full Text OpenURL

  20. Srivastava, HM, Choi, J: Series Associated with the Zeta and Related Functions, Kluwer Academic, Dordrecht (2001)

  21. Srivastava, HM, Pinter, A: Remarks on some relationships between the Bernoulli and Euler polynomials. Appl. Math. Lett.. 17, 375–380 (2004). Publisher Full Text OpenURL

  22. Kim, M-S, Hu, S: Sums of products of Apostol-Bernoulli numbers. Ramanujan J.. 28, 113–123 (2012). Publisher Full Text OpenURL

  23. Miki, H: A relation between Bernoulli numbers. J. Number Theory. 10, 297–302 (1978). Publisher Full Text OpenURL

  24. Matiyasevich, Y: Identities with Bernoulli numbers. http://logic.pdmi.ras.ru/~yumat/Journal/Bernoulli/bernulli.htm

  25. Dunne, GV, Schubert, C: Bernoulli number identities from quantum field theory. Preprint (2004). arXiv:math.NT/0406610

  26. Crabb, MC: The Miki-Gessel Bernoulli number identity. Glasg. Math. J.. 47, 327–328 (2005). Publisher Full Text OpenURL

  27. Gessel, IM: On Miki’s identity for Bernoulli numbers. J. Number Theory. 110, 75–82 (2005). Publisher Full Text OpenURL

  28. Sun, Z-W, Pan, H: Identities concerning Bernoulli and Euler polynomials. Acta Arith.. 125, 21–39 (2006). Publisher Full Text OpenURL

  29. Cohen, H: Number Theory, Vol. II: Analytic and Modern Tools, Springer, New York (2007)

  30. Robert, AM: A Course in p-Adic Analysis, Springer, New York (2000)

  31. Shiratani, K, Yamamoto, S: On a p-adic interpolation function for the Euler numbers and its derivatives. Mem. Fac. Sci., Kyushu Univ., Ser. A, Math.. 39, 113–125 (1985)

  32. Osipov, JV: p-adic zeta functions. Usp. Mat. Nauk. 34, 209–210 (1979) (in Russian)

  33. Maïga, H: Some identities and congruences concerning Euler numbers and polynomials. J. Number Theory. 130, 1590–1601 (2010). Publisher Full Text OpenURL

  34. Kim, M-S, Hu, S: On p-adic Hurwitz-type Euler zeta functions. J. Number Theory. 132, 2977–3015 (2012). Publisher Full Text OpenURL

  35. Sun, Z-W: Introduction to Bernoulli and Euler polynomials. A Lecture Given in Taiwan on June 6, 2002. http://math.nju.edu.cn/~zwsun/BerE.pdf