Prime number theorem

Let be the prime-counting function defined to be the number of primes less than or equal to , for any real number . For example, because there are four prime numbers (2, 3, 5 and 7) less than or equal to 10. The prime number theorem then states that is a good approximation to (where log here means the natural logarithm), in the sense that the limit of the quotient of the two functions and as increases without bound is 1: : \lim_{x\to\infty}\frac{\;\pi(x)\;}{\;\left[ \frac{x}{\log(x)}\right]\;} = 1, known as the asymptotic law of distribution of prime numbers. Using asymptotic notation this result can be restated as : \pi(x)\sim \frac{x}{\log x}. This notation (and the theorem) does not say anything about the limit of the difference of the two functions as increases without bound. Instead, the theorem states that approximates in the sense that the relative error of this approximation approaches 0 as increases without bound. The prime number theorem is equivalent to the statement that the th prime number satisfies : p_n \sim n\log(n), the asymptotic notation meaning, again, that the relative error of this approximation approaches 0 as increases without bound. For example, the th prime number is , and ()log() rounds to , a relative error of about 6.4%. On the other hand, the following asymptotic relations are logically equivalent: : \begin{align} \lim_{x\rightarrow \infty}\frac{\pi(x)\log x}{x}&=1,\text{ and}\\ \lim_{x\rightarrow \infty}\frac{\pi(x)\log \pi(x)}{x}\,&=1. \end{align} As outlined below, the prime number theorem is also equivalent to : \lim_{x\to\infty} \frac{\vartheta (x)}x = \lim_{x\to\infty} \frac{\psi(x)}x=1, where and are the first and the second Chebyshev functions respectively, and to : \lim_{x \to \infty} \frac{M(x)}{x}=0, where M(x)=\sum_{n \leq x} \mu(n) is the Mertens function. == History of the proof of the asymptotic law of prime numbers ==

Based on the tables by Anton Felkel and Jurij Vega, Adrien-Marie Legendre conjectured in 1797 or 1798 that is approximated by the function where and are unspecified constants. In the second edition of his book on number theory (1808) he then made a more precise conjecture, with and Carl Friedrich Gauss considered the same question at age 15 or 16 "in the year 1792 or 1793", according to his own recollection in 1849. In 1838 Peter Gustav Lejeune Dirichlet came up with his own approximating function, the logarithmic integral (under the slightly different form of a series, which he communicated to Gauss). Both Legendre's and Dirichlet's formulas imply the same conjectured asymptotic equivalence of and stated above, although it turned out that Dirichlet's approximation is considerably better if one considers the differences instead of quotients. In two papers from 1848 and 1850, the Russian mathematician Pafnuty Chebyshev attempted to prove the asymptotic law of distribution of prime numbers. His work is notable for the use of the zeta function , for real values of the argument "", as in works of Leonhard Euler, as early as 1737. Chebyshev's papers predated Riemann's celebrated memoir of 1859, and he succeeded in proving a slightly weaker form of the asymptotic law, namely, that if the limit as goes to infinity of exists at all, then it is necessarily equal to one. He was able to prove unconditionally that this ratio is bounded above and below by 0.92129 and 1.10555, for all sufficiently large . During the 20th century, the theorem of Hadamard and de la Vallée Poussin also became known as the Prime Number Theorem. Several different proofs of it were found, including the "elementary" proofs of Atle Selberg (1949) Hadamard's and original proofs are long and elaborate; later proofs introduced various simplifications through the use of Tauberian theorems but remained difficult to digest. A short proof was discovered in 1980 by the American mathematician Donald J. Newman. Newman's proof is arguably the simplest known proof of the theorem, although it is not "elementary" since it uses Cauchy's integral theorem from complex analysis. == Proof sketch ==

