# Sums of Products Involving Power Sums ofφ

Authors: Jitender Singh
Journal: Journal of Numbers
 (2014)

---

## 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Ψ
k
(
x
,
n
)
:
=
∑
d
|
n
‍
μ
(
d
)
d
k
S
k
x
/
d
,
n
∈
ℤ
+ which are defined via the Möbius function μ and the usual power sum S
k
(
x
) of a real or complex variable x. The power sum S
k
(
x
) is expressible in terms of the well-known Bernoulli polynomials by S
k
(
x
)
:
=
(
B
k
+
1
(
x
+
1
)
-
B
k
+
1
(
1
)
)
/
(
k
+
1
).


---

## Body

## 1. Introduction

Singh [1] introduced the power sum Ψ
k
(
x
,
n
) of real or complex variable x and positive integer n defined by the generating function
(1)
∑
d
∣
n
‍
μ
(
d
)
e
(
1
+
(
x
/
d
)
)
t
d
-
e
t
d
e
t
d
-
1
=
∑
k
=
0
∞
‍
Ψ
k
(
x
,
n
)
t
k
k
!
from which he derived the following closed form formula for these power sums:
(2)
Ψ
k
(
x
,
n
)
=
1
k
+
1
∑
m
=
0
[
k
/
2
]
‍
(
k
+
1
2
m
)
B
2
m
x
k
+
1
-
2
m
∏
p
∣
n
‍
(
1
-
p
2
m
-
1
)
for all k
=
0,1
,
… where B
m are the Bernoulli numbers and p runs over all prime divisors of n. In particular, Ψ
k
(
n
,
n
) gives the sum of kth power of those positive integers which are less than n and relatively prime to n. We will call Ψ
k
(
x
,
n
) as Möbius-Bernoulli power sums. Present work is aimed at describing sums of products of the power sums Ψ
k
(
x
,
n
) via introducing yet another sequence of rational numbers which we will call as the sequence of Möbius-Bernoulli numbers. The rational sequence {
B
k
} that appears in (2) is defined via the generating function t
/
(
e
t
-
1
)
=
∑
k
=
0
∞
‍
B
k
(
t
k
/
k
!
), |
t
|
<
2
π, 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]:
(3)
B
k
=
∑
m
=
1
k
1
m
+
1
∑
n
=
1
m
‍
(
-
1
)
n
(
m
n
)
n
k
,
k
=
0,1
,
…
.
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 φ
(
n
) 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 numbersM
k
(
n
), k
=
0,1
,
…, via the generating function
(4)
∑
d
∣
n
‍
t
μ
(
d
)
e
t
d
-
1
=
∑
k
=
0
∞
‍
M
k
(
n
)
t
k
k
!
,
|
t
|
<
2
π
n
,
∀
n
∈
ℤ
+
.We immediately notice from (4) that the Möbius-Bernoulli numbers are given by
(5)
M
k
(
1
)
=
B
k
;
M
k
(
n
)
=
B
k
∏
p
∣
n
‍
(
1
-
p
k
-
1
)
,
hhhhhhhhhhhhhhhhhhhh
∀
n
≥
2
,
k
=
0,1
,
…
.Note that, for a fixedk, the Möbius Bernoulli number M
k
(
n
) is a multiplicative function of n. Singh [1] has obtained the following identity relating the function Ψ
k
(
x
,
n
) to the Möbius-Bernoulli numbers. (
d
/
d
x
)
Ψ
k
(
x
,
n
)
=
k
Ψ
k
-
1
(
x
,
n
)
+
(
-
1
)
k
M
k
(
n
), from which we observe that Ψ
k
(
x
,
1
)
=
S
k
(
x
), Ψ
0
(
x
,
n
)
=
φ
(
n
), and Ψ
k
(
0
,
n
)
=
0 for all k
=
0,1
,
…, where φ
(
n
) 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 z and integer k by c
k
(
z
)
:
=
∑
q
∈
Φ
(
k
)
‍
e
2
π
ι
q
z
/
k, where Φ
(
k
)
=
{
q
∣
1
≤
q
≤
k
,
gcd
(
q
,
k
)
=
1
}. Möbius Bernoulli numbers are also related to Jordan’s totient (a generalization of Euler’s totient) J
k
(
n
)
:
=
n
k
∏
p
∣
n
‍
(
1
-
p
-
k
) by n
2
k
-
1
M
2
k
(
n
)
=
-
B
2
k
c
(
n
)
2
k
-
1
J
2
k
-
1
(
n
), where c
(
n
) is the square free part of n. The notion of Bernoulli polynomials generalizes to Möbius Bernoulli polynomials which we define next.Definition 2.
We sayM
k
(
x
,
n
) is the kth Möbius Bernoulli polynomial defined by the generating function
(6)
∑
d
∣
n
‍
μ
(
d
)
t
e
x
t
e
t
d
-
1
=
∑
k
=
0
∞
‍
M
k
(
x
,
n
)
t
k
k
!
,
|
t
|
<
2
π
n
,
hhhhhhhhhhhhhhhhhhhhhhhhhhh
x
∈
ℂ
,
n
∈
ℤ
+
.From (6), we see that M
0
(
x
,
n
)
=
M
0
(
n
), M
k
(
0
,
n
)
=
M
k
(
n
), and (
d
/
d
x
)
M
k
(
x
,
n
)
=
k
M
k
-
1
(
x
,
n
) for all k
∈
ℤ
+ and that M
k
(
x
,
1
)
=
B
k
(
x
), where B
k
(
x
) is the kth Bernoulli polynomial. It also follows from (6) that Möbius Bernoulli polynomials are given by
(7)
M
k
(
x
,
n
)
=
∑
d
∣
n
‍
μ
(
d
)
B
k
(
x
d
)
d
k
-
1
=
∑
m
=
0
k
‍
(
k
m
)
M
k
-
m
(
n
)
x
m
.IfΨ
k
(
x
,
n
) is as before, we have Ψ
0
(
x
,
n
)
=
M
1
(
x
,
n
)
-
M
1
(
0
,
n
)
=
x
φ
(
n
)
/
n and
(8)
Ψ
k
(
x
,
n
)
=
x
k
∑
d
∣
n
‍
μ
(
d
)
+
M
k
+
1
(
x
,
n
)
-
M
k
+
1
(
0
,
n
)
k
+
1
,
hhhhhhhhhhihhhhhhhhhhhhhhhhhhhh
k
=
1,2
,
…
.Definition 3.
Letn, N be positive integers, and let k be a nonnegative integer. Define higher order Möbius-Bernoulli numbers by
(9)
M
k
N
(
n
)
:
=
∑
∑
i
=
1
N
‍
k
i
=
k
‍
(
k
k
1
,
…
,
k
N
)
M
k
1
(
n
)
⋯
M
k
N
(
n
)
,
which are described by the generating function
(10)
H
(
t
)
=
(
∑
d
∣
n
‍
t
μ
(
d
)
e
t
d
-
1
)
N
=
∑
k
=
0
∞
‍
M
k
N
(
n
)
t
k
k
!
,
|
t
|
<
2
π
n
.Note thatM
k
N
(
1
) is the higher order Bernoulli number [4]. Also, M
k
1
(
n
)
=
M
k
(
n
) for all k
=
0,1
,
…. Some of the first few higher order Möbius-Bernoulli numbers are given by the following:
(11)
M
0
N
(
n
)
=
(
φ
(
n
)
)
N
,
M
2
N
(
n
)
=
N
(
φ
(
n
)
)
N
-
1
M
2
(
n
)
,
M
4
N
(
n
)
=
N
(
φ
(
n
)
)
N
-
2
{
3
(
N
-
1
)
(
M
2
(
n
)
)
2
+
φ
(
n
)
M
4
(
n
)
}
,
M
6
N
(
n
)
=
N
(
φ
(
n
)
)
N
-
3
×
{
15
(
N
-
1
)
(
N
-
2
)
(
M
2
(
n
)
)
2
hih
+
15
(
N
-
1
)
φ
(
n
)
M
2
(
n
)
M
4
(
n
)
hih
+
φ
(
n
)
2
M
6
(
n
)
(
M
2
(
n
)
)
2
}
etc
.Note thatM
k
N
(
n
)
=
lim
⁡
t
→
0
(
d
k
H
(
t
)
/
d
t
k
). 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.
LetN be a positive integer, and let f
:
ℝ
→
ℝ be a function of class C
k
,
k
≥
1. Then
(12)
D
k
(
f
(
x
)
)
N
=
N
!
∑
j
=
1
k
(
f
(
x
)
)
N
-
j
(
N
-
j
)
!
∑
∑
i
=
1
j
‍
k
i
=
k
(
k
k
1
,
…
,
k
j
)
hhhhhhhhhhhhhhhhhhhhhh
×
D
k
1
f
(
x
)
⋯
D
k
j
f
(
x
)
λ
(
k
1
)
!
⋯
λ
(
k
j
)
!
,
where D
k
=
d
k
/
d
x
k, k
=
1,2
,
…; λ
(
k
i
) is the multiplicity of occurrence of k
i in the partition {
k
1
,
…
,
k
j
} of n of length j; and λ
(
k
i
)
! contributes only once in the above product.Proof.
We use induction onk in proving the result. For k
=
1, we see that j
=
1 in the RHS of (12) and it reduces to N
(
f
(
x
)
)
N
-
1
D
1
(
f
(
x
)
)
=
D
1
(
f
(
x
)
N
). This proves that the result is true for k
=
1. Let us assume that formula (12) holds for all positive integers
≤
k. Now assume that f is of class C
k
+
1 and consider
(13)
D
k
+
1
(
(
f
(
x
)
)
N
)
=
N
!
∑
j
=
1
k
-
1
‍
(
f
(
x
)
)
N
-
(
j
+
1
)
(
N
-
(
j
+
1
)
)
!
hhhhhhhhh
×
∑
∑
i
=
1
j
k
i
=
k
(
k
k
1
,
…
,
k
j
)
hhhhhhhhhh
×
D
1
f
(
x
)
D
k
1
f
(
x
)
⋯
D
k
j
f
(
x
)
λ
(
k
1
)
!
⋯
λ
(
k
j
)
!
hhhhhhhhhhh
+
N
!
(
f
(
x
)
)
N
-
k
-
1
(
N
-
k
-
1
)
!
(
D
f
(
x
)
)
k
+
1
hhhhhhhhhhh
+
N
!
∑
j
=
1
k
(
f
(
x
)
)
N
-
j
(
N
-
j
)
!
hhhhhhhhhhhhhh
×
∑
∑
i
=
1
j
‍
k
i
=
k
(
k
k
1
,
…
,
k
j
)
hhhhhhhhhhhhhhhhhhhh
×
D
(
D
k
1
f
(
x
)
⋯
D
k
j
f
(
x
)
)
λ
(
k
1
)
!
⋯
λ
(
k
j
)
!
.
At this point, observe that any partitionπ
′ of k
+
1 can be obtained from a partition π of k by adjoining 1, and let us denote the set of all such partitions of k
+
1 by S. Denote by T the set of remaining all partitions of k
+
1 where each π
′ is obtained simply by adding 1 to exactly one member of π. In each of these cases one has k
+
1 choices of doing so for a fixed π. In the former case for each π
′
∈
S, |
π
′
|
=
|
π
|
+
1 which happens in the first summation above in (13) which reduces to the following:
(14)
(
k
+
1
)
N
!
∑
j
=
1
k
‍
(
f
(
x
)
)
N
-
j
(
N
-
j
)
!
hhhhhhhhhh
×
∑
{
k
1
′
,
…
,
k
j
+
1
′
}
∈
S
‍
(
k
k
1
′
,
…
,
k
j
+
1
′
)
hhhhhhhhhhhhhh
×
D
k
1
′
f
(
x
)
⋯
D
k
j
+
1
′
f
(
x
)
λ
(
k
1
′
)
!
⋯
λ
(
k
j
′
)
!
.
In the latter case for eachπ
′
∈
T, |
π
′
|
=
|
π
|, and the terms after first summation in (13) reduce to
(15)
(
k
+
1
)
N
!
∑
j
=
1
k
+
1
‍
(
f
(
x
)
)
N
-
j
(
N
-
j
)
!
hhhhhhhhh
×
∑
{
k
1
′
,
…
,
k
j
+
1
′
}
∈
T
‍
(
k
k
1
′
,
…
,
k
j
+
1
′
)
hhhhhhhhhhhhhh
×
D
k
1
′
f
(
x
)
⋯
D
k
j
+
1
′
f
(
x
)
λ
(
k
1
′
)
!
⋯
λ
(
k
j
′
)
!
,
where the term N
!
(
(
f
(
x
)
)
N
-
k
-
1
/
(
N
-
k
-
1
)
!
)
(
D
f
(
x
)
)
k
+
1 corresponds to j
=
k
+
1. The result follows by substituting (14) and (15) into (13). This completes the final step of induction.Theorem 5.
For each positive integersN and n, the higher order Möbius-Bernoulli numbers are given by
(16)
M
k
N
(
n
)
=
N
!
∑
j
=
1
k
‍
(
φ
(
n
)
)
N
-
j
(
N
-
j
)
!
hhhh
×
∑
∑
i
=
1
j
‍
k
i
=
k
(
k
k
1
,
…
,
k
j
)
hhihhhhhh
×
M
k
1
(
n
)
⋯
M
k
j
(
n
)
λ
(
k
1
)
!
⋯
λ
(
k
j
)
!
.Proof.
First note from definition thatM
k
N
(
n
)
=
lim
⁡
t
→
0
D
k
(
H
(
0
)
), the result follows at once by applying Lemma 4 to the function f
(
t
)
=
∑
d
∣
n
‍
(
t
μ
(
d
)
/
(
e
t
d
-
1
)
), 0
<
t
<
2
π
/
n, and then taking limit t
→
0 throughout and using φ
(
n
)
=
n
∑
d
∣
n
‍
(
μ
(
d
)
/
d
) therein.Proposition 6.
M
2
k
-
1
N
(
n
)
=
0 for all positive integers k and N and n
>
1.Proof.
Observe from (10) that, for a positive integer N, the following holds:
(17)
H
(
-
t
)
=
(
t
∑
d
∣
n
‍
μ
(
d
)
+
∑
d
∣
n
‍
t
μ
(
d
)
e
t
d
-
1
)
N
=
(
t
δ
1
n
+
∑
d
∣
n
‍
t
μ
(
d
)
e
t
d
-
1
)
N
=
H
(
t
)
for all n
>
1 where the arithmetic function ∑
d
∣
n
‍
μ
(
d
)
=
δ
1
n is the Kronecker delta. We have proved that H is an even function of t for n
>
1. Thus the coefficient of t
2
k
-
1 in the RHS of (10) (which is precisely M
2
k
-
1
N
(
n
)) vanishes for each k
=
1,2
,
….Remark 7.
If we extend the definition of higher Möbius-Bernoulli numbers to complexN
≠
0, the formula (18) for M
k
N
(
n
) is still valid just on replacing N
!
/
(
N
-
j
)
! by N
(
N
-
1
)
⋯
(
N
-
j
+
1
) in it. In this regard we note that M
2
k
-
1
N
(
n
>
1
)
=
0 holds for all k
=
1,2
,
… and N
∈
ℂ.In view of the Theorem5 and the Proposition 6, we have for all positive integers n
>
1 and k
=
1,2
,
…,
(18)
M
2
k
N
(
n
)
=
N
!
∑
j
=
1
k
(
φ
(
n
)
)
N
-
2
j
(
N
-
2
j
)
!
hihh
×
∑
∑
i
=
1
j
k
i
=
k
‍
(
2
k
2
k
1
,
…
,
2
k
j
)
hhhhhhhh
×
M
2
k
1
(
n
)
⋯
M
2
k
j
(
n
)
λ
(
2
k
1
)
!
⋯
λ
(
2
k
j
)
!
.
As an example, let k
=
4. There are five partitions of 4 which are given by 4
1, 3
1
1
1, 2
2, 2
1
1
2, and 1
4, and therefore from (18), we obtain
(19)
M
8
N
(
n
)
=
N
!
{
(
M
2
(
n
)
)
4
4
!
φ
(
n
)
N
-
2
(
N
-
2
)
!
M
8
(
n
)
+
φ
(
n
)
N
-
4
(
N
-
4
)
!
hihhhhh
×
(
8
!
6
!
2
!
M
6
(
n
)
M
2
(
n
)
+
8
!
4
!
4
!
(
M
4
(
n
)
)
2
2
!
)
hihhhhh
+
φ
(
n
)
N
-
6
(
N
-
6
)
!
(
8
!
4
!
2
!
2
!
M
4
(
n
)
(
M
2
(
n
)
)
2
2
!
)
hihhhhh
+
φ
(
n
)
N
-
8
(
N
-
8
)
!
8
!
2
!
2
!
2
!
2
!
(
M
2
(
n
)
)
4
4
!
}
=
N
!
{
φ
(
n
)
N
-
2
(
N
-
2
)
!
M
8
(
n
)
+
φ
(
n
)
N
-
4
(
N
-
4
)
!
hihhhhh
×
(
28
M
6
(
n
)
M
2
(
n
)
+
35
(
M
4
(
n
)
(
M
4
(
n
)
)
2
)
hihhhhh
+
φ
(
n
)
N
-
6
(
N
-
6
)
!
(
210
M
4
(
n
)
(
M
2
(
n
)
)
2
)
hihhhhh
+
φ
(
n
)
N
-
8
(
N
-
8
)
!
105
(
M
2
(
n
)
)
4
}
.Remark 8.
The formula (18) is not suitable for explicit evaluation of M
k
N for large k. Because number of partitions of k increases at a faster rate than k. For example the number of partitions of 10 is 42, which is the number of terms in the expression for M
20
N. 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.
Ifn
=
p
s for some positive integer s and prime p, then the simplest possible formula (18) for the higher order Möbius-Bernoulli numbers can be found as follows:
(20)
H
(
t
)
|
n
=
p
s
=
(
t
e
t
-
1
-
t
e
t
p
-
1
)
N
=
∑
m
=
0
N
(
N
m
)
(
t
e
t
-
1
)
m
(
t
p
e
t
p
-
1
)
N
-
m
(
-
p
)
m
-
N
.
Therefore
(21)
M
k
N
(
p
s
)
=
lim
⁡
t
→
0
d
k
H
(
t
)
d
t
k
=
∑
m
=
0
N
‍
(
N
m
)
(
-
p
)
m
-
N
∑
j
=
0
k
‍
(
k
j
)
B
j
m
B
k
-
j
N
-
m
p
k
-
j
,
where we have utilized the Leibniz product rule for higher order derivatives and B
j
m is the higher order Bernoulli number given by (see for more details Srivastava and Todorov [4])
(22)
B
j
m
=
∑
ℓ
=
0
j
(
j
+
m
j
-
ℓ
)
(
j
+
m
-
1
ℓ
)
×
j
!
(
j
+
ℓ
)
!
∑
h
=
0
ℓ
(
-
1
)
h
(
ℓ
h
)
h
j
+
ℓ
.
Similarly, if we take n
=
p
1
s
1
p
2
s
2 for some positive integers s
1, s
2 and distinct primes p
1, p
2, then
(23)
H
(
t
)
|
n
=
p
1
s
1
p
2
s
2
=
(
t
e
t
-
1
-
t
e
t
p
1
-
1
-
t
e
t
p
2
-
1
+
t
e
t
p
1
p
2
-
1
)
N
=
∑
∑
i
=
0
3
‍
m
i
=
N
(
N
m
0
,
…
,
m
3
)
(
t
e
t
-
1
)
m
0
×
(
t
p
1
e
t
p
1
-
1
)
m
1
(
t
p
2
e
t
p
2
-
1
)
m
2
×
(
t
p
1
p
2
e
t
p
1
p
2
-
1
)
m
3
(
-
p
1
)
-
m
1
(
-
p
2
)
-
m
2
(
p
1
p
2
)
-
m
3
which gives
(24)
M
k
N
(
p
1
s
1
p
2
s
2
)
=
∑
∑
i
=
0
3
‍
m
i
=
N
(
N
m
0
,
⋯
,
m
3
)
lim
⁡
t
→
0
⁡
D
k
hhhhihhhhh
×
{
(
t
e
t
-
1
)
m
0
(
t
p
1
e
t
p
1
-
1
)
m
1
hhhhhhhhhhhh
×
(
t
p
2
e
t
p
2
-
1
)
m
2
(
t
p
1
p
2
e
t
p
1
p
2
-
1
)
m
3
}
hhhhhhhhh
×
(
-
p
1
)
-
m
1
(
-
p
2
)
-
m
2
(
p
1
p
2
)
-
m
3
=
∑
∑
i
=
0
3
‍
m
i
=
N
(
N
m
0
,
⋯
,
m
3
)
×
∑
∑
j
=
0
3
‍
k
i
=
k
(
k
k
0
,
⋯
,
k
3
)
hihhhhhhhhhhh
×
B
k
0
m
0
B
k
1
m
1
B
k
2
m
2
B
k
3
m
3
p
1
k
1
-
m
1
p
2
k
2
-
m
2
hihhhhhhhhhhh
×
(
p
1
p
2
)
k
3
-
m
3
(
-
1
)
m
1
+
m
2
.
These formulas involve products of higher order Bernoulli numbers. So in general, the formulas for M
k
N
(
n
) 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:
(25)
M
k
N
(
∏
α
=
1
β
‍
p
α
e
α
)
=
N
!
k
!
(
N
+
k
)
!
×
∑
∑
i
=
0
2
β
-
1
‍
m
i
=
N
‍
∑
∑
j
=
0
2
β
-
1
‍
k
i
=
k
‍
(
N
+
k
m
0
,
…
,
m
2
β
-
1
,
k
0
,
…
,
k
2
β
-
1
)
hhihhhhhhhhhhhhhhh
×
B
k
0
m
0
B
k
1
m
1
⋯
B
k
2
β
-
1
m
2
β
-
1
χ
hhihhhhhhhhhhhhhhh
×
(
k
1
,
…
,
k
2
β
-
1
,
m
1
,
…
,
m
2
β
-
1
)
,
where
(26)
χ
(
k
1
,
…
,
k
2
β
-
1
,
m
1
,
…
,
m
2
β
-
1
)
=
∏
s
0
=
1
β
‍
(
-
p
s
0
)
k
s
0
-
m
s
0
hhihhhhhh
⋯
∏
s
β
-
2
<
s
β
-
1
=
β
β
(
(
-
1
)
β
p
s
0
p
s
1
⋯
p
s
β
-
1
)
k
2
β
-
1
-
m
2
β
-
1
.


## 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Ψ
k
(
x
,
n
).Definition 10.
We define sums of products of the Möbius Bernoulli power sums asΨ
k
N
(
x
,
n
)
:
=
∑
∑
i
=
1
N
‍
k
i
=
k
‍
(
k
k
1
⋯
k
N
)
Ψ
k
1
(
x
,
n
)
⋯
Ψ
k
N
(
x
,
n
) for nonnegative integers k and N which are described by the generating function
(27)
(
∑
d
∣
n
‍
μ
(
d
)
e
(
1
+
(
x
/
d
)
)
t
d
-
e
t
d
e
t
d
-
1
)
N
=
∑
k
=
0
∞
‍
Ψ
k
N
(
x
,
n
)
t
k
k
!
.The next result evaluates the sums of productsΨ
k
N
(
x
,
n
).Theorem 11.
For a positive integerN and nonnegative integer k,
(28)
Ψ
k
N
(
x
,
n
)
=
k
!
N
!
(
k
+
N
)
!
∑
j
=
0
k
(
k
+
N
j
)
M
j
N
(
n
)
S
(
k
+
N
-
j
,
N
)
x
k
+
N
-
j
for all n
=
2,3
,
… where S
(
ℓ
,
m
) are the Stirling numbers of second kind.Proof.
Observe from the generating function forΨ
k
N
(
x
) that
(29)
(
∑
d
∣
n
‍
μ
(
d
)
e
(
1
+
(
x
/
d
)
)
t
d
-
e
t
d
e
t
d
-
1
)
N
=
(
∑
d
∣
n
‍
μ
(
d
)
t
e
t
d
e
t
d
-
1
)
N
(
e
x
t
-
1
t
)
N
=
(
t
∑
d
∣
n
‍
μ
(
d
)
+
∑
d
∣
n
‍
t
μ
(
d
)
e
t
d
-
1
)
N
(
e
x
t
-
1
t
)
N
=
(
t
δ
1
n
+
∑
d
∣
n
‍
t
μ
(
d
)
e
t
d
-
1
)
N
(
e
x
t
-
1
t
)
N
=
∑
m
=
0
N
‍
(
N
m
)
t
N
-
m
δ
1
n
N
-
m
(
∑
d
∣
n
‍
t
μ
(
d
)
e
t
d
-
1
)
m
(
e
x
t
-
1
t
)
N
which on further simplification gives
(30)
∑
k
=
0
∞
‍
Ψ
k
N
(
x
,
n
)
t
k
k
!
=
∑
k
=
0
∞
‍
∑
m
=
0
N
‍
(
N
m
)
∑
j
=
0
k
‍
(
k
j
)
δ
1
n
M
k
-
j
m
(
n
)
hhhhhhhhhhhhhhh
×
j
!
N
!
(
j
+
N
)
!
S
(
j
+
N
,
N
)
x
j
+
N
t
k
k
!
,
where we have used the identity (
(
e
x
t
-
1
)
/
t
)
N
=
∑
j
=
0
∞
‍
(
(
k
!
N
!
/
(
j
+
N
)
!
)
S
(
j
+
N
,
N
)
x
j
+
N
(
t
j
/
j
!
)
). On comparing like powers of t in (30) we obtain
(31)
Ψ
k
N
(
x
,
n
)
=
∑
m
=
0
N
(
N
m
)
∑
j
=
0
k
(
k
j
)
(
δ
1
n
)
N
-
m
M
k
-
j
m
(
n
)
j
!
N
!
(
j
+
N
)
!
hhhhhhhhhhhhhh
×
S
(
j
+
N
,
N
)
x
j
+
N
for all k
=
0,1
,
…; n
,
N
∈
ℤ
+, where we define δ
1
n
0
:
=
1 for all n
∈
ℤ
+. The result follows now.Note that, from Theorem11, we recover for N
=
1,
(32)
Ψ
k
1
(
x
,
n
)
=
1
k
+
1
∑
j
=
0
k
‍
(
k
+
1
j
)
x
k
+
1
-
j
B
j
∏
p
∣
n
‍
(
1
-
p
j
-
1
)
,
where we have used M
j
1
(
n
)
=
B
j
∏
p
∣
n
‍
(
1
-
p
j
-
1
) and S
(
k
+
1
-
j
,
1
)
=
1 for all j
=
0,1
,
…
,
k.



---

*Source: 158351-2014-02-27.xml*