OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
EXAMPLE
4 is in the sequence because 4^4+1 = 257 and 6^4+1 = 1297 are both prime.
MAPLE
for n from 0 by 2 to 3500 do: if type(n^4+1, prime)=true and type((n+2)^4+1, prime)=true then printf(`%d, `, n):else fi:od:
MATHEMATICA
lst={}; Do[p=n^4+1; q=(n+2)^4+1; If[PrimeQ[p] && PrimeQ[q], AppendTo[lst, n]], {n, 0, 3000}]; lst
Select[Range[3500], AllTrue[{#^4+1, (#+2)^4+1}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Jul 29 2015 *)
PROG
(Magma) [n: n in [0..3300] | IsPrime(n^4 + 1) and IsPrime((n + 2)^4 + 1)]; // Vincenzo Librandi, Oct 13 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michel Lagneau, Oct 12 2012
STATUS
approved