Dirichlet convolution explained

In mathematics, the Dirichlet convolution (or divisor convolution) is a binary operation defined for arithmetic functions; it is important in number theory. It was developed by Peter Gustav Lejeune Dirichlet.

Definition

If

f,g:N\toC

are two arithmetic functions from the positive integers to the complex numbers, the Dirichlet convolution is a new arithmetic function defined by:

(f*g)(n) = \sumd\midf(d)g\left(

n
d

\right) = \sumab=nf(a)g(b)

where the sum extends over all positive divisors d of n, or equivalently over all distinct pairs of positive integers whose product is n.

This product occurs naturally in the study of Dirichlet series such as the Riemann zeta function. It describes the multiplication of two Dirichlet series in terms of their coefficients:

\left(\sumn\geq

f(n)
ns

\right) \left(\sumn\geq

g(n)
ns

\right)  = \left(\sumn\geq

(f*g)(n)
ns

\right).

Properties

The set of arithmetic functions forms a commutative ring, the , under pointwise addition, where is defined by, and Dirichlet convolution. The multiplicative identity is the unit function ε defined by if and if . The units (invertible elements) of this ring are the arithmetic functions f with .

Specifically,[1] Dirichlet convolution is associative,

(f*g)*h=f*(g*h),

distributive over addition

f*(g+h)=f*g+f*h

,commutative,

f*g=g*f

,and has an identity element,

f*\varepsilon

=

\varepsilon*f=f

.Furthermore, for each

f

having

f(1)0

, there exists an arithmetic function

f-1

with

f*f-1=\varepsilon

, called the of

f

.

The Dirichlet convolution of two multiplicative functions is again multiplicative, and every not constantly zero multiplicative function has a Dirichlet inverse which is also multiplicative. In other words, multiplicative functions form a subgroup of the group of invertible elements of the Dirichlet ring. Beware however that the sum of two multiplicative functions is not multiplicative (since

(f+g)(1)=f(1)+g(1)=21

), so the subset of multiplicative functions is not a subring of the Dirichlet ring. The article on multiplicative functions lists several convolution relations among important multiplicative functions.

h

, pointwise multiplication by

h

distributes over Dirichlet convolution:

(f*g)h=(fh)*(gh)

.[2] The convolution of two completely multiplicative functions is multiplicative, but not necessarily completely multiplicative.

Properties and Examples

In these formulas, we use the following arithmetical functions:

\varepsilon

is the multiplicative identity:

\varepsilon(1)=1

, otherwise 0 (

\varepsilon(n)=\lfloor\tfrac1n\rfloor

).

1

is the constant function with value 1:

1(n)=1

for all

n

. Keep in mind that

1

is not the identity. (Some authors denote this as

\zeta

because the associated Dirichlet series is the Riemann zeta function.)

1C

for

C\subsetN

is a set indicator function:

1C(n)=1

iff

n\inC

, otherwise 0.

Id

is the identity function with value n:

Id(n)=n

.

Idk

is the kth power function:
k
Id
k(n)=n
.

The following relations hold:

1*\mu=\varepsilon

, the Dirichlet inverse of the constant function

1

is the Möbius function (see proof). Hence:

g=f*1

if and only if

f=g*\mu

, the Möbius inversion formula

\sigmak=Idk*1

, the kth-power-of-divisors sum function σk

\sigma=Id*1

, the sum-of-divisors function

\tau=1*1

, the number-of-divisors function

Idk=\sigmak*\mu

,  by Möbius inversion of the formulas for σk, σ, and τ

Id=\sigma*\mu

1=\tau*\mu

\phi*1=Id

, proved under Euler's totient function

\phi=Id*\mu

, by Möbius inversion

\sigma=\phi*\tau

 , from convolving 1 on both sides of

\phi*1=Id

λ*|\mu|=\varepsilon

 where λ is Liouville's function

λ*1=1Sq

 where Sq = is the set of squares

Idk*(Idk\mu)=\varepsilon

\tau3*1=(\tau*1)2

Jk*1=Idk

, Jordan's totient function

(IdsJr)*Js=Js

Λ*1=log

, where

Λ

is von Mangoldt's function

|\mu|\ast1=2\omega,

where

\omega(n)

is the prime omega function counting distinct prime factors of n

\Omega\ast\mu=1l{P

}, the characteristic function of the prime powers.

\omega\ast\mu=1P

where

1P(n)\mapsto\{0,1\}

is the characteristic function of the primes.

This last identity shows that the prime-counting function is given by the summatory function

\pi(x)=\sumn(\omega\ast\mu)(n)=

x
\sum
d=1

\omega(d)M\left(\left\lfloor

x
d

\right\rfloor\right)

where

M(x)

is the Mertens function and

\omega

is the distinct prime factor counting function from above. This expansion follows from the identity for the sums over Dirichlet convolutions given on the divisor sum identities page (a standard trick for these sums).[3]

