Abstract

A sequence of rational numbers as a generalization of the sequence of Bernoulli numbers is introduced. Sums of products involving the terms of this generalized sequence are then obtained using an application of Faà di Bruno's formula. These sums of products are analogous to the higher order Bernoulli numbers and are used to develop the closed form expressions for the sums of products involving the power sums which are defined via the Möbius function μ and the usual power sum of a real or complex variable . The power sum is expressible in terms of the well-known Bernoulli polynomials by .

1. Introduction

Singh [1] introduced the power sum of real or complex variable and positive integer defined by the generating function from which he derived the following closed form formula for these power sums: for all where are the Bernoulli numbers and runs over all prime divisors of . In particular, gives the sum of th power of those positive integers which are less than and relatively prime to . We will call as Möbius-Bernoulli power sums. Present work is aimed at describing sums of products of the power sums via introducing yet another sequence of rational numbers which we will call as the sequence of Möbius-Bernoulli numbers. The rational sequence that appears in (2) is defined via the generating function , , and was known to Faulhaber and Bernoulli. Many explicit formulas for the Bernoulli numbers are also well known in the literature. One such formula is as follows [2]: The rest of the paper is organized as follows. Möbius Bernoulli numbers are introduced in Section 2 and their sums of products are discussed via Faà di Bruno’s formula. In Section 3 sums of products of power sums of integers are obtained in closed form using sums of products of Möbius Bernoulli numbers.

2. Möbius-Bernoulli Numbers

Definition 1. We define Möbius-Bernoulli numbers , , via the generating function

We immediately notice from (4) that the Möbius-Bernoulli numbers are given by

Note that, for a fixed , the Möbius Bernoulli number is a multiplicative function of . Singh [1] has obtained the following identity relating the function to the Möbius-Bernoulli numbers. , from which we observe that , , and for all , where is Euler’s totient. Use of Möbius Bernoulli numbers is inherent in studies recently done by Alkan [3] on averages of Ramanujan sums which are defined for any complex number and integer by , where . Möbius Bernoulli numbers are also related to Jordan’s totient (a generalization of Euler’s totient) by , where is the square free part of . The notion of Bernoulli polynomials generalizes to Möbius Bernoulli polynomials which we define next.

Definition 2. We say is the th Möbius Bernoulli polynomial defined by the generating function

From (6), we see that , , and for all and that , where is the th Bernoulli polynomial. It also follows from (6) that Möbius Bernoulli polynomials are given by

If is as before, we have and

Definition 3. Let , be positive integers, and let be a nonnegative integer. Define higher order Möbius-Bernoulli numbers by which are described by the generating function

Note that is the higher order Bernoulli number [4]. Also, for all . Some of the first few higher order Möbius-Bernoulli numbers are given by the following:

Note that . In this regard, a formula for the higher order Möbius-Bernoulli numbers can be obtained from the following version of the well-known Faà di Bruno’s formula [5].

Lemma 4. Let be a positive integer, and let be a function of class . Then where , ; is the multiplicity of occurrence of in the partition of of length ; and contributes only once in the above product.

Proof. We use induction on in proving the result. For , we see that in the RHS of (12) and it reduces to . This proves that the result is true for . Let us assume that formula (12) holds for all positive . Now assume that is of class and consider
At this point, observe that any partition of can be obtained from a partition of by adjoining , and let us denote the set of all such partitions of by . Denote by the set of remaining all partitions of where each is obtained simply by adding to exactly one member of . In each of these cases one has choices of doing so for a fixed . In the former case for each , which happens in the first summation above in (13) which reduces to the following:
In the latter case for each , , and the terms after first summation in (13) reduce to where the term corresponds to . The result follows by substituting (14) and (15) into (13). This completes the final step of induction.

Theorem 5. For each positive integers and , the higher order Möbius-Bernoulli numbers are given by

Proof. First note from definition that , the result follows at once by applying Lemma 4 to the function , , and then taking limit throughout and using therein.

Proposition 6. for all positive integers and and .

Proof. Observe from (10) that, for a positive integer , the following holds: for all where the arithmetic function is the Kronecker delta. We have proved that is an even function of for . Thus the coefficient of in the RHS of (10) (which is precisely ) vanishes for each .

Remark 7. If we extend the definition of higher Möbius-Bernoulli numbers to complex , the formula (18) for is still valid just on replacing by in it. In this regard we note that holds for all and .

In view of the Theorem 5 and the Proposition 6, we have for all positive integers and , As an example, let . There are five partitions of which are given by , , , , and , and therefore from (18), we obtain

Remark 8. The formula (18) is not suitable for explicit evaluation of for large . Because number of partitions of increases at a faster rate than . For example the number of partitions of is , which is the number of terms in the expression for . In this regard, it will be good to see a formula for the higher order Möbius Bernoulli numbers which can describe them better than the one we have given above.

Remark 9. If for some positive integer and prime , then the simplest possible formula (18) for the higher order Möbius-Bernoulli numbers can be found as follows: Therefore where we have utilized the Leibniz product rule for higher order derivatives and is the higher order Bernoulli number given by (see for more details Srivastava and Todorov [4]) Similarly, if we take for some positive integers , and distinct primes , , then which gives These formulas involve products of higher order Bernoulli numbers. So in general, the formulas for involve sums containing product of several higher order Bernoulli numbers and such a formula in the above sense would be complicated and will take the following form: where

3. Sums of Products

Having developed the expressions for the Möbius Bernoulli numbers in the previous section, we will now use them in expressing the sums of products of the Möbius-Bernoulli power sums .

Definition 10. We define sums of products of the Möbius Bernoulli power sums as for nonnegative integers and which are described by the generating function

The next result evaluates the sums of products .

Theorem 11. For a positive integer and nonnegative integer , for all where are the Stirling numbers of second kind.

Proof. Observe from the generating function for that which on further simplification gives where we have used the identity . On comparing like powers of in (30) we obtain for all ; , where we define for all . The result follows now.

Note that, from Theorem 11, we recover for , where we have used and for all .

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Many suggestions regarding presentation of the paper by Professor László Tóth are gratefully acknowledged. The author is also thankful to the anonymous referees for their suggestions.