A339920 Primes p such that p^2 - p*q + q^2 is prime, where q is the next prime after p.
2, 3, 53, 151, 167, 263, 373, 443, 467, 509, 523, 571, 1063, 1103, 1117, 1217, 1493, 1553, 1901, 1973, 2161, 2207, 2281, 2399, 2713, 2837, 2963, 3259, 3347, 3511, 4073, 4297, 4373, 4463, 4523, 4673, 4691, 4877, 5147, 5237, 5303, 5419, 5471, 5851, 6211, 6311, 6367
Offset: 1
Keywords
Links
- Robert Israel, Table of n, a(n) for n = 1..10000
Crossrefs
Cf. A339698.
Programs
-
Maple
q:= 2: count:= 0: R:= NULL: while count < 100 do p:= q; q:= nextprime(q); if isprime(p^2-p*q+q^2) then count:= count+1; R:= R, p; fi od: R; # Robert Israel, Dec 24 2020
-
PARI
forprime(p=1, 1e4, my(q=nextprime(p+1)); if(ispseudoprime(p^2-p*q+q^2), print1(p, ", ")));
Comments