OFFSET
1,1
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MATHEMATICA
Select[Prime[Range[2000]], PrimeQ[#^6 - #^3 + 1]&] (* Vincenzo Librandi, Apr 16 2013 *)
PROG
(PARI) forprime(p=2, 10000, if(isprime(p^6-p^3+1), print1(p, ", ")))
(Magma) [p: p in PrimesUpTo(9000) | IsPrime(p^6-p^3+1)]; // Vincenzo Librandi, Apr 16 2013
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Dec 03 2010
STATUS
approved