OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Seiichi Manyama)
FORMULA
a(k) = A108814(k) - 1. - Jeppe Stig Nielsen, Feb 26 2016
MATHEMATICA
Select[Range[5000], AllTrue[{#^2+1, (#+2)^2+1}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Jul 23 2014 *)
Select[Range[5000], PrimeQ[#^2 + 1] && PrimeQ[(# + 2)^2 + 1] &] (* Vincenzo Librandi, Feb 27 2016 *)
PROG
(Magma) [n: n in [1..5000] | IsPrime(n^2+1) and IsPrime((n+2)^2+1)]; // Vincenzo Librandi, Feb 27 2016
(PARI) isok(n) = isprime(n^2+1) && isprime((n+2)^2+1); \\ Michel Marcus, Feb 27 2016
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Jason Earls, Jul 20 2004
STATUS
approved