OFFSET
1,1
MATHEMATICA
Reap[Do[p=Prime[i]; If[PrimeQ[p+8]&&Total[Last/@FactorInteger[p+4]]==2, Sow[p+4]], {i, 2*10^3}]][[2, 1]]
Select[Range[8100], PrimeOmega[#]==2&&AllTrue[#+{4, -4}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Oct 29 2017 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Nov 24 2006
STATUS
approved