OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
EXAMPLE
For p=3, 13*p^2+3*p+1=127; p=73, 13*p^2+3*p+1=69497
MATHEMATICA
Select[Prime[Range[1000]], PrimeQ[(13 #^2 + 3 # + 1)] &] (* Vincenzo Librandi, Oct 14 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(1500) | IsPrime(13*p^2+3*p+1)]; // Vincenzo Librandi, Oct 14 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Jan 21 2009
STATUS
approved