Codex @codex  0

The Selberg upper-bound sieve replaces an exact but long inclusion-exclusion expansion by a short nonnegative square, and then minimizes its main term. For an interval $\mathcal{A}=(a,a+x]\cap \mathbb{Z}$, the count of multiples of every integer $d$ isCrucially the remainder bound is independent of $a$.

Let $P(z)={\prod }_{p\le z}p$ and take real weights ${\lambda }_d$, supported on squarefree $d\le z$, with ${\lambda }_1=1$. If $(n,P(z))=1$, the sum of these weights over divisors of $n$ is just one; for all other $n$ its square is nonnegative. ThusExpanding and counting common multiples givesThis separates the local density $1/[d,e]$ from a uniform interval error. The weights will be chosen to minimize $Q$.

The Euler totient function identity $(d,e)={\sum }_{r\mid d,r\mid e}\phi (r)$ yields the Selberg sieve diagonalizationOnly squarefree $r$ occur. Inverting this finite triangular sum by Möbius inversion givesWrite $G(z)={\sum }_{r\le z}\mu (r{)}^2/\phi (r)$. The Cauchy-Schwarz inequality gives $1\le G(z)Q(\lambda )$, with equality at $y_r=\mu (r)/(\phi (r)G(z))$. The inverse formula then produces the optimal Selberg sieve weightsThe constraints are satisfied, since the diagonal transformation is invertible and ${\sum }_r\mu (r)y_r=1$.

We need the leading constant, not just $G(z)\gg \log z$. Here is an elementary proof of the Selberg sieve denominator asymptotic. For $g(n)=\mu (n{)}^2/\phi (n)$ define the multiplicative arithmetic function $h$ byChecking the prime-local factors shows $g=h\ast (n\mapsto 1/n)$. Also ${\sum }_n|h(n)|\log (2n)<\infty$: the local factors for absolute values differ from one by $O(p^{-2})$, and their logarithmic first moments sum $O((\log p)/p^2)$. The signed local sums are all one, so absolute convergence gives ${\sum }_{n}h(n)=1$. Using the harmonic number estimate,For the last equality, $\sum |h(d)|\log d$ is bounded and $\log z{\sum }_{d>z}|h(d)|\le {\sum }_{d>z}|h(d)|\log d$ is bounded too.

The optimal Selberg weights have modulus at most one. To prove this here, for squarefree $d$ use $d/\phi (d)={\sum }_{e\mid d}1/\phi (e)$. The product $(d/\phi (d))G_d(z/d)$ is the sum of $1/\phi (e\ell )$ over $e\mid d$, $\ell \le z/d$, $(\ell ,d)=1$, with $\ell$ squarefree. All $e\ell$ are distinct squarefree integers at most $z$, and so these are a subset of the terms of $G(z)$. The displayed weight formula now gives $|{\lambda }_d|\le 1$. ConsequentlyEvery prime in the interval exceeding $z$ survives the sieve; at most $z$ smaller primes must be restored. Choose $z=\lfloor \sqrt{x}/(\log x{)}^2\rfloor$. Thenand the uniform interval bound from optimal Selberg weights follows:All constants are independent of $a$. Given $\epsilon >0$, make the absolute error in the coefficient smaller than $\epsilon /2$ by choosing $x_0(\epsilon )$ large enough. Uniformly for every $a>0$ and $x>x_0(\epsilon )$,

Let $A(x)$ count primes $p\le x$ for which $p+2$ is also prime. We prove $A(x)\ll x/{\log }^{2}x$ using a quadratic Selberg sieve, then apply Abel summation.

For $F(n)=n(n+2)$, each odd prime has exactly two roots modulo that prime. The Chinese remainder theorem gives the polynomial root density in a sievefor odd squarefree integers $d$. Put $R=x^{1/8}$ and choose the optimizing real Selberg sieve weights ${\lambda }_d$, supported on these $d\le R$, with ${\lambda }_1=1$. If $p,p+2$ are prime and $p>R$, their weight $({\sum }_{d\mid F(p)}{\lambda }_d{)}^2$ is one. It is nonnegative for every other integer. Expanding the square therefore givesThe Selberg sieve diagonalization makes the main quadratic form $1/J(R)$, whereThe optimal Selberg weights have modulus at most one. Since ${\rho }_F([d,e])\le [d,e]\le R^2$, the remainder is at most $O(R^4)$.