Here is a sketch of the proof referred to in one of Terence Tao's lectures. Like most proofs of the PNT, it starts out by reformulating the problem in terms of a less intuitive, but better-behaved, prime-counting function. The idea is to count the primes (or a related set such as the set of prime powers) with weights to arrive at a function with smoother asymptotic behavior. The most common such generalized counting function is the Chebyshev function , defined by : \psi(x) = \sum_{k \geq 1} \sum_\overset{p^k \le x,}{\!\!\!\!p \text{ is prime}\!\!\!\!} \log p \; . is easily confused. After some experimentation, it's more robust to take the space out of the limit (symmetrically so it remains centred), so the spacing determined by the other symbols is not affected.--> This is sometimes written as : \psi(x) = \sum_{n\le x} \Lambda(n) \; , where is the von Mangoldt function, namely : \Lambda(n) = \begin{cases} \log p & \text{ if } n = p^k \text{ for some prime } p \text{ and integer } k \ge 1, \\ 0 & \text{otherwise.} \end{cases} It is now relatively easy to check that the PNT is equivalent to the claim that : \lim_{x\to\infty} \frac{\psi(x)}{x} = 1 \; . Indeed, this follows from the easy estimates : \psi(x) = \sum_\overset{p\le x}{\!\!\!\! p \text{ is prime}\!\!\!\!} \log p \left\lfloor \frac{\log x}{\log p} \right\rfloor \le \sum_\overset{p\le x}{\!\!\!\! p \text{ is prime}\!\!\!\!} \log x = \pi(x)\log x and (using big notation) for any , : \psi(x) \ge \sum_{\!\!\!\!\overset{x^{1-\varepsilon}\le p\le x}{p \text{ is prime}}\!\!\!\!} \log p \ge \sum_{\!\!\!\!\overset{x^{1-\varepsilon}\le p\le x}{p \text{ is prime}}\!\!\!\!} (1-\varepsilon)\log x=(1-\varepsilon)\left(\pi(x)+O\left(x^{1-\varepsilon}\right)\right)\log x \; . The next step is to find a useful representation for . Let be the Riemann zeta function. It can be shown that is related to the von Mangoldt function , and hence to , via the relation : -\frac{\zeta'(s)}{\zeta(s)} = \sum_{n = 1}^\infty \Lambda(n) \, n^{-s} \; . A delicate analysis of this equation and related properties of the zeta function, using the Mellin transform and Perron's formula, shows that for non-integer the equation : \psi(x) = x \; - \; \log(2\pi) \; - \!\!\!\! \sum\limits_{\rho :\, \zeta(\rho) = 0} \frac{x^\rho}{\rho} holds, where the sum is over all zeros (trivial and nontrivial) of the zeta function. This striking formula is one of the so-called explicit formulas of number theory, and is already suggestive of the result we wish to prove, since the term (claimed to be the correct asymptotic order of ) appears on the right-hand side, followed by (presumably) lower-order asymptotic terms. The next step in the proof involves a study of the zeros of the zeta function. The trivial zeros −2, −4, −6, −8, ... can be handled separately: : \sum_{n=1}^\infty \frac{1}{2n\,x^{2n}} = -\frac{1}{2}\log\left(1-\frac{1}{x^2}\right), which vanishes for large . The nontrivial zeros, namely those on the critical strip , can potentially be of an asymptotic order comparable to the main term if , so we need to show that all zeros have real part strictly less than 1. === Non-vanishing on Re(s) = 1 === To do this, we take for granted that is meromorphic in the half-plane , and is analytic there except for a simple pole at , and that there is a product formula : \zeta(s)=\prod_p\frac{1}{1-p^{-s}} for . This product formula follows from the existence of unique prime factorization of integers, and shows that is never zero in this region, so that its logarithm is defined there and : \log\zeta(s)=-\sum_p\log \left(1-p^{-s} \right)=\sum_{p,n}\frac{p^{-ns}}{n} \; . Write ; then : \big| \zeta(x+iy) \big| = \exp\left( \sum_{n,p} \frac{\cos ny\log p}{np^{nx}} \right) \; . Now observe the identity : 3 + 4 \cos \phi+ \cos 2 \phi = 2 ( 1 + \cos \phi )^2\ge 0 \; , so that : \left| \zeta(x)^3 \zeta(x+iy)^4 \zeta(x+2iy) \right| = \exp\left( \sum_{n,p} \frac{3 + 4 \cos(ny\log p) + \cos( 2 n y \log p )}{np^{nx}} \right) \ge 1 for all . Suppose now that . Certainly is not zero, since has a simple pole at . Suppose that and let tend to 1 from above. Since \zeta(s) has a simple pole at and stays analytic, the left hand side in the previous inequality tends to 0, a contradiction. Finally, we can conclude that the PNT is heuristically true. To rigorously complete the proof there are still serious technicalities to overcome, due to the fact that the summation over zeta zeros in the explicit formula for does not converge absolutely but only conditionally and in a "principal value" sense. There are several ways around this problem but many of them require rather delicate complex-analytic estimates. Edwards's book provides the details. Another method is to use Ikehara's Tauberian theorem, though this theorem is itself quite hard to prove. D.J. Newman observed that the full strength of Ikehara's theorem is not needed for the prime number theorem, and one can get away with a special case that is much easier to prove. == Newman's proof of the prime number theorem ==

D.J. Newman gives a quick proof of the prime number theorem (PNT). The proof is "non-elementary" by virtue of relying on complex analysis, but uses only elementary techniques from a first course in the subject: Cauchy's integral formula, Cauchy's integral theorem and estimates of complex integrals. Here is a brief sketch of this proof. See for the complete details. The proof uses the same preliminaries as in the previous section except instead of the function \ \psi\ , the Chebyshev function \ \vartheta(x) = \sum_{p\le x} \log p\ is used, which is obtained by dropping some of the terms from the series for \ \psi~. Similar to the argument in the previous proof based on Tao's lecture, we can show that and for any Thus, the PNT is equivalent to \ \lim _{x \to \infty} \tfrac{\ \vartheta(x)\ }{ x } = 1 ~. Likewise instead of \ - \tfrac{\ \zeta '(s)\ }{ \zeta(s) }\ the function \ \Phi(s) = \sum_{p\le x} \tfrac{\ \log p\ }{\; p^s\ }\ is used, which is obtained by dropping some terms in the series for \ - \tfrac{\ \zeta '(s)\ }{ \zeta(s) } ~. The functions \ \Phi(s)\ and \ -\tfrac{\ \zeta'(s)\ }{ \zeta(s) }\ differ by a function holomorphic on \ \Re(s) = 1 ~. Since, as was shown in the previous section, \ \zeta(s)\ has no zeroes on the line \ \Re = 1\ , and \ \Phi(s) - \tfrac{ 1 }{\ s - 1\ }\ has no singularities on \ \Re(s) = 1 ~. One further piece of information needed in Newman's proof, and which is the key to the estimates in his simple method, is that \ \tfrac{\ \vartheta(x)\ }{ x }\ is bounded. This is proved using an ingenious and easy method due to Chebyshev. Integration by parts shows how \ \vartheta(x)\ and \ \Phi(s)\ are related: For \ \Re(s) > 1\ , : \Phi(s) ~=~ \int _1^\infty \frac{ \operatorname{d}\vartheta(x) }{\; x^s\ } ~=~ s \int_1^\infty \frac{ \vartheta(x) }{\ x^{s+1}\ } \operatorname{d}x ~=~ s \int_0^\infty \frac{\; \vartheta(e^t)\ }{\; e^{st}\ } \operatorname{d}t ~. Newman's method proves the PNT by showing the integral : I \equiv \int_0 ^\infty \left( \frac{\ \vartheta(e^t)\ }{ e^t } -1 \right) \operatorname{d}t ~. converges, and therefore the integrand goes to zero as \ t \to \infty\ , which is the PNT. In general, the convergence of the improper integral does not imply that the integrand goes to zero at infinity, since it may oscillate, but since \ \vartheta\ is increasing, it is easy to show in this case. To show the convergence of \ I\ , for \ \Re(z) > 0\ let : g_T(z) \equiv \int_0^T f(t)\ e^{-zt} \operatorname{d}t \quad and \quad g(z) \equiv \int_0^\infty f(t)\ e^{-zt} \operatorname{d}t \quad where \quad f(t) \equiv \frac{\ \vartheta(e^t)\ }{ e^t } - 1\ then : \lim_{T \to \infty} g_T(z) ~=~ g(z) ~=~ \frac{\ \Phi(s)\ }{ s } - \frac{ 1 }{\ s - 1\ } \qquad \text{where} \quad z \equiv s - 1\ which is equal to a function holomorphic on the line \ \Re(z) = 0 ~. The convergence of the integral \ I\ , and thus the PNT, is proved by showing that \ \lim_{T \to \infty} g_T(0) ~=~ g(0) ~. This involves change of order of limits since it can be written \ \lim_{T \to \infty}\ \lim_{z \to 0} g_T(z) ~=~ \lim_{z \to 0}\ \lim_{T \to \infty}g_T(z)\ and therefore classified as a Tauberian theorem. The difference \ g(0) - g_T(0)\ is expressed using Cauchy's integral formula and then shown to be small for large \ T\ by estimating the integrand: Fix \ R > 0\ and \ \delta > 0\ so that \ g(z)\ is holomorphic in the region where \ |z| \le R ~ and ~ \Re(z) \ge - \delta\ , and let \ C\ be the boundary of that region. Since is in the interior of the region, Cauchy's integral formula gives : g(0) - g_T(0) ~=~ \frac{ 1 }{\ 2 \pi i\ }\int_C \bigl(\ g(z) - g_T(z)\ \bigr)\ \frac{\ \operatorname{d}z\ }{ z } ~=~ \frac{ 1 }{2 \pi i }\int_C \bigl(\ g(z) - g_T(z)\ \bigr)\ F(z)\ \frac{\ \operatorname{d}z\ }{ z }\ where \ F(z) \equiv e^{zT} \left( 1 + \frac{\; z^2\ }{\; R^2\ } \right)\ is the factor introduced by Newman, which does not change the integral since \ F\ is entire and \ F(0) = 1 ~. To estimate the integral, break the contour \ C\ into two parts, \ C = C_+ + C_-\ where \ C_+ \equiv C \cap \left\{ z\ \vert\ \Re(z) > 0 \right\}\ and \ C_- \equiv C \cap \left\{ z\ \vert\ \Re(z) \le 0 \right\} ~. Then : \ g(0) - g_T(0) ~~=~~ \int_{C_+} \int_T^\infty H(t,z) \operatorname{d}t\ \operatorname{d}z ~~-~ \int_{C_-}\int_0^T H(t,z) \operatorname{d}t \operatorname{d}z ~~+~ \int_{C_-} g(z)\ F(z) \frac{ \operatorname{d}z }{\ 2\pi i z\ }\ , where \ H(t,z) \equiv f(t)\ e^{-tz}\frac{ F(z) }{\ 2 \pi i\ } ~. Note that \ \frac{\ \vartheta(x)\ }{ x }\ , and hence \ f(t)\ , are bounded; so let \ B\ be some upper bound: \ B \ge \bigl| f(t) \bigr| ~. This bound, combined with the estimate \ \left| F \right|\ \le\ 2\ \exp\!\Bigl( T\ \Re(z) \Bigr) \frac{\ \left| \Re(z) \right|\ }{ R }\ for \ |z| = R\ , together give that the absolute value of the first integral must be \ \le \frac{\ B\ }{ R } ~. The integrand over \ C_-\ in the second integral is entire, so by Cauchy's integral theorem, the contour \ C_-\ can be modified to a semicircle of radius \ R\ in the left half-plane without changing the integral, and the same argument as for the first integral gives the absolute value of the second integral must be \ \le \frac{\ B\ }{ R } ~. Finally, letting \ T \to \infty\ , the third integral goes to zero since \ e^{zT}\ and hence \ F\ goes to zero on the contour. Combining the two estimates and the limit get : \limsup_{T \to \infty }\ \bigl|\ g(0) - g_T(0)\ \bigr|\ \le\ \frac{\ 2 B\ }{ R }\ ~. This holds for any \ R\ so \ \lim_{T \to \infty} g_T(0) = g(0)\ , and the PNT follows. == Prime-counting function in terms of the logarithmic integral ==