Dirichlet inverse

Examples

Given an arithmetic function

f

its Dirichlet inverse

g=f-1

may be calculated recursively: the value of

g(n)

is in terms of

g(m)

for

m<n

.

For

n=1

:

(f*g)(1)=f(1)g(1)=\varepsilon(1)=1

, so

g(1)=1/f(1)

. This implies that

f

does not have a Dirichlet inverse if

f(1)=0

.

For

n=2

:

(f*g)(2)=f(1)g(2)+f(2)g(1)=\varepsilon(2)=0

,

g(2)=-(f(2)g(1))/f(1)

,

For

n=3

:

(f*g)(3)=f(1)g(3)+f(3)g(1)=\varepsilon(3)=0

,

g(3)=-(f(3)g(1))/f(1)

,

For

n=4

:

(f*g)(4)=f(1)g(4)+f(2)g(2)+f(4)g(1)=\varepsilon(4)=0

,

g(4)=-(f(4)g(1)+f(2)g(2))/f(1)

,

and in general for

n>1

,

g(n) = 

-1
f(1)
\sumd\midf\left(
d<n
n
d

\right)g(d).

Properties

The following properties of the Dirichlet inverse hold:[4]

(f\astg)-1=f-1\astg-1

.

f-1(n)=\mu(n)f(n)

.

(fg)-1=fg-1

whenever

g(1)0

and where

denotes pointwise multiplication of functions.

Other formulas

Arithmetic function Dirichlet inverse:[5]
Constant function with value 1 Möbius function μ

n\alpha

\mu(n)n\alpha

Liouville's function λ Absolute value of Möbius function
Euler's totient function

\varphi

\sumd|nd\mu(d)

The generalized sum-of-divisors function

\sigma\alpha

\sumd|nd\alpha\mu(d)\mu\left(

n
d

\right)

An exact, non-recursive formula for the Dirichlet inverse of any arithmetic function f is given in Divisor sum identities. A more partition theoretic expression for the Dirichlet inverse of f is given by

f-1(n)=

\Omega(n)
\sum
k=1

\left\{

\sum
{λ1+2λ2+ … +kλk=n

\atop{λ1,λ2,\ldots,λk|n}}

(λ1+λ2+ … +λk)!
1!2!k!

(-1)kf(λ1)

2
f(λ
2)

k\right\}.
f(λ
k)
The following formula provides a compact way of expressing the Dirichlet inverse of an invertible arithmetic function f :

f-1

+infty
=\sum
k=0
(f(1)\varepsilon-f)*k
f(1)k+1

where the expression

(f(1)\varepsilon-f)*k

stands for the arithmetic function

f(1)\varepsilon-f

convoluted with itself k times. Notice that, for a fixed positive integer

n

, if

k>\Omega(n)

then

(f(1)\varepsilon-f)*k(n)=0

, this is because

f(1)\varepsilon(1)-f(1)=0

and every way of expressing n as a product of k positive integers must include a 1, so the series on the right hand side converges for every fixed positive integer n.

Dirichlet series

If f is an arithmetic function, the Dirichlet series generating function is defined by

DG(f;s)=

infty
\sum
n=1
f(n)
ns

for those complex arguments s for which the series converges (if there are any). The multiplication of Dirichlet series is compatible with Dirichlet convolution in the following sense:

DG(f;s)DG(g;s)=DG(f*g;s)

for all s for which both series of the left hand side converge, one of them at least converging absolutely (note that simple convergence of both series of the left hand side does not imply convergence of the right hand side!). This is akin to the convolution theorem if one thinks of Dirichlet series as a Fourier transform.

Related concepts

The restriction of the divisors in the convolution to unitary, bi-unitary or infinitary divisors defines similar commutative operations which share many features with the Dirichlet convolution (existence of a Möbius inversion, persistence of multiplicativity, definitions of totients, Euler-type product formulas over associated primes, etc.).

Dirichlet convolution is a special case of the convolution multiplication for the incidence algebra of a poset, in this case the poset of positive integers ordered by divisibility.

The Dirichlet hyperbola method computes the summation of a convolution in terms of its functions and their summation functions.

See also

References

. Hugh L. Montgomery . Hugh Montgomery (mathematician) . Robert C. Vaughan . Robert Charles Vaughan (mathematician) . Multiplicative number theory I. Classical theory . Cambridge tracts in advanced mathematics . 97 . 2007 . 978-0-521-84903-6 . 38 . Cambridge Univ. Press . Cambridge .

Notes and References

  1. Proofs are in Chan, ch. 2
  2. A proof is in the article Completely multiplicative function#Proof of distributive property.
  3. Book: Apostol's Introduction to Analytic Number Theory. Schmidt . Maxie. This identity is a little special something I call "croutons". It follows from several chapters worth of exercises in Apostol's classic book.
  4. Again see Apostol Chapter 2 and the exercises at the end of the chapter.
  5. See Apostol Chapter 2.