OFFSET
1,1
MATHEMATICA
lst={}; Do[p=Prime[n]; r=Sqrt[p-1]-1; If[PrimeQ[r], AppendTo[lst, p]], {n, 4*8!}]; lst
PROG
(Magma) [a: p in PrimesUpTo(950) | IsPrime(a) where a is p^2+2*p+2]; // Vincenzo Librandi, Dec 20 2010
(PARI) forprime(p=2, 1000, m=p^2+2*p+2; if(isprime(m), print1(m, ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Mar 01 2009
STATUS
approved