Product rule: Difference between revisions
en>InverseHypercube m →A Brief Proof: uncap. section |
en>Monkbot |
||
Line 1: | Line 1: | ||
[[File:Sum of reciprocals of primes.svg|thumb|300px|The sum of the reciprocal of the primes increasing without bound. The x axis is in log scale, showing that the divergence is very slow. The purple function is a lower bound that also diverges.]] | |||
The '''sum of the [[multiplicative inverse|reciprocal]]s of all [[prime number]]s [[Divergent series|diverges]]'''; that is: | |||
:<math>\sum_{p\text{ prime }}\frac1p = \frac12 + \frac13 + \frac15 + \frac17 + \frac1{11} + \frac1{13} + \frac1{17} + \cdots = \infty</math> | |||
This was proved by [[Leonhard Euler]] in 1737, and strengthens [[Euclid]]'s 3rd-century-BC result that [[Euclid's theorem|there are infinitely many prime number]]s. | |||
There are a variety of proofs of Euler's result, including a lower bound for the partial sums stating that | |||
:<math>\sum_{\scriptstyle p\text{ prime }\atop \scriptstyle p\le n}\frac1p \ge \log \log (n+1) - \log\frac{\pi^2}6</math> | |||
for all natural numbers ''n''. The double [[natural logarithm]] indicates that the divergence might be very slow, which is indeed the case, see [[Meissel–Mertens constant]]. | |||
==The harmonic series== | |||
First, we describe how Euler originally discovered the result. He was considering the [[harmonic series (mathematics)|harmonic series]] | |||
: <math> | |||
\sum_{n=1}^\infty \frac{1}{n} = | |||
1 + \frac{1}{2} + \frac{1}{3} + \frac{1}{4} + \cdots | |||
</math> | |||
He had already used the following "[[Riemann zeta function#Euler product formula|product formula]]" to show the existence of infinitely many primes. | |||
: <math> | |||
\sum_{n=1}^\infty \frac{1}{n} = \prod_{p} \frac{1}{1-p^{-1}} | |||
= \prod_{p} \left( 1+\frac{1}{p}+\frac{1}{p^2}+\cdots \right) | |||
</math> | |||
(Here, the product is taken over all primes ''p''; in the following, a sum or product taken over ''p'' always represents a sum or product taken over a specified set of primes, unless noted otherwise.) | |||
Such infinite products are today called [[Euler product]]s. The product above is a reflection of the [[fundamental theorem of arithmetic]]. Of course, the above "equation" is not necessary because the harmonic series is known (by other means) to diverge. This type of formal manipulation was common at the time, when mathematicians were still experimenting with the new tools of [[calculus]].{{Citation needed|date=October 2010}} | |||
Euler noted that if there were only a finite number of primes, then the product on the right would clearly converge, contradicting the divergence of the harmonic series. (In modern language, we now say that the existence of infinitely many primes is reflected by the fact that the [[Riemann zeta function]] has a [[pole (complex analysis)|simple pole]] at ''s'' = 1.) | |||
==Proofs== | |||
===First=== | |||
Euler took the above product formula and proceeded to make a sequence of audacious leaps of logic. First, he took the natural logarithm of each side, then he used the Taylor series expansion for ln(''x'') as well as the sum of a geometric series: | |||
: <math> | |||
\begin{align} | |||
\ln \left( \sum_{n=1}^\infty \frac{1}{n}\right) & {} = \ln\left( \prod_p \frac{1}{1-p^{-1}}\right) | |||
= -\sum_p \ln \left( 1-\frac{1}{p}\right) \\ | |||
& {} = \sum_p \left( \frac{1}{p} + \frac{1}{2p^2} + \frac{1}{3p^3} + \cdots \right) \\ | |||
& {} = \left( \sum_{p}\frac{1}{p} \right) + \sum_p \frac{1}{p^2} \left( \frac{1}{2} + \frac{1}{3p} + \frac{1}{4p^2} + \cdots \right) \\ | |||
& {} < \left( \sum_p \frac{1}{p} \right) + \sum_p \frac{1}{p^2} \left( 1 + \frac{1}{p} + \frac{1}{p^2} + \cdots \right) \\ | |||
& {} = \left( \sum_p \frac{1}{p} \right) + \left( \sum_p \frac{1}{p(p-1)} \right) \\ | |||
& {} = \left( \sum_p \frac{1}{p} \right) + C | |||
\end{align} | |||
</math> | |||
for a fixed constant ''C'' < 1. Since the sum of the reciprocals of the first ''n'' positive integers is asymptotic to ln(''n''), (i.e. their ratio approaches one as ''n'' approaches infinity), Euler then concluded | |||
:<math>\frac{1}{2} + \frac{1}{3} + \frac{1}{5} + \frac{1}{7} + \frac{1}{11} + \cdots = \ln \ln (+ \infty)</math> | |||
It is almost certain that Euler meant that the sum of the reciprocals of the primes less than ''n'' is asymptotic to ln(ln(''n'')) as ''n'' approaches infinity. It turns out this is indeed the case; Euler had reached a correct result by questionable means.{{Citation needed|date=January 2012}} | |||
====A variation==== | |||
: <math> | |||
\begin{align} | |||
& {} \quad \log \left( \sum_{n=1}^\infty \frac{1}{n}\right) = \log \left( \prod_p \frac{1}{1-p^{-1}}\right) = \sum_p \log \left( \frac{p}{p-1}\right) = \sum_p \log\left(1+\frac{1}{p-1}\right) | |||
\end{align} | |||
</math> | |||
Since | |||
: <math> e^x = 1 + x + \frac{x^2}{2!} + \frac{x^3}{3!} + \cdots</math> | |||
Shows that <math>\scriptstyle e^x \,>\, 1 \,+\, x</math> therefore <math>\scriptstyle \log(e^x) \,>\, \log(1 \,+\, x)</math> , so <math> \scriptstyle x \,>\, \log(1 \,+\, x)</math>. So | |||
: <math> \sum_p \frac{1}{p - 1}>\sum_p \log\left(1+\frac{1}{p-1}\right)=\quad \log \left( \sum_{n=1}^\infty \frac{1}{n}\right)</math> | |||
Hence <math>\scriptstyle \sum_p \frac{1}{p - 1}</math> diverges. But <math>\scriptstyle \frac{1}{p_i - 1} \,<\, \frac{1}{p_{i - 1}}</math> (consider the '''''i''''' from 3). Where <math>\scriptstyle p_i</math> is the '''''i'''''<sup> th</sup> prime, (because <math>\scriptstyle {p_i-1} \,>\, {p_{i-1}}</math>). | |||
Hence <math>\scriptstyle \sum_p \frac{1}{p}</math> diverges. | |||
===Second=== | |||
The following [[proof by contradiction]] is due to [[Paul Erdős]]. | |||
Let ''p''<sub>''i''</sub> denote the ''i''<sup>''th''</sup> prime number. Assume that the [[series (mathematics)|sum]] of the reciprocals of the primes [[convergent series|converges]]; i.e., | |||
:<math>\sum_{i=1}^\infty {1\over p_i} < \infty</math> | |||
Then there exists a smallest [[Positive number|positive]] [[integer]] ''k'' such that | |||
:<math>\sum_{i=k+1}^\infty {1\over p_i} < {1\over 2} \qquad(1)</math> | |||
For a positive integer ''x'' let ''M<sub>x</sub>'' denote the set of those ''n'' in {1, 2, . . ., ''x''} which are not [[divisible]] by any prime greater than ''p<sub>k</sub>''. We will now derive an upper and a lower estimate for the [[cardinality|number]] |''M<sub>x</sub>''| of elements in ''M<sub>x</sub>''. For large ''x'', these bounds will turn out to be contradictory. | |||
====Upper estimate==== | |||
Every ''n'' in ''M<sub>x</sub>'' can be written as ''n'' = ''r m''<sup>2</sup> with positive integers ''m'' and ''r'', where ''r'' is [[square-free integer|square-free]]. Since only the ''k'' primes ''p''<sub>1</sub>, …, ''p<sub>k</sub>'' can show up (with exponent 1) in the [[Fundamental theorem of arithmetic|prime factorization]] of ''r'', there are at most 2<sup>''k''</sup> different possibilities for ''r''. Furthermore, there are at most √''x'' possible values for ''m''. This gives us the upper estimate | |||
:<math>|M_x| \le 2^k\sqrt{x} \qquad(2)</math> | |||
====Lower estimate==== | |||
The remaining ''x'' − |''M<sub>x</sub>''| numbers in the [[set difference]] {1, 2, . . ., ''x''} \ ''M<sub>x</sub>'' are all divisible by a prime greater than ''p''<sub>''k''</sub>. Let ''N<sub>i,x</sub>'' denote the set of those ''n'' in {1, 2, . . ., x} which are divisible by the ''i''<sup>th</sup> prime ''p<sub>i</sub>''. Then | |||
:<math>\{1,2,\ldots,x\}\setminus M_{x}=\bigcup_{i=k+1}^\infty N_{i,x}</math> | |||
Since the number of integers in ''N<sub>i,x</sub>'' is at most ''x''/''p<sub>i</sub>'' (actually zero for ''p<sub>i</sub>'' > ''x''), we get | |||
:<math>x-|M_x| \le \sum_{i=k+1}^\infty |N_{i,x}|< \sum_{i=k+1}^\infty {x\over p_i}</math> | |||
Using (1), this implies | |||
:<math>{x\over 2} < |M_x| \qquad(3)</math> | |||
====Contradiction==== | |||
When ''x'' ≥ 2<sup>2''k'' + 2</sup>, the estimates (2) and (3) cannot both hold, because <math>\tfrac{x}{2}\ge 2^k\sqrt{x}</math>. | |||
===Third === | |||
Here is another proof that actually gives a lower estimate for the partial sums; in particular, it shows that these sums grow at least as fast as log(log(''n'')). The proof is an adaptation of the product expansion idea of [[Euler]]. In the following, a sum or product taken over ''p'' always represents a sum or product taken over a specified set of primes. | |||
The proof rests upon the following four inequalities: | |||
* Every positive integer ''i'' can be uniquely expressed as the product of a square-free integer and a square. This gives the inequality | |||
::<math> \sum_{i=1}^n{\frac{1}{i}} \le \prod_{p \le n}{\left(1 + \frac{1}{p}\right)}\sum_{k=1}^n{\frac{1}{k^2}}</math> | |||
:where for every ''i'' between 1 and ''n'' the (expanded) product corresponds to the [[Radical of an integer|square-free part]] of ''i'' and the sum corresponds to the square part of ''i'' (see [[fundamental theorem of arithmetic]]). | |||
* The upper estimate for the [[natural logarithm]] | |||
::<math> | |||
\log(n+1) | |||
= \int_1^{n+1}\frac{dx}x | |||
= \sum_{i=1}^n\underbrace{\int_i^{i+1}\frac{dx}x}_{{} \,<\, 1/i} | |||
< \sum_{i=1}^n{\frac{1}{i}} | |||
</math> | |||
* The lower estimate 1 + ''x'' < exp(''x'') for the [[exponential function]], which holds for all ''x'' > 0. | |||
* Let n ≥ 2. The upper bound (using a [[telescoping sum]]) for the partial sums (convergence is all we really need) | |||
::<math> | |||
\sum_{k=1}^n{\frac{1}{k^2}} | |||
< 1 + \sum_{k=2}^n\underbrace{\left(\frac1{k - \frac{1}{2}} - \frac1{k + \frac{1}{2}}\right)}_{=\, 1/(k^2 - 1/4) \,>\, 1/k^2} | |||
= 1 + \frac23 - \frac1{n + \frac{1}{2}} < \frac53 | |||
</math> | |||
Combining all these inequalities, we see that | |||
:<math>\begin{align} | |||
{} & {} \log(n+1) \\ | |||
< &\sum_{i=1}^n\frac{1}{i} \\ | |||
\le &\prod_{p \le n}{\left(1 + \frac{1}{p}\right)}\sum_{k=1}^n{\frac{1}{k^2}} \\ | |||
< &\frac53\prod_{p \le n}{\exp\left(\frac{1}{p}\right)} \\ | |||
= &\frac53\exp\left(\sum_{p \le n}{\frac{1}{p}}\right) | |||
\end{align}</math> | |||
Dividing through by {{sfrac|5|3}} and taking the natural logarithm of both sides gives | |||
:<math>\log \log(n + 1) - \log\frac53 < \sum_{p \le n}{\frac{1}{p}}</math> | |||
as desired. [[Q.E.D.|∎]] | |||
Using | |||
:<math>\sum_{k=1}^\infty{\frac{1}{k^2}} = \frac{\pi^2}6</math> | |||
(see [[Basel problem]]), the above constant ln ({{sfrac|5|3}}) = 0.51082… can be improved to ln({{sfrac|π<sup>2</sup>|6}}) = 0.4977…; in fact it turns out that | |||
:<math> | |||
\lim_{n \to \infty } \left( | |||
\sum_{p \leq n} \frac{1}{p} - \log \log(n) | |||
\right) = M | |||
</math> | |||
where ''M'' = 0.261497… is the [[Meissel–Mertens constant]] (somewhat analogous to the much more famous [[Euler–Mascheroni constant]]). | |||
===Fourth=== | |||
From [[Dusart's inequality]], we get | |||
:<math> p_n < n \log n + n \log \log n \quad\mbox{for } n \ge 6</math> | |||
Then | |||
:<math>\begin{align} | |||
\sum_{n=1}^\infty \frac1{ p_n} | |||
&\ge \sum_{n=6}^\infty \frac1{ p_n} \\ | |||
&\ge \sum_{n=6}^\infty \frac1{ n \log n + n \log \log n} \\ | |||
&\ge \sum_{n=6}^\infty \frac1{2n \log n} \\ | |||
&= \infty | |||
\end{align}</math> | |||
by the [[integral test for convergence]]. This shows that the series on the left diverges. | |||
==See also== | |||
*[[Euclid's theorem]] that there are infinitely many primes | |||
*[[Small set (combinatorics)]] | |||
*[[Brun's theorem]] | |||
==References== | |||
* {{cite book | author=William Dunham | authorlink=William Dunham (mathematician) | title=Euler The Master of Us All | publisher=[[Mathematical Association of America|MAA]] | year=1999 | isbn=0-88385-328-0 | pages=61–79 }} | |||
==External links== | |||
* [http://www.utm.edu/research/primes/infinity.shtml Chris K. Caldwell: There are infinitely many primes, but, how big of an infinity?] | |||
{{DEFAULTSORT:Sum Of The Reciprocals Of The Primes Diverges}} | |||
[[Category:Mathematical series]] | |||
[[Category:Articles containing proofs]] | |||
[[Category:Theorems about prime numbers]] |
Revision as of 23:22, 28 January 2014
The sum of the reciprocals of all prime numbers diverges; that is:
This was proved by Leonhard Euler in 1737, and strengthens Euclid's 3rd-century-BC result that there are infinitely many prime numbers.
There are a variety of proofs of Euler's result, including a lower bound for the partial sums stating that
for all natural numbers n. The double natural logarithm indicates that the divergence might be very slow, which is indeed the case, see Meissel–Mertens constant.
The harmonic series
First, we describe how Euler originally discovered the result. He was considering the harmonic series
He had already used the following "product formula" to show the existence of infinitely many primes.
(Here, the product is taken over all primes p; in the following, a sum or product taken over p always represents a sum or product taken over a specified set of primes, unless noted otherwise.)
Such infinite products are today called Euler products. The product above is a reflection of the fundamental theorem of arithmetic. Of course, the above "equation" is not necessary because the harmonic series is known (by other means) to diverge. This type of formal manipulation was common at the time, when mathematicians were still experimenting with the new tools of calculus.Potter or Ceramic Artist Truman Bedell from Rexton, has interests which include ceramics, best property developers in singapore developers in singapore and scrabble. Was especially enthused after visiting Alejandro de Humboldt National Park.
Euler noted that if there were only a finite number of primes, then the product on the right would clearly converge, contradicting the divergence of the harmonic series. (In modern language, we now say that the existence of infinitely many primes is reflected by the fact that the Riemann zeta function has a simple pole at s = 1.)
Proofs
First
Euler took the above product formula and proceeded to make a sequence of audacious leaps of logic. First, he took the natural logarithm of each side, then he used the Taylor series expansion for ln(x) as well as the sum of a geometric series:
for a fixed constant C < 1. Since the sum of the reciprocals of the first n positive integers is asymptotic to ln(n), (i.e. their ratio approaches one as n approaches infinity), Euler then concluded
It is almost certain that Euler meant that the sum of the reciprocals of the primes less than n is asymptotic to ln(ln(n)) as n approaches infinity. It turns out this is indeed the case; Euler had reached a correct result by questionable means.Potter or Ceramic Artist Truman Bedell from Rexton, has interests which include ceramics, best property developers in singapore developers in singapore and scrabble. Was especially enthused after visiting Alejandro de Humboldt National Park.
A variation
Since
Shows that therefore , so . So
Hence diverges. But (consider the i from 3). Where is the i th prime, (because ).
Second
The following proof by contradiction is due to Paul Erdős.
Let pi denote the ith prime number. Assume that the sum of the reciprocals of the primes converges; i.e.,
Then there exists a smallest positive integer k such that
For a positive integer x let Mx denote the set of those n in {1, 2, . . ., x} which are not divisible by any prime greater than pk. We will now derive an upper and a lower estimate for the number |Mx| of elements in Mx. For large x, these bounds will turn out to be contradictory.
Upper estimate
Every n in Mx can be written as n = r m2 with positive integers m and r, where r is square-free. Since only the k primes p1, …, pk can show up (with exponent 1) in the prime factorization of r, there are at most 2k different possibilities for r. Furthermore, there are at most √x possible values for m. This gives us the upper estimate
Lower estimate
The remaining x − |Mx| numbers in the set difference {1, 2, . . ., x} \ Mx are all divisible by a prime greater than pk. Let Ni,x denote the set of those n in {1, 2, . . ., x} which are divisible by the ith prime pi. Then
Since the number of integers in Ni,x is at most x/pi (actually zero for pi > x), we get
Using (1), this implies
Contradiction
When x ≥ 22k + 2, the estimates (2) and (3) cannot both hold, because .
Third
Here is another proof that actually gives a lower estimate for the partial sums; in particular, it shows that these sums grow at least as fast as log(log(n)). The proof is an adaptation of the product expansion idea of Euler. In the following, a sum or product taken over p always represents a sum or product taken over a specified set of primes.
The proof rests upon the following four inequalities:
- Every positive integer i can be uniquely expressed as the product of a square-free integer and a square. This gives the inequality
- where for every i between 1 and n the (expanded) product corresponds to the square-free part of i and the sum corresponds to the square part of i (see fundamental theorem of arithmetic).
- The upper estimate for the natural logarithm
- The lower estimate 1 + x < exp(x) for the exponential function, which holds for all x > 0.
- Let n ≥ 2. The upper bound (using a telescoping sum) for the partial sums (convergence is all we really need)
Combining all these inequalities, we see that
Dividing through by Electrical Engineer Stephan from Merrickville-Wolford, likes to spend time body building, property developers ec in singapore singapore and storytelling. Gets inspiration through travel and just spent 6 days at Würzburg Residence with the Court Gardens. and taking the natural logarithm of both sides gives
as desired. ∎
Using
(see Basel problem), the above constant ln (Electrical Engineer Stephan from Merrickville-Wolford, likes to spend time body building, property developers ec in singapore singapore and storytelling. Gets inspiration through travel and just spent 6 days at Würzburg Residence with the Court Gardens.) = 0.51082… can be improved to ln(Electrical Engineer Stephan from Merrickville-Wolford, likes to spend time body building, property developers ec in singapore singapore and storytelling. Gets inspiration through travel and just spent 6 days at Würzburg Residence with the Court Gardens.) = 0.4977…; in fact it turns out that
where M = 0.261497… is the Meissel–Mertens constant (somewhat analogous to the much more famous Euler–Mascheroni constant).
Fourth
From Dusart's inequality, we get
Then
by the integral test for convergence. This shows that the series on the left diverges.
See also
- Euclid's theorem that there are infinitely many primes
- Small set (combinatorics)
- Brun's theorem
References
- 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
My blog: http://www.primaboinca.com/view_profile.php?userid=5889534