OFFSET
1,1
MATHEMATICA
Select[Range[10^6/2], PrimeQ[10^6# + 1] && PrimeQ[10^6# + 3] && PrimeQ[10^6# + 7] && PrimeQ[10^6# + 9] &]
Select[Range[400000], AllTrue[10^6*#+{1, 3, 7, 9}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Mar 11 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Oct 29 2001
STATUS
approved