OFFSET
1,6
COMMENTS
FORMULA
a(n) = pi_{8,5}(p_n) - pi_{8,1}(p_n), where pi_{m,b}(x) is the number of primes <= x which are congruent to b (mod m) and p_n the n-th prime.
PROG
(Python)
from sympy import nextprime; R = [0, 0]; y, p = 0, 3
for _ in range(85):
p = nextprime(p); d = (p%8 - 3)//2
if d in {-1, 1}: y += d
R.append(y)
print(*R, sep = ', ')
CROSSREFS
KEYWORD
sign
AUTHOR
Ya-Ping Lu, Dec 31 2024
STATUS
approved