 Solution

In the level-$D$ convention used here, the Selberg upper-bound sieve states that a nonnegative sequence with sieve distribution $A_d=Xg(d)+r(d)$ satisfiesIt follows by minimizing the quadratic form with the Selberg sieve weights and using the Selberg least-common-multiple weights bound $|{\lambda }_d|\le 3^{\omega (d)}$.

Sieve the interval $x<n\le x+z$ by two and by the primes congruent to three modulo four, up to $L=\sqrt{D}$. The interval divisor counts giveThus $X=z$, $g(p)=1/p$, and $r(d)=O(1)$ uniformly in the location of the interval. Take $D=z^{1/2}$ and $L=z^{1/4}$.

To bound $J$, let $y=L^{1/4}$ and form the Euler product over the forbidden primes at most $y$. The given reciprocal-prime sum in residue class one modulo four, subtracted from the Mertens second theorem, yieldsMoreover, the Mertens first theorem bounds the logarithmic mean by ${\sum }_{p\le y}(\log p)/p=\log y+O(1)\le \frac{1}{2}\log L$ for large $L$. The truncated Euler-product lower bound gives $J\ge W/2\gg \sqrt{\log z}$.

The full-range summatory bound for three to the prime omega bounds the error by $O(D{\log }^{2}D)=O(\sqrt{z}{\log }^{2}z)$, which is $O(z/\sqrt{\log z})$. Every integer whose prime factors are all congruent to one modulo four survives this finite sieve. Therefore the half-dimensional interval sieve givesThe constants are independent of $x$. Enlarging them covers the bounded range of $z$ before the asymptotic product estimates apply, so the result holds throughout $1000\le z\le x$.