OFFSET
1,1
LINKS
Felix Fröhlich, Table of n, a(n) for n = 1..10000
EXAMPLE
733 - 666 = 67, 733 + 666 = 1399 and 67, 733, 1399 are all prime, so 733 is a term of the sequence. - Felix Fröhlich, May 26 2022
MATHEMATICA
lst={}; b=666; Do[p=Prime[n]; If[PrimeQ[p+b]&&PrimeQ[p-b], AppendTo[lst, p]], {n, 5!+2, 7!}]; lst
Select[Prime[Range[122, 600]], AllTrue[#+{666, -666}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Apr 08 2018 *)
PROG
(PARI) forprime(p=1, 3900, if(ispseudoprime(p+666) && ispseudoprime(p-666), print1(p, ", "))) \\ Felix Fröhlich, May 26 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Aug 30 2008
EXTENSIONS
Name edited by Felix Fröhlich, May 26 2022
STATUS
approved