In a handwritten note on a reprint of his 1838 paper "", which he mailed to Gauss, Dirichlet conjectured (under a slightly different form appealing to a series rather than an integral) that an even better approximation to is given by the offset logarithmic integral function , defined by : \operatorname{Li}(x) = \int_2^x \frac{dt}{\log t} = \operatorname{li}(x) - \operatorname{li}(2). Indeed, this integral is strongly suggestive of the notion that the "density" of primes around should be . This function is related to the logarithm by the asymptotic expansion : \operatorname{Li}(x) \sim \frac{x}{\log x} \sum_{k=0}^\infty \frac{k!}{(\log x)^k} = \frac{x}{\log x} + \frac{x}{(\log x)^2} + \frac{2x}{(\log x)^3} + \cdots So, the prime number theorem can also be written as . In fact, in another paper in 1899 de la Vallée Poussin proved that : \pi(x) = \operatorname{Li} (x) + O \left(x e^{-a\sqrt{\log x}}\right) \quad\text{as } x \to \infty for some positive constant , where is the big notation. This has been improved to : \pi(x) = \operatorname{li} (x) + O \left(x \exp \left( -\frac{A(\log x)^\frac35}{(\log \log x)^\frac15} \right) \right) where A = 0.2098. In 2016, Timothy Trudgian proved an explicit upper bound for the difference between \pi(x) and \operatorname{li}(x): : \big| \pi(x) - \operatorname{li}(x) \big| \le 0.2795 \frac{x}{(\log x)^{3/4}} \exp \left( -\sqrt{ \frac{\log x}{6.455} } \right) for x \ge 229. The connection between the Riemann zeta function and is one reason the Riemann hypothesis has considerable importance in number theory: if established, it would yield a far better estimate of the error involved in the prime number theorem than is available today. More specifically, Helge von Koch showed in 1901 that if the Riemann hypothesis is true, the error term in the above relation can be improved to : \pi(x) = \operatorname{Li} (x) + O\left(\sqrt x \log x\right) (this last estimate is in fact equivalent to the Riemann hypothesis). The constant involved in the big notation was estimated in 1976 by Lowell Schoenfeld, assuming the Riemann hypothesis: : \big|\pi(x) - \operatorname{li}(x)\big| for all . He also derived a similar bound for the Chebyshev prime-counting function : : \big|\psi(x) - x\big| for all  . This latter bound has been shown to express a variance to mean power law (when regarded as a random function over the integers) and  noise and to also correspond to the Tweedie compound Poisson distribution. (The Tweedie distributions represent a family of scale invariant distributions that serve as foci of convergence for a generalization of the central limit theorem.) A lower bound is also derived by J. E. Littlewood, assuming the Riemann hypothesis: : \big|\pi(x) - \operatorname{li}(x)\big| = \Omega \left(\sqrt x\frac{\log\log\log x}{\log x} \right) The logarithmic integral is larger than for "small" values of . This is because it is (in some sense) counting not primes, but prime powers, where a power of a prime is counted as of a prime. This suggests that should usually be larger than by roughly \ \tfrac{1}{2} \operatorname{li}(\sqrt{x})\ , and in particular should always be larger than . However, in 1914, Littlewood proved that \ \pi(x) - \operatorname{li}(x)\ changes sign infinitely often. The first value of where exceeds is probably around ; see the article on Skewes' number for more details. (On the other hand, the offset logarithmic integral is smaller than already for ; indeed, , while .) == Elementary proofs ==

