OFFSET
1,1
COMMENTS
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[n_]:=PrimeQ[n^2+(n+1)^2]; lst={}; Do[p=Prime[n]; If[f[p], AppendTo[lst, p]], {n, 6!}]; lst (* Vladimir Joseph Stephan Orlovsky, Oct 27 2009 *)
PROG
(PARI) forprime(p=2, 1e4, if(isprime(2*p*(p+1)+1), print1(p", "))) \\ Charles R Greathouse IV, Sep 09 2009
(Magma) [p: p in PrimesUpTo(2000) | IsPrime(2*p^2+2*p+1)]; // Vincenzo Librandi, Jun 22 2014
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Zak Seidov, Mar 20 2009
STATUS
approved