OFFSET
1,1
COMMENTS
The difference d(x) = pi(x,4,3) - pi(x,4,1) changes sign infinitely often, see link "Prime Quadratic Effect". But this does not say anything about the amplitudes of these oscillations. For diagrams, see A349518, "Oscillations of d(x)". If d(x) has no lower limit, the current sequence is infinite. Regarding the upper limit, see A349518.
Note the gaps between 2, 26861 and 616897, 623681 and 12315529, 12366589 and 951821281.
LINKS
Eric Weisstein's World of Mathematics, Prime Quadratic Effect.
EXAMPLE
primes 4*j+1: 5, 13, 17, ...
4*j+3: 3, 7, 11, ...
d(x) = pi(x,4,3) - pi(x,4,1)
.
n x pi(x,4,3) pi(x,4,1) d(x)=1-n?
- ----- --------- --------- -----------
1 2 0 0 0=0 true a(1) = 2
2 3 1 0 1=-1 false a(2) != 3
2 5 1 1 2=-1 false a(2) != 5
...........................
2 26861 1472 1473 -1=-1 true a(3) = 26861
PROG
(Maxima) block(w:[2], su:0, sum:0, n:1, p:2, nmax: 25,
/* returns nmax terms */
while n<nmax do(
p: next_prime(p), su:su+mod(p, 4)-2,
if su<sum then(n:n+1, sum:su, w: append(w, [p]) ) ) ,
w);
CROSSREFS
KEYWORD
nonn
AUTHOR
Gerhard Kirchner, Nov 20 2021
STATUS
approved