OFFSET
1,1
COMMENTS
If p == -q (mod 3) then p^2 - p*q + q^2 is divisible by 3.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
a(5) = 19 is a term because 19 and 23 are consecutive primes and (19^2 - 19*23 + 23^2)/3 = 151 is prime.
MAPLE
R:= NULL: q:= 2: count:= 0:
while count < 100 do
p:= q; q:= nextprime(p);
r:= (p^2-p*q+q^2)/3;
if r::integer and isprime(r) then
count:= count+1; R:= R, p;
fi;
od:
R;
PROG
(Python)
from sympy import isprime, nextprime
def aupto(limit):
p, q, num, alst = 2, 3, 7, []
while p <= limit:
if num%3 == 0 and isprime(num//3): alst.append(p)
p, q = q, nextprime(q)
num = p**2 - p*q + q**2
return alst
print(aupto(2819)) # Michael S. Branicky, Mar 18 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Mar 18 2021
STATUS
approved