In the first half of the twentieth century, some mathematicians (notably G. H. Hardy) believed that there exists a hierarchy of proof methods in mathematics depending on what sorts of numbers (integers, reals, complex) a proof requires, and that the prime number theorem (PNT) is a "deep" theorem by virtue of requiring complex analysis. This belief was somewhat shaken by a proof of the PNT based on Wiener's tauberian theorem, though Wiener's proof ultimately relies on properties of the Riemann zeta function on the line Re(s)=1, where complex analysis must be used. In March 1948, Atle Selberg established, by "elementary" means, the asymptotic formula : \vartheta ( x )\log ( x ) + \sum\limits_{p \le x} {\log ( p )}\ \vartheta \left( {\frac{x}{p}} \right) = 2x\log ( x ) + O( x ) where : \vartheta ( x ) = \sum\limits_{p \le x} {\log ( p )} for primes . By July of that year, Selberg and Paul Erdős These proofs effectively laid to rest the notion that the PNT was "deep" in that sense, and showed that technically "elementary" methods were more powerful than had been believed to be the case. On the history of the elementary proofs of the PNT, including the Erdős–Selberg priority dispute, see an article by Dorian Goldfeld. However, this does not address the question of whether or not the standard proof of PNT can be formalized in PA. A more recent "elementary" proof of the prime number theorem uses ergodic theory, due to Florian Richter. The prime number theorem is obtained there in an equivalent form that the Cesàro sum of the values of the Liouville function is zero. The Liouville function is (-1)^{\omega(n)} where \omega(n) is the number of prime factors, with multiplicity, of the integer n. Bergelson and Richter (2022) then obtain this form of the prime number theorem from an ergodic theorem which they prove: : Let X be a compact metric space, T a continuous self-map of X, and \mu a T-invariant Borel probability measure for which T is uniquely ergodic. Then, for every f\in C(X), \tfrac1N\sum_{n=1}^Nf(T^{\omega(n)}x)\to \int_Xf\,d\mu,\quad\forall x\in X. This ergodic theorem can also be used to give "soft" proofs of results related to the prime number theorem, such as the Pillai–Selberg theorem and Erdős–Delange theorem. == Computer verifications ==

In 2005, Avigad et al. employed the Isabelle theorem prover to devise a computer-verified variant of the Erdős–Selberg proof of the PNT. This was the first machine-verified proof of the PNT. Avigad chose to formalize the Erdős–Selberg proof rather than an analytic one because while Isabelle's library at the time could implement the notions of limit, derivative, and transcendental function, it had almost no theory of integration to speak of. By developing the necessary analytic machinery, including the Cauchy integral formula, Harrison was able to formalize "a direct, modern and elegant proof instead of the more involved 'elementary' Erdős–Selberg argument". == Prime number theorem for arithmetic progressions ==

Let denote the number of primes in the arithmetic progression that are less than . Dirichlet and Legendre conjectured, and de la Vallée Poussin proved, that if and are coprime, then : \pi_{d,a}(x) \sim \frac{ \operatorname{Li}(x) }{ \varphi(d) } \ , where is Euler's totient function. In other words, the primes are distributed evenly among the residue classes modulo with  . This is stronger than Dirichlet's theorem on arithmetic progressions (which only states that there is an infinity of primes in each class) and can be proved using similar methods used by Newman for his proof of the prime number theorem. The Siegel–Walfisz theorem gives a good estimate for the distribution of primes in residue classes. Bennett et al. proved the following estimate that has explicit constants and (Theorem 1.3): Let \ge 3 be an integer and let be an integer that is coprime to . Then there are positive constants and such that : \left | \pi_{d,a}(x) - \frac{\ \operatorname{Li}(x)\ }{\ \varphi(d)\ } \right | where : A = \frac{1}{\ 840\ } \quad \text{ if } \quad 3 \leq d \leq 10^4 \quad \text{ and } \quad A = \frac{1}{\ 160\ } \quad \text{ if } \quad d > 10^4 ~, and : B = 8 \cdot 10^9 \quad \text{ if } \quad 3 \leq d \leq 10^5 \quad \text{ and } \quad B = \exp(\ 0.03\ \sqrt{d\ }\ (\log{d})^3 \ ) \quad \text{ if } \quad d > 10^5\ . Prime number race Although we have in particular : \pi_{4,1}(x) \sim \pi_{4,3}(x) \ , empirically the primes congruent to 3 are more numerous and are nearly always ahead in this "prime number race"; the first reversal occurs at . However Littlewood showed in 1914 Granville and Martin give a thorough exposition and survey. This follows that 1 and 9 are quadratic residues modulo 10, and 3 and 7 are quadratic nonresidues modulo 10. == Non-asymptotic bounds on the prime-counting function ==

The prime number theorem is an asymptotic result. It gives an ineffective bound on as a direct consequence of the definition of the limit: for all , there is an such that for all , : (1-\varepsilon)\frac {x}{\log x} \; However, better bounds on are known, for instance Pierre Dusart's : \frac{x}{\log x}\left(1+\frac{1}{\log x}\right) \; The first inequality holds for all and the second one for . The proof by de la Vallée Poussin implies the following bound: For every , there is an such that for all , : \frac {x}{\log x - (1 - \varepsilon)} \; The value gives a weak but sometimes useful bound for : : \frac {x}{\log x + 2} \; In Pierre Dusart's thesis there are stronger versions of this type of inequality that are valid for larger . Later in 2010, Dusart proved: : \begin{align} \frac {x}{\log x - 1} \; & Note that the first of these obsoletes the condition on the lower bound. == Approximations for the nth prime number ==

As a consequence of the prime number theorem, one gets an asymptotic expression for the th prime number, denoted by : : p_n \sim n \log n. A better approximation is by Cesàro (1894): : p_n = n B_2(\log n),\text{ where} : B_2(x) = x + \log x - 1 + \frac{\log x - 2}{x} - \frac{(\log x)^2 - 6 \log x + 11}{2x^2} + o\left(\frac 1{x^2}\right). Again considering the th prime number , assuming the trailing error term is zero gives an estimate of ; the first 5 digits match and relative error is about 0.46 parts per million. Cipolla (1902) showed that these are the leading terms of an infinite series which may be truncated at arbitrary degree, with : B_k(x) = x + \log x - 1 - \sum_{i=1}^k (-1)^i \frac{P_i(\log x)}{ix^i} + O\left(\frac{(\log x)^{k+1}}{x^{k+1}}\right), where each is a degree- monic polynomial. (, , , and so on. found tighter bounds using the form of the Cesàro/Cipolla approximations but varying the lowest-order constant term. is the same function as above, but with the lowest-order constant term replaced by a parameter : : \begin{align} p_n \; &> \; n B_0(\log n; 1) && \text{for } n \ge 2,\text{ and}\\ p_n \; & \; n B_1(\log n; 2.25) && \text{for } n \ge 2,\text{ and}\\ p_n \; & The upper bounds can be extended to smaller by loosening the parameter. For example, for all . Axler (2019) extended this to higher order, showing: : \begin{align} p_n \; &> \; n B_2(\log n;11.321) \quad \text{for } n \ge 2, \text{ and }\\ p_n \; & Again, the bound on may be decreased by loosening the parameter. For example, for . == Table of π(x), x / log x, and li(x) ==

The table compares exact values of to the two approximations and . The approximation difference columns are rounded to the nearest integer, but the "% error" columns are computed based on the unrounded approximations. The last column, , is the average prime gap below . : The value for was originally computed assuming the Riemann hypothesis; it has since been verified unconditionally. == Analogue for irreducible polynomials over a finite field ==

There is an analogue of the prime number theorem that describes the "distribution" of irreducible polynomials over a finite field; the form it takes is strikingly similar to the case of the classical prime number theorem. To state it precisely, let be the finite field with elements, for some fixed , and let be the number of monic irreducible polynomials over whose degree is equal to . That is, we are looking at polynomials with coefficients chosen from , which cannot be written as products of polynomials of smaller degree. In this setting, these polynomials play the role of the prime numbers, since all other monic polynomials are built up of products of them. One can then prove that : N_n \sim \frac{q^n}{n}. If we make the substitution , then the right hand side is just : \frac{x}{\log_q x}, which makes the analogy clearer. Since there are precisely monic polynomials of degree (including the reducible ones), this can be rephrased as follows: if a monic polynomial of degree is selected randomly, then the probability of it being irreducible is about . One can even prove an analogue of the Riemann hypothesis, namely that : N_n = \frac{q^n}n + O\left(\frac{q^\frac{n}{2}}{n}\right). The proofs of these statements are far simpler than in the classical case. It involves a short, combinatorial argument, summarised as follows: every element of the degree extension of is a root of some irreducible polynomial whose degree divides ; by counting these roots in two different ways one establishes that : q^n = \sum_{d\mid n} d N_d, where the sum is over all divisors of . Möbius inversion then yields : N_n = \frac{1}{n} \sum_{d\mid n} \mu\left(\frac{n}{d}\right) q^d, where is the Möbius function. (This formula was known to Gauss.) The main term occurs for , and it is not difficult to bound the remaining terms. The "Riemann hypothesis" statement depends on the fact that the largest proper divisor of can be no larger than . == See also ==