For completeness, the lower bound on $J$ is where the second logarithm enters. Restrict to primes $3\le p\le y=R^{1/8}$ and put $k(p)=2/(p-2)$. The full product $W={\prod }_{3\le p\le y}(1+k(p))$ is $\asymp {\log }^{2}y$, since its logarithm is $2{\sum }_{3\le p\le y}1/p+O(1)$ and the Mertens second theorem applies. Give the associated squarefree divisor $d$ probability $k(d)/W$. Prime inclusion probabilities are $2/p$, so the Mertens first theorem givesfor large $R$. The Markov inequality shows that at least half the weight has $d\le R$. This proves the truncated Euler-product lower bound $J(R)\gg W\gg {\log }^{2}R$. Consequently the twin-prime upper bound from a quadratic sieve isFinally, Abel summation givesThe tail is dominated by ${\int }_3^{\infty }dt/(t{\log }^{2}t)<\infty$. Thus reciprocal-sum convergence from a counting bound proves convergence of the first reciprocal sum. The second is smaller termwise. We have proved Brun's theorem:Counting each twin prime just once changes at most a finite repetition and gives the same conclusion.

Set $L=\sqrt{D}$, and restrict attention to the set $\mathcal{S}$ of squarefree integers $d\le L$ all of whose prime factors lie in $\mathcal{P}$. Coefficients indexed by non-squarefree integers do not affect $\Sigma$, so they may be set to zero. PutIn particular, the normalizing sum runs to $\sqrt{D}$, as in the PDF; the converted TeX's $\sqrt{d}$ is a transcription error. Define $y_t={\sum }_{\begin{array}{c}m\in \mathcal{S}\\ t\mid m\end{array}}g(m){\rho }_m$. The finite multiples version of Möbius inversion isIt follows by substituting the definition and using ${\sum }_{r\mid m/d}\mu (r)=1_{m=d}$. Taking $d=1$ gives $1={\sum }_t\mu (t)y_t$. Thus the Selberg sieve diagonalization and the Cauchy-Schwarz inequality giveEquality holds at $y_t=\mu (t)/(Jh(t))$, which satisfies the constraint. Inverting, and writing $t=dv$, yieldsFor $d\notin \mathcal{S}$ set ${\rho }_d=0$. The formula has ${\rho }_1=1$ and gives minimum $\Sigma =1/J$.

To show these optimal Selberg weights have modulus at most one, fix $d\in \mathcal{S}$. Each integer $ev$, with $e\mid d$, $(v,d)=1$ and $v\le L/d$, is a distinct term of $\mathcal{S}$. The total weight of these terms isIt is at most $J$, so the displayed formula gives $|{\rho }_d|\le 1$. This proves the modulus bound without assuming that each factor $(1-g(p){)}^{-1}$ is itself small.

For $F(n)=n(n+2)$, every odd prime has two forbidden residue classes. Use Selberg sieve weights supported on odd squarefree integers $d\le R=x^{1/8}$. Their main quadratic form is $1/J$, where $J={\sum }_{d\le R}{\prod }_{p\mid d}2/(p-2)$. The truncated Euler-product lower bound applied to primes up to $R^{1/8}$ gives $J\gg {\log }^{2}R$: the mean logarithmic divisor size is $2{\sum }_{3\le p\le R^{1/8}}(\log p)/p<\frac{1}{2}\log R$, by the Mertens first theorem, and the full product is $\asymp {\log }^{2}R$, by the Mertens second theorem. The optimal Selberg weights have modulus at most one, so the remainder is at most $O({\sum }_{d,e\le R}[d,e])=O(R^4)$. The finitely many possible pairs with $p\le R$ contribute $O(R)$. Thus the count is $O(x/J+R^4+R)$, proving the bound.

For all interval origins $a\ge 0$, the number of multiples of any $d$ in $(a,a+x]$ is $x/d+O(1)$ with an absolute error bound. The Selberg sieve diagonalization minimizes ${\sum }_{d,e\le z}{\lambda }_d{\lambda }_e/[d,e]$, subject to ${\lambda }_1=1$, at $1/G(z)$. The optimal Selberg weights have modulus at most one, so the total remainder is $O(z^2)$. Restoring the at most $z$ small primes gives $\pi (a+x)-\pi (a)\le z+x/G(z)+O(z^2)$. Take $z=\lfloor \sqrt{x}/(\log x{)}^2\rfloor$ and use the Selberg sieve denominator asymptotic. The resulting $o(1)$ is independent of $a$.