OFFSET
1,1
EXAMPLE
74 is in the sequence because 74^16+1 = 808551180810136214718004658177 and 76^16+1 = 1238846438084943599707227160577 are both prime.
MATHEMATICA
lst={}; Do[p=n^16+1; q=(n+2)^16+1; If[PrimeQ[p] && PrimeQ[q], AppendTo[lst, n]], {n, 0, 34000}]; lst
Select[Range[34000], AllTrue[{#^16+1, (#+2)^16+1}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Jun 24 2015 *)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michel Lagneau, Oct 17 2012
STATUS
approved