Multiplicative Number Theory I-Montgomery
Multiplicative Number Theory I-Montgomery
Prime numbers are the multiplicative building blocks of natural numbers. Un-
derstanding their overall influence and especially their distribution gives rise
to central questions in mathematics and physics. In particular their finer distri-
bution is closely connected with the Riemann hypothesis, the most important
unsolved problem in the mathematical world. Assuming only subjects covered
in a standard degree in mathematics, the authors comprehensively cover all the
topics met in first courses on multiplicative number theory and the distribution
of prime numbers. They bring their extensive and distinguished research exper-
tise to bear in preparing the student for intelligent reading of the more advanced
research literature. The text, which is based on courses taught successfully over
many years at Michigan, Imperial College and Pennsylvania State, is enriched
by comprehensive historical notes and references as well as over 500 exercises.
HUGH L. MONTGOMERY
University of Michigan, Ann Arbor
ROBERT C. VAUGHAN
Pennsylvania State University, University Park
cambridge university press
Cambridge, New York, Melbourne, Madrid, Cape Town, Singapore, São Paulo
Cambridge University Press has no responsibility for the persistence or accuracy of urls
for external or third-party internet websites referred to in this publication, and does not
guarantee that any content on such websites is, or will remain, accurate or appropriate.
Dedicated to our teachers:
P. T. Bateman
J. H. H. Chalk
H. Davenport
T. Estermann
H. Halberstam
A. E. Ingham
Talet är tänkandets början och slut.
Med tanken föddes talet.
Utöfver talet når tanken icke.
Preface page xi
List of notation xiii
1 Dirichlet series: I 1
1.1 Generating functions and asymptotics 1
1.2 Analytic properties of Dirichlet series 11
1.3 Euler products and the zeta function 19
1.4 Notes 31
1.5 References 33
2 The elementary theory of arithmetic functions 35
2.1 Mean values 35
2.2 The prime number estimates of Chebyshev and of Mertens 46
2.3 Applications to arithmetic functions 54
2.4 The distribution of (n) − ω(n) 65
2.5 Notes 68
2.6 References 71
3 Principles and first examples of sieve methods 76
3.1 Initiation 76
3.2 The Selberg lambda-squared method 82
3.3 Sifting an arithmetic progression 89
3.4 Twin primes 91
3.5 Notes 101
3.6 References 104
4 Primes in arithmetic progressions: I 108
4.1 Additive characters 108
4.2 Dirichlet characters 115
4.3 Dirichlet L-functions 120
vii
viii Contents
14 Zeros 452
14.1 General distribution of the zeros 452
14.2 Zeros on the critical line 456
14.3 Notes 460
14.4 References 461
APPENDICES
A The Riemann–Stieltjes integral 486
A.1 Notes 492
A.2 References 493
Our object is to introduce the interested student to the techniques, results, and
terminology of multiplicative number theory. It is not intended that our discus-
sion will always reach the research frontier. Rather, it is hoped that the material
here will prepare the student for intelligent reading of the more advanced re-
search literature.
Analytic number theorists are not very uniformly distributed around the
world and it possible that a student may be working without the guidance of an
experienced mentor in the area. With this in mind, we have tried to make this
volume as self-contained as possible.
We assume that the reader has some acquaintance with the fundamentals of
elementary number theory, abstract algebra, measure theory, complex analysis,
and classical harmonic analysis. More specialized or advanced background
material in analysis is provided in the appendices.
The relationship of exercises to the material developed in a given section
varies widely. Some exercises are designed to illustrate the theory directly
whilst others are intended to give some idea of the ways in which the theory can
be extended, or developed, or paralleled in other areas. The reader is cautioned
that papers cited in exercises do not necessarily contain a solution.
This volume is the first instalment of a larger project. We are preparing a
second volume, which will cover such topics as uniform distribution, bounds for
exponential sums, a wider zero-free region for the Riemann zeta function, mean
and large values of Dirichlet polynomials, approximate functional equations,
moments of the zeta function and L functions on the line σ = 1/2, the large
sieve, Vinogradov’s method of prime number sums, zero density estimates,
primes in arithmetic progressions on average, sums of primes, sieve methods,
the distribution of additive functions and mean values of multiplicative func-
tions, and the least prime in an arithmetic progression. The present volume was
xi
xii Preface
xiii
xiv List of notation
Symbol Meaning
∞ Found on page
(s, a) = a e−w w s−1 dw; the incomplete 327
Gamma function.
γ The imaginary part of a zero of the 172
zeta function or of an L-function.
N −1
N (θ ) = 1 + 2 n=1 (1 − n/N ) cos 2π nθ; 174
known as the Fejér
kernel.
κ 1/2
ε(χ ) = τ (χ )/ i q . 332
ζ (s) = ∞ n=1 n −s
for σ > 1, known as the 2
Riemann zeta function.
ζ (s, α) = ∞ n=0 (n + α)
−s
for σ > 1; known 30
as the Hurwitz zeta function.
−s
ζ K (s) a N (a) ; known as the Dedekind 343
zeta function of the algebraic number
field K .
= sup ρ 430, 463
ϑ(x) = p≤x log p. 46
∞
= n=−∞ e−π n z for z > 0.
2
ϑ(z) 329
ϑ(x; q, a) The sum of log p over primes p ≤ x 128, 377ff
for which p ≡ a (mod q).
ϑ(x, χ) = p≤x χ ( p) log p. 377ff
κ = (1 − χ (−1))/2. 332
(n) = log p if n = p k , = 0 otherwise; 23
known as the von Mangoldt Lambda
function.
2 (n) = (n) log n + bc=n (b)(c). 251
(x; q, a) The sum of λ(n) over those n ≤ x 383
such that n ≡ a (mod q).
(x, χ ) = n≤x χ (n)λ(n). 383
λ(n) = (−1)(n) ; known as the Liouville 21
lambda function.
µ(n) = (−1)ω(n) for square-free n, = 0 21
otherwise. Known as the Möbius mu
function.
µ(σ ) the Lindelöf mu function 330
ξ (s) = 12 s(s − 1)ζ (s) (s/2)π −s/2 . 328
ξ (s, χ ) = L(s, χ ) ((s + κ)/2)(q/π )(s+κ)/2 333
where χ is a primitive character
modulo q, q > 1.
xvi List of notation
for |z| < 1, then the n th power series coefficient of f (z)s is the number rk,s (n)
of representations of n as a sum of s positive k th powers,
n = m k1 + m k2 + · · · + m ks .
We can recover rk,s (n) from f (z)s by means of Cauchy’s coefficient formula:
1 f (z)s
rk,s (n) = dz.
2πi z n+1
By choosing an appropriate contour, and estimating the integrand, we can de-
termine the asymptotic size of rk,s (n) as n → ∞, provided that s is sufficiently
large, say s > s0 (k). This is the germ of the Hardy–Littlewood circle method,
but considerable effort is required to construct the required estimates.
To appreciate why power series are useful in dealing with additive prob-
lems, note that if A(z) = ak z k and B(z) = bm z m then the power series
1
2 Dirichlet series: I
The terms are grouped according to the sum of the indices, because
z k z m = z k+m .
A Dirichlet series is a series of the form α(s) = ∞ n=1 an n
−s
where s is
∞ −s
a complex variable. If β(s) = m=1 bm m is a second Dirichlet series and
γ (s) = α(s)β(s), then (ignoring questions relating to the rearrangement of terms
of infinite series)
∞
∞ ∞ ∞ ∞
γ (s) = ak k −s bm m −s = ak bm (km)−s = ak bm n −s .
k=1 m=1 k=1 m=1 n=1 km=n
(1.2)
∞ −s
That is, we expect that γ (s) is a Dirichlet series, γ (s) = n=1 cn n , whose
coefficients are
cn = a k bm . (1.3)
km=n
This corresponds to (1.1), but the terms are now grouped according to the
product of the indices, since k −s m −s = (km)−s .
Since we shall employ the complex variable s extensively, it is useful to have
names for its real and complex parts. In this regard we follow the rather peculiar
notation that has become traditional: s = σ + it.
Among the Dirichlet series we shall consider is the Riemann zeta function,
which for σ > 1 is defined by the absolutely convergent series
∞
ζ (s) = n −s . (1.4)
n=1
As a first application of (1.3), we note that if α(s) = β(s) = ζ (s) then the
manipulations in (1.3) are justified by absolute convergence, and hence we see
that
∞
d(n)n −s = ζ (s)2 (1.5)
n=1
for σ > 1. Here d(n) is the divisor function, d(n) = d|n 1.
From the rate of growth or analytic behaviour of generating functions we
glean information concerning the sequence of coefficients. In expressing our
findings we employ a special system of notation. For example, we say, ‘f (x) is
asymptotic to g(x)’ as x tends to some limiting value (say x → ∞), and write
1.1 Generating functions and asymptotics 3
0000
0000
0000
0000
Figure 1.1 Graph of π(x) (solid) and x/ log x (dotted) for 2 ≤ x ≤ 106 .
proved that π (x) x/ log x. This is of course weaker than the Prime Number
Theorem, but it was derived much earlier, in 1852. Chebyshev also showed
that π(x) x/ log x. In general, we say that f (x) g(x) if there is a positive
constant c such that f (x) ≥ cg(x) and g is non-negative. In this situation both
f and g take only positive values. If both f g and f g then we say that f
and g have the same order of magnitude, and write f g. Thus Chebyshev’s
estimates can be expressed as a single relation,
x
π (x) .
log x
The estimate (1.6) is best possible to the extent that the error term is not
o(x(log x)−2 ). We have also a special notation to express this:
x x
π(x) − = .
log x (log x)2
In general, if lim supx→∞ | f (x)|/g(x) > 0 then we say that f (x) is ‘Omega of
g(x)’, and write f (x) = (g(x)). This is precisely the negation of the statement
‘ f (x) = o(g(x))’. When studying numerical values, as in Figure 1.1, we find
that the fit of x/ log x to π (x) is not very compelling. This is because the error
term in the approximation is only one logarithm smaller than the main term.
This error term is not oscillatory – rather there is a second main term of this
1.1 Generating functions and asymptotics 5
size:
x x x
π (x) = + + O .
log x (log x)2 (log x)3
This is also best possible, but the main term can be made still more elaborate to
give a smaller error term. Gauss was the first to propose a better approximation to
π(x). Numerical studies led him to observe that the density of prime numbers in
the neighbourhood of x is approximately 1/ log x. This suggests that the number
of primes not exceeding x might be approximately equal to the logarithmic
integral,
x
1
li(x) = du.
2 log u
(Orally, ‘li’ rhymes with ‘pi’.) By repeated integration by parts we can show
that
K −1
(k − 1)! x
li(x) = x + OK
k=1
(log x)k (log x) K
for any positive integer K ; thus the secondary main terms of the approximation
to π(x) are contained in li(x).
In Chapter 6 we shall prove the Prime Number Theorem in the sharper
quantitative form
x
π(x) = li(x) + O √
exp(c log x)
√
for some suitable positive constant c. Note that exp(c log x) tends to infinity
faster than any power of log x. The error term above seems to fall far from
what seems to be the truth. Numerical evidence, such as that in Table 1.1,
√
suggests that the error term in the Prime Number Theorem is closer to x in
size. Gauss noted the good fit, and also that π(x) < li(x) for all x in the range of
his extensive computations. He proposed that this might continue indefinitely,
but the numerical evidence is misleading, for in 1914 Littlewood showed that
1/2
x log log log x
π(x) − li(x) = ± .
log x
Here the subscript ± indicates that the error term achieves the stated or-
der of magnitude infinitely often, and in both signs. In particular, the dif-
ference π − li has infinitely many sign changes. More generally, we write
f (x) = + (g(x)) if lim supx→∞ f (x)/g(x) > 0, we write f (x) = − (g(x))
if lim infx→∞ f (x)/g(x) < 0, and we write f (x) = ± (g(x)) if both these re-
lations hold.
6 Dirichlet series: I
1.1.1 Exercises
1. Let r (n) be the number of ways that n cents of postage can be made, using
only 1 cent, 2 cent, and 3 cent stamps. That is, r (n) is the number of ordered
triples (x1 , x2 , x3 ) of non-negative integers such that x1 + 2x2 + 3x3 = n.
(a) Show that
∞
1
r (n)z n =
n=0
(1 − z)(1 − z 2 )(1 − z 3 )
a b c d e f
+ + + + +
(z − 1) 3 (z − 1) 2 z−1 z+1 z−ω z−ω
where ω = e2πi/3 and ω = e−2πi/3 are the primitive cube roots of unity.
(c) Show that r (n) is the integer nearest (n + 3)2 /12.
(d) Show that r (n) is the number of ways of writing n = y1 + y2 + y3 with
y1 ≥ y2 ≥ y3 ≥ 0.
2. Explain why
∞
k
1 + z2 = 1 + z + z2 + · · ·
k=0
K
z ak 1
=
k=1
1 − zmk 1−z
e2πia K /m K
∼
m K (1 − r )
for all n ≥ 2.
(b) Let P(z) = ∞ n
n=0 A(n)z . Show that
P(z)2 = P(z) − z.
Deduce that
√
1− 1 − 4z ∞
1/2 2n−1
P(z) = = 2 (−1)n−1 z n .
2 n=1
n
2n−2
(c) Conclude that A(n) = n−1 /n for all n ≥ 1. These are called the Cata-
lan numbers.
1.1 Generating functions and asymptotics 9
(d) What needs to be said concerning the convergence of the series used
above?
6. (a) Let n k denote the total number of monic polynomials of degree k in
F p [x]. Show that n k = p k .
(b) Let P1 , P2 , . . . be the irreducible monic polynomials in F p [x], listed in
some (arbitrary) order. Show that
∞
(1 + z deg Pr + z 2 deg Pr + z 3 deg Pr + · · · ) = 1 + pz + p 2 z 2
r =1
+ p3 z 3 + · · ·
(j) Show that gn > 0 for all p and all n ≥ 1. (If P ∈ F p [x] is irreducible and
has degree n, then the quotient ring F p [x]/(P) is a field of p n elements.
Thus we have proved that there is such a field, for each prime p and
integer n ≥ 1. It may be further shown that the order of a finite field
is necessarily a prime power, and that any two finite fields of the same
order are isomorphic. Hence the field of order p n , whose existence we
have proved, is essentially unique.)
7. (E. Berlekamp) Let p be a prime number. We recall that polynomials in a
single variable (mod p) factor uniquely into irreducible polynomials. Thus
a monic polynomial f (x) can be expressed uniquely (mod p) in the form
g(x)h(x)2 where g(x) is square-free (mod p) and both g and h are monic. Let
sn denote the number of monic square-free polynomials (mod p) of degree
n. Show that
∞
∞ ∞
sk z k p m z 2m = pn z n
k=0 m=0 n=0
∞
1 − pz 2
sk z k = ,
k=0
1 − pz
induction that
∞
k
A(z) − B(z) = 1 − z2
k=0
since R(u) is constant in the interval [n − 1, n). The integrals combine to give
(1.7).
If |R(u)| ≤ ε for all u ≥ M and if σ > σ0 , then from (1.7) we see that
N ∞
|s − s0 |
an n −s ≤ 2ε + ε|s − s0 | u σ0 −σ −1 du ≤ 2 + ε.
n=M+1 M σ − σ0
For s in the prescribed region we see that
|s − s0 | ≤ σ − σ0 + |t − t0 | ≤ (H + 1)(σ − σ0 ),
N
so that the sum M+1 an n −s is uniformly small, and the result follows by the
uniform version of Cauchy’s principle.
Weierstrass that α(s) is analytic for σ > σc , and that the differentiated series is
locally uniformly convergent to α (s):
∞
α (s) = − an (log n)n −s (1.8)
n=1
need not exist. However, α(s) is continuous in the sector S of Theorem 1.1, in
view of the uniform convergence there. That is,
lim α(s) = α(s0 ),
s→s0
(1.9)
s∈S
Let φ denote the left-hand side of (1.11). If θ > φ then A(x) x θ where the
implicit constant may depend on the an and on θ . Thus if σ > θ, then the integral
in (1.10) is absolutely convergent. Thus we obtain (1.10) by letting N → ∞,
since the first term above tends to 0 as N → ∞.
Suppose that σc < 0. By Corollary 1.2 we know that A(x) tends to a finite
limit as x → ∞, and hence φ ≤ 0, so that (1.10) holds for all σ > 0.
14 Dirichlet series: I
Now suppose that σc ≥ 0. By Corollary 1.2 we know that the series in (1.10)
diverges when σ < σc . Hence φ ≥ σc . To complete the proof it suffices to show
that φ ≤ σc . Choose σ0 > σc . By (1.7) with s = 0 and M = 0 we see that
N
A(N ) = −R(N )N σ0 + σ0 R(u)u σ0 −1 du.
0
Since R(u) is a bounded function, it follows that A(N ) N σ0 where the implicit
constant may depend on the an and on σ0 . Hence φ ≤ σ0 . Since this holds for
any σ0 > σc , we conclude that φ ≤ σc .
converges for σ > 0, but it is absolutely convergent only for σ > 1. In general
we let σa denote the infimum of those σ for which ∞ n=1 |an |n
−σ
< ∞. Then σa ,
the abscissa of absolute convergence, is the abscissa of convergence of the series
∞ −s
n=1 |an |n , and we see that an n −s is absolutely convergent if σ > σa ,
but not if σ < σa . We now show that the strip σc ≤ σ ≤ σa of conditional
convergence is never wider than in the example (1.12).
Proof The first inequality is obvious. To prove the second, suppose that ε > 0.
Since the series an n −σc −ε is convergent, the summands tend to 0, and hence
an n σc +ε where the implicit constant may depend on the an and on ε. Hence
the series an n −σc −1−2ε is absolutely convergent by comparison with the series
−1−ε
n .
Theorem 1.5 Suppose that α(s) = an n −s has abscissa of convergence σc .
If δ and ε are fixed, 0 < ε < δ < 1, then
α(s) τ 1−δ+ε
By the example found in Exercise 8 at the end of this section, we see that
the bound above is reasonably sharp.
Since the series α(σc + ε) converges, we know that an n σc +ε , and also that
R(u) 1. Thus the above is
M
|σc + ε − s| σc +ε−σ
n −δ+ε + M −δ+ε + M .
n=1
σ − σc − ε
By Theorem 1.4 this sum is absolutely convergent for σ > σ0 + 1. Since each
term tends to 0 as σ → ∞, we see that the right-hand side tends to 0, by
the principle of dominated convergence. Hence c N = 0, and by induction we
deduce that this holds for all N .
16 Dirichlet series: I
for σ > 1. That this is an entire function follows from Theorem 10.2.) Since a
Dirichlet series does not in general have a singularity on its line of convergence,
it is noteworthy that a Dirichlet series with non-negative coefficients not only
has a singularity on the line σc + it, but actually at the point σc .
Theorem 1.7 (Landau) Let α(s) = an n −s be a Dirichlet series whose ab-
scissa of convergence σc is finite. If an ≥ 0 for all n then the point σc is a
singularity of the function α(s).
It is enough to assume that an ≥ 0 for all sufficiently large n, since any finite
N
sum n=1 an n −s is an entire function.
Proof By replacing an by an n −σc , we may assume that σc = 0. Suppose that
α(s) is analytic at s = 0, so that α(s) is analytic in the domain D = {s : σ >
0} ∪ {|s| < δ} if δ > 0 is sufficiently small. We expand α(s) as a power series
at s = 1:
∞
α(s) = ck (s − 1)k . (1.13)
k=0
for |s − 1| < 1 + δ . If s < 1 then all terms above are non-negative. Since
series of non-negative numbers may be arbitrarily rearranged, for −δ < s < 1
we may interchange the summations over k and n to see that
∞
∞
(1 − s)k (log n)k
α(s) = an n −1
n=1 k=0
k!
∞
∞
= an n −1 exp (1 − s) log n = an n −s .
n=1 n=1
Hence this last series converges at s = −δ /2, contrary to the assumption that
σc = 0. Thus α(s) is not analytic at s = 0.
1.2.1 Exercises
1. Suppose that α(s) is a Dirichlet series, and that the series α(s0 ) is boundedly
oscillating. Show that σc = σ0 .
2. Suppose that α(s) = ∞ n=1 an n
−s
is a Dirichlet series with abscissa of con-
vergence σc . Suppose that α(0) converges, and put R(x) = n>x an . Show
that σc is the infimum of those numbers θ such that R(x) xθ .
3. Let Ak (x) = n≤x an (log n) .k
marks following the proof of Theorem 1.1 imply only that σc ≤ σc .)
4. (Landau 1909b) Let α(s) = an n −s be a Dirichlet series with abscissa of
convergence σc and abscissa of absolute convergence σa > σc . Let C(x) =
−σc
n≤x an n and A(x) = n≤x |an |n −σc .
(a) By a suitable application of Theorem 1.3, or otherwise, show that
C(x) x ε and that A(x) x σa −σc +ε for any ε > 0, where the implicit
constants may depend on ε and on the sequence {an }.
(b) Show that if σ > σc then
∞
an n −s = −C(N )N σc −s + (s − σc ) C(u)u σc −s−1 du.
n>N N
18 Dirichlet series: I
α(s) τ θ (σ )+ε
lim α(σ ) = a1 .
σ →∞
(b) Show that ζ (s) = − ∞ n=1 (log n)n
−s
for σ > 1.
(c) Show that limσ →∞ ζ (σ ) = 0.
(d) Show that there is no half-plane in which 1/ζ (s) can be written as a
convergent Dirichlet series.
6. Let α(s) = an n −s be a Dirichlet series with an ≥ 0 for all n. Show that
σc = σa , and that
2tr
(a) Show that tr an = 0.
1.3 Euler products and the zeta function 19
(b) Show that if tr ≤ x < 2tr for some r , then A(x) = [x]itr where A(x) =
n≤x an .
(c) Show that A(x) 1 uniformly for x ≥ 1.
(d) Deduce that α(s) converges for σ > 0.
(e) Show that α(it) does not converge; conclude that σc = 0.
(f) Show that if σ > 0, then
R
2tr ∞
α(s) = an n −s + s A(x)x −s−1 d x .
r =1 n=tr t R+1
(i) Show that if n ≤ x < n + 1, then (n it R x −it R ) ≥ 1/2. Deduce that
2t R
[x]it R x −σ −it R −1 d x t R−σ .
tR
(j) Suppose that δ > 0 is fixed. Conclude that if R ≥ R0 (δ), then |α(σ +
it R )| t R1−σ uniformly for δ ≤ σ ≤ 1 − δ.
(k) Show that |an |n −σ < ∞ when σ > 1. Deduce that σa = 1.
The mere convergence of α(s) and β(s) is not sufficient to justify (1.2).
Indeed, the square of the series (1.12) can be shown to have abscissa of conver-
gence ≥ 1/4.
20 Dirichlet series: I
(1 + f ( p) p −s + f ( p 2 ) p −2s + f ( p 3 ) p −3s + · · · )
p
Set n = p1k1 p2k2 · · · prkr . Since f is multiplicative, the above is f (n)n −s . More-
over, this correspondence between products of prime powers and positive inte-
gers n is one-to-one, in view of the fundamental theorem of arithmetic. Hence
after rearranging the terms, we obtain the sum f (n)n −s . That is, we expect
that
∞
f (n)n −s = (1 + f ( p) p −s + f ( p 2 ) p −2s + · · · ). (1.14)
n=1 p
The product on the right-hand side is called the Euler product of the Dirichlet
series. The mere convergence of the series on the left does not imply that the
product converges; as in the case of the identity (1.2), we justify (1.14) only
under the stronger assumption of absolute convergence.
Theorem 1.9 If f is multiplicative and | f (n)|n −σ < ∞, then (1.14) holds.
Hence
∞
y − f (n)n −s ≤ | f (n)|n −σ .
n=1 n ∈N
/
Let ω(n) denote the number of distinct primes dividing n, and let (n) be
the number of distinct prime powers dividing n. That is,
ω(n) = 1, (n) = 1= k. (1.16)
p|n p k |n p k n
It is easy to distinguish these functions, since ω(n) ≤ (n) for all n, with equal-
ity if and only if n is square-free. These functions are examples of additive
functions because they satisfy the functional relation f (mn) = f (m) + f (n)
whenever (m, n) = 1. Moreover, (n) is totally additive because this func-
tional relation holds for all pairs m, n. An exponential of an additive function is
a multiplicative function. In particular, the Liouville lambda function is the to-
tally multiplicative function λ(n) = (−1)(n) . Closely related is the Möbius mu
function, which is defined to be µ(n) = (−1)ω(n) if n is square-free, µ(n) = 0
otherwise. By the fundamental theorem of arithmetic we know that a multi-
plicative (or additive) function is uniquely determined by its values at prime
powers, and similarly that a totally multiplicative (or totally additive) function
is uniquely determined by its values at the primes. Thus µ(n) is the unique
multiplicative function that takes the value −1 at every prime, and the value 0
at every higher power of a prime, while λ(n) is the unique totally multiplicative
function that takes the value −1 at every prime. By using Theorem 1.9 we can
22 Dirichlet series: I
determine the Dirichlet series generating functions of λ(n) and of µ(n) in terms
of the Riemann zeta function.
Corollary 1.10 For σ > 1,
∞
n −s = ζ (s) = (1 − p −s )−1 , (1.17)
n=1 p
∞
1
µ(n)n −s = = (1 − p −s ), (1.18)
n=1
ζ (s) p
and
∞
ζ (2s)
λ(n)n −s = = (1 + p −s )−1 . (1.19)
n=1
ζ (s) p
Proof All three series are absolutely convergent, since n −σ < ∞ for σ >
1, by the integral test. Since the coefficients are multiplicative, the Euler product
formulae follow by Theorem 1.9. In the first and third cases use the variant
(1.15). On comparing the Euler products in (1.17) and (1.18), it is immediate
that the second of these Dirichlet series is 1/ζ (s). As for (1.19), from the identity
1 + z = (1 − z 2 )/(1 − z) we deduce that
−2s
−s p (1 − p ) ζ (s)
(1 + p ) = −s
= .
p p (1 − p ) ζ (2s)
∞
log ζ (s) = log(1 − p −s )−1 = k −1 p −ks .
p p k=1
for σ > 1. This is a Dirichlet series, whose n th coefficient is the von Mangoldt
lambda function: (n) = log p if n is a power of p, (n) = 0 otherwise.
and
ζ (s) ∞
− = (n)n −s .
ζ (s) n=1
for σ > 1.
We now continue the zeta function beyond the half-plane in which it was
initially defined.
Here {u} denotes the fractional part of u, so that {u} = u − [u] where [u]
denotes the integral part of u.
We evaluate the first integral on the right-hand side, and integrate the second
one by parts. Thus the above is
∞
x 1−s
= + {x}x −s + {u} du −s .
s−1 x
Since (u −s ) = −su −s−1 , the desired formula now follows by Theorem A.3.
The integral in (1.23) is convergent in the half-plane σ > 0, and uniformly so
for σ ≥ δ > 0. Since the integrand is an analytic function of s, it follows that the
integral is itself an analytic function for σ > 0. By the uniqueness of analytic
continuation the formula (1.23) holds in this larger half-plane.
1.3 Euler products and the zeta function 25
10
0 1 5
–10
In addition,
1
= log x + C0 + O(1/x) (1.26)
n≤x n
Proof The first estimate follows by crudely estimating the integral in (1.23):
∞ ∞
x −σ
{u}u −s−1 du u −σ −1 du = .
x x σ
As for the second estimate, we note that the sum is
x x x
u −1 d[u] = u −1 du − u −1 d{u}
1− 1− 1−
x
= log x + 1 − {x}/x − {u}u −2 du.
1
x ∞ ∞
The result now follows by writing 1 = 1 − x , and noting that
∞ ∞
{u}u −2 du u −2 du = 1/x.
x x
Proof The first assertion is clear from (1.24). When |t| is larger, we obtain
a bound for |ζ (s)| by estimating the sum in (1.25). Assume that x ≥ 2. We
observe that
x
n −s n −σ 1+ u −σ du
n≤x n≤x 1
1.3.1 Exercises
1. Suppose that f (mn) = f (m) f (n) whenever (m, n) = 1, and that f is not
identically 0. Deduce that f (1) = 1, and hence that f is multiplicative.
2. (Stieltjes 1887) Suppose that an converges, that |b | < ∞, and that
n
cn is given by (1.3). Show that cn converges to ( an )( bn ). (Hint:
Write n≤x cn = n≤x bn A(x/n) where A(y) = n≤y an .)
3. Determine ϕ(n)n −s , σ (n)n −s , and |µ(n)|n −s in terms of the zeta
function. Here ϕ(n) is Euler’s ‘totient function’, which is the number of a,
1 ≤ a ≤ n, such that (a, n) = 1.
4. Let q be a positive integer. Show that if σ > 1, then
∞
n −s = ζ (s) (1 − p −s ).
n=1 p|q
(n,q)=1
6. Let σa (n) = d|n d a . Show that
∞
σa (n)σb (n)n −s = ζ (s)ζ (s − a)ζ (s − b)ζ (s − a − b)/ζ (2s − a − b)
n=1
for σ > σc .
(b) Show that
∞
(−1)n−1 n −s = (1 − 21−s )ζ (s)
n=1
for σ > 0.
1.3 Euler products and the zeta function 29
Thus if
Nk
εk log 1 (1.31)
Nk−1
then the series is not uniformly convergent in D.
(b) By using Corollary 1.15, or otherwise, show that if (a, b] ⊆ (Nk−1 , Nk ],
then
εk
an n −1 .
a<n≤b
tk
Hence if
∞
εk
< ∞, (1.32)
k=1
tk
then the series α(1) converges.
30 Dirichlet series: I
(c) Show that the parameters can be chosen so that (1.30)–(1.32) hold, say
1/2
by taking Nk = exp(1/εk ) and tk = εk with εk tending rapidly to 0.
14. Let t(n) = (−1)(n)−ω(n) p|n ( p − 1)−1 , and put T (s) = n t(n)n −s .
(a) Show that for σ > 0, T (s) has the absolutely convergent Euler product
1
T (s) = 1+ .
p ( p − 1)( p s + 1)
(c) Deduce that ζ (s, α) is an analytic function of s for σ > 0 apart from a
simple pole at s = 1 with residue 1.
(d) Show that
∞
1 {u}
lim ζ (s, α) − = 1/α − log α − du.
s→1 s−1 0 (u + α)2
(e) Show that
1 1 {x}
lim ζ (s, α) − = − log(x + α) +
s→1 s−1 0≤n≤x
n + α x +α
∞
{u}
− du.
x (u + α)2
(f) Let x → ∞ in the above, and use (C.2), (C.10) to show that
1
lim ζ (s, α) − =− (α).
s→1 s−1
(This is consistent with Corollary 1.16, in view of (C.11).)
1.4 Notes 31
1.4 Notes
Section 1.1. For a brief introduction to the Hardy–Littlewood circle method,
including its application to Waring’s problem, see Davenport (2005). For a
comprehensive account of the method, see Vaughan (1997). Other examples
of the fruitful use of generating functions are found in many sources, such as
Andrews (1976) and Wilf (1994).
Algorithms for the efficient computation of π(x) have been developed
by Meissel (Lehmer, 1959), Mapes (1963), Lagarias, Miller & Odlyzko
(1985), Deléglise & Rivat (1996), and by X. Gourdon. For discussion
of these methods, see Chapter 1 of Riesel (1994) and the web page of
Gourdon & Sebah at http://numbers.computation.free.fr/Constants/Primes/
countingPrimes.html.
The ‘big oh’ notation was introduced by Paul Bachmann (1894, p. 401). The
‘little oh’ was introduced by Edmund Landau (1909a, p. 61). The notation
was introduced by Hardy (1910, p. 2). Our notation f ∼ g also follows Hardy
(1910). The Omega notation was introduced by G. H. Hardy and J. E. Littlewood
(1914, p. 225). Ingham (1932) replaced the R and L of Hardy and Littlewood
by + and − . The notation is due to I. M. Vinogradov.
Section 1.2. The series an n −s is called an ordinary Dirichlet series,
to distinguish it from a generalized Dirichlet series, which is a sum of the
form an e−λn s where 0 < λ1 < λ2 < · · · , λn → ∞. We see that generalized
Dirichlet series include both ordinary Dirichlet series (λn = log n) and power
series (λn = n). Theorems 1.1, 1.3, 1.6, and 1.7 extend naturallyto generalized
∞
Dirichlet series, and even to the more general class of functions 0 e−us d A(u)
where A(u) is assumed to have finite variation on each finite interval [0, U ].
The proof of the general form of Theorem 1.6 must be modified to depend on
uniform, rather than absolute, convergence, since a generalized Dirichlet series
may be never more than conditionally convergent (e.g., (−1)n (log n)−s ).
If we put a = lim sup(log n)/λn , then the general form of Theorem 1.4
reads σc ≤ σa ≤ σc + a. Hardy & Riesz (1915) have given a detailed ac-
count of this subject, with historical attributions. See also Bohr & Cramér
(1923).
Jensen (1884) showed that the domain of convergence of a generalized
Dirichlet series is always a half-plane. The more precise information provided
by Theorem 1.1 is due to Cahen (1894) who proved it not only for ordinary
Dirichlet series but also for generalized Dirichlet series.
The construction in Exercise 1.2.8 would succeed with the simpler choice
an = n itr for tr ≤ n ≤ 2tr , an = 0 otherwise, but then to complete the argu-
ment one would need a further tool, such as the Kusmin–Landau inequality
32 Dirichlet series: I
(cf. Mordell 1958). The square of the Dirichlet series in Exercise 1.2.8 has ab-
scissa of convergence 1/2; this bears on the result of Exercise 2.1.9. Information
concerning the convergence of the product of two Dirichlet series is found in
Exercises 1.3.2, 2.1.9, 5.2.16, and in Hardy & Riesz (1915).
Theorem 1.7 originates in Landau (1905). The analogue for power series had
been proved earlier by Vivanti (1893) and Pringsheim (1894). Landau’s proof
extends to generalized Dirichlet series (including power series).
Section 1.3. The hypothesis | f (n)|n −σ < ∞ of Theorem 1.9 is equivalent
to the assertion that
which is slightly stronger than merely asserting that the Euler product converges
absolutely. We recall that a product n (1 + an ) is said to be absolutely con-
vergent if n (1 + |an |) < ∞. To see that the hypothesis p (1 + | f ( p) p −s +
· · · |) < ∞ is not sufficient, consider the following example due to Ingham:
For every prime p we take f ( p) = 1, f ( p 2 ) = −1, and f ( p k ) = 0 for k > 2.
Then the product is absolutely convergent at s = 0, but the terms f (n) do not
tend to 0, and hence the series f (n) diverges. Indeed, it can be shown that
−2 −3
n≤x f (n) ∼ cx as x → ∞ where c = p 1 − 2 p + p > 0.
Euler (1735) defined the constant C0 , which he denoted C.
Mascheroni (1790) called the constant γ , which is in common use, but
we wish to reserve this symbol for the imaginary part of a zero of the
zeta function or an L-function. It is conjectured that Euler’s constant C0
is irrational. The early history of the determination of the initial digits of
C0 has been recounted by Nielsen (1906, pp. 8–9). More recently, Wrench
(1952) computed 328 digits, Knuth (1963) computed 1,271 digits, Sweeney
(1963) computed 3,566 digits, Beyer & Waterman (1974) computed 4,879
digits, Brent (1977) computed 20,700 digits, Brent & McMillan (1980)
computed 30,100 digits. At this time, it seems that more than 108 digits
have been computed – see the web page of X. Gourdon & P. Sebah at
http://numbers.computation.free.fr/Constants/Gamma/gamma.html. To 50
places, Euler’s constant is
C0 = 0.57721 56649 01532 86060 65120 90082 40243 10421 59335 93992.
work on Dirichlet series with natural boundaries see Estermann (1928a,b) and
Kurokawa (1987).
1.5 References
Andrews, G. E. (1976). The Theory of Partitions, Reprint. Cambridge: Cambridge Uni-
versity Press (1998).
Bachmann, P. (1894). Zahlentheorie, II, Die analytische Zahlentheorie, Leipzig:
Teubner.
Beyer, W. A. & Waterman, M. S. (1974). Error analysis of a computation of Euler’s
constant and ln 2, Math. Comp. 28, 599–604.
Bohr, H. (1910). Bidrag til de Dirichlet’ske Rækkers theori, København: G. E. C. Gad;
Collected Mathematical Works, Vol. I, København: Danske Mat. Forening, 1952.
A3.
Bohr, H. & Cramér, H. (1923). Die neuere Entwicklung der analytischen Zahlentheo-
rie, Enzyklopädie der Mathematischen Wissenschaften, 2, C8, 722–849; H. Bohr,
Collected Mathematical Works, Vol. III, København: Dansk Mat. Forening, 1952,
H; H. Cramér, Collected Works, Vol. 1, Berlin: Springer-Verlag, 1952, pp. 289–
416.
Brent, R. P. (1977). Computation of the regular continued fraction of Euler’s constant,
Math. Comp. 31, 771–777.
Brent, R. P. & McMillan, E. M. (1980). Some new algorithms for high-speed computation
of Euler’s constant, Math. Comp. 34, 305–312.
Cahen, E. (1894). Sur la fonction ζ (s) de Riemann et sur des fonctions analogues, Ann.
de l’École Normale (3) 11, 75–164.
Davenport, H. (2005). Analytic Methods for Diophantine Equations and Diophantine
Inequalities. Second edition, Cambridge: Cambridge University Press.
Deléglise, M. & Rivat, J. (1996). Computing π (x): the Meissel, Lehmer, Lagarias, Miller,
Odlyzko method, Math. Comp. 65, 235–245.
Estermann, T. (1928a). On certain functions represented by Dirichlet series, Proc. Lon-
don Math. Soc. (2) 27, 435–448.
(1928b). On a problem of analytic continuation, Proc. London Math. Soc. (2) 27,
471–482.
Euler, L. (1735). De Progressionibus harmonicus observationes, Comm. Acad. Sci. Imper.
Petropol. 7, 157; Opera Omnia, ser. 1, vol. 14, Teubner, 1914, pp. 93–95.
Hardy, G. H. (1910). Orders of Infinity. Cambridge Tract 12, Cambridge: Cambridge
University Press.
Hardy, G. H. & Littlewood, J. E. (1914). Some problems of Diophantine approximation
(II), Acta Math. 37, 193–238; Collected Papers, Vol I. Oxford: Oxford University
Press. 1966, pp. 67–112.
Hardy, G. H. & Riesz, M. (1915). The General Theory of Dirichlet’s Series, Cambridge
Tract No. 18. Cambridge: Cambridge University Press. Reprint: Stechert–Hafner
(1964).
Ingham, A. E. (1932). The Distribution of Prime Numbers, Cambridge Tract 30. Cam-
bridge: Cambridge University Press.
34 Dirichlet series: I
1 N
lim F(n) = c.
N →∞ N n=1
In this section we develop a simple method by which mean values can be shown
to exist in many interesting cases.
If two arithmetic functions f and F are related by the identity
F(n) = f (d), (2.1)
d|n
This is the Möbius inversion formula. Conversely, if (2.2) holds for all n then
so also does (2.1). If f is generally small then F has an asymptotic mean value.
To see this, observe that
F(n) = f (d).
n≤x n≤x d|n
By iterating the sums in the reverse order, we see that the above is
= f (d) 1= f (d)[x/d].
d≤x n≤x d≤x
d|n
35
36 The elementary theory of arithmetic functions
by Corollary 1.10. From Corollary B.3 we know that ζ (2) = π 2 /6; hence the
proof is complete.
Let Q(x) denote the number of square-free integers not exceeding x, Q(x) =
n≤x µ(n) . We now calculate the asymptotic density of these numbers.
2
√
This is a relation of the shape (2.1) where f (d) = µ( d) if d is a perfect square,
and f (d) = 0 otherwise. Hence by (2.3),
µ(d)
Q(x) = x +O 1 .
d 2 ≤x
d2 d 2 ≤x
The error term is x 1/2 , and the sum in the main term is treated as in the
preceding proof.
We note that the argument above is routine once the appropriate identity
(2.4) is established. This relation can be discovered by considering (2.2), or by
using Dirichlet series: Let Q denote the class of square-free numbers. Then for
σ > 1,
1 − p −2s ζ (s)
n −s = (1 + p −s ) = −s
= .
n∈Q p p 1− p ζ (2s)
Now 1/ζ (2s) can be written as a Dirichlet series in s, with coefficients f (n) =
µ(d) if n = d 2 , f (n) = 0 otherwise. Hence the convolution equation (2.4) gives
the coefficients of the product Dirichlet series ζ (s) · 1/ζ (2s).
Suppose that ak , bm , cn are joined by the convolution relation
cn = ak bm , (2.5)
km=n
and that A(x), B(x), C(x) are their respective summatory functions. Then
C(x) = ak bm , (2.6)
km≤x
and it is useful to note that this double sum can be iterated in various ways. On
one hand we see that
C(x) = ak B(x/k); (2.7)
k≤x
this is the line of reasoning that led to (2.3) (take ak = f (k), bm = 1). At the
opposite extreme,
C(x) = bm A(x/m), (2.8)
m≤x
for 0 < y ≤ x. This is obvious once it is observed that the first term on the right
sums those terms ak bm for which km ≤ x, k ≤ y, the second sum includes the
38 The elementary theory of arithmetic functions
pairs (k, m) for which km ≤ x, m ≤ x/y, and the third term subtracts those ak bm
for which k ≤ y, m ≤ x/y, since these (k, m) were included in both the previous
terms. The advantage of (2.9) over (2.7) is that the number of terms is reduced
( y + x/y instead of x), and at the same time A and B are evaluated only
at large values of the argument, so that asymptotic formulæ for these quantities
may be expected to be more accurate. For example, if we wish to estimate the
average size of d(n) we take ak = bm = 1, and then from (2.3) we see that
d(n) = x log x + O(x).
n≤x
To obtain a more accurate estimate we observe that the first term on the
right-hand side of (2.9) is
[x/k] = x 1/k + O(y).
k≤y k≤y
Here the error term is minimized by taking y = x 1/2 . The second term
on the right in (2.9) is then identical to the first, and the third term is
[x 1/2 ]2 = x + O(x 1/2 ), and we have
We often construct estimates with one or more parameters, and then choose
values of the parameters to optimize the result. The instance above is typical –
we minimized x/y + y by taking y = x 1/2 . Suppose, more generally, that we
wish to minimize T1 (y) + T2 (y) where T1 is a decreasing function, and T2 is
an increasing function. We could differentiate and solve for a root of T1 (y) +
T2 (y) = 0, but there is a quicker method: Find y0 so that T1 (y0 ) = T2 (y0 ). This
does not necessarily yield the exact minimum value of T1 (y) + T2 (y), but it is
easy to see that
so the bound obtained in this way is at most twice the optimal bound.
Despite the great power of analytic techniques, the ‘method of the hyperbola’
used above is a valuable tool. The sequence cn given by (2.5) is called the
Dirichlet convolution of ak and bm ; in symbols, c = a ∗ b. Arithmetic functions
form a ring when equipped with pointwise addition, (a + b)n = an + bn , and
2.1 Mean values 39
Dirichlet convolution for multiplication. This ring is called the ring of formal
Dirichlet series. Manipulations of arithmetic functions in this way correspond
to manipulations of Dirichlet series without regard to convergence. This is
analogous to the ring of formal power series, in which multiplication is provided
by Cauchy convolution, cn = k+m=n ak bm .
In the ring of formal Dirichlet series we let O denote the arithmetic function
that is identically 0; this is the additive identity. The multiplicative identity is i
where i 1 = 1, i n = 0 for n > 1. The arithmetic function that is identically 1 we
denote by 1, and we similarly abbreviate µ(n), (n), and log n by µ, Λ, and
L. In this notation, the characteristic property of µ(n) is that µ ∗ 1 = i, which
is to say that µ and 1 are convolution inverses of each other, and the Möbius
inversion formula takes the compact form
a∗1=b ⇐⇒ a = b ∗ µ.
In the elementary study of prime numbers the relations Λ ∗ 1 = L, L ∗ µ = Λ
are fundamental.
2.1.1 Exercises
1. (de la Vallée Poussin 1898; cf. Landau 1911) Show that
{x/n} = (1 − C0 )x + O x 1/2
n≤x
for σ > 1.
40 The elementary theory of arithmetic functions
for x ≥ 1.
4. (cf. Evelyn & Linfoot 1930) Let N be a positive integer, and suppose that
P is square-free.
(a) Show that the number of residue classes n (mod P 2 ) for which (n, P 2 )
is square-free and (N − n, P 2 ) is square-free is
1 2
P 2
1− 2 1− 2 .
p|P
p p|P
p
p 2 |N p 2 N
(b) Show that the number of integers n, 0 < n < N , for which (n, P 2 ) is
square-free and (N − n, P 2 ) is square-free is
1 2
N 1− 2 1 − 2 + O(P 2 ).
p|P
p p|P
p
p 2 |N p 2 N
(c) Show that the number of n, 0 < n < N , such that n is divisible by the
square of a prime > y is N /y.
(d) Take P to be the product of all primes not exceeding y. By letting y
tend to infinity slowly, show that the number of ways of writing N as
a sum of two square-free integers is ∼ c(N )N where
1 2
c(N ) = a 1+ 2 , a= 1− 2 .
p 2 |N
p −2 p p
5. (cf. Hille 1937) Suppose that f (x) and F(x) are complex-valued functions
defined on [1, ∞). Show that
F(x) = f (x/n)
n≤x
for all x.
6. (cf. Hartman & Wintner 1947) Suppose that | f (n)|d(n) < ∞, and that
|F(n)|d(n) < ∞. Show that
F(n) = f (m)
m
n|m
2.1 Mean values 41
7. (Jarnı́k 1926; cf. Bombieri & Pila 1989) Let C be a simple closed curve in
the plane, of arc length L. Show that the number of ‘lattice points’ (m, n),
m, n ∈ Z, lying on C is at most L + 1. Show that if C is strictly convex
then the number of lattice points on C is 1 + L 2/3 , and that this estimate
is best possible.
8. Let C be a simple closed curve in the plane, of arc length L that encloses
a region of area A. Let N be the number of lattice points inside C. Show
that |N − A| ≤ 3(L + 1).
9. Let r (n) be the number of pairs ( j, k) of integers such that j 2 + k 2 = n.
Show that
r (n) = π x + O x 1/2 .
n≤x
10. (Stieltjes 1887) Suppose that an , bn are convergent series, and that
−1/2
cn = km=n ak bm . Show that cn n converges. (Hence if two Dirichlet
series have abscissa of convergence ≤ σ then the product series γ (s) =
α(s)β(s) has abscissa of convergence σc ≤ σ + 1/2.)
11. (a) Show that n≤x ϕ(n) = (3/π 2 )x 2 + O(x log x) for x ≥ 2.
(b) Show that
1 = −1 + 2 ϕ(n)
m≤x n≤x
n≤x
(m,n)=1
for x ≥ 1. Deduce that the expression above is (6/π 2 )x 2 + O(x log x).
12. Let σ (n) = d|n d. Show that
π2 2
σ (n) = x + O(x log x)
n≤x 12
for x ≥ 2.
13. (Landau 1900, 1936; cf. Sitaramachandrarao 1982, 1985, Nowak 1989)
(a) Show that n/ϕ(n) = d|n µ(d)2 /ϕ(d).
(b) Show that
n ζ (2)ζ (3)
= x + O(log x)
n≤x ϕ(n) ζ (6)
for x ≥ 2.
42 The elementary theory of arithmetic functions
where
1 1
C= 1+ .
8 log 2 p>2
p( p − 2)
(b) Show that for any real number x ≥ 1 and any positive integer q,
1 log p ϕ(q)
= log x + C0 + + O 2ω(q) /x .
m≤x m p|q
p−1 q
(m,q)=1
2.1 Mean values 43
(c) Show that for any real number x ≥ 2 and any positive integer q,
log p
1 ζ (2)ζ (3) p
= 1− 2 log x + C0 +
n≤x ϕ(n) ζ (6) p|q p − p+1 p|q
p−1
(n,q)=1
log p log x
− + O 2ω(q) .
pq
p2 − p + 1 x
18. Let dk (n) be the number of ordered k-tuples (d1 , . . . , dk ) of positive integers
such that d1 d2 · · · dk = n.
(a) Show that dk (n) = d|n dk−1 (d).
∞
(b) Show that n=1 dk (n)n −s = ζ (s)k for σ > 1.
(c) Show that for every fixed positive integer k,
dk (n) = x Pk (log x) + O x 1−1/k (log x)k−2
n≤x
embarrassing that this is the best-known upper bound for gaps between
sums of two squares.)
22. (Feller & Tornier 1932) Let f (n) denote the multiplicative function such
that f ( p) = 1 for all p, and f ( p k ) = −1 whenever k > 1.
(a) Show that
∞
f (n) 2
= ζ (s) 1 −
n=1
ns p p 2s
for σ > 1.
(b) Deduce that
f (n) = µ(d)2ω(d) .
d 2 |n
Let (x) be as in Exercise 23(b). Show that (x) x 1/3 (log x)2 .
27. Let R(x) be as in Exercise 24(c). Show that R(x) x 1/3 log x.
The proof we give below establishes only that there is an x0 such that
ψ(x) x uniformly for x ≥ x0 . However, both ψ(x) and x are bounded away
from 0 and from ∞ in the interval [2, x0 ], and hence the implicit constants can
be adjusted so that ψ(x) x uniformly for x ≥ 2. In subsequent situations of
2.2 Estimates of Chebyshev and of Mertens 47
this sort, we shall assume without comment that the reader understands that it
suffices to prove the result for all sufficiently large x.
which arise in the main terms. To avoid this problem we introduce an idea that
is fundamental to much of prime number theory, namely we replace µ(d) by
an arithmetic function ad that in some way forms a truncated approximation to
µ(d). Suppose that D is a finite set of numbers, and that ad = 0 when d ∈ / D.
Then by (2.11) we see that
ad log d
ad T (x/d) = (x log x − x) ad /d − x + O(log 2x).
d∈D d∈D d∈D
d
(2.12)
Here the implicit constant depends on the choice of ad , which we shall consider
to be fixed. Since we want the above to approximate the relation (2.10), and
since we are hoping that ψ(x) x, we restrict our attention to ad that satisfy
the condition
ad
= 0, (2.13)
d∈D
d
48 The elementary theory of arithmetic functions
in view of (1.20). Thus E(y) will be near 1 for y not too large if ad is near µ(d)
for small d. Moreover, by (2.13) we see that E(y) = − d∈D ad {y/d}, so that
E(y) is periodic with period dividing lcmd∈D d. Hence for a given choice of
the ad , the behaviour of E(y) can be determined by a finite calculation.
The simplest realization of this approach involves taking a1 = 1, a2 = −2,
ad = 0 for d > 2. Then (2.13) holds, the expression (2.14) is log 2, E(y) has
period 2 and E(y) = 0 for 0 ≤ y < 1, E(y) = 1 for 1 ≤ y < 2. Hence for this
choice of the ad the sum in (2.15) satisfies the inequalities
ψ(x) − ψ(x/2) = (k) ≤ (k)E(x/k) ≤ (k) = ψ(x).
x/2<k≤x k≤x k≤x
Thus ψ(x) ≥ (log 2)x + O(log x), which is a lower bound of the desired shape.
In addition,
and
By computing the implicit constants one can use this method to determine a
constant x0 such that ψ(2x) − ψ(x) > x/2 for all x > x0 . Since the contribution
of the proper prime powers is small, it follows that there is at least one prime
in the interval (x, 2x], when x > x0 . After separate consideration of x ≤ x0 ,
one obtains Bertrand’s postulate: For each real number x > 1, there is a prime
number in the interval (x, 2x).
Chebyshev said it, but I’ll say it again:
There’s always a prime between n and 2n.
N. J. Fine
and
ψ(x) x
π(x) = +O .
log x (log x)2
Proof Clearly
∞
ψ(x) = log p = ϑ x 1/k .
p k ≤x k=1
and that
1
∼ log log x.
p≤x p
However, these assertions are weaker than PNT, as we can derive them from
Theorem 2.4.
By Theorem 2.4 the error term is x. Thus (2.11) gives (a). The sum in (b)
differs from that in (a) by the amount
log p log p
≤ 1.
p k ≤x
pk p p( p − 1)
k≥2
by Theorem 2.4. We now prove (d) without determining the value of the con-
stant b. We express (b) in the form L(x) = log x + R(x) where R(x) 1.
Then
1 x x
1 x
d R(u)
= (log u)−1 d L(u) = d log u +
p≤x p 2− 2− log u 2− log u
x
du R(u) x x
= + − R(u) d(log u)−1
2− u log u log u 2− 2−
x
R(u) R(x)
= log log x − log log 2 + 1 + du. +
2 u(log u)2 log x
∞ ∞
The
∞ penultimate term is 1/ log x, and the integral is 2 − x =
2 +O(1/ log x), so we have (d) with
∞
R(u)
b = 1 − log log 2 + du.
2 u(log u)2
As for (e), we note that
−1 −1
1 1 1 1
log 1 − = + log 1 − − .
p≤x p p≤x p p≤x p p
−1
1
log 1 − = log log x + c + O(1/ log x) (2.16)
p≤x p
where c = b + p k≥2 (kp k )−1 . Since e z = 1 + O(|z|) for |z| ≤ 1, on expo-
nentiating we deduce that
−1
1
1− = ec log x + O(1).
p≤x p
To complete the proof it suffices to show that c = C0 . To this end we first note
that if p ≤ x and p k > x, then k ≥ (log x)/ log p. Hence
1 log p log p 1 log p 1
p −k ,
p≤x kp k p≤x (log x) p k p log x k≥2
log x p p 2 log x
p k >x p k >x
52 The elementary theory of arithmetic functions
Since this is trivial when 1 ≤ x < 2, the above holds for all x ≥ 1. We
express this briefly as T1 = T2 + T3 + T4 , and estimate the quantities Ii =
∞ −1−δ
δ 1 x Ti (x) d x. On comparing the results as δ → 0+ we shall deduce
that c = C0 . By Theorem 1.3, Corollary 1.11, and Corollary 1.13 we see that
1
I1 = log ζ (1 + δ) = log + O(δ)
δ
as δ → 0+ . Secondly,
∞
1 ∞ ∞
1 −δn
I2 = δ x −1−δ d x = e = log(1 − e−δ )−1
n=1
n en n=1
n
= log(δ + O(δ 2 ))−1 = log 1/δ + O(δ).
Thirdly,
I3 = c − C 0 ,
and finally
∞ e1/δ ∞
dx dx
I4 δ x −1−δ δ+δ + δ2 x −1−δ d x δ log 1/δ.
1 log 2x 2 x log x e1/δ
Then there is an x0 such that ψ(x) ≤ (a + ε)x for all x ≥ x0 , and hence
x x0 x
ψ(u)u −2 du ≤ ψ(u)u −2 du +(a + ε) u −1 du ≤ (a + ε) log x + Oε (1).
1 1 x0
x
Since this holds for arbitrary ε > 0, it follows that 1 ψ(u)u −2 du ≤ (a +
o(1)) log x. Thus by Theorem 2.7(c) we have a ≥ 1. Similarly lim inf ψ(u)/u
≤ 1.
2.2.1 Exercises
1. (a) Let dn = [1, 2, . . . , n]. Show that dn = eψ(n) .
1
(b) Let P ∈ Z[x], deg P ≤ n. Put I = I (P) = 0 P(x) d x. Show that
I dn+1 ∈ Z, and hence that dn+1 ≥ 1/|I | if I = 0.
(c) Show that there is a polynomial P as above so that I dn+1 = 1.
20≤x≤1 |x2 (1 − x) (2x
(d) Verify that max 2 2
− 1)| = 5−5/2 .
2n
(e) For P(x) = x (1 − x) (2x − 1) , verify that 0 < I < 5−5n .
(f) Show that ψ(10n + 1) ≥ ( 12 log 5) · 10n.
2. Let A be the set of integers composed entirely of primes p ≤ A1 , and
let B be the set of integers composed entirely of primes p > A1 . Then n
is uniquely of the form n = ab, a ∈ A, b ∈ B. Let δ(A1 , A2 ) denote the
density of those n such that a ≤ A2 .
(a) Give a formula for δ(A1 , A2 ).
(b) Show that δ(A1 , A2 ) (log A2 )/ log A1 for 2 ≤ A2 ≤ A1 .
3. Let an = 1 + cos log n, and note that an ≥ 0 for all n.
(a) Show that
∞
1 1
an n −s = ζ (s) + ζ (s + i) + ζ (s − i)
n=1
2 2
for σ > 1.
(b) By Corollary 1.15, or otherwise, show that
an
= log x + O(1).
n≤x n
5. (Chebyshev 1850) From Corollaries 2.5 and 2.8 we see that if there is a
number a such that ψ(x) = (a + o(1))x as x → ∞, then we must have
a = 1. We now take this a step further.
(a) Suppose that there is a number a such that
as x → ∞. Deduce that
x
ψ(u)
du = log x + (a + o(1)) log log x
2 u2
as x → ∞.
(b) By comparing the above with Theorem 2.7(c), deduce that if (2.17)
holds, then necessarily a = 0.
(c) Suppose that there is a constant A such that
x x
π (x) = +o (2.18)
log x − A (log x)2
x
as x → ∞. By writing ϑ(x) = 2− log u dπ (u), integrating by parts,
and estimating the expressions that arise, show that if (2.18) holds,
then
as x → ∞.
(d) Deduce that if (2.18) holds, then A = 1.
Proof Let R be the set of those n for which ϕ(n)/n < ϕ(m)/m for all m < n.
We first prove the inequality for these ‘record-breaking’ n ∈ R. Suppose that
ω(n) = k, and let n ∗ be the product of the first k primes. If n = n ∗ then n ∗ < n
and ϕ(n ∗ )/n ∗ < ϕ(n)/n. Hence R is the set of n of the form
n= p. (2.19)
p≤y
which gives the desired result for n ∈ R. If n ∈/ R then there is an m < n such
that m ∈ R, ϕ(m)/m < ϕ(n)/n. Hence
ϕ(n) ϕ(m) 1 1
> = e−C0 + O
n m log log m log log m
1 1
≥ e−C0 + O .
log log n log log n
We note that equality holds for n of the type (2.19), so the proof is complete.
We now consider the maximum order of d(n). From the pairing d ↔ n/d
√
of divisors, and the fact that at least one of these is ≤ n, it is immediate that
√
d(n) ≤ 2 n. On the other hand, if n is square-free then d(n) = 2ω(n) , which
56 The elementary theory of arithmetic functions
√
can be large, but not nearly as large as n. Indeed, for each ε > 0 there is a
constant C(ε) such that
d(n) ≤ C(ε)n ε (2.20)
for all n ≥ 1. To see this we express n in terms of its canonical factorization,
n = p pa , so that
d(n) a+1
= = f p (a),
nε p paε p
say. Let α p be an integral value of a for which f p (a) is maximized. From the
inequalities f p (α p ) ≥ f p (α p ± 1) we see that
( p ε − 1)−1 − 1 ≤ α p ≤ ( p ε − 1)−1 ,
so that we may take α p = [( p ε − 1)−1 ]. Hence (2.20) holds with
C(ε) = f p (α p ).
p
This constant is best possible, since equality holds when n = p p α p . By
analysing the rate at which C(ε) grows as ε → 0+ , we derive
Theorem 2.11 For all n ≥ 3
log n
log d(n) ≤ (log 2 + O(1/ log log n)) .
log log n
We note that this bound is sharp for n of the form in (2.19).
Proof It suffices to show that there is an absolute constant K such that
C(ε) ≤ exp K ε 2 21/ε , (2.21)
since the stated bound then follows by taking ε = (log 2)/ log log n. We observe
that α p = 0 if p > 21/ε , that α p = 1 if (3/2)1/ε < p ≤ 21/ε , and that α p 1/ε
when p ≤ (3/2)1/ε . Hence
log C(ε) log(2/ p ε ) + log(1/ε).
p≤21/ε p≤(3/2)1/ε
Here the second sum is π (3/2) 1/ε
log 1/ε ε 2 21/ε . The first sum is
(log 2)π (2 ) − εϑ(2 ), and by Corollary 2.5 this is ε 2 21/ε . Thus we have
1/ε 1/ε
so that ω(n) = p X p (n). If we were to treat the X p as though they
were independent random variables then we would have E(X p ) = 1/ p,
Var(X p ) = (1 − 1/ p)/ p. Hence we expect that the average of ω(n) should be
approximately
1
E Xp = E(X p ) = = log log n + O(1),
p≤n p≤n p≤n p
1
ω(n) = x + O (π (x)) .
n≤x p≤x p
and
(ω(n) − log log n)2 x log log x. (2.24)
1<n≤x
so we have
2.3 Applications to arithmetic functions 59
Note that in analytic number theory we say ‘almost all’ when the excep-
tional set has asymptotic density 0; this conflicts with the usage in some
parts of algebra, where the term means that there are at most finitely many
exceptions.
Proof of Theorem 2.12 To prove (2.23) we first multiply out the square on the
left, and write the sum as
x 1
2
1
≤x ≤x = x(log log x)2 + O(x log log x)
p1 = p2
p 1 p 2 p1 p2 ≤x p 1 p 2 p≤x p
p1 = p2
(2.26)
The estimate (2.23) now follows by inserting this and (2.22) in (2.25).
We derive (2.24) from (2.23) by applying the triangle inequality x −
y ≤ x − y for vectors. This gives
1/2 1/2
(ω(n) − log log n)2 − (ω(n) − log log x)2
1<n≤x 1<n≤x
1/2
≤ (log log x − log log n)2 .
1<n≤x
60 The elementary theory of arithmetic functions
and (2.24) follows by squaring both sides and applying (2.23). We omit the
similar argument for (n).
Since 2ω(n) ≤ d(n) ≤ 2(n) for all n, Corollary 2.13 carries an interesting
piece of information for d(n):
for almost all n. Since this is smaller than the average size of d(n), we see that
the average is determined not by the usual size of d(n) but by a sparse set of n for
which d(n) is disproportionately large. Since the first moment (i.e., average) of
d(n) is inflated by the ‘tail’ in its distribution, it is not surprising that this effect
is more pronounced for the higher moments. As was originally suggested by
Ramanujan, it can be shown that for any fixed real number κ there is a positive
constant c(κ) such that
κ
d(n)κ ∼ c(κ)x(log x)2 −1 (2.27)
n≤x
as x → ∞.
In order to handle the error terms that arise in our arguments we are frequently
led to estimate the mean value of multiplicative functions. In most such cases
the method of the hyperbola or the simpler identity (2.3) will suffice, but the
labour involved quickly becomes tiresome. It will therefore be convenient to
have the following result on record, as it is very readily applied.
Then for x ≥ 2,
x f (n)
f (n) (A + 1) .
n≤x log x n≤x n
We note that this is sharper than the trivial estimate
f (n) ≤ x f (n)/n (2.30)
n≤x n≤x
for any fixed κ. Though weaker than (2.27), this is all that is needed in many
cases. We can similarly show that for any fixed real κ,
n κ
x. (2.32)
n≤x ϕ(n)
62 The elementary theory of arithmetic functions
2.3.1 Exercises
1. Let σ (n) = d|n d.
(a) Show that σ (n)ϕ(n) ≤ n 2 for all n ≥ 1 .
(b) Deduce that n + 1 ≤ σ (n) ≤ eC0 n log log n + O(1) for all n ≥ 3.
2.3 Applications to arithmetic functions 63
√
2. Show that d(n) ≤ 3n with equality if and only if n = 12.
3. Let f (n) = p|n (1 + p −1/2 ).
(a) Show that there is a constant a such that if n ≥ 3, then
f (n) < exp a(log n)1/2 (log log n)−1 .
(b) Show that n≤x f (n) = cx + O x 1/2 where c = p (1 + p −3/2 ).
4. Let dk (n) be as in Exercise 2.1.18. Show that if k and κ are fixed, then
κ
dk (n)κ x(log x)k −1 .
n≤x
for x ≥ 2.
5. (Davenport 1932) Let
µ(d) log d
f (n) = − .
d|n
d
(d) Show that the right-hand side above is (log y)/ log x.
2.4 The distribution of (n) − ω(n) 65
(e) Deduce that the second and third terms in (2.35) are 1.
(f) Conclude that
2 = x(log log x)2 + (2b + 1) log log x + O(x)
where b is the constant in Theorem 2.7(d).
(g) Show that the left-hand side of (2.23) is = x log log x + O(x).
(h) Show that the left-hand side of (2.24) is = x log log x + O(x).
9. (cf. Pomerance 1977, Shan 1985) Note that ϕ(n)|(n − 1) when n is prime. An
old – and still unsolved – problem of D. H. Lehmer asks whether there exists
a composite integer n such that ϕ(n)|(n − 1). Let S denote the (presumably
empty) set of such numbers.
(a) Show that if n ∈ S, then n is square-free.
(b) Suppose that mp ∈ S. Show that m ≡ 1 (mod p − 1).
(c) Let p be given. Show that the number of m such that mp ≤ x and mp ∈ S
is x/ p 2 .
(d) Show that the number of n ∈ S, n ≤ x, such that n has a prime factor
> y is x/(y log y).
(e) Suppose that x/y < n ≤ x and that n is composed entirely of primes
p ≤ y. Show that ω(n) ≥ (log x)/(log y) − 1.
(f) By Exercise 4, or otherwise, show that the number of n ≤ x such that
ω(n) ≥ z is x(log x)2 /3z .
(g) Conclude that the number of n ≤ x such that n ∈ S is
√
x/ exp( log x).
uniform in k. This is indeed the case, as we see from the following quantitative
form of Rényi’s theorem.
∞
ζ (s) ζ (s)
µ(n)2 n −s = (1 + p −s ) = (1 + p −s )−1 = λ(d)d −s ,
n=1 p f
ζ (2s) p| f
ζ (2s) d∈D
(n, f )=1
by Theorem 2.2. But d∈D λ(d)/d = p| f (1 + 1/ p)−1 and d∈D d −1/2 =
−1/2 −1
p| f (1 − p ) , so that the proof is complete.
Proof of Theorem 2.16 Let Q denote the set of square-free numbers and F
denote the set of ‘power-full’ numbers (i.e., those f such that p| f ⇒ p 2 | f ).
Every number is uniquely expressible in the form n = q f , q ∈ Q, f ∈ F,
(q, f ) = 1. Hence
Nk = 1.
f ≤x q≤x/ f
f ∈F q∈Q
( f )−ω( f )=k (q, f )=1
2.4 The distribution of (n) − ω(n) 67
1 ⎜ ⎟
6 ⎜ 1/2 ⎟
−1/2 −1 ⎟
x (1 + p ) −1 −1
+O⎜
⎜x f −1/2
1− p ⎟.
π2 f ≤x
f p| f ⎝ f ≤x p| f ⎠
f ∈F f ∈F
( f )−ω( f )=k ( f )−ω( f )=k
In order to appreciate the nature of these sums it is helpful to observe that each
member of F is uniquely of the form a 2 b3 with b square-free, so that there are
x 1/2 members of F not exceeding x. Suppose that z ≥ 1. Then the sum in
the error term is
−1
≤ z −k z ( f )−ω( f ) f −1/2 1 − p −1/2 .
f ≤x p| f
f ∈F
To see that (2.36) holds, it suffices to multiply this by z k and sum over k.
2.4.1 Exercise
1. Let dk be as in (2.36). Show that
dk = c2−k + O(5−k )
where
1 1 −1
c= 1− .
4 p>2
( p − 1)2
2.5 Notes
Section 2.1. Mertens (1874 a) showed that n≤x ϕ(n) = 3x 2 /π 2 + O(x log x).
This refines an earlier estimate of Dirichlet, and is equivalent to Theorem 2.1,
by partial summation. Let R(x) denote the error term in Theorem 2.1. Chowla
(1932) showed that
x
x
R(u)2 du ∼
1 2π 2
as x → ∞, and Walfisz (1963, p. 144) showed that
R(x) (log x)2/3 (log log x)4/3 .
In the opposite direction, Pillai & Chowla (1930) showed (cf. Exercise
7.3.6) that R(x) = (log log log x). That the error term changes sign in-
finitely often was first proved by Erdős & Shapiro (1951), who showed that
R(x) = ± (log log log log x). More recently, Montgomery (1987) showed that
√
R(x) = ± ( log log x). It may be speculated that R(x) log log x and that
R(x) = ± (log log x).
Theorem 2.2 is due to Gegenbauer (1885).
Theorem 2.3 is due to Dirichlet (1849). The problem of improving the error
term in this theorem is known as the Dirichlet divisor problem. Let (x) denote
the error term. Voronoı̈ (1903) showed that (x) x 1/3 log x (see Exercises
2.1.23, 2.1.25, 2.1.26). van der Corput (1922) used estimates of exponential
sums to show that (x) x 33/100+ε . This exponent has since been reduced
2.5 Notes 69
by van der Corput (1928), Chih (1950), Richert (1953), Kolesnik (1969, 1973,
1982, 1985), Iwaniec & Mozzochi (1988), and by Huxley (1993), who showed
that (x) x 23/73+ε . In the opposite direction, Hardy (1916) showed that
(x) = ± (x 1/4 ). Soundararajan (2003) showed that
(x) = x 1/4 (log x)1/4 (log log x)b (log log log x)−5/8
with b = 34 (24/3 − 1), and it is plausible that the first three exponents above are
optimal.
The result of Exercise 2.1.12 generalizes to Rn : A lattice point
(a1 , a, . . . , an ∈ Zn ) is said to be primitive if gcd(a1 , a2 , . . . , an ) = 1. The
asymptotic density of primitive lattice points is easily shown to be 1/ζ (n).
In addition, Cai & Bach (2003) have shown that the density of lattice points
a ∈ Zn such that gcd(ai , a j ) = 1 for all pairs with 1 ≤ i < j ≤ n is
1 n n 1 n−1
1− + 1− .
p p p p
pp. 285–288). Polynomials can be found that produce better constants, but
Gorshkov (1956) showed that the supremum of such constants is < 1, so
the Prime Number Theorem cannot be established by this method. For more
on this subject, see Montgomery (1994, Chapter 10), Pritsker (1999), and
Borwein (2002, Chapter 10).
Theorem 2.7(b)–(e) is due to Mertens (1874a, b). Our determination of the
constant in Theorem 2.7(e) incorporates an expository finesse due to Heath-
Brown.
Section 2.3. Theorem 2.9 is due to Landau (1903). Runge (1885) proved
(2.20), and Wigert (1906/7) showed that d(n) < n (log 2+ε)/ log log n for n > n 0 (ε).
Ramanujan (1915a, b) established the upper bound of Theorem 2.11, first with
an extra log log log n in the error term, and then without. Ramanujan (1915b)
also proved that
log d(n)
< li(n) + O n exp − c log n
log 2
for all n ≥ 2, and that
log d(n)
> li(n) + O n exp − c log n
log 2
for infinitely many n. For a survey of extreme value estimates of arithmetic
functions, see Nicolas (1988).
Theorem 2.12 is due to Turán (1934), although Corollary 2.13 and the es-
timate (2.22) used in the proof of Theorem 2.12 were established earlier by
Hardy & Ramanujan (1917). Kubilius (1956) generalized Turán’s inequality to
arbitrary additive functions. See Tenenbaum (1995, pp. 302–304) for a proof,
and discussion of the sharpest constants.
Theorem 2.14 is due to Hall & Tenenbaum (1988, pp. 2, 11). It represents
a weakening of sharper estimates that can be derived with more work. For
example, Wirsing (1961) showed that if f is a multiplicative function such that
f (n) ≥ 0 for all n, if there is a constant C < 2 such that f ( p k ) C k for all
k ≥ 2, and if
f ( p) ∼ κ x/ log x
p≤x
(1984, 1986, 1987). For a comprehensive account of the mean values of (not
necessarily non-negative) multiplicative functions, see Tenenbaum (1995, pp.
48–50, 308–310, 325–357). The two sides of (2.31) are of the same order of
magnitude, and with more work one can derive a more precise asymptotic
estimate; see Wilson (1922).
Section 2.4. Rényi (1955) gave a qualitative form of Theorem 2.16. Robinson
(1966) gave formulæ for the densities dk . Kac (1959, pp. 64–71) gave a proof
by probabilistic techniques. Generalizations have been given by Cohen (1964)
and Kubilius (1964). Sharper estimates for the error term have been derived
by Delange (1965, 1967/68, 1973), Kátai (1966), Saffari (1970), and Schwarz
(1970).
For a much more detailed historical account of the development of prime
number theory, see Narkiewicz (2000).
2.6 References
Bateman, P. T. (1949). Note on the coefficients of the cyclotomic polynomial, Bull. Amer.
Math. Soc. 55, 1180–1181.
Bateman, P. T. & Grosswald, E. (1958). On a theorem of Erdős and Szekeres, Illinois J.
Math. 2, 88–98.
Bombieri, E. & Pila, J. (1989). The number of integral points on arcs and ovals, Duke
Math. J. 59, 337–357.
Borwein, P. (2002). Computational excursions in analysis and number theory. Canadian
Math. Soc., New York: Springer.
Cai, J.-Y. & Bach, E. (2003). On testing for zero polynomials by a set of points with
bounded precision, Theoret. Comp. Sci. 296, 15–25.
Chebyshev, P. L. (1848). Sur la fonction qui détermine la totalité des nombres premiers
inférieurs à une limite donné, Mem. Acad. Sci. St. Petersburg 6, 1–19.
(1850). Mémoire sur nombres premiers, Mem. Acad. Sci. St. Petersburg 7, 17–33.
(1946). Collected works of P. L. Chebyshev, Vol. 1, Akad. Nauk SSSR, Moscow–
Leningrad.
Chih, T.-T. (1950). A divisor problem, Acta Sinica Sci. Record 3, 177–182.
Chowla, S. (1932). Contributions to the analytic theory of numbers, Math. Zeit. 35,
279–299.
Cohen, E. (1964). Some asymptotic formulas in the theory of numbers, Trans. Amer.
Math. Soc. 112, 214–227.
van der Corput, J. G. (1922). Vereschärfung der Abschätzung beim Teilerproblem, Math.
Ann. 87, 39–65.
(1928). Zum Teilerproblem, Math. Ann. 98, 697–716.
Costa Pereira, N. (1989). Elementary estimates for the Chebyshev function ψ(x) and
for the Möbius function M(x), Acta Arith. 52, 307–337.
Davenport, H. (1932). On a generalization of Euler’s function φ(n), J. London Math. Soc.
7, 290–296; Collected Works, Vol. IV. London: Academic Press, pp. 1827–1833.
72 The elementary theory of arithmetic functions
Huxley, M. N. (1993). Exponential sums and lattice points II. Proc. London Math. Soc.
(3) 66, 279–301.
Iwaniec, H. & Mozzochi, C. J. (1988). On the divisor and circle problems, J. Number
Theory 29, 60–93.
Jarnı́k, V. (1926). Über die Gitterpunkte auf konvexen Curven, Math. Z. 24, 500–
518.
Kac, M. (1959). Statistical Independence in Probability, Analysis and Number Theory,
Carus Monograph 12. Washington: Math. Assoc. Amer.
Kátai, I. (1966). A remark on H. Delange’s paper “Sur un théorème de Rényi”, Magyar
Tud. Akad. Mat. Fiz. Oszt. Közl. 16, 269–273.
Kolesnik, G. (1969). The improvement of the error term in the divisor problem, Mat.
Zametki 6, 545–554.
(1973). On the estimation of the error term in the divisor problem, Acta Arith. 25,
7–30.
(1982). On the order of ζ ( 12 + it) and (R), Pacific J. Math. 82, 107–122.
(1985). On the method of exponent pairs, Acta Arith. 45, 115–143.
Kubilius, J. (1956). Probabilistic methods in the theory of numbers (in Russian), Uspehi
Mat. Nauk 11, 31–66; Amer. Math. Soc. Transl. (2) 19 (1962), 47–85.
(1964). Probabilistic Methods in the Theory of Numbers, Translations of Mathematical
Monographs, Vol. 11. Providence: American Mathematical Society.
Landau, E. (1900). Ueber die zahlentheoretische Function ϕ(n) und ihre Beziehung zum
Goldbachschen Satz, Nachr. Akad. Wiss. Göttingen, 177–186; Collected Works,
Vol. 1. Essen: Thales Verlag, 1985, pp. 106–115.
(1903). Über den Verlauf der zahlentheoretischen Funktion ϕ(x), Arch. Math. Phys.
(3) 5, 86–91; Collected Works, Vol. 1. Essen: Thales Verlag, 1985, pp. 378–383.
(1911). Sur les valeurs moyennes de certaines fonctions arithmétiques, Bull. Acad.
Royale Belgique, 443–472; Collected Works, Vol. 4. Essen: Thales Verlag, 1986,
pp. 377–406.
(1936). On a Titchmarsh–Estermann sum, J. London Math. Soc. 11, 242–245;
Collected Works, Vol. 9. Essen: Thales Verlag, 1987, pp. 393–396.
Linfoot, E. H. & Evelyn, C. J. A. (1929). On a problem in the additive theory of numbers,
I, J. Reine Angew. Math. 164, 131–140.
Ma̧kowski, A. (1960). Partitions into unequal primes, Bull. Acad. Pol. Sci. 8, 125–126.
Massias, J.-P. & Robin, G. (1996). Bornes effectives pour certaines fonctions concernant
les nombres premiers, J. Théor. Nombres Bordeaux 8, 215–242.
Mertens, F. (1874a). Ueber einige asymptotische Gesetze der Zahlentheorie, J. Reine
Angew. Math. 77, 289–338.
(1874b). Ein Beitrag zur analytischen Zahlentheorie, J. Reine Angew. Math. 78,
46–62.
Montgomery, H. L. (1987). Fluctuations in the mean of Euler’s phi function, Proc. Indian
Acad. Sci. (Math. Sci.) 97, 239–245.
(1994). Ten Lectures on the Interface of Analytic Number Theory and Harmonic
Analysis, CBMS 84. Providence: Amer. Math. Soc.
Narkiewicz, W. (2000). The Development of Prime Number Theory. Berlin: Springer-
Verlag.
Nicolas, J.-L. (1988). On Highly Composite Numbers. Ramanujan Revisited (G. E.
Andrews, R. A. Askey, B. C. Berndt, K. G. Ramanathan, R. A. Rankin, eds.). New
York: Academic Press, pp. 215–244.
74 The elementary theory of arithmetic functions
Shan, Z. (1985). On composite n for which ϕ(n)|(n − 1), J. China Univ. Sci. Tech. 15,
109–112.
Sitaramachandrarao, R. (1982). On an error term of Landau, Indian J. Pure Appl. Math.
13, 882–885.
(1985). On an error term of Landau, II, Rocky Mountain J. Math. 15, 579–588.
Soundararajan, K. (2003). Omega results for the divisor and circle problems, Int. Math.
Res. Not., 1987–1998.
Stieltjes, T. J. (1887). Note sur la multiplication de deux séries, Nouvelles Annales (3)
6, 210–215.
Sylvester, J. J. (1881). On Tchebycheff’s theory of the totality of the prime numbers
comprised within given limits, Amer. J. Math. 4, 230–247.
Tenenbaum, G. (1995). Introduction to Analytic and Probabilistic Number Theory, Cam-
bridge Studies 46, Cambridge: Cambridge University Press.
Turán, P. (1934). On a theorem of Hardy and Ramanujan, J. London Math. Soc. 9,
274–276.
de la Vallée Poussin, C. J. (1898). Sur les valeurs moyennes de certaines fonctions
arithmétiques, Ann. Soc. Sci. Bruxelles 22, 84–90.
Voronoı̈, G. (1903). Sur un problème du calcul des fonctions asymptotiques, J. Reine
Angew. Math. 126, 241–282.
Walfisz, A. (1963). Weylsche Exponentialsummen in der neueren Zahlentheorie, Math-
ematische Forschungsberichte 15, Berlin: VEB Deutscher Verlag Wiss.
Ward, D. R. (1927). Some series involving Euler’s function, J. London Math. Soc. 2,
210–214.
Wigert, S. (1906/7). Sur l’ordre de grandeur du nombre des diviseurs d’un entier, Ark.
Mat. 3, 1–9.
Wilson, B. M. (1922). Proofs of some formulæ enunciated by Ramanujan, Proc. London
Math. Soc. 21, 235–255.
Wintner, A. (1944). The Theory of Measure in Arithmetic Semigroups. Baltimore:
Waverly Press.
Wirsing, E. (1961). Das asymptotische Verhalten von Summen über multiplikative Funk-
tionen, Math. Ann. 143, 75–102.
(1967). Das asymptotische Verhalten von Summen über multiplikative Funktionen,
II, Acta Math. Acad. Sci. Hungar. 18, 411–467.
3
Principles and first examples of sieve methods
3.1 Initiation
The aim of sieve theory is to construct estimates for the number of integers
remaining in a set after members of certain arithmetic progressions have been
discarded. If P is given, then the asymptotic density of the set of integers
relatively prime to P is ϕ(P)/P; with the aid of sieves we can estimate how
quickly this asymptotic behaviour is approached. Throughout this chapter we
let S(x, y; P) denote the numbers of integers n in the interval x < n ≤ x + y
for which (n, P) = 1. A first (weak) result is provided by
Proof From the characteristic property (1.20) of the Möbius µ-function, and
the fact that d|(n, P) if and only if d|n and d|P, we see that
S(x, y; P) = µ(d)
x<n≤x+y d|n
d|P
= µ(d) 1
d|P x<n≤x+y
d|n
x + y x
= µ(d) − . (3.1)
d|P
d d
76
3.1 Initiation 77
where
(
s
s = card Sr j .
1≤r1 <···<rs ≤R j=1
where ε(y) → 0 as y → ∞. This bound is very weak, but has the interesting
property of being uniform in x. Since the bound for the error term in Theorem 3.1
is very crude, we might expect that more is true, so that perhaps
ϕ(P)
S(x, y; P) ∼ y
P
even when z is fairly large. However, as we have already noted in our remarks
following Theorem 2.11, this asymptotic formula fails when z = y 1/2 .
In order to derive a sharper estimate for S(x, y; P), we replace µ(d) by a more
general arithmetic function λd that in some sense is a truncated approximation
to µ(d). This is reminiscent of our derivation of the Chebyshev bounds, but in
fact the specific properties required of the λd are now rather different. Suppose
that we seek an upper bound for S(x, y; P). Let λ+ n be a function such that
1 if n = 1,
λ+
d ≥ (3.5)
d|n
0 otherwise.
Such a λ+d we call an ‘upper bound sifting function’, and by arguing as in the
proof of Theorem 3.1 we see that
S(x, y; P) ≤ λ+
d = y λ+d /d + O |λ+
d| . (3.6)
x<n≤x+y d|n d|P d|P
d|P
This will be useful if d|P λ+ d /d is not much larger than ϕ(P)/P, and if
+
d|P |λ d | is much smaller than 2ω(P) . Brun (1915) was the first to succeed
with an argument of this kind. He took his λ+n to be of the form
µ(n) if n ∈ D+ ,
λ+
n =
0 otherwise,
where D+ is a judiciously chosen set of integers. A sieve of this kind is called
‘combinatorial’. With Brun’s choice of D+ it is easy to verify (3.5), and it
is not hard to bound d|P |λ+ d |, but the determination of the asymptotic size
of the main term d|P λ+ d /d presents some technical difficulties. We do not
develop a detailed account of Brun’s method, but the spirit of the approach can
be appreciated by considering the following simple choice of D+ : Let r be an
integer at our disposal, and put
D+ = {n : ω(n) ≤ 2r }.
We observe that
2r
2r
ω(P)
λ+
d = µ(d) = (−1) j
.
d|P j=0 d|P j=0
j
ω(d)= j
3.1 Initiation 79
3.1.1 Exercises
1. (Charles Dodgson) In a very hotly fought battle, at least 70% of the combat-
ants lost an eye, at least 75% an ear, at least 80% an arm, and at least 85% a
leg. What can you say about the percentage that lost all four members?
80 Principles and first examples of sieve methods
2. (P. T. Bateman) Would you believe a market investigator who reports that of
1000 people, 816 like candy, 723 like ice cream, 645 like cake, while 562
like both candy and ice cream, 463 like both candy and cake, 470 like both
ice cream and cake, while 310 like all three?
3. (Erdős 1946) For x > 0 write
ϕ(k)
1= x + E k (x).
1≤n≤x
k
(n,k)=1
(c) By using the result of Exercise B.10, or otherwise, show that if d|k and
e|k, then
k
(d, e)2
B1 ({x/d})B1 ({x/e}) d x = k.
0 12de
(d) Show that if k > 1, then
k
1 ω(k)
E k (x)2 d x = 2 ϕ(k).
0 12
(e) Deduce that if k > 1, then
1/2
ω(k)/2 ϕ(k)
max |E k (x)| 2 .
x k
4. (Lehmer 1955; cf. Vijayaraghavan 1951) Let E k (x) be defined as above.
(a) Show that |E k (x)| ≤ 2ω(k)−1 for all k > 1.
(b) Suppose that k is composed of distinct primes p ≡ 3 (mod 4), and that
ω(k) is even. Show that if d|k, then µ(d)B1 ({k/(4d)}) = −1/4.
(c) Show that there exist infinitely many numbers k for which
5. (Behrend 1948; cf. Heilbronn 1937, Rohrbach 1937, Chung 1941, van der
Corput 1958) Let a1 , . . . , a J be positive integers, and let T (a1 , . . . , a J ) de-
note the asymptotic density of the set of those positive integers that are not
divisible by any of the ai .
(a) Show that T (a1 , . . . , a J ) = Jj=0 (−1) j j where
1
j = .
1≤i 1 <···<i j ≤J
[ai1 , . . . , ai j ]
This simple observation can be used to obtain an upper bound for S(x, y; P);
namely
⎛ ⎞2
⎜ ⎟
S(x, y; P) ≤ ⎝ d ⎠
x<n≤x+y d|n
d|P
= d e 1
d|P x<n≤x+y
e|P d|n,e|n
x+y x
= d e −
d|P
[d, e] [d, e]
e|P
⎛ 2 ⎞
d e
=y +O⎝ |d | ⎠ . (3.10)
d|P
[d, e] d|P
e|P
Theorem 3.2 Let x, y, and z be real numbers such that y > 0 and z ≥ 1. For
any positive integer P we have
y
S(x, y; P) ≤ + O(z 2 L P (z)−2 )
L P (z)
3.2 The Selberg lambda-squared method 83
where
µ(n)2
L P (z) = .
n≤z ϕ(n)
n|P
Hence
d e d e
= ϕ( f )
d|P,e|P
[d, e] f |P d
d e e
f |d|P f |e|P
= ϕ( f )yf2
f |P
where
d
yf = . (3.11)
d
d
f |d|P
Moreover, from these formulæ we see that d = 0 for all d > z if and only if
yf = 0 for all f > z. Thus we have diagonalized the quadratic form in (3.10),
and by (3.12) we see that the constraint 1 = 1 is equivalent to the linear
condition
yf µ( f ) = 1. (3.13)
f |P
2
µ( f ) 1
ϕ( f )yf2 = ϕ( f ) yf − + . (3.14)
f |P f |P
ϕ( f )L P (z) L P (z)
f ≤z
84 Principles and first examples of sieve methods
In order to apply Theorem 3.2, we require a lower bound for the sum L P (z).
To this end we show that
µ(n)2
> log z (3.18)
n≤z ϕ(n)
for all z ≥ 1. Let s(n) denote the largest square-free number dividing n (some-
times called the ‘square-free kernel of n’). Then for square-free n,
1
1 1 1 1
= 1 + + 2 + ··· = ,
ϕ(n) n p|n p p m m
s(m)=n
Here the last inequality is obtained by the integral test. With more work one can
derive an asymptotic formula for the the sum in (3.18) (recall Exercise 2.1.17).
By taking z = y 1/2 in Theorem 3.2, and appealing to (3.18), we obtain
Theorem 3.3 Let P = p≤√ y p. Then for any x and any y ≥ 2,
2y 1
S(x, y; P) ≤ 1+O .
log y log y
By combining the above with (3.3) we obtain an immediate application to
the distribution of prime numbers.
Corollary 3.4 For any x ≥ 0 and any y ≥ 2,
2y 1
π (x + y) − π (x) ≤ 1+O .
log y log y
In Theorem 3.3 we consider only a very special sort of P, but the following
lemma enables us to obtain corresponding results for more general P.
Lemma 3.5 Put M(y; P) = maxx S(x, y; P). If (P, q) = 1, then
q
M(y; P) ≤ M(y; q P).
ϕ(q)
Proof It suffices to show that
q
ϕ(q)S(x, y; P) = S(x + Pm, y; q P), (3.19)
m=1
since the right-hand side is bounded above by q M(y; q P). Suppose that x +
Pm < n ≤ x + Pm + y and that (n, q P) = 1. Put r = n − Pm. Then x <
r ≤ x + y, (r, P) = 1, and (r + Pm, q) = 1. Thus the right-hand side above is
1 = 1.
m r x<r ≤x+y 1≤m≤q
(r,P)=1 (r +Pm,q)=1
Since (P, q) = 1, the map m → r + Pm permutes the residue classes (mod q).
Hence the inner sum above is ϕ(q), and we have (3.19).
Proof Let
P1 = p, q1 = p.
p|P
√ pP
p≤ y √
p≤ y
Theorem 3.3 provides an upper bound for M(y; q1 P1 ), and hence by Lemma
3.5 we have an upper bound for M(y; P1 ). To complete the argument it suffices
to note that S(x, y; P) ≤ S(x, y; P1 ) ≤ M(y; P1 ), and to appeal to Mertens’
formula (Theorem 2.7(e)).
We note that Theorem 3.3 is a special case of Theorem 3.6. Although we have
taken great care to derive uniform estimates, for many purposes it is enough to
know that
1
S(x, y; P) y 1− . (3.20)
p|P
p
p≤y
This follows from Theorem 3.6 since √ y< p≤y (1 − 1/ p)−1 1 by Mertens’
formula. To obtain an estimate in the opposite direction, write P = P1 q1 where
P1 is composed entirely of primes > y, and q1 is composed entirely of primes
≤ y. Since the integers in the interval (0, y] have no prime factor > y, we see
that M(y; P1 ) ≥ [y] . Hence by Lemma 3.5,
1
M(y; P) ≥ [y] 1− . (3.21)
p|P
p
p≤y
Then by Exercise 2.1.17 and Mertens’ estimates (Theorem 2.7) it follows that
this is 14 (3 − 2 log 2) log y + O(1).
3.2.1 Exercises
1. Let d be defined as in the proof of Theorem 3.2.
(a) Show that
d 2z
d log
L P (z)ϕ(d) d
for d ≤ z.
(b) Use the above to give a second proof of (3.17).
2. Show that for y ≥ 2 the number of prime powers p k in the interval
(x, x + y] is
2y 1
≤ 1+O .
log y log y
3. (Chowla 1932) Let f (n) be an arithmetic function, put
g(n) = f (d) f (e),
[d,e]=n
5. (Hensley 1978)
(a) Let P = p≤√ y p. Show that the number of n, x < n ≤ x + y, such
that (n) = 2, is
x + y
x
≤ S(x, y; P) + π −π .
√
p≤ y
p p
(b) By using Theorem 3.3 and Corollary 3.4, show that for y ≥ 2,
2y log log y 1
1≤ 1+O .
x<n≤x+y log y log log y
(n)=2
(c) Let λ−
be a lower bound sifting function such that λ−
d d = 0 for d > z.
Show that for any q,
ϕ(q) λ− λ−
d
≥ d
.
q d
d d
d
(d,q)=1
Theorem 3.8 Suppose that (a, q) = 1, that (P, q) = 1, and that x and y are
real numbers with y ≥ 2q. The number of n, x < n ≤ x + y, such that n ≡ a
(mod q) and (n, P) = 1 is
⎛ ⎞
⎜
C0 y ⎜ 1 ⎟ 1
≤e 1− ⎟ 1+O .
q ⎝ p|P p ⎠ log y/q
√
p≤ y/q
Thus the n that remain after sifting are precisely the n for which (a(n), P) = 1.
By the sieve we obtain upper and lower bounds for the number of remaining n
of the form
λm = λm 1. (3.24)
x<n≤x+y m|(a(n),P) m|P x<n≤x+y
m|a(n)
Now p|a(n) if and only if n ∈ B( p) (mod p). By the Chinese remainder theo-
rem, this will be the case for all p|m when n lies in one of precisely p|m b( p)
residue classes modulo m. The b( p) are defined only for primes, but it is con-
venient now to extend the definition to all positive integers by putting
b(m) = b( p)α .
p α m
Thus b(m) is the totally multiplicative function generated by the b( p). For
square-free m, b(m) represents the number of deleted residue classes modulo
m. We are now in a position to estimate the inner sum above. We partition the
interval (x, x + y] into [y/m] intervals of length m, and one interval of length
{y/m}m. In each interval of length m there are precisely b(m) values of n for
which m|a(n). In the final shorter interval, the number of such n lies between
0 and b(m). Thus the inner sum on the right above is = yb(m)/m + O(b(m)),
and hence the expression (3.24) is
b(m)λm
=y + O b(m)|λm | . (3.25)
m|P
m m|P
92 Principles and first examples of sieve methods
To continue from this point, one should specify the choice of λm , and then
estimate the main term and error term. In the context of Selberg’s 2 method,
we have real d with 1 and d = 0 for d > z. The number of n ∈ (x, x + y]
that survive sifting is
2
≤ d = d e 1
x<n≤x+y d|(a(n),P) d|P e|P x<n≤x+y
[d,e]|a(n)
b([d, e])
=y d e + O g([d, e])|d e | . (3.26)
d|P e|P
[d, e] d|P e|P
This is (3.25) with λm = [d,e]=m d e .
We consider first the main term above. Clearly [d, e] = de/(d, e) and
b([d, e]) = b(d)b(e)/b((d, e)). For square-free m put
b( p)
g(m) = . (3.27)
p|m
p − b( p)
By applying this with m = (d, e) we see that the first sum in (3.26) is
b(d)d b(e)e (d, e) b(d)d b(e)e 1
· · = ·
d|P
d e b((d, e)) d|P
d e f |d
g( f )
e|P e|P f |e
1 b(d) b(e)
= d e
f |P
g( f ) d d e e
f |d|P f |e|P
1 2
= y (3.28)
f |P
g( f ) f
where
b(d)
yf = d . (3.29)
d
d
f |d|P
3.4 Twin primes 93
By the above formulæ we see that the condition that d = 0 for d > z is
equivalent to the condition that yf = 0 for f > z. Also, the condition that
1 = 1 is equivalent to
yf µ( f ) = 1. (3.31)
f |P
where
L= µ( f )2 g( f ) . (3.33)
f ≤z
f |P
Hence
L= µ(m)2 g(m) = c(k) d(n)/n.
m≤z k≤z n≤z/k
The Euler product in (3.35) is absolutely convergent for σ > −1/2. Hence
|c(k)|k −σ < ∞ for σ > −1/2. Thus the two sums in the error terms above
are convergent. Also,
1 ∞
1
|c(k)| ≤ |c(k)| log k .
k>z
log z k=1
log z
Thus by taking s = 0 in (3.35) we find that
1
L = (log z)2 + O(log z). (3.36)
2c
3.4 Twin primes 95
It remains to bound the error term in (3.26). Since 0 ≤ b([d, e]) ≤ b(d)b(e),
the error term is
2
b(d)|d | .
d≤z
Hence
1
b(d)|d | µ(d)2 dg(d) µ(m)2 g(m)
d≤z
L d≤z m≤z/d
1
= µ(m)2 g(m) µ(d)2 dg(d) .
L m≤z d≤z/m
We note that
p
p
1 = 1. (3.37)
b=1 x<n≤x+y x<n≤x+y b=1
/ 1 ( p )
b∈B n ∈B
/ 1 n ∈B
/ 1 / 1 ( p )
b∈B
n≡b ( p ) b≡n ( p )
Consider the inner sum on the right. Since n ∈ / B1 ( p ), the variable b is restricted
to lie in one of p − b1 ( p ) − 1 residue classes. Hence the right-hand side above
is
= ( p − b1 ( p ) − 1)M(x, y; b1 ).
Since there are p − b1 ( p ) values of b in the outer sum on the left-hand side of
(3.37), it follows that there is a choice of b such that b ∈ / B1 ( p ) and
p − b1 ( p ) − 1
1 ≥ M(x, y; b1 ) .
x<n≤x+y p − b1 ( p )
n ∈B
/ 1
n≡b ( p )
3.4 Twin primes 97
Thus
p − b1 ( p)
M(x, y; b1 ) ≤ M(x, y; b2 ) ,
p|P
p − b2 ( p)
and that
r (n)2 x. (3.39)
n≤x
The first of these estimates is easy: Put y = [(log x)/ log 2]. If 0 ≤ k ≤ y − 1,
then 2k ≤ x/2, and if also p ≤ x/2, then p + 2k ≤ x. Thus the sum in (3.38)
is
x
≥ π (x/2)y log x x
log x
for x ≥ 4.
To prove (3.39), we first observe that the sum on the left-hand side is
= 1.
p1 , p2 , j,k
p1 +2 j ≤x
p2 +2k ≤x
p1 +2 j = p2 +2k
3.4.1 Exercises
1. For each prime p let B( p) be the union of b( p) ‘bad’ arithmetic progressions
)
with common difference p. Put B = p|P B( p), and let
m(x, y; b) = min 1
B
x<n≤x+y
n ∈B
/
where the minimum is over all choices of the B( p) with b( p) fixed. Show
that if b1 ( p) ≤ b2 ( p) for all p, then
b1 ( p) −1 b2 ( p) −1
m(x, y; b1 ) 1− ≥ m(x, y; b2 ) 1− .
p p p p
100 Principles and first examples of sieve methods
(d) Show that if a and b are fixed real numbers with a < b, then
(b log p − d( p)) 4(b − a)2 x .
p≤x
a log p≤d( p)≤b log p
(i) Take b = a + 1/8, and suppose that d( p) ≥ a log p for all p > p0 .
Show that the estimates of (f) and (h) are inconsistent if a > 15/16.
Thus conclude that
d( p) 15
lim inf ≤ .
p→∞ log p 16
4. Let r (n) be defined as in the proof of Theorem 3.15. Show that
x
r (n) ∼ .
n≤x log 2
5. Let r (n) be defined as in the proof of Theorem 3.15. Show that
x
r (n) .
n≤x log x
2|n
6. (Erdős 1950)
(a) Show that if n ≡ 1 (mod 3) and k ≡ 0 (mod 2), then 3|(n − 2k ).
(b) Show that if n ≡ 1 (mod 7) and k ≡ 0 (mod 3), then 7|(n − 2k ).
(c) Show that if n ≡ 2 (mod 5) and k ≡ 1 (mod 4), then 5|(n − 2k ).
(d) Show that if n ≡ 8 (mod 17) and k ≡ 3 (mod 8), then 17|(n − 2k ).
(e) Show that if n ≡ 11 (mod 13) and k ≡ 7 (mod 12), then 13|(n − 2k ).
(f) Show that if n ≡ 121 (mod 241) and k ≡ 23 (mod 24), then 241|
(n − 2k ).
(g) Show that every integer k satisfies at least one of the congruences
k ≡ 0 (mod 2), k ≡ 0 (mod 3), k ≡ 1 (mod 4), k ≡ 3 (mod 8), k ≡
7 (mod 12), k ≡ 23 (mod 24).
(h) Show that if n satisfies all the congruences n ≡ 1 (mod 3), n ≡ 1
(mod 7), n ≡ 2 (mod 5), n ≡ 8 (mod 17), n ≡ 11 (mod 13), n ≡
121 (mod 241), then n − 2k is divisible by at least one of the primes
3, 7, 5, 17, 13, 241.
(i) Show that these congruential conditions are equivalent to the single
condition n ≡ 172677 (mod 3728270).
(j) An integer n satisfying the above might still be representable in the
form p + 2k , but if it is, then the prime in question must be one of the
six primes listed. Show that if in addition, n ≡ 9 or 11 or 15 (mod 16),
then n cannot be expressed as a sum of a prime and a power of 2.
3.5 Notes
Sections 3.1, 3.2. The modern era of sieve methods began with the work
of Brun (1915, 1919). Hardy & Littlewood (1922) used Brun’s method to
establish the estimate (3.9). The sharp form of this in Corollary 3.4 is due
102 Principles and first examples of sieve methods
µ(q)2 b( p)
L= . (3.43)
q≤Q 1 + 32 q Q/N p|q
p − b( p)
Here b( p) is the number of residue classes modulo p that are deleted. This is
both a generalization and a sharpening of Theorem 3.2.
Section 3.3. Titchmarsh (1930) used Brun’s method to obtain Theorem 3.9,
but with a larger constant instead of 2. Montgomery & Vaughan (1973) have
shown that Corollary 3.4 and Theorem 3.9 are still valid when the error terms are
omitted. See also Selberg (1991, Section 22). The first significant improvement
of Theorem 3.9 was obtained by Motohashi (1973). Other improvements of
various kinds have been derived by Motohashi (1974), Hooley (1972, 1975),
Goldfeld (1975), Iwaniec (1982), and Friedlander & Iwaniec (1997).
In Lemmas 3.5 and 3.12, and in Exercises 3.2.7, 3.2.9, 3.2.10, 3.4.1 we see
evidence of a monotonicity principle that permeates sieve theory; cf. Selberg
(1991, pp. 72–73).
3.5 Notes 103
Hooley (1994) has shown that quite sharp sieve bounds can be derived using
the interrupted inclusion–exclusion idea that Brun started with. This approach
has been developed further by Ford & Halberstam (2000). An exposition of
sieves based on these ideas is given by Bateman & Diamond (2004, Chapters 12,
13). Still more extensive accounts of sieve methods have been given by Greaves
(2001), Halberstam & Richert (1974), Iwaniec & Kowalski (2004, Chapter
6), Motohashi (1983), and Selberg (1971, 1991). In addition, a collection of
applications of sieves to arithmetic problems has been given by Hooley (1976),
and additional sieve ideas are found in Bombieri (1977), Bombieri, Friedlander
& Iwaniec (1986, 1987, 1989), Fouvry & Iwaniec (1997), Friedlander & Iwaniec
(1998a, b), and Iwaniec (1978, 1980a, b, 1981).
Section 3.4. The twin prime conjecture is a special case of the prime k-tuple
conjecture. Suppose that d1 , . . . , dk are distinct integers, and let b( p) denote
the number of distinct residue classes modulo p found among the di . The prime
k-tuple conjecture asserts that if b( p) < p for every prime number p, then there
exist infinitely many positive integers n such that the k numbers n + di are all
prime. Hardy & Littlewood (1922) put this in a quantitative form: If b( p) < p
for all p, then the number of n ≤ N for which the k numbers n + di are all
prime is conjectured to be
N
∼ S(d) (3.44)
(log N )k
as N → ∞ where
−k
b( p) 1
S(d) = 1− 1− . (3.45)
p p p
This product is absolutely convergent, since b( p) = k for all sufficiently large
primes p. Although this remains unproved, by sifting we can obtain an upper
bound of the expected order of magnitude. In particular, from (3.43) it can be
shown that the number of n, M + 1 ≤ n ≤ M + N , for which the numbers
n + di are all prime is
N
2k k!S(d) . (3.46)
(log N )k
Corollarys 3.4 and 3.14 are special cases of this.
Theorem 3.15 is due to Romanoff (1934). Once the bound for the number
of twin primes is in place, the hardest part of the proof is to establish the
estimate (3.41). Romanoff’s original proof of this was rather difficult. Erdős
& Turán (1935) gave a simpler proof, but the clever proof we have given is
due to Erdős (1951). Let r (n) be defined as in the proof of Theorem 3.15.
Erdős (1950) showed that r (n) = (log log n), and that n≤x r (n)k k x for
104 Principles and first examples of sieve methods
any positive k. Presumably r (n) = o(log n), but for all we know there could be,
although it seems unlikely, infinitely many n such that n − 2k is prime whenever
0 < 2k < n. The number n = 105 has this property, and is probably the largest
such number. The best upper bound we have for the number of such n not
exceeding X is (Vaughan 1973),
c log X log log log X
X exp − .
log log X
3.6 References
Ankeny, N. C. & Onishi, H. (1964/1965). The general sieve, Acta Arith. 10, 31–62.
Bateman, P. T. & Diamond, H. (2004). Analytic Number Theory, Hackensack: World
Scientific.
Behrend, F. A. (1948). Generalization of an inequality of Heilbronn and Rohrbach, Bull.
Amer. Math. Soc. 54, 681–684.
Bombieri, E. (1977). The asymptotic sieve, Rend. Accad. Naz. XL (5) 1/2 (1975/76),
243–269.
Bombieri, E., Friedlander, J. B., & Iwaniec, H. (1986). Primes in arithmetic progressions
to large moduli, Acta Math. 156, 203–251.
(1987). Primes in arithmetic progressions to large moduli, II, Math. Ann. 277, 361–
393.
(1989). Primes in arithmetic progressions to large moduli, III, J. Amer. Math. Soc. 2,
215–224.
Brun, V. (1915). Über das Goldbachsche Gesetz und die Anzahl der Primzahlpaare,
Archiv for Math. og Naturvid. B 34, no. 8, 19 pp.
(1919). La série 1/5 + 1/7 + 1/11 + 1/13 + 1/17 + 1/19 + 1/29 + 1/31 +
1/41 + 1/43 + 1/59 + 1/61 + · · · où les dénominateurs sont “nombres premiers
jumeaus” est convergente ou finie, Bull. Sci. Math. (2) 43, 100–104; 124–128.
(1967). Reflections on the sieve of Eratosthenes, Norske Vid. Selsk. Skr. Trondheim,
no. 1, 9 pp.
Buchstab, A. A. (1938). New improvements in the method of the sieve of Eratosthenes,
Mat. Sb. (N. S.) 4 (46), 375–387.
Chowla, S. (1932). Contributions to the analytic theory of numbers, Math. Z. 35, 279–
299.
Chung, K.-L. (1941). A generalization of an inequality in the elementary theory of
numbers, J. Reine Angew. Math. 183, 193–196.
van der Corput, J. G. (1958). Inequalities involving least common multiple and other
arithmetical functions, Nederl. Akad. Wetensch. Proc. Ser. A 61 (= Indag. Math.
20), 5–15.
Erdős, P. (1940). The difference of consecutive primes, Duke Math. J. 6, 438–441.
(1946). On the coefficients of the cyclotomic polynomial, Bull. Amer. Math. Soc. 52,
179–184.
3.6 References 105
(1950). On integers of the form 2k + p and some related problems, Summa Brasil.
Math. 2, 113–123.
(1951). On some problems of Bellman and a theorem of Romanoff, J. Chinese Math.
Soc. (N. S.) 1, 409–421.
Erdős, P. & Turán, P. (1935). Ein zahlentheoretischer Satz, Mitt. Forsch. Inst. Math.
Mech. Univ. Tomsk 1, 101–103.
Ford, K. & Halberstam, H. (2000). The Brun–Hooley sieve, J. Number Theory 81,
335–350.
Fouvry, E. & Iwaniec, H. (1997). Gaussian primes, Acta Arith. 79 (1997), 249–287.
Friedlander, J. B. & Iwaniec, H. (1997). The Brun–Titchmarsh theorem, Analytic Number
Theory (Kyoto, 1996). London Math. Soc. Lecture Note Ser. 247, Cambridge:
Cambridge University Press, pp. 85–93.
(1998a). The polynomial X 2 + Y 4 captures its primes, Ann. of Math. (2) 148, 945–
1040.
(1998b). Asymptotic sieve for primes, Ann. of Math. (2) 148, 1041–1065.
Goldfeld, D. M. (1975). A further improvement of the Brun–Titchmarsh theorem, J.
London Math. Soc. (2) 11, 434–444.
Greaves, G. (2001). Sieves in Number Theory. Berlin: Springer.
Halberstam, H. (1985). Lectures on the linear sieve, Topics in Analytic Number Theory
(Austin, 1982). Austin: University of Texas Press, pp. 165–220.
Halberstam, H. & Richert, H.-E. (1973). Brun’s method and the fundamental lemma,
Acta Arith. 24, 113–133.
(1974). Sieve Methods. London: Academic Press.
(1975). Brun’s method and the fundamental lemma. II, Acta Arith. 27, 51–59.
Hardy, G. H. & Littlewood, J. E. (1922). Some problems of ‘Partitio Numerorum’: III.
On the expression of a number as a sum of primes, Acta Math. 44, 1–70; Collected
Papers, Vol. I, London: Oxford University Press, 1966, pp. 561–630.
Heilbronn, H. (1937). On an inequality in the elementary theory of numbers, Proc.
Cambridge Philos. Soc. 33, 207–209.
Hensley, D. (1978). An almost-prime sieve, J. Number Theory 10, 250–262; Corrigen-
dum, 12, (1980), 437.
Hooley, C. (1972). On the Brun–Titchmarsh theorem, J. Reine Angew. Math. 255,
60–79.
(1975). On the Brun–Titchmarsh theorem, II, Proc. London Math. Soc. (3) 30, 114–
128.
(1976). Applications of Sieve Methods to the Theory of Prime Numbers, Cambridge
Tract 70. Cambridge: Cambridge University Press.
(1994). An almost pure sieve, Acta Arith. 66, 359–368.
Iwaniec, H. (1978). Almost-primes represented by quadratic polynomials, Invent. Math.
47, 171–188.
(1980a). Rosser’s sieve, Acta Arith. 36, 171–202.
(1980b). A new form of the error term in the linear sieve, Acta Arith. 37, 307–320.
(1981). Rosser’s sieve – bilinear forms of the remainder terms – some applications.
Recent Progress in Analytic Number Theory, Vol. 1. New York: Academic Press,
pp. 203–230.
(1982). On the Brun–Titchmarsh theorem, J. Math. Soc. Japan 34, 95–123.
Iwaniec, H. & Kowalski, E. (2004). Analytic Number Theory, Colloquium Publications
53. Providence: Amer. Math. Soc.
106 Principles and first examples of sieve methods
Boston: Academic Press, pp. 467–484; Collected Papers, Vol. 1. Berlin: Springer-
Verlag, 1989, pp. 675–69.
(1991). Lectures on Sieves, Collected Papers, Vol. 2. Berlin: Springer-Verlag,
pp. 65–247.
Titchmarsh, E. C. (1930). A divisor problem, Rend. Circ. Math. Palermo 54, 414–429.
Tsang, K. M. (1989). Remarks on the sieving limit of the Buchstab–Rosser sieve, Number
Theory, Trace Formulas and Discrete Groups (Oslo, 1987). Boston: Academic
Press, pp. 485–502.
Vaughan, R. C. (1973). Some applications of Montgomery’s sieve, J. Number Theory 5,
64–79.
Vijayaraghavan, T. (1951). On a problem in elementary number theory, J. Indian Math.
Soc. (N.S.) 15, 51–56.
4
Primes in arithmetic progressions: I
or
1 1 ∞
f (z) − f (−z) = cn z n .
2 2 n=0
n≡1 (2)
108
4.1 Additive characters 109
unless ζ = 1. Hence
1
q
1 if n ≡ a (mod q),
e(−ka/q)e(kn/q) = (4.1)
q k=1 0 otherwise,
Here the exact values that k runs through are immaterial, as long as the set of
these values forms a complete residue system modulo q. Hence we may replace
k by −k in the above, and so we see that
q
f (n) = *
f (k)e(kn/q). (4.3)
k=1
and Fourier expansion of a function f ∈ L 1 (T), but the situation here is simpler
because our sums have only finitely many terms.
Let v(h) be the vector v(h) = [e(h/q), e(2h/q), . . . , e((q − 1)h/q), 1].
From (4.1) we see that two such vectors v(h 1 ) and v(h 2 ) are orthogonal un-
less h 1 ≡ h 2 (mod q). These vectors are not normalized, but they all have the
√
same length q, so apart from some rescaling, the transformation from f to * f
is an isometry. More precisely, if f has period q and * f is given by (4.2), then
by (4.3),
q
q
q 2
| f (n)|2 = *
f (k)e(kn/q) .
n=1 n=1 k=1
By expanding and taking the sum over n inside, we see that this is
q q
q
= *
f ( j) *
f (k) e( jn/q)e(−kn/q).
j=1 k=1 n=1
Finally,
µ(q/(q, n))
cq (n) = dµ(q/d) = ϕ(q). (4.7)
d|(q,n)
ϕ(q/(q, n))
Proof The first assertion is evident, as each term in the sum (4.5) has period
q. As for the second, suppose that q = q1 q2 where (q1 , q2 ) = 1. By the Chinese
Remainder Theorem, for each a (mod q) there is a unique pair a1 , a2 with ai
determined (mod qi ), so that a ≡ a1 q2 + a2 q1 (mod q). Moreover, under this
correspondence we see that (a, q) = 1 if and only if (ai , qi ) = 1 for i = 1, 2.
Then
q1
q2
cq (n) = e((a1 q2 + a2 q1 )n/(q1 q2 ))
a1 =1 a2 =1
(a1 ,q1 )=1 (a2 ,q2 )=1
⎛ ⎞⎛ ⎞
⎜ ⎟⎜
q1 q2
⎟
=⎝ e(a1 n/q1 )⎠ ⎝ e(a2 n/q2 )⎠
a1 =1 a2 =1
(a1 ,q1 )=1 (a2 ,q2 )=1
q
q
e(na/q) = e(na/q)
a=1 d|q a=1
(a,q)=d
q/d
= e(nb/(q/d))
d|q b=1
(b,q/d)=1
= cq/d (n).
d|q
By (4.1), the left-hand side above is q when q|n, and is 0 otherwise. Thus we
have (4.6).
The first formula in (4.7) is merely the Möbius inverse of (4.6). To obtain
the second formula in (4.7), we begin by considering the special case in which
q is a prime power, q = p k .
k
p
c pk (n) = e(na/ p k )
a=1
pa
k k−1
p p
= e(na/ p ) −k
e(na/ p k−1 ).
a=1 a=1
112 Primes in arithmetic progressions: I
Here the first sum is p k if p k |n, and is 0 otherwise. Similarly, the second
sum is p k−1 if p k−1 |n, and is 0 otherwise. Hence the above is
⎧
⎨0 if p k−1 n,
= −p k−1
if p k−1 n,
⎩ k
p − p k−1 if p k |n
k
µ p /(n, p k )
= k ϕ( p k ).
ϕ p /(n, p k )
The general case of (4.7) now follows because cq (n) is a multiplicative function
of q.
4.1.1 Exercises
√
1. Let U = [u kn ] be the q × q matrix with elements u kn = e(kn/q)/ q. Show
that UU ∗ = U ∗ U = I , i.e., that U is unitary.
2. (Friedman 1957; cf. Reznick 1995)
(a) Show that
1
2r 2r
ue(θ/2) + ve(−θ/2) dθ = u r vr
0 r
for any non-negative integer r and arbitrary complex numbers u, v.
(b) Show that if u = (x − i y)/2, v = (x + i y)/2, then
x cos π θ + y sin π θ = ue(θ/2) + ve(−θ/2)
for all θ.
(c) Show that
1
2r 2r
x cos π θ + y sin π θ dθ = 2−2r (x 2 + y 2 )r
0 r
for any non-negative integer r and arbitrary real or complex numbers
x, y.
(d) Show that
q
πia/q
−πia/q 2r 2r
ue + ve =q u r vr
a=1
r
if r is an integer, 0 ≤ r < q.
(e) Show that
q
2r
(x cos πa/q + y sin πa/q)2r = q 2−2r (x 2 + y 2 )r
a=1
r
if r is an integer, 0 ≤ r < q.
4.1 Additive characters 113
where
−1
6µ(q) 1
aq = 2 2 1− 2 .
π q p|q
p
by taking
1 1
aq = lim F(n)cq (n) . (4.9)
ϕ(q) x→∞ x n≤x
In the following, suppose that f (r ) is chosen so that F(n) = r |n f (r ) for
all n.
(a) Suppose that
∞
| f (r )|
< ∞. (4.10)
r =1
r
114 Primes in arithmetic progressions: I
as x → ∞.
(b) Suppose that (4.10) holds. Show that
1 ∞
f (r )
lim F(n)cq (n) = ϕ(q) .
x→∞ x r
n≤x r =1
q|r
(c) Put
∞
f (r )
aq = .
r =1
r
q|r
Show that if
∞
| f (r )|d(r )
<∞ (4.11)
r =1
r
∞
then (4.8) and (4.9) hold, and moreover that q=1 |aq cq (n)| < ∞.
∞
8. (Ramanujan 1918) Show that if q > 1, then n=1 cq (n)/n = −(q). (See
also Exercise 8.3.4.)
9. Let q (z) denote the q th cyclotomic polynomial, i.e., the monic polynomial
whose roots are precisely the primitive q th roots of unity, so that
q
q (z) = (z − e(n/q)).
n=1
(n,q)=1
and that (z d − 1)µ(q/d) has a power series expansion, valid when |z| < 1,
with integer coefficients. Deduce that q (z) ∈ Z[z].
(b) Suppose that z ∈ Z and p | q (z) and let e denote the order of z modulo
p. Show that e | q and that if p | (z d − 1) then e | d.
(c) Choose t so that p t (z e − 1). Show that for m ∈ N with p m one has
p t (z me − 1).
(d) Show that if p q, then p ht q (z) where h = µ(q/d). Deduce that
e|d|q
e = q and that q | ( p − 1).
(e) By taking z to be a suitable multiple of q, or otherwise, show that there
are infinitely many primes p with p ≡ 1 (mod q).
4.2 Dirichlet characters 115
Lemma 4.2 Suppose that G is cyclic of order n, say G = (a). Then there are
exactly n characters of G, namely χk (a m ) = e(km/n) for 1 ≤ k ≤ n. Moreover,
n if χ = χ0 ,
χ (g) = (4.12)
g∈G
0 otherwise,
and
n if g = e,
χ (g) = (4.13)
*
0 otherwise.
χ ∈G
* is cyclic, G
In this situation, G * = (χ1 ).
Since the characters are now known explicitly, the remaining assertions are
easily verified.
Next we describe the characters of the direct product of two groups in terms
of the characters of the factors.
Lemma 4.3 Suppose that G 1 and G 2 are finite abelian groups, and that G =
G 1 ⊗ G 2 . If χi is a character of G i , i = 1, 2, and g ∈ G is written g = (g1 , g2 ),
gi ∈ G i , then χ (g) = χ1 (g1 )χ2 (g2 ) is a character of G. Conversely, if χ ∈ G, *
then there exist unique χi ∈ G i such that χ (g) = χ1 (g1 )χ2 (g2 ). The identities
(4.12) and (4.13) hold for G if they hold for both G 1 and G 2 .
* corresponds to a pair (χ1 , χ2 ) ∈ G
We see here that each χ ∈ G *1 × G
* 2 . Thus
∼ * *
G = G1 ⊗ G2.
Proof The first assertion is clear. As for the second, put χ1 (g1 ) = χ ((g1 , e2 )),
*i for i = 1, 2, and χ1 (g1 )χ2 (g2 ) = χ (g). The
χ2 (g2 ) = χ ((e1 , g2 )). Then χi ∈ G
χi are unique, for if g = (g1 , e2 ), then
χ (g) = χ ((g1 , e2 )) = χ1 (g1 )χ2 (e2 ) = χ1 (g1 ),
and similarly for χ2 . If χ (g) = χ1 (g1 )χ2 (g2 ), then
χ (g) = χ1 (g1 ) χ2 (g2 ) ,
g∈G g1 ∈G 1 g2 ∈G 2
so that (4.12) holds for G if it holds for G 1 and for G 2 . Similarly, if g = (g1 , g2 ),
then
⎛ ⎞⎛ ⎞
χ (g) = ⎝ χ1 (g1 )⎠ ⎝ χ2 (g2 )⎠ ,
*
χ ∈G *1
χ1 ∈G *2
χ1 ∈ G
Though G and G * are isomorphic, the isomorphism is not canonical. That is,
no particular one-to-one correspondence between the elements of G and those
* is naturally distinguished.
of G
4.2 Dirichlet characters 117
If (n, q) = 1, then
ϕ(q) if n ≡ 1 (mod q),
χ (n) = (4.15)
χ
0 otherwise,
where the sum is extended over the ϕ(q) Dirichlet characters χ (mod q).
Proof The first assertion follows immediately from the observations that
χ1 (n)χ2 (n) is totally multiplicative, that it vanishes if (n, [q1 , q2 ]) > 1, and
that it has period [q1 , q2 ]. As for the second assertion, we may suppose that
(n, q) = 1. By the Chinese Remainder Theorem we see that
(Z/qZ)× ∼
= (Z/q1 Z)× ⊗ (Z/q2 Z)×
Our proof of Theorem 4.4 depends on Abel’s theorem that any finite abelian
group is isomorphic to the direct product of cyclic groups, but we can prove
Corollary 4.5 without appealing to this result, as follows. By the Chinese Re-
mainder Theorem we see that
+
(Z/qZ)× ∼ = (Z/ p α Z)× .
p α q
If p is odd, then the reduced residue classes (mod p α ) form a cyclic group; in
classical language we say there is a primitive root g. Thus if (n, p) = 1, then
there is a unique ν (mod ϕ( p α )) such that g ν ≡ n (mod p α ). The number ν is
118 Primes in arithmetic progressions: I
called the index of n, and is denoted ν = indg n. From Lemma 4.2 it follows
that the characters (mod p α ), p > 2, are given by
k indg n
χk (n) = e (4.16)
ϕ( p α )
for (n, p) = 1. We obtain ϕ( p α ) different characters by allowing k to assume
integral values in the range 1 ≤ k ≤ ϕ( p α ). By Lemma 4.3 it follows that if q
is odd, then the general character (mod q) is given by
k indg n
χ (n) = e (4.17)
p α q
ϕ( p α )
when (n, q) = 1.
By definition, if f (n) is totally multiplicative, f (n) = 0 whenever (n, q) > 1,
and f (n) has period q, then f is a Dirichlet character (mod q). It is useful to
note that the first condition can be relaxed.
Theorem 4.7 If f is multiplicative, f (n) = 0 whenever (n, q) > 1, and f has
period q, then f is a Dirichlet character modulo q.
Proof It suffices to show that f is totally multiplicative. If (mn, q) > 1, then
f (mn) = f (m) f (n) since 0 = 0. Suppose that (mn, q) = 1. Hence in partic-
ular (m, q) = 1, so that the map k → n + kq (mod m) permutes the residue
classes (mod m). Thus there is a k for which n + kq ≡ 1 (mod m), and
4.2 Dirichlet characters 119
4.2.1 Exercises
1. Let G be a finite abelian group of order n. Let g1 , g2 , . . . , gn denote the
elements of G, and let χ1 (g), χ2 (g), . . . , χn (g) denote the characters of G.
√
Let U = [u i j ] be the n × n matrix with elements u i j = χi (g j )/ n. Show
that UU ∗ = U ∗ U = I , i.e., that U is unitary.
2. Show that for arbitrary real or complex numbers c1 , . . . , cq ,
q 2
q
cn χ (n) = ϕ(q) |cn |2
χ n=1 n=1
(n,q)=1
where the sum on the left-hand side runs over all Dirichlet characters
χ (mod q).
3. Show that for arbitrary real or complex numbers cχ ,
q 2
cχ χ (n) = ϕ(q) |cχ |2
n=1 χ χ
where the sum over χ is extended over all Dirichlet characters (mod q).
4. Let (a, q) = 1, and suppose that k is the order of a in the multiplicative group
of reduced residue classes (mod q).
(a) Show that if χ is a Dirichlet character (mod q), then χ (a) is a k th root
of unity.
(b) Show that if z is a k th root of unity, then
k if z = 1,
1 + z + ··· + z k−1
=
0 otherwise.
(b) Show that each k th root of unity occurs precisely ϕ(q)/k times among the
numbers χ (a) as a runs over the ϕ(q) reduced residue classes (mod q).
6. Let χ be a character (mod q) such that χ (a) = ±1 whenever (a, q) = 1, and
q
put S(χ ) = n=1 nχ (n). Thus S(χ ) is an integer.
(a) Show that if (a, q) = 1 then aχ (a)S(χ ) ≡ S(χ ) (mod q).
(b) Show that there is an a such that (a, q) = 1 and (aχ (a) − 1, q)|12.
(c) Deduce that 12S(χ ) ≡ 0 (mod q).
In algebraic number fields we encounter not only Dirichlet characters, but
also characters of ideal class groups and of Galois groups. In addition, algebraic
number fields possessing one or more complex embeddings also have a further
kind of character, Hecke’s Grössencharaktere. In a sequence of exercises, be-
ginning with the one below, we develop the basic properties of these characters
√
for the Gaussian field Q( −1).
7. Let K be the Gaussian field,
√
K =Q −1 = {a + bi : a, b ∈ Q},
O K = {a + bi : a, b ∈ Z}.
Elements α = a + bi ∈ K have a norm, N (α) = a 2 + b2 , and we observe
that N (αβ) = N (α)N (β). An element α of a ring is a unit if α has an inverse
in the ring. The ring O K has precisely four units, namely i k for k = 0, 1, 2, 3.
Two elements α, β ∈ O K are associates if α = uβ for some unit u. For each
integer m we define the Hecke Grössencharakter
4mi arg α
e if α = 0,
χm (α) =
0 if α = 0.
for k = 1, 2, 3, . . . . Hence
χ (n) ≤ q (4.23)
n≤x
for any x, so that by Theorem 1.3, the series (4.20) converges for σ > 0. This
result is best possible since the terms in (4.20) do not tend to 0 when σ = 0. On
the other hand, we shall show in Chapter 10 that the function L(s, χ) is entire
if χ = χ0 . For σ > 1 we can take logarithms in (4.21), and differentiate, as in
Corollary 1.11, and thus we obtain
In these last formulæ we see how relations for L-functions parallel those
for the zeta functions. Indeed, when manipulating Dirichlet series formally, the
only property of n −s that is used is that it is totally multiplicative. Hence all
such calculations can be made with n −s replaced by χ (n)n −s . For example, we
know that µ(n)2 n −s = ζ (s)/ζ (2s) for σ > 1. Hence formally
∞
µ(n)2 χ (n)n −s = L(s, χ )/L(2s, χ 2 ). (4.26)
n=1
where the sum is extended over all characters χ (mod q). This is the multiplica-
tive analogue of (4.1). Hence if (a, q) = 1 then
∞
1 ∞
(n)n −s = (n)n −s χ (a)χ (n)
n=1
ϕ(q) n=1 χ
n≡a (q)
−1 L
= χ (a) (s, χ) (4.28)
ϕ(q) χ L
for σ > 1. As L(s, χ0 ) has a simple pole at s = 1, the function LL (s, χ) has a
simple pole at 1 with residue −1. Thus the term arising from χ0 on the right-hand
side above is
1
+ Oq (1) (4.29)
ϕ(q)(s − 1)
Suppose that (a, q) = 1. Then the above, with (4.28), (4.29), and (4.30) give
the estimate
∞
1
(n)n −s = + Oq (1)
n=1
ϕ(q)(s − 1)
n≡a (q)
4.3 Dirichlet L-functions 123
as s → 1+ . Consequently
∞
(n)
= ∞.
n=1
n
n≡a (q)
We call a character real if all its values are real (i.e., χ (n) = 0 or ±1 for all
n). Otherwise a character is complex. A character is quadratic if it has order
2 in the character group: χ 2 = χ0 but χ = χ0 . Thus a quadratic character is
real, and a real character is either principal or quadratic. In Chapter
9 we shall
express quadratic characters in terms of the Kronecker symbol dn .
If we take s = σ > 1, then the sum above is a non-negative real number, and
hence we see that
L(σ, χ) ≥ 1 (4.32)
χ
for σ > 1. Now L(s, χ0 ) has a simple pole at s = 1, but the other L(s, χ)
are analytic at s = 1. Thus L(1, χ ) = 0 can hold for at most one χ , since
otherwise the product in (4.32) would tend to 0 as σ → 1+ . If χ is a character
(mod q), then χ is a character (mod q), and χ = χ if χ is complex. Moreover
124 Primes in arithmetic progressions: I
Hence r (n) ≥ 0 for all n, and r (n 2 ) ≥ 1 for all n. Suppose that L(1, χ ) = 0.
Then ζ (s)L(s, χ) is analytic for σ > 0, and by Landau’s theorem (Theorem
1.7) the series r (n)n −s converges for σ > 0. But this is false, since
∞
∞
∞
r (n)n −1/2 ≥ r (n 2 )n −1 ≥ n −1 = +∞.
n=1 n=1 n=1
Hence L(1, χ ) = 0. Since L(σ, χ ) > 0 for σ > 1 when χ is quadratic, we see
in fact that L(1, χ ) > 0 in this case.
By using the techniques of Chapter 2 we can prove more than the mere
divergence of the series in Corollary 4.10.
This last error term is χ 1, and then (a) follows from (4.33) and the fact that
L(1, χ) = 0. The derivation of (b) from (a), and of (c) from (b) proceeds as in
the proof of Theorem 2.7. Continuing as in that proof, we see from (c) that
(n)χ (n)
1
= c(χ ) + Oχ
1<n≤x
n log n log x
where
χ ( pk )
c(χ ) = b(χ ) + .
k
kp k
p
k>1
We let s → 1+ in (4.24), and deduce by Theorem 1.1 that c(χ) = log L(1, χ ).
To complete the derivation of (d) it suffices to argue as in the proof of
Theorem 2.7.
log p 1
(b) = log x + Oq (1),
p≤x p ϕ(q)
1
n≡a (q)
1 1
(c) = log log x + b(q, a) + Oq ,
p≤x p ϕ(q) log x
n≡a (q)
−1
1 1
(d) 1− = c(q, a)(log x)1/ϕ(q) 1 + Oq
p≤x p log x
n≡a (q)
where
1 1 1
b(q, a) = C0 + log 1 − + χ (a) log L(1, χ) −
ϕ(q) p χ =χ0
kp k
p|q k
p ≡a (q)
k>1
and
1/ϕ(q)
ϕ(q) 1 −χ ( p) χ ( p)
c(q, a) = e C0
L(1, χ )χ (a) 1− 1− .
q χ =χ0 p p p
Proof To derive (a) from Theorem 4.11(a) we use (4.27) and the estimate
(n)χ0 (n)
= log x + Oq (1),
n≤x n
which follows from Theorem 2.7(a) since
log p log p
= q 1.
k
pk p|q
p−1
p
p|q
We derive (b) and (c) similarly from the corresponding parts of Theorem 4.11.
In the latter case we use the estimate
χ0 ( p)
1
= log log x + b(χ0 ) + Oq
p≤x p log x
where
1 χ0 ( p k )
b(χ0 ) = C0 + log 1 − − .
p|q
p k
kp k
p
k>1
To derive (d) we observe first that
−1
χ0 ( p) −1 1 1
1− = 1− 1− ,
p≤x p p≤x p p≤x p
p|q
which by Theorem 2.7(e) is
⎛ ⎞−1
ϕ(q) ⎜ 1 ⎟ −C0 1
= ⎝ 1− ⎠ e (log x) 1 + O .
q p|q
p log x
p>x
4.3 Dirichlet L-functions 127
Here each term in the product is 1 + O(1/x), and the number of factors is
≤ ω(q), so the product is 1 + Oq (1/x), and hence the above is
ϕ(q) 1
= eC0 (log x) 1 + Oq .
q log x
To complete the proof it suffices to combine this with Theorem 4.11(d)
in (4.27).
4.3.1 Exercises
1. Let χ be a Dirichlet character (mod q). Show that if σ > 1, then
∞
(a) (−1)n−1 χ (n)n −s = (1 − χ (2)21−s )L(s, χ );
n=1
∞
L(s, χ)4
(b) d(n)2 χ (n)n −s = .
L(2s, χ 2 )
n=1
2. (Mertens 1895a,b) Let r (n) = d|n χ (d).
(a) Show that if χ is a non-principal character (mod q), then
χ (n) 1
√ χ √ .
n>x n x
(b) Show that if χ is a non-principal character (mod q), then
r (n)
1/2
= 2x 1/2 L(1, χ ) + Oχ (1).
n≤x n
(c) Recall that if χ is quadratic then r (n) ≥ 0 for all n, and that r (n 2 ) ≥ 1.
Deduce that if χ is a quadratic character, then the left-hand side above
is log x.
(d) Conclude that if χ is a quadratic character, then L(1, χ ) > 0.
3. (Mertens 1897, 1899) For u ≥ 0, put f (u) = m≤u (1 − m/u).
(a) Show that f (u) ≥ 0, that f (u) is continuous, and that if u is not an
integer, then
[u]([u] + 1)
f (u) = ;
2u 2
deduce that f is increasing.
(b) Show also that
u 1 u u 1
f (u) = − {v} dv = − + O(1/u) .
2 u 0 2 2
(c) Let r (n) = d|n χ (d), and assume that χ is non-principal. Show that
r (n)(1 − n/x) = χ (d) f (x/d) .
n≤x d≤x
128 Primes in arithmetic progressions: I
(d) Write d≤x = d≤y + y<d≤x = S1 + S2 where 1 ≤ y ≤ x. Use
part (b) to show that S1 = 2 x L(1, χ ) + Oχ (x/y) + O(y 2 /x).
1
(d) Show that for any positive integer q there is a small number cq and a
large number Cq such that if x ≥ 2Cq and (a, q) = 1, then
log p
> cq .
x/Cq < p≤x
p
p≡a (q)
(e) Show that for any positive integer q there is a Cq such that if (a, q) = 1,
then
x
π (x; q, a) q
log x
uniformly for x ≥ Cq .
(f) Show that if (a, q) = 1, then
π(x; q, a) 1 π (x; q, a) 1
lim inf ≤ , lim sup ≥ .
x→∞ x/ log x ϕ(q) x→∞ x/ log x ϕ(q)
6. (a) Show that
x
ϑ(x) ≤ π (x) log x ≤ ϑ(x) + O
log x
for x ≥ 2.
(b) Let P denote a set of prime numbers, and put
πP (x) = 1, ϑP (x) = log p.
p≤x p≤x
p∈P p∈P
Show that
x
ϑP (x) = πP (x) log x + O
log x
for x ≥ 2, where the implicit constant is absolute.
(c) Let
n= p.
p≤y
p∈P
as y → ∞.
7. Let R(n) denote the number of ordered pairs a, b such that a 2 + b2 = n
with a ≥ 0 and b > 0. Also, let r (n) denote
the number of such pairs for
−4
which (a, b) = 1. Finally, let χ−4 = n be the non-principal character
(mod 4). We recall that if the prime factorization of n is written in the form
n = 2α pβ qγ ,
γ
β
p n q n
p≡1 (4) q≡3 (4)
then r (n) > 0 if and only if γ = 0 for all primes q and α ≤ 1. We also
recall that
p (β + 1) if 2|γ for all q,
R(n) = r (n/d 2 ) = χ−4 (d) =
d 2 |n d|n
0 otherwise.
(a) Show that ∞ R(n)n −s = ζ (s)L(s, χ−4 ) for σ > 1.
n=1
∞
(b) Show that n=1 r (n)n −s = ζ (s)L(s, χ−4 )/ζ (2s) for σ > 1.
(c) Show that if x ≥ 0 and y ≥ 2, then
y
card{n ∈ (x, x + y] : r (n) > 0} √ .
log y
(d) Show that
x
card{n ≤ x : R(n) > 0} √
log x
for x ≥ 2.
(e) Suppose that n is of the form
n= p.
p≤y
p≡1 (4)
In the above it is noteworthy that although R(n) ≤ d(n) for all n, that
R(n) is usually 0 and has a smaller average value (cf. Exercise 2.1.9)
than d(n) (cf. Theorem 2.3), the maximum order of magnitude of R(n)
is the same as for d(n).
4.3 Dirichlet L-functions 131
√
8. Let K = Q( −1) be the Gaussian field, O K = {a + ib : a, b ∈ Z} the ring
of integers in K . Ideals a in O K are principal, a = (a + ib), and have norm
N (a) = a 2 + b2 .
(a) Explain why the number of ideals a with N (a) ≤ x is π4 x + O(x 1/2 ).
(b) For σ > 1, let ζ K (s) = a N (a)−s be the Dedekind zeta function of
K . Show that ζ K (s) = ζ (s)L s, χ−4 .
(c) For the Gaussian field K , show that N (ab) = N (a)N (b). (This is true
in any algebraic number field.)
(d) Assume that ideals in K factor uniquely into prime ideals. (This is true
in any algebraic number field, and is particularly easy to establish for
the Gaussian field since it has a division algorithm.) Deduce that if
σ > 1, then
1 −1
ζ K (s) = 1−
p N (p)
for σ > 1.
(f) Let a and b be given ideals. Show that
1 if gcd(a, b) = 1,
µ(d) =
d|a
0 otherwise.
d|b
(g) Among pairs a, b of ideals with N (a) ≤ x, N (b) ≤ x, show that the
probability that gcd(a, b) = 1 is
1 6
+ O x −1/2 = 2 + O x 1/2 .
ζ K (2) π L 2, χ−4
9. (Erdős 1946, 1949, 1957, Vaughan 1974, Saffari, unpublished, but see
Bateman, Pomerance & Vaughan 1981; cf. Exercise 2.3.7) Let q (z) =
µ(q/d)
d|q (z − 1)
d
denote the q th cyclotomic polynomial. Suppose that
q= p
p≤y
p≡±2 (5)
(d) Deduce that q (z) has a coefficient whose absolute value is at least
exp q (log 2−ε)/ log log q
if y > y0 (ε).
√
10. Grössencharaktere for Q( −1), continued from Exercise 4.2.7.
(a) For σ > 1 put
1
L(s, χm ) = χm (α)N (α)−s = χm (a + bi)(a 2 + b2 )−s
α∈O
4 a,b∈Z
K
(a,b)=(0,0)
where α denotes a sum over unassociated members of O K . Show
that the above sum is absolutely convergent in this half-plane.
(b) We recall that members of O K factor uniquely into Gaussian primes.
Also, the Gaussian primes are obtained by factoring the rational primes:
The prime 2 ramifies, 2 = i 3 (1 + i)2 , the rational primes p ≡ 1 (mod 4)
split into two distinct Gaussian primes, p = (a + bi)(a − bi), and the
rational primes q ≡ 3 (mod 4) are inert. Show that
L(s, χm ) = (1 − χm (p)N (p)−s )−1
p
(f) Show that if m = 0, then the Dirichlet series L(s, χm ) is convergent for
σ > 1/2.
(g) Show that L(s, χm ) and L(s, χ−m ) are identically equal, and hence that
L(σ, χm ) ∈ R for σ > 1/2.
4.4 Notes 133
4.4 Notes
Section 4.1. Ramanujan’s sum was introduced by Ramanujan (1918). Incredi-
bly, both Hardy and Ramanujan missed the fact that cq (n) be written in closed
form: The formula on the extreme right of (4.7) is due to Hölder (1936). Nor-
mally one would say that a function f is even if f (x) = f (−x). However, in
the present context, an arithmetic function f with period q is said to be even
if f (n) is a function only of (n, q). Thus cq (n) is an even function. The space
of almost-even functions is rather small, but includes several arithmetic func-
tions of interest. For such functions one may hope for a representation in the
∞
form f (n) = q=1 aq cq (n), called a Ramanujan expansion. For a survey of the
theory of such expansions, see Schwarz (1988). Hildebrand (1984) established
definitive results concerning the pointwise convergence of Ramanujan expan-
sions. An appropriate Parseval identity has been established for mean-square
summable almost-even functions; see Hildebrand, Schwarz & Spilker (1988).
Section 4.2. The first instance of characters of a non-cyclic group occurs in
Gauss’s analysis of the genus structure of the class group of binary quadratic
forms. The quotient of the class group by the principal genus is isomorphic to
C2 ⊗ C2 ⊗ · · · ⊗ C2 , and the associated characters are given by Kronecker’s
symbol. Dirichlet (1839) defined the Dirichlet characters for the multiplicative
group (Z/qZ)× of reduced residues modulo q, and the same technique suffices
to construct the characters for any finite Abelian group. More generally, if
G is a group, then a homomorphism h : G −→ G L(n, C) is called a group
representation, and the trace of h(g) is a group character. Note that if a and
b are conjugate elements of G, say a = gbg −1 , then h(a) and h(b) are similar
matrices. Hence they have the same eigenvalues, and in particular tr h(a) =
tr h(b). Thus a group character is constant on conjugacy classes. In the case of a
finite Abelian group it suffices to take n = 1, and in this case the representation
and its trace are essentially the same. For an introduction to characters in a
wider setting, see Serre (1977).
Section 4.3. Dirichlet (1837a,b,c) first proved Corollary 4.10 in the case that
q is prime. The definition of the Dirichlet characters is not difficult in that case,
since the multiplicative group (Z/ pZ)× of reduced residues is cyclic. The most
challenging part of the proof is to show that L(1, χ ) when χ is the Legendre
symbol (mod p). If p ≡ 3 (mod 4), then
p−1
p−1
a p( p − 1)
a ≡ a= ≡ 1 (mod 2),
a=1
p a=1
2
and hence the sum on the left is non-zero. It follows by (9.9) that L(1, χ p ) = 0
in this case. If p ≡ 1 (mod 4), then one has the identity of Exercise 9.3.7(c),
134 Primes in arithmetic progressions: I
4.5 References
Baker, A., Birch, B. J., & Wirsing, E. A. (1973). On a problem of Chowla, J. Number
Theory 5, 224–236.
Bateman, P. T. (1959). Theorems implying the non-vanishing of χ(m)m −1 for real
residue-characters, J. Indian
Math. Soc. 23, 101–115.
(1966). Lower bounds for h(m)/m for arithmetical function h similar to real
residue characters, J. Math. Anal. Appl. 15, 2–20.
4.5 References 135
The interchange of limits here is difficult to justify, since α(s) may not be
uniformly convergent, and because the integral in (5.3) is neither uniformly nor
absolutely convergent. Moreover, if x is an integer, then the term n = x in (5.4)
gives rise to the integral (5.3) with y = 1, and this integral does not converge,
although its Cauchy principal value exists:
σ0 +i T
1 ds 1
lim = (5.5)
T →∞ 2πi σ −i T s 2
0
137
138 Dirichlet series: II
Here indicates that if x is an integer, then the last term is to be counted with
weight 1/2.
Proof Choose N so large that N > 2x + 2, and write
α(s) = an n −s + an n −s = α1 (s) + α2 (s),
n≤N n>N
here the justification is trivial since there are only finitely many terms. As for
α2 (s), we observe that
∞ ∞
α2 (s) = u −s d(A(u) − A(N )) = s (A(u) − A(N ))u −s−1 du.
N N
We have now established a precise relationship between (5.1) and (5.2), but
Theorem 5.1 is not sufficiently quantitative to be useful in practice. We express
the error term more explicitly in terms of the sine integral
∞
sin u
si(x) = − du.
x u
By integration by parts we see that si(x) 1/x for x ≥ 1, and hence that
si(x) min(1, 1/x) (5.6)
for x > 0. We also note that
+∞
sin u
si(x) + si(−x) = − du = −π. (5.7)
−∞ u
Theorem 5.2 If σ0 > max(0, σa ) and x > 0, then
1 σ0 +i T
xs
an = α(s) ds + R (5.8)
n≤x 2πi σ0 −i T s
where
1 x
R= an si T log
π x/2<n<x
n
1 n 4σ0 + x σ0 |an |
− an si T log +O σ0
.
π x<n<2x x T n n
Proof Since the series α(s) is absolutely convergent on the interval [σ0 −
i T, σ0 + i T ], we see that
1 σ0 +i T
xs 1 σ0 +i T
x s ds
α(s) ds = an .
2πi σ0 −i T s n 2πi σ0 −i T n s
Thus it suffices to show that
⎧
⎪
⎪ 1 + O(y σ0 /T ) if y ≥ 2,
σ0 +i T ⎨ σ0
1 ds 1 + 1
si(T log y) + O(2 /T ) if 1 ≤ y ≤ 2,
ys = π
σ0
2πi σ0 −i T s ⎪
⎪ − si(T log 1/y) + O(2 /T ) if 1/2 ≤ y ≤ 1,
1
⎩ π σ0
O(y /T ) if y ≤ 1/2
(5.9)
for σ0 > 0.
To establish the first part of this formula, suppose that y ≥ 2, and let C be
the piecewise linear path from −∞ − i T to σ0 − i T to σ0 + i T to −∞ + i T .
Then by the calculus of residues we see that
1 ds
ys = 1,
2πi C s
140 Dirichlet series: II
In classical
1 harmonic analysis, for f ∈ L1 (T) we define Fourier coefficients
*
f (k) = 0 f (x)e(−kα) dα, and we expect that the Fourier series *f (k)e(kα)
provides a useful formula for f (α). As it happens, the Fourier series may
diverge, or converge to a value other than f (α), but for most f a satisfactory
alternative can be found. For example, if f is of bounded variation, then
f (α − ) + f (α + ) K
= lim *
f (k)e(kα).
2 K →∞
−K
As in the case of Fourier series, this may fail, but it is not difficult to show that
if f is of bounded variation on [−A, A] for every A, then
f (α − ) + f (α + ) T
*
= lim f (t)e(t x) dt. (5.12)
2 T →∞ −T
The relationship between (5.1) and (5.2) is precisely the same as between
(5.10) and (5.11). Indeed, if we take f (x) = A(e2π x )e−2π σ x , then f ∈ L 1 (R) by
Theorem 1.3, and by changing variables in (5.1) we find that
* α(σ + it)
f (t) = .
2π (σ + it)
Thus (5.2) is equivalent to (5.11), and an appeal to (5.12) provides a second
(real variable) proof of Theorem 5.1.
In general, if
∞
F(s) = f (x)x s−1 d x, (5.13)
0
then we say that F(s) is the Mellin transform of f (x). By (5.10) and (5.11) we
expect that
σ0 +i∞
1
f (x) = F(s)x −s ds, (5.14)
2πi σ0 −i∞
and when this latter formula holds we say that f is the inverse Mellin transform
of F. Thus if A(x) is the summatory function of a Dirichlet series α(s), then
α(s)/s is the Mellin transform of A(1/x) for σ > max(0, σc ), and Perron’s
formula (Theorem 5.1) asserts that if σ0 > max(0, σc ), then A(1/x) is the inverse
142 Dirichlet series: II
and that
σ0 +i∞
1
Aw (x) = α(s)K (s)x s ds. (5.16)
2πi σ0 −i∞
Alternatively, we may start with a kernel K (s), and define the weight w(x)
to be its inverse Mellin transform. The precise conditions under which these
identities hold depends on the weight or kernel; we mention several important
examples.
1. Cesàro weights. For a positive integer k, put
1
Ck (x) = an (x − n)k . (5.17)
k! n≤x
x
Then Ck (x) = 0 Ck−1 (u) du for k ≥ 1 where C0 (x) = A(x), and hence
Ck (x) x θ for θ > k + max(0, σc ). (The implicit constant here may depend
on k, on θ, and on the an .) By integrating (5.1) by parts repeatedly, we see
that
∞
α(s) = s(s + 1) · · · (s + k) Ck (x)x −s−k−1 d x (5.18)
1
for σ > max(0, σc ). By following the method used to prove Theorem 5.1, it
may also be shown that
σ0 +i∞
1 x s+k
Ck (x) = α(s) ds (5.19)
2πi σ0 −i∞ s(s + 1) · · · (s + k)
when x > 0 and σ0 > max(0, σc ). Here the critical step is to show that if y ≥ 1
and σ0 > 0, then
σ0 +i∞ k
1 ys ys
ds = Res
2πi σ0 −i∞ s(s + 1) · · · (s + k) j=0
s(s + 1) · · · (s + k) s=− j
5.1 The inverse Mellin transform 143
for σ > max(0, σc ). By following the method used to prove Theorem 5.1 we
also find that
σ0 +i∞
1 xs
Rk (x) = α(s) ds (5.22)
2πi σ0 −i∞ s k+1
when x > 0 and σ0 > max(0, σc ). Here the critical observation is that if y ≥ 1
and σ0 > 0, then
σ0 +i∞ s
1 ys y 1
ds = Res k+1 = (log y)k .
2πi σ0 −i∞ s k+1 s s=0 k!
3. Abelian weights. For σ > 0 we have
∞ ∞
(s) = e−u u s−1 du = n s e−nx x s−1 d x.
0 0
where
∞
P(x) = an e−nx . (5.24)
n=1
These operations are valid for σ > max(0, σa ), but by partial summation
P(x) x −θ as x → 0+ for θ > max(0, σc ), so that the integral in (5.23) is
absolutely convergent in the half-plane σ > max(0, σc ). Hence the integral is
an analytic function in this half-plane, so that by the principle of uniqueness
144 Dirichlet series: II
of analytic continuation it follows that (5.23) holds for σ > max(0, σc ). In the
opposite direction,
σ0 +i∞
1
P(x) = α(s) (s)x −s ds (5.25)
2πi σ0 −i∞
for x > 0, σ > max(0, σc ). To prove this we recall from Theorem 1.5 that
α(s) τ uniformly for σ ≥ ε + max(0, σc ), and from Stirling’s formula
π
(Theorem C.1) we see that | (s)| e− 2 |t| |t|σ −1/2 as |t| → ∞ with σ bounded.
Thus the value of the integral is independent of σ0 , and in particular we may
assume that σ0 > max(0, σa ). Consequently the terms in α(s) can be integrated
individually, and it suffices to appeal to Theorem C.4.
The formulæ (5.23) and (5.25) provide an important link between the Dirich-
let series α(s) and the power series generating function P(x). Indeed, these
formulæ hold for complex x, provided that x > 0. In particular, by taking
x = δ − 2πiα we find that
∞
1 σ0 +i∞
an e(nα)e−nδ = α(s) (s)(δ − 2πiα)−s ds.
n=1
2πi σ0 −i∞
It may be noted in the above examples that smoother weights w(x) give rise
to kernels K (s) that tend to 0 rapidly as |t| → ∞. Further useful kernels can
be constructed as linear combinations of the above kernels.
Since the Mellin transform is a Fourier transform with altered variables, all
results pertaining to Fourier transforms can be reformulated in terms of Mellin
transforms. Particularly useful is Plancherel’s identity, which asserts that if f ∈
L 1 (R) ∩ L 2 (R), then f 2 = *
f 2 . This is the analogue for Fourier transforms
of Parseval’s identity for Fourier series, which asserts that k | * f (k)|2 = f 22 .
By the changes of variables we noted before, we obtain
∞
Theorem 5.4 (Plancherel’s identity) Suppose that 0 |w(x)|x −σ −1 d x < ∞,
∞ ∞
and also that 0 |w(x)|2 x −2σ −1 d x < ∞. Put K (s) = 0 w(x)x −s−1 d x. Then
∞ +∞
2π |w(x)|2 x −2σ −1 d x = |K (σ + it)|2 dt.
0 −∞
5.1.1 Exercises
1. Show that if σc < σ0 < 0, then
1 σ0 +i T
xs
lim α(s) ds = a .
n>x n
T →∞ 2πi σ0 −i T s
2. (a) Show that if y ≥ 0, then
π
− = si(0) ≤ si(y) ≤ si(π ) = 0.28114 . . . .
2
(b) Let β > 0 be fixed. Show that if x > 0 and σ0 > max(0, σc ), then
1 σ0 +i∞ ∞
β
α(s) (s/β)x s ds = β an e−(n/x) .
2πi σ0 −i∞ n=1
(b) Explain why the values of the integrals above are independent of the
value of σ0 . Hence show that if σ0 = −b/a 2 , then the above is
+∞
e−b /(2a
2 2
)
1
e−a t /2
e−b /a .
2 2 2 2
= dt = √
2π −∞ 2π a
(c) Show that if a > 0, x > 0 and σ0 > σc , then
1 σ0 +i∞
a 2 s 2/2 s 1 ∞
(log x/n)2
α(s)e x ds = √ an exp − .
2πi σ0 −i∞ 2π a n=1 2a 2
for σ > σw .
(a) Show that Aw (x) = ∞
n=1 an w(x/n) satisfies Aw (x) x θ for θ >
max(σw , σc ).
(b) Show that
∞
K (s)α(s) = Aw (x)x −s−1 d x
0
10. Suppose that F is strictly increasing, and that for i = 1, 2 the functions f i
are real-valued with f i ∈ L 1 (R) ∩ L 2 (R) and F( f i ) ∈ L 1 (R) ∩ L 2 (R).
(a) Show that
+∞
( f 1 (x) − f 2 (x))(F( f 1 (x)) − F( f 2 (x))) d x
−∞
+∞
= f*1 (t) − f*2 (t) F(
f 1 )(t) − F( f 2 )(t) dt.
−∞
5.2 Summability
We say that an infinite series an is Abel summable to a, and write an = a
(A) if
∞
lim− an r n = a.
r →1
n=0
Abel proved that if a series converges, then it is A-summable to the same value.
Because of this historical antecedent, we call a theorem ‘Abelian’ if it states
that one kind of summability implies another. Perhaps the simplest Abelian
theorem asserts that if ∞ n=1 an converges to a, then
N
n
lim 1− an = a. (5.27)
N →∞ N
n=1
1 N
lim sn = a. (5.28)
N →∞ N n=1
∞
lim tmn = 1. (5.31)
m→∞
n=1
We now show that regular transformations preserve limits, and relegate the
verification of the converse to exercises.
The important special case (5.28) is obtained by noting that the (semi-infinite)
matrix [tmn ] with
1/m if 1 ≤ n ≤ m,
tmn =
0 if n > m
To establish the second assertion, suppose that ε > 0 and that |an | < ε for
n > N = N (ε). Now
N
|bm | ≤ |tmn an | + |tmn an | = 1 + 2 ,
n=1 n>N
say. From (5.29) and the argument above with A = ε we see that 2 ≤ Cε.
From (5.30) we see that limm→∞ 1 = 0. Hence lim supm→∞ |bm | ≤ Cε, and
we have the desired conclusion since ε is arbitrary. Finally, suppose that T is
regular and that limn→∞ an = a. We write an = a + αn , so that
∞ ∞
bm = a tmn + tmn αn .
n=1 n=1
To see how this may also be derived from Theorem 5.5, let {sm } be an arbitrary
sequence of points of S for which limm→∞ sm = s0 . It suffices to show that
limm→∞ α(sm ) = α(s0 ). Take
tmn = n s0 −sm − (n + 1)s0 −sm ,
so that
∞
n
−s0
α(sm ) = tmn ak k .
n=1 k=1
In view of Theorem 5.5, it suffices to show that [tmn ] is regular. The conditions
(5.30) and (5.31) are clearly satisfied, and (5.29) follows on observing that if
s ∈ S, then s − s0 H σ − σ0 , so that
n+1
n s0 −s − (n + 1)s0 −s = (s − s0 ) u s0 −s−1 du
n
n+1
H
(σ − σ0 ) u σ0 −σ −1 du
n
= n σ0 −σ − (n + 1)σ0 −σ .
Thus we have the result. Abel’s analogous theorem on the convergence of power
series can be derived similarly from Theorem 5.5.
150 Dirichlet series: II
The converse of Abel’s theorem on power series is false, but Tauber (1897)
proved a partial converse: If an = o(1/n) and an = a (A), then an = a.
Following Hardy and Littlewood, we call a theorem ‘Tauberian’ if it provides
a partial converse of an Abelian theorem. The qualifying hypothesis (‘an =
o(1/n)’ in the above) is the ’Tauberian hypothesis’. For simplicity we begin
with partial converses of (5.27).
Theorem 5.6 If ∞ n=1 an = a (C, 1), then an = a provided that one of the
following hypotheses holds:
(a) an ≥ 0 for n ≥ 1;
(b) an = O(1/n) for n ≥ 1;
(c) There is a constant A such that an ≥ −A/n for all n ≥ 1.
Proof Clearly (a) implies (c). If (b) holds, then both an and an satisfy (c).
Thus it suffices to prove that an = a when (c) holds. We observe that if H
is a positive integer, then
N
N + H N +H
n N N
n
an = an 1 − − an 1 −
n=1
H n=1
N+H H n=1 N
1
− an (N + H − n) (5.33)
H N <n<N +H
= T1 − T2 − T3 ,
say. Take H = [εN ] for some ε > 0. By hypothesis, lim N →∞ T1 = a(1 + ε)/ε,
and lim N →∞ T2 = a/ε. From (c) we see that
1 AH
T3 ≥ −A ≥− ≥ −Aε.
N <n<N +H
n N
Hence on combining these estimates in (5.33) we see that
N
lim sup an ≤ a + Aε.
N →∞ n=1
so that
N
lim inf an ≥ a,
N →∞
n=1
If we had argued from (a) or (b), then the treatment of the term T3 above
would have been simpler, since from (a) it follows that T3 ≥ 0, while from
(b) we have T3 ε.
Our next objective is to generalize and strengthen Theorem 5.6. The type of
generalization we have in mind is exhibited in the following result, which can
be established by adapting the above proof: Let β be fixed, β ≥ 0. If
N
n
an 1 − = (a + o(1))N β ,
n=1
N
β (β) = (β + 1) (5.38)
when β > 0.
The amount of unsmoothing required in deriving (5.37) from (5.35) is now
much greater than it was in the proof of Theorem 5.6. Nevertheless we follow
the same line of attack. To obtain the proper perspective we review the preceding
proof. Let J = [0, 1], let χJ (u) be its characteristic function, and put K (u) =
N N
max(0, 1 − u) for u ≥ 0. Thus n=1 an = n an χJ (n/N ), and n=1 an (1 −
n/N ) = n an K (n/N ). Our strategy was to approximate to χJ (u) by linear
combinations of K (κu) for various values of κ, κ > 0. The relation underlying
(5.33) and (5.34) is both simple and explicit:
1 1
K (u) − (1− ε)K (u/(1 − ε)) ≤ χJ (u) ≤ ((1+ ε)K (u/(1+ ε)) − K (u));
ε ε
(5.39)
we took ε = H/N . In the present situation we wish to approximate to χJ (u) by
linear combinations of e−κu , κ > 0. We make the change of variable x = e−u ,
so that 0 ≤ x ≤ 1, and we put J = [1/e, 1]. Then we want to approximate to
χJ (x) by a linear combination P(x) of the functions x κ , κ > 0. In fact it suffices
to use only integral values of κ, so that P(x) is a polynomial that vanishes at
the origin. In place of (5.33), (5.34) and (5.39) we shall substitute
Lemma 5.8 Let ε be given, 0 < ε < 1/4, and put J = [1/e, 1], K =
[e−1−ε , e−1+ε ]. There exist polynomials P± (x) such that for 0 ≤ x ≤ 1 we have
and
Proof Let g(x) = (χJ (x) − x)/(x(1 − x)). Then g is continuous in [0, 1]
apart from a jump discontinuity at x = 1/e of height e2 /(e − 1) < 5. Hence
by Weierstrass’s theorem on the uniform approximation of continuous func-
tions by polynomials we see that there are polynomials Q ± (x) such that
Q − (x) ≤ g(x) ≤ Q + (x) for 0 ≤ x ≤ 1, and for which
for 0 ≤ x ≤ 1. Then the polynomials P± (x) = x + x(1 − x)Q ± (x) have the
desired properties.
and that
∞ ∞
v β−1 e−vδ dv = δ −β w β−1 e−w dw = δ −β (β).
0 0
Hence if b(u) = a(u) − α(u + 1)β−1 / (β), then b(u) ≥ −B(u + 1)β−1 , and
∞
b(u)e−uδ du = o(δ −β ).
0
U
Thus 0 b(u) du = o(U β ), so that
U
α
a(u) du = U β + o(U β ),
0 β (β)
and we have (5.37), in view of (5.38).
For the remaining case, β = 0, it suffices to consider b(u) = a(u) −
αχ[0,1] (u).
∞
Corollary 5.9 Suppose that p(z) = n=0 an z n converges for |z| < 1, and
that β ≥ 0. If p(x) = (α + o(1))(1 − x)−β as x → 1− , and if an ≥ −An β−1
for n ≥ 1, then
N
α
an = + o(1) N β .
n=0
(β + 1)
∞
Proof Take β = 1, p(z) = ∞ n=0 sn z = (1 − z)
n −1 n
n=0 an z in Corollary
N
5.9. Then n=0 sn = (α + o(1))N , which is the desired result.
so that
an
α
≥ + o(1) (log N )β − A1 (log N )β−1 .
n≤N
n (β + 1)
then
N
1 2
an (N − n) = N + O N 3/2 . (5.45)
n=1
2
156 Dirichlet series: II
This is best possible (take an = 1 + n −1/2 ), but if the error term is oscilla-
√
tory, then smoothing may reduce its size (consider an = cos n). Conversely if
(5.45) holds and if the sequence an is bounded, then the method used to prove
Theorem 5.6 can be used to show that
N
an = N + O N 3/4 . (5.46)
n=1
This error term, though weak, is best possible (take an = 1 + cos(log n)2 ).
For Dirichlet series it can be shown that if
∞
1
α(s) = an n −s = + O(1)
n=1
s−1
This is also best possible (take an = 1 + cos(log log n)2 ), but we can obtain a
sharper result by strengthening our analytic hypothesis. For example, it can be
shown that if α(s) is analytic in a neighbourhood of 1 and if the sequence an is
bounded, then
N
an
= O(1).
n=1
n
However, even this stronger assumption does not allow us to deduce that
N
an = o(N ),
n=1
5.2.1 Exercises
1. Let T be a regular matrix such that tmn ≥ 0 for all m, n. Show that if
limn→∞ an = +∞, then limm→∞ bm = +∞.
2. Show that if T = [tmn ] and U = [u mn ] are regular matrices, then so is
T U = V = [vmn ] where
∞
vmn = tmk u kn .
k=1
3. Show that if b = T a and limm→∞ bm = a whenever limn→∞ an = a, then
T is regular.
4. For n = 0, 1, 2, . . . let tn (x) be defined on [0, 1), and suppose that the tn
satisfy the following conditions:
(i) There is a constant C such that if x ∈ [0, 1), then ∞ n=0 |tn (x)| ≤ C.
(ii) For all n, limx→1− tn (x) = 0.
(iii) limx→1− ∞ n=0 tn (x) = 1.
Show that if limn→∞ an = a and if b(x) = ∞ n=0 an tn (x), then
limx→1− b(x) = a.
5. (Kojima 1917) Suppose that the numbers tmn satisfy the following
conditions:
(i) There is a constant C such that ∞ n=1 |tmn | ≤ C for all m.
(ii) For all n, limm→∞ tmn exists.
(iii) limm→∞ ∞ n=1 tmn exists.
Show that if limn→∞ an exists and if bm = ∞ n=1 tmn an , then limm→∞ bm
exists.
6. For positive
∞ integers n let K n (x) be a function defined on [0, ∞) such that
(i) 0 K n (x) d x → 1 as n → ∞;
∞
(ii) 0 |K n (x)| d x ≤ C for all n;
(iii) limn→∞ K n (x) = 0 uniformly for 0 ≤ x ≤ X . ∞
Suppose that a(x) is a bounded function, and that bn = 0 a(x)K n (x) d x.
Show that if limx→∞ a(x) = a, then limn→∞ bn = a.
7. Let rm be a sequence of positive real numbers with rm → 1− as m → ∞ .
For m ≥ 1, n ≥ 1, put tmn = nrmn−1 (1 − rm )2 .
(a) Show that [tmn ] is regular.
(b) Show that if an = n−1 k=0 ck (1 − k/n) and bm is defined by (5.32), then
∞
bm = k=0 ck rmk .
(c) Show that if cn = c (C, 1), then cn = c (A).
8. Suppose that T = [tmn ] is given by
⎧
⎪
⎪ 0 if n = 0,
⎨ m!n
tmn = if m ≥ n > 0,
⎪
⎪ m n+1 (m − n)!
⎩
0 if m < n.
158 Dirichlet series: II
for 1 ≤ k ≤ m .
(b) Verify that T is regular.
(c) Show that if an = nk=0 x k /k! for n ≥ 0, then bm = (1 + x/m)m for
m ≥ 1.
9. (Mercer’s theorem) Suppose that
1 1 a1 + a2 + · · · + am
bm = am + ·
2 2 m
for m ≥ 1. Show that
2n 2
n−1
an = bn − mbm .
n+1 n(n + 1) m=1
but
N
lim an (1 − n/N ) = 0.
N →∞
n=1
N
(c) Let B(N ) = n=1 nan . Show that if an converges, then B(N ) =
o(N ) as N → ∞.
(d) Show that if P(δ) converges for δ > 0, then
B(N ) N
1 e−u/N e−u/N
s N − P(1/N ) = + B(u) − − du
N 1 u2 u2 uN
∞
u du
+ B(u)e−u/N −1 .
N N u2
(e) Show that if B(N ) = o(N ), then s N − P(1/N ) = o(1).
(f) Show that if an = a (A), then an = a if and only if B(N ) = o(N ).
∞
21. (a) Using Ramanujan’s identity n=1 d(n)2 n −s = ζ (s)4 /ζ (2s) and Theo-
rem 5.11, show that n≤x d(n)2 /n ∼ (4π 2 )−1 (log x)4 .
(b) Show that if n≤x d(n)2 ∼ cx(log x)3 as x → ∞, then c = 1/π 2 .
22. Show that ∞ n=1 1/(d(n)n ) ∼ c(s − 1)
s −1/2
as s → 1+ where
p
c= ( p 2 − p)1/2 log .
p p−1
Deduce that
1 2c
∼ √ (log x)1/2
n≤x nd(n) π
as x → ∞.
23. Show that if n≤N an /n = O(1) and lims→1+ ∞ n=1 an n
−s
= a, then
an
log n
lim 1− = a.
n≤x n log x
x→∞
27. Suppose that for every ε > 0 there is an η > 0 such that
|an | < ε whenever N > 1/η. Show that if an = a (A),
N <n≤(1+η)N
then an = a.
28. Show that if an = a (C, 1) and if an+1 − an = O(|an |/n), then an = a.
29. (Hardy & Littlewood 1913, Theorem 27) Show that if an = a (A) and if
an+1 − an = O(|an |/n), then an = a.
30. (Hardy 1907) Show that
∞
k
lim− (−1)k x 2
x→1
k=0
5.3 Notes
Section 5.1. Theorem 5.1 and the more general (5.22) were first proved rig-
orously by Perron (1908). Although the Mellin transform had been used by
Riemann and Cahen, it was Mellin (1902) who first described a general class
of functions for which the inversion succeeds. Hjalmar Mellin was Finnish, but
his family name is of Swedish origin, so it is properly pronounced mĕ · lēn .
However, in English-speaking countries the uncultured pronunciation mĕl · ı̆n
is universal.
In connection with Theorem 5.4, it should be noted that Plancherel’s formula
f 2 = * f 2 holds not just for all f ∈ L 1 (R) ∩ L 2 (R) but actually for all
f ∈ L (R). However, in this wider setting one must adopt a new definition for
2
*f , since the definition we have taken is valid only for f ∈ L 1 (R). See Goldberg
(1961, pp. 46–47) for a resolution of this issue.
For further material concerning properties of Dirichlet series, one should
consult Hardy & Riesz (1915), Titchmarsh (1939, Chapter 9), or Widder (1971,
Chapter 2). Beyond the theory developed in these sources, we call attention to
two further topics of importance in number theory. Wiener (1932, p. 91) proved
that if the Fourier series of f ∈ L 1 (T) is absolutely convergent and is never zero,
then the Fourier series of 1/ f is also absolutely convergent. Wiener’s proof was
rather difficult, but Gel’fand (1941) devised a simpler proof depending on his
theory of normed rings. Lévy (1934) proved more generally that the Fourier
series of F( f ) is absolutely convergent provided that F is analytic at all points
in the range of f . Elementary proofs of these theorems have been given by
Zygmund (1968, pp. 245–246) and Newman (1975). These theorems were
generalized to absolutely convergent Dirichlet series by Hewitt & Williamson
(1957), who showed that if α(s) = an n −s is absolutely convergent for σ ≥
σ0 , then 1/α(s) is represented by an absolutely convergent Dirichlet series
5.3 Notes 163
in the same half-plane, if and only if the values taken by α(s) in this half-
plane are bounded away from 0. Ingham (1962) noted a fallacy in Zygmund’s
account of Lévy’s theorem, corrected it, and gave an elementary proof of the
generalization to absolutely convergent Dirichlet series. See also Goodman &
Newman (1984). Secondly, Bohr (1919) developed a theory concerning the
values taken on by an absolutely convergent Dirichlet series. This is described
by Titchmarsh (1986, Chapter 11), and in greater detail by Apostol (1976,
Chapter 8). For a small footnote to this theory, see Montgomery & Schinzel
(1977).
Section 5.2. That conditions (5.29)–(5.31) are necessary and sufficient for
the transformation T to preserve limits was proved by Toeplitz (1911) for upper
triangular matrices, and by Steinhaus (1911) in general. See also Kojima (1917)
and Schur (1921). For more on the Toeplitz matrix theorem and various aspects
of Tauberian theorems, see Peyerimhoff (1969).
Theorem 5.6 under the hypothesis (a) is trivial by dominated convergence.
Theorem 5.6(b) is a special case of a theorem of Hardy (1910), who considered
the more general (C,k) convergence, and Theorem 5.6(c) is similarly a special
case of a theorem of Landau (1910, pp. 103–113).
Tauber (1897) proved two theorems, the second of which is found in Exer-
cise 5.2.18. Littlewood (1911) derived his strengthening of Tauber’s first theo-
rem by using high-order derivatives. Subsequently Hardy & Littlewood (1913,
1914a, b, 1926, 1930) used the same technique to obtain Theorem 5.8 and
its corollaries. Karamata (1930, 1931a, b) introduced the use of Weierstrass’s
approximation theorem. Karamata also considered a more general situation,
in which the right-hand sides of (5.35) and (5.36) are multiplied by a slowly
oscillating function L(1/δ), and the right-hand side of (5.37) is multiplied by
L(U ). Our exposition employs a further simplification due to Wielandt (1952).
Other proofs of Littlewood’s theorem have been given by Delange (1952) and
by Eggleston (1951). Ingham (1965) observed that a peak function similar
to Littlewood’s can be constructed by using high-order differencing instead
of differentiation. Since many proofs of the Weierstrass theorem involve con-
structing a peak function, the two methods are not materially different. Sharp
quantitative Tauberian theorems have been given by Postnikov (1951), Kore-
vaar (1951, 1953, 1954a–d), Freud (1952, 1953, 1954), Ingham (1965), and
Ganelius (1971).
For other accounts of the Hardy–Littlewood theorem, see Hardy (1949) or
Widder (1946, 1971). For a brief survey of applications of summability to
classical analysis, see Rubel (1989).
Wiener (1932, 1933) invented a general Tauberian theory that contains the
Hardy–Littlewood theorems for power series (Theorem 5.8 and its corollaries)
164 Dirichlet series: II
as a special case. Wiener’s theory is discussed by Hardy (1949), Pitt (1958), and
Widder (1946). Among the longer expositions of Tauberian theory, the recent
accounts of Korevaar (2002, 2004) are especially recommended.
5.4 References
Apostol, T. (1976). Modular Functions and Dirichlet Series in Number Theory, Graduate
Texts Math. 41. New York: Springer-Verlag.
Bohr, H. (1909). Über die Summabilität Dirichletscher Reihen, Nachr. König. Gesell.
Wiss. Göttingen Math.-Phys. Kl., 247–262; Collected Mathematical Works, Vol. I.
København: Dansk Mat. Forening, 1952, A2.
(1919). Zur Theorie algemeinen Dirichletschen Reihen, Math. Ann. 79, 136–156;
Collected Mathematical Works, Vol. I. København: Dansk Mat. Forening, 1952,
A13.
Delange, H. (1952). Encore une nouvelle démonstration du théorème taubérien de Lit-
tlewood, Bull. Sci. Math. (2) 76, 179–189.
Edwards, D. A. (1957). On absolutely convergent Dirichlet series, Proc. Amer. Math.
Soc. 8, 1067–1074.
Eggleston, H. G. (1951). A Tauberian lemma, Proc. London Math. Soc. (3) 1, 28–45.
Freud, G. (1952). Restglied eines Tauberschen Satzes, I, Acta Math. Acad. Sci. Hungar.
2, 299–308.
(1953). Restglied eines Tauberschen Satzes, II, Acta Math. Acad. Sci. Hungar. 3,
299–307.
(1954). Restglied eines Tauberschen Satzes, III, Acta Math. Acad. Sci. Hungar. 5,
275–289.
Ganelius, T. (1971). Tauberian Remainder Theorems, Lecture Notes Math. 232. Berlin:
Springer-Verlag.
Gel’fand, I. M. (1941). Über absolut konvergente trigonometrische Reihen und Integrale,
Mat. Sb. N. S. 9, 51–66.
Goldberg R. R. (1961). Fourier Transforms, Cambridge Tract 52. Cambridge: Cambridge
University Press.
Goodman, A. & Newman, D. J. (1984). A Wiener type theorem for Dirichlet series,
Proc. Amer. Math. Soc. 92, 521–527.
Hardy, G. H. (1907). On certain oscillating series, Quart. J. Math. 38, 269–288; Collected
Papers, Vol. 6. Oxford: Clarendon Press, 1974, pp. 146–167.
(1910). Theorems relating to the summability and convergence of slowly oscillating
series, Proc. London Math. Soc. (2) 8, 301–320; Collected Papers, Vol. 6. Oxford:
Clarendon Press, 1974, pp. 291–310.
(1949). Divergent Series, Oxford: Oxford University Press.
Hardy, G. H. & Littlewood, J. E. (1913). Contributions to the arithmetic theory of
series, Proc. London Math. Soc. (2) 11, 411–478; Collected Papers, Vol. 6. Oxford:
Clarendon Press, 1974, pp. 428–495.
(1914a). Tauberian theorems concerning power series and Dirichlet series whose co-
efficients are positive, Proc. London Math. Soc. (2) 13, 174–191; Collected Papers,
Vol. 6. Oxford: Clarendon Press, 1974, pp. 510–527.
5.4 References 165
(1914b). Some theorems concerning Dirichlet’s series, Messenger Math. 43, 134–147;
Collected Papers, Vol. 6. Oxford: Clarendon Press, 1974, pp. 542–555.
(1926). A further note on the converse of Abel’s theorem, Proc. London Math.
Soc. (2) 25, 219–236; Collected Papers, Vol. 6. Oxford: Clarendon Press, 1974,
pp. 699–716.
(1930). Notes on the theory of series XI: On Tauberian theorems, Proc. London
Math. Soc. (2) 30, 23–37; Collected Papers, Vol. 6. Oxford: Clarendon Press, 1974,
pp. 745–759.
Hardy, G. H. & Riesz, M. (1915). The General Theory of Dirichlet’s Series, Cambridge
Tract No. 18. Cambridge: Cambridge University Press. Reprint: Stechert–Hafner
(1964).
Hewitt, E. & Williamson, H. (1957). Note on absolutely convergent Dirichlet series,
Proc. Amer. Math. Soc. 8, 863–868.
Ingham, A. E. (1962). On absolutely convergent Dirichlet series. Studies in Mathemati-
cal Analysis and Related Topics. Stanford: Stanford University Press, pp. 156–164.
(1965). On tauberian theorems, Proc. London Math. Soc. (3) 14A, 157–173.
Karamata, J. (1930). Über die Hardy–Littlewoodschen Umkehrungen des Abelschen
Stetigkeitssatzes, Math. Z. 32, 319–320.
(1931a). Neuer Beweis und Verallgemeinerung einiger Tauberian-Sätze, Math. Z. 33,
294–300.
(1931b). Neuer Beweis und Verallgemeinerung der Tauberschen Sätze, welche die
Laplacesche und Stieltjessche Transformation betreffen, J. Reine Angew. Math.
164, 27–40.
Kojima, T. (1917). On generalized Toeplitz’s theorems on limit and their application,
Tôhoku Math. J. 12, 291–326.
Korevaar, J. (1951). An estimate of the error in Tauberian theorems for power series,
Duke Math. J. 18, 723–734.
(1953). Best L 1 approximation and the remainder in Littlewood’s theorem, Proc.
Nederl. Akad. Wetensch. Ser. A 56 (= Indagationes Math. 15), 281–293.
(1954a). A very general form of Littlewood’s theorem, Proc. Nederl. Akad. Wetensch.
Ser. A 57 (= Indagationes Math. 16), 36–45.
(1954b). Another numerical Tauberian theorem for power series, Proc. Nederl. Akad.
Wetensch. Ser. A 57 (= Indagationes Math. 16), 46–56.
(1954c). Numerical Tauberian theorems for Dirichlet and Lambert series, Proc.
Nederl. Akad. Wetensch. Ser. A 57 (= Indagationes Math. 16), 152–160.
(1954d). Numerical Tauberian theorems for power series and Dirichlet series, I, II,
Proc. Nederl. Akad. Wetensch. Ser. A 57 (= Indagationes Math. 16), 432–443,
444–455.
(2001). Tauberian theory, approximation, and lacunary series of powers, Trends in
approximation theory (Nashville, 2000), Innov. Appl. Math. Nashville: Vanderbilt
University Press, pp. 169–189.
(2002). A century of complex Tauberian theory, Bull. Amer. Math. Soc. (N.S.) 39,
475–531.
(2004). Tauberian Theory. A Century of Developments. Grundl. Math. Wiss. 329.
Berlin: Springer-Verlag.
Landau, E. (1907). Über die Multiplikation Dirichletscher Reihen, Rend. Circ. Mat.
Palermo 24, 81–160.
166 Dirichlet series: II
(1908). Zwei neue Herleitungen für die asymptotische Anzahl der Primzahlen unter
einer gegebenen Grenze, Sitzungsberichte Akad. Wiss. Berlin 746–764; Collected
Works, Vol.4. Essen: Thales Verlag, 1986, pp. 21–39.
(1909). Handbuch der Lehre von der Verteilung der Primzahlen, Leipzig: Teubner.
Reprint: Chelsea (New York), 1953.
(1910). Über die Bedeutung einiger neuerer Grenzwertsätze der Herren Hardy und
Axer, Prace mat.-fiz. (Warsaw) 21, 97–177; Collected Works, Vol. 4. Essen: Thales
Verlag, 1986, pp. 267–347.
(1913). Einige Ungleichungen für zweimal differentiierbare Funktionen, Proc. Lon-
don Math. Soc. (2) 13, 43–49; Collected Works, Vol. 6. Essen: Thales Verlag, 1986,
pp. 49–55.
Lévy, P. (1934). Sur la convergence absolue des séries de Fourier, Compositio Math. 1,
1–14.
Littlewood, J. E. (1911). The converse of Abel’s theorem on power series, Proc. London
Math. Soc. (2) 9, 434–448; Collected Papers, Vol. 1. Oxford: Oxford University
Press, 1982, pp. 757–773.
(1986). Littlewood’s Miscellany, Bollobas, B. Ed., Cambridge: Cambridge University
Press.
van de Lune, J. (1986). An Introduction to Tauberian Theory: From Tauber to Wiener.
CWI Syllabus 12. Amsterdam: Mathematisch Centrum.
Mellin, H. (1902). Über den Zusammenhang zwischen den linearen Differential- und
Differenzengleichungen, Acta Math. 25, 139–164.
Montgomery, H. L. & Schinzel, A. (1977). Some arithmetic properties of polynomials in
several variables. Transcendence Theory: Advances and Applications (Cambridge,
1976). London: Academic Press, pp. 195–203.
Newman, D. J. (1975). A simple proof of Wiener’s 1/ f theorem, Proc. Amer. Math. Soc.
48, 264–265.
Perron, O. (1908). Zur Theorie der Dirichletschen Reihen, J. Reine Angew. Math. 134,
95–143.
Peyerimhoff, A. (1969). Lectures on summability, Lecture Notes Math. 107. Berlin:
Springer-Verlag.
Pitt, H. R. (1958). Tauberian Theorems. Tata Monographs. London: Oxford University
Press.
Postnikov, A. G. (1951). The remainder term in the Tauberian theorem of Hardy and
Littlewood, Dokl. Akad. Nauk SSSR N. S. 77, 193–196.
Riesz, M. (1909). Sur la sommation des séries de Dirichlet, C. R. Acad. Sci. Paris 149,
18–21.
Rubel, L. (1989). Summability theory: a neglected tool of analysis, Amer. Math. Monthly
96, 421–423.
Schoenberg, I. J. (1973). The elementary cases of Landau’s problem of inequalities
between derivatives, Amer. Math. Monthly 80, 121–158.
Schur, I. (1921). Über lineare Transformationen in der Theorie der unendlichen Reihen,
J. Reine Angew. Math. 151, 79–111.
Steinhaus, H. (1911). Kilka slów o uogólnieniu pojȩcia granicy, Warsaw: Prace mat-fiz
22, 121–134.
Tauber, A. (1897). Ein Satz aus der Theorie der unendlichen Reihen, Monat. Math. 8,
273–277.
5.4 References 167
168
6.1 A zero-free region 169
put
K
R 2 − zz k
g(z) = f (z) .
k=1
R(z − z k )
th
The k factor of the product has been constructed so that it has a pole at z k , and
so that it has modulus 1 on the circle |z| = R. Hence g is an analytic function
in the disc |z| ≤ R, and if |z| = R, then |g(z)| = | f (z)| ≤ M. Hence by the
maximum modulus principle, |g(0)| ≤ M. But
K
R
|g(0)| = | f (0)| .
|z |
k=1 k
We now show that a bound for the modulus of an analytic function can be
derived from a one-sided bound for its real part in a slightly larger region.
Lemma 6.2 (The Borel–Carathéodory Lemma) Suppose that h(z) is analytic
in a domain containing the disc |z| ≤ R, that h(0) = 0, and that h(z) ≤ M
for |z| ≤ R. If |z| ≤ r < R, then
2Mr
|h(z)| ≤
R −r
and
2M R
|h (z)| ≤ .
(R − r )2
Proof It suffices to show that
h (k) (0) 2M
≤ k (6.1)
k! R
for all k ≥ 1, for then
∞
h (k) (0) k ∞
r k 2Mr
|h(z)| ≤ r ≤ 2M = ,
k=1
k! k=1
R R −r
and
∞
|h (k) (0)|kr k−1 2M ∞
r k−1 2M R
|h (z)| ≤ ≤ k = .
k=1
k! R k=1 R (R − r )2
To prove (6.1) we first note that
1
1 dz
h(Re(θ )) dθ = h(z) = h(0) = 0.
0 2πi |z|=R z
170 The Prime Number Theorem
and
1
Rk R k h (k) (0)
h(Re(θ ))e(−kθ) dθ = h(z)z −k−1 dz = .
0 2πi |z|=R k!
By forming a linear combination of these identities we see that if k > 0, then
1
R k e(−φ)h (k) (0)
h(Re(θ ))(1 + cos 2π(kθ + φ)) dθ = .
0 2 · k!
By taking real parts it follows that
1
1 k
R e(−φ)h (k) (0)/k! ≤ M (1 + cos 2π (kθ + φ)) dθ = M
2 0
for k > 0. Since this holds for any real φ, we are free to choose φ so that
e(−φ)h (k) (0) = |h (k) (0)|. Then the above inequality gives (6.1), and the proof
is complete.
K
If P(z) = c k=1 (z − z k ), then
P K
1
(z) = .
P k=1
z − zk
We now generalize this to analytic functions f (z), to the extent that f / f can
be approximated by a sum over its nearby zeros.
Lemma 6.3 Suppose that f (z) is analytic in a domain containing the disc
|z| ≤ 1, that | f (z)| ≤ M in this disc, and that f (0) = 0. Let r and R be fixed,
0 < r < R < 1. Then for |z| ≤ r we have
f K
1 M
(z) = + O log
f k=1
z − zk | f (0)|
where the sum is extended over all zeros z k of f for which |z k | ≤ R. (The implicit
constant depends on r and R, but is otherwise absolute.)
Proof If f (z) has zeros on the circle |z| = R, then we replace R by a very
slightly larger value. Thus we may assume that f (z) = 0 for |z| = R. Set
K
R 2 − zz k
g(z) = f (z) .
k=1
R(z − z k )
6.1 A zero-free region 171
where τ = |t| + 4 and the sum is extended over all zeros ρ of ζ (s) for which
|ρ − (3/2 + it)| ≤ 5/6.
Proof We apply Lemma 6.3 to the function f (z) = ζ (z + (3/2 + it)), with
R = 5/6 and r = 2/3. To complete the proof it suffices to note that | f (0)| 1
by the (absolutely convergent) Euler product formula (1.17), and that f (z) τ
for |z| ≤ 1 by Corollary 1.17.
172 The Prime Number Theorem
We now use Lemmas 6.4 and 6.5 to establish the existence of a zero-free
region for the zeta function.
Theorem 6.6 There is an absolute constant c > 0 such that ζ (s) = 0 for
σ ≥ 1 − c/ log τ .
Proof Since ζ (s) is given by the absolutely convergent product (1.17) for
σ > 1, it suffices to consider σ ≤ 1. From (1.24) we see that
∞
s |s|
ζ (s) − ≤ |s| u −σ −1 du = (6.5)
s−1 1 σ
for σ > 0. From this we see that ζ (s) = 0 when σ > |s − 1|, i.e., in the parabolic
region σ > (1 + t 2 )/2. In particular, ζ (s) = 0 in the rectangle 8/9 ≤ σ ≤ 1,
|t| ≤ 7/8. Now suppose that ρ0 = β0 + iγ0 is a zero of the zeta function with
5/6 ≤ β0 ≤ 1, |γ0 | ≥ 7/8. Since ρ ≤ 1 for all zeros ρ of ζ (s), it follows that
1/(s − ρ) > 0 whenever σ > 1. Hence by Lemma 6.4 with s = 1 + δ + iγ0
we see that
ζ 1
− (1 + δ + iγ0 ) ≤ − + c1 log(|γ0 | + 4).
ζ 1 + δ − β0
Similarly, by Lemma 6.4 with s = 1 + δ + 2iγ0 we find that
ζ
− (1 + δ + 2iγ0 ) ≤ c1 log(|2γ0 | + 4).
ζ
From Corollary 1.13 we see that
ζ 1
− (1 + δ) = + O(1).
ζ δ
On combining these estimates in Lemma 6.5 we conclude that
3 4
− + c2 log(|γ0 | + 4) ≥ 0.
δ 1 + δ − β0
We take δ = 1/(2c2 log(|γ0 | + 4)). Thus the above gives
4
7c2 log(|γ0 | + 4) ≥ ,
1 + δ − β0
which is to say that
1 4
1+ − β0 ≥ .
2c2 log(|γ0 | + 4) 7c2 log(|γ0 | + 4)
Hence
1
1 − β0 ≥ ,
14c2 log(|γ0 | + 4)
so the proof is complete.
It is useful to have bounds for the zeta function and its logarithmic derivative
in the zero-free region.
Theorem 6.7 Let c be the constant in Theorem 6.6. If σ > 1 − c/(2 log τ )
and |t| ≥ 7/8, then
ζ
(s) log τ , (6.6)
ζ
| log ζ (s)| ≤ log log τ + O(1) , (6.7)
and
1
log τ . (6.8)
ζ (s)
ζ
On the other hand, if 1 − c/(2 log τ ) < σ ≤ 2 and |t| ≤ 7/8, then (s) =
ζ
−1/(s − 1) + O(1), log ζ (s)(s − 1) 1, and 1/ζ (s) |s − 1|.
Proof If σ > 1, then by Corollary 1.11 and the triangle inequality we see that
ζ ∞
ζ 1
(s) ≤ (n)n −σ = − (σ ) .
ζ n=1
ζ σ −1
Hence (6.6) is obvious if σ ≥ 1 + 1/ log τ . Let s1 = 1 + 1/ log τ + it. In par-
ticular we have
ζ
(s1 ) log τ. (6.9)
ζ
From this estimate and Lemma 6.4 we deduce that
1
log τ (6.10)
ρ s1 − ρ
where the sum is over those zeros ρ for which |ρ − (3/2 + it)| ≤ 5/6. Suppose
that 1 − c/(2 log τ ) ≤ σ ≤ 1 + 1/ log τ . Then by Lemma 6.4 we see that
ζ ζ 1 1
(s) − (s1 ) = − + O(log τ ). (6.11)
ζ ζ ρ s−ρ s1 − ρ
6.1 A zero-free region 175
But by Theorem 1.14 we know that ζ (σ ) < 1 + 1/(σ − 1), so that (6.7)
holds when σ ≥ 1 + 1/ log τ . In particular (6.7) holds at the point s1 =
1 + 1/ log τ + it, so that to treat the remaining s it suffices to bound the
difference
s
ζ
log ζ (s) − log ζ (s1 ) = (w) dw.
s1 ζ
We take the path of integration to be the line segment joining the endpoints.
Then the length of this interval multiplied by the bound (6.6) gives the error
term O(1) in (6.7).
The estimate (6.8) follows directly from (6.7), since log 1/|ζ | = − log ζ .
The remaining estimates follow trivially from (6.5).
The ideas we have used enable us not only to derive a zero-free region but
also to place a bound on the number of zeros ρ that might lie near the point
1 + it.
Theorem 6.8 Let n(r ; t) denote the number of zeros ρ of ζ (s) in the disc
|ρ − (1 + it)| ≤ r . Then n(r ; t) r log τ , uniformly for r ≤ 3/4.
6.1.1 Exercises
1. (a) Show that if |z| < R, |w| ≤ R, and z = w, then
zw − R 2
≥ 1.
(z − w)R
(b) Show that if |w| ≤ ρ < R, |z| = r < R, and z = w, then
zw − R 2 rρ + R 2
≥ .
(z − w)R (r + ρ)R
(c) Suppose that f is analytic in the disc |z| ≤ R. For r ≤ R put M(r ) =
max|z|≤r | f (z)|. Show that if 0 < r < R and 0 < ρ < R, then the num-
ber of zeros of f in the disc |z| ≤ ρ does not exceed
M(R)
log
M(r )
.
rρ + R 2
log (r + ρ)R
2. Suppose that R, M, and ε are positive real numbers, and set h(z) =
2M z/(z + R + ε).
(a) Show that h(0) = 0, that h(z) is analytic for |z| < R + ε, and that
h(z) ≤ M for |z| ≤ R + ε.
(b) Show that if 0 < r < R, then
2Mr
max |h(z)| = −h(−r ) = .
|z|≤r R +ε−r
(c) Show that if 0 < r < R, then
2M(R + ε)
max |h (z)| = h (−r ) = .
|z|≤r (R + ε − r )2
3. Show that, in the situation of the Borel–Carathéodory lemma (Lemma 6.2),
if |z| ≤ r < R, then
4M R
|h (z)| ≤ .
(R − r )3
4. (Mertens 1898) Use the Dirichlet series expansion of log ζ (s) to show that
if σ > 1, then
|ζ (σ )3 ζ (σ + it)4 ζ (σ + 2it)| ≥ 1.
The method used to establish a zero-free region for the zeta function can be
applied to any particular Dirichlet L-function, though the constants involved
may depend on the function. We shall pursue this systematically in Chapter 11,
but in the exercise below we treat one interesting example.
6.1 A zero-free region 177
5. Let χ0 denote the principal character (mod 4), and χ1 the non-principal
character (mod 4).
(a) Show that L(1, χ1 ) = π/4, and hence that there is a neighbourhood of
1 in which L(s, χ1 ) = 0.
(b) Show that if σ > 1, then
L L L
−3 (σ, χ0 ) − 4 (σ + it, χ1 ) − (σ + 2it, χ0 ) ≥ 0.
L L L
(c) Show that there is a constant c > 0 such that L(s, χ1 ) = 0 for σ >
1 − c/ log τ .
(d) Show that there is a constant c > 0 such that if σ > 1 − c/ log τ , then
L
(s, χ1 ) log τ,
L
| log L(s, χ1 )| ≤ log log τ + O(1),
1
log τ.
L(s, χ1 )
6. (a) Show that if 1 < σ1 ≤ σ2 , then
ζ (σ2 ) ζ (σ2 + it) ζ (σ1 )
≤ ≤
ζ (σ1 ) ζ (σ1 + it) ζ (σ2 )
for all real t.
(b) Show that if 1 < σ1 ≤ σ2 ≤ 2, then
σ1 − 1 ζ (σ2 + it) σ2 − 1
σ2 − 1 ζ (σ1 + it) σ1 − 1
uniformly in t.
7. (Montgomery & Vaughan 2001)
(a) Show that if σ > 1, then
ζ (σ + i(t + 1)) ∞
(n) 1
≤ exp 2 sin 2 log n
ζ (σ + it) n=1
n σ log n
for σ ≥ 1 − θ(t + 1)/3 where the sum is over zeros ρ for which |ρ −
(1 + θ (t + 1) + it)| ≤ 5θ (t + 1)/3.
(b) Show that there is an absolute constant c > 0 such that ζ (s) = 0 for
θ (2t + 1)
σ ≥1−c .
φ(2t + 1)
(c) Show that the zero-free region (6.26) follows from the estimate (6.25).
6.2 The Prime Number Theorem 179
and then use partial summation to derive an estimate for π (x). It would be more
direct to apply Perron’s formula to log ζ (s), but our approach is technically
simpler since log ζ (s) has a logarithmic singularity at s = 1 while ζζ (s) has
only a simple pole there.
since we may suppose that 0 < c < 1. Thus the proof of (6.12) is complete.
To derive (6.13) it suffices to combine (6.12) with the first estimate of Corol-
lary 2.5. As for (6.14), we note that
x x
1 1
π(x) = dϑ(u) = li(x) + d(ϑ(u) − u).
2− log u 2− log u
By integrating by parts we see that this last integral is
ϑ(u) − u x x
ϑ(u) − u
+ 2
du,
log u 2− 2 u(log u)
√
and by (6.13) it follows that this is x exp(−c log x). Thus we have (6.14),
and the proof is complete.
182 The Prime Number Theorem
The method we used to derive Theorem 6.9 is very flexible, and can be
applied to many other situations. For example, the summatory function
M(x) = µ(n)
n≤x
6.2.1 Exercises
1. (Landau 1901b; cf. Rosser & Schoenfeld 1962) Use Theorem 6.9 to show
that
π(2x) − 2π (x) = −2(log 2)x(log x)−2 + O(x(log x)−3 ).
Deduce that for all large x, the interval (x, 2x] contains fewer prime num-
bers than the interval (0, x].
2. Use Theorem 6.9 to show that if n is of the form n = p≤y p where y is
sufficiently large, then d(n) > n (log 2)/ log log n .
3. (a) Use Theorem 6.9 to show that
1 log y
= log + O exp − c log x .
x< p≤y p log x
(b) Use the above and Theorem 2.7 to show that
1
= log log x + b + O exp − c log x
p≤x p
where b = C0 − p ∞ k
k=2 1/(kp ) .
4. Show that for x ≥ 2,
(n)
= log x − C0 + O exp − c log x .
n≤x n
6.2 The Prime Number Theorem 183
5. (cf. Cipolla 1902; Rosser 1939) Let p1 < p2 < · · · denote the prime num-
bers. Show that
log log n 2 (log log n)2
pn = n log n + log log n − 1 + − +O .
log n log n (log n)2
6. (Landau 1900) Let πk (x) denote the number of integers not exceeding x
that are composed of exactly k distinct primes.
(a) Show that
π2 (x) = π(x/ p) + O x(log x)−2 .
√
p≤ x
(c) Using Theorem 6.9 and integration by parts, show that the sum above
is
√
x
du
x + O(x/ log x).
2 u(log x/u) log u
(d) Conclude that π2 (x) = x(log log x)/ log x + O(x/ log x).
7. (D. E. Knutson) Let dn denote the least common multiple of the numbers
1, 2, . . . , n.
(a) Show that dn = exp(ψ(n)).
(b) Let E(z) = ∞ n=1 z /dn . Show that this power series has radius of
n
convergence e.
(c) Show that E(1) is irrational.
8. (Landau 1905) Let Q(x) denote the number of square-free integers not
exceeding x, and define R(x) by the relation Q(x) = (6/π 2 )x + R(x).
(a) Show that
R(x) = M(y){x/y 2 } − µ(d){x/d 2 }
d≤y
∞
+ M x/m − 2x M(u)u −3 du.
m≤x/y 2 y
√
(b) Taking y = x 1/2 exp(−c log x) where c is sufficiently small, show
√
that R(x) x 1/2 exp(−c log x).
9. Let N = N (Q) = 1 + q≤Q ϕ(q) be the number of Farey points of order
Q, and for 0 ≤ α ≤ 1 write
for 0 ≤ α ≤ 1.
(d) Show that R Q uniformly for 0 ≤ α ≤ 1.
10. (Landau 1903b; Massias, Nicolas & Robin 1988, 1989) Let f (n) denote
the maximal order of any element of the symmetric group Sn .
(a) Show that f (n) = max lcm(n 1 , n 2 , . . . , n k ) where the maximum is ex-
tended over all sets {n 1 , n 2 , . . . , n k ) of natural numbers for which
n 1 + n 2 + · · · + n k ≤ n.
(b) Choose y as large as possible so that p≤y p ≤ n. Show that
log f (n) ≥ log p = (1 + o(1))(n log n)1/2 .
p≤y
16. (Landau 1899b, 1901a, 1903c) Use the method of proof of Theorem 6.9 to
show that
∞
µ(n) log n
(a) = −1;
n=1
n
∞
µ(n)(log n)2
(b) = −2C0 ;
n=1
n
∞
λ(n) log n
(c) = −ζ (2).
n=1
n
17. Taking (6.18) and a quantitative form of the first part of the preceding
exercise for granted, use elementary reasoning to show that if q ≤ x then
µ(n)
(a) exp − c log x ,
n≤x n
(n,q)=1
µ(n) log n q
(b) =− + O exp − c log x .
n≤x n ϕ(q)
(n,q)=1
18. (Hardy 1921) Use the method of proof of Theorem 6.9 to show that
∞
µ(n)
(a) = 0;
n=1
ϕ(n)
∞
µ(n) log n
(b) = 0;
n=1
ϕ(n)
186 The Prime Number Theorem
∞
µ(n)(log n)2
(c) = 4A log 2
n=1
ϕ(n)
where A = p>2 1 − ( p−1)
1
2 .
19. Let Q(x) denote the number of square-free integers not exceeding x, and
recall Theorem 2.2.
(a) Show that
6 µ(n)
Q(x) = x−x − µ(n){x/n 2 }
π 2 √
n> x
n 2 √
n≤ x
as x → ∞.
6.2 The Prime Number Theorem 187
(d) Use the estimate d≤y µ(d)/d (log 2y)−2 to show that
x µ(d) dv
1.
1 d≤x/v
d v
(e) Mimic the proof of Theorem 5.5, or use Exercise 5.2.6 to show that if
(i) holds, then
f (n) n
lim 1− = c.
x→∞
n≤x n x
(f) Use Theorem 5.6 to show that if (i) holds and f (n) = O(1), then (ii)
follows.
(g) Take f (n) = µ(n) to deduce that ∞ n=1 µ(n)/n = 0. (Of course we
used much more above in (d). For a result in the converse direction, see
Exercise 8.1.5.)
21. (Landau 1908b) Let R be the set of positive integers that can be expressed
as a sum of two squares, let R(x) denote the number of such integers not
exceeding x, and let χ1 denote the non-principal character (mod 4), as in
Exercise 6.1.5.
(a) Show that
n −s = (1 − 2−s )−1 (1 − p −s )−1 (1 − p −2s )−1
n∈R p≡1 (4) p≡3 (4)
for σ > 1.
√
(b) Show that the Dirichlet series above is f (s) ζ (s)L(s, χ1 ) where
f (s) = (1 − 2−s )−1/2 (1 − p −2s )−1/2
p≡3 (4)
where
f (s)
g(s) = (s − 1)ζ (s)L(s, χ1 )
s
is analytic in a neighbourhood of 1.
(f) Show that
,
π
g(1) = (1 − p −2 )−1/2 .
2 p≡3 (4)
b = 2−1/2 (1 − p −2 )−1/2 .
p≡3 (4)
22. Let A denote the set of those positive integers that are composed entirely
of the prime 2 and primes ≡ 1 (mod 4), and let B be the the set of those
positive integers that are composed entirely of primes ≡ 3 (mod 4).
(a) Explain why any positive integer n has a unique representation in the
form n = a(n)b(n) where a(n) ∈ A and b(n) ∈ B.
(b) Let A(x) denote the number of a ∈ A, a ≤ x. Show that
αx x
A(x) = √ +O
log x (log x)3/2
√
where α = 1/ 2.
(c) Let B(x) denote the number of b ∈ B, b ≤ x. Show that
βx x
B(x) = √ +O
log x (log x)3/2
√
where β = 2/π .
6.2 The Prime Number Theorem 189
(d) For 0 ≤ κ ≤ 1 let Nκ (x) denote the number of n ≤ x such that a(n) ≤
n κ . Show that
Nκ (x) = 1.
a≤x κ a 1/κ−1 ≤b≤x/a
a∈A b∈B