OFFSET
1,1
COMMENTS
The primes p^2 - p + 1 are in A074268.
LINKS
T. D. Noe, Table of n, a(n) for n = 1..1000
MATHEMATICA
Select[Prime[Range[500]], PrimeQ[#^2-#+1]&] (* Harvey P. Dale, Oct 06 2015 *)
PROG
(PARI) { n=0; for (m=1, 10^9, p=prime(m); if (isprime(p^2 - p + 1), write("b065508.txt", n++, " ", p); if (n==1000, return)) ) } \\ Harry J. Smith, Oct 20 2009
(Magma) [n: n in [0..10000]| IsPrime(n) and IsPrime(n^2 - n + 1)] // Vincenzo Librandi, Aug 07 2010
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Vladeta Jovovic, Nov 26 2001
STATUS
approved