OFFSET
1,1
COMMENTS
7*41=287+-6->primes,..
MATHEMATICA
q=7; lst={}; Do[p=Prime[n]; If[PrimeQ[q*p-Mod[p, q]]&&PrimeQ[q*p+Mod[p, q]], AppendTo[lst, p]], {n, 8!}]; lst
Select[Prime[Range[3000]], AllTrue[7#+{Mod[#, 7], -Mod[#, 7]}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Feb 13 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, May 26 2010
STATUS
approved