OFFSET
1,1
MATHEMATICA
tpQ[n_]:=Module[{x=n^6+1091}, PrimeQ[x]&&PrimeQ[x+2]]
Select[Range[5000000], tpQ] (* Harvey P. Dale, Feb 17 2011 *)
PROG
(Magma) IsA181114:=func< n | IsPrime(k+1091) and IsPrime(k+1093) where k is n^6 >; [ n: n in [1..5000000] | IsA181114(n) ]; // Klaus Brockhaus, Feb 17 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Jason Earls, Oct 04 2010
STATUS
approved