A065508 Primes p such that p^2 - p + 1 is prime.
2, 3, 7, 13, 67, 79, 139, 151, 163, 193, 337, 349, 379, 457, 541, 613, 643, 727, 769, 919, 991, 1021, 1093, 1117, 1201, 1231, 1381, 1423, 1549, 1567, 1597, 1621, 1693, 1747, 1789, 1801, 1933, 1987, 2011, 2017, 2113, 2137, 2143, 2239, 2281, 2557, 2647
Offset: 1
Links
- T. D. Noe, Table of n, a(n) for n = 1..1000
Programs
-
Magma
[n: n in [0..10000]| IsPrime(n) and IsPrime(n^2 - n + 1)] // Vincenzo Librandi, Aug 07 2010
-
Mathematica
Select[Prime[Range[500]],PrimeQ[#^2-#+1]&] (* Harvey P. Dale, Oct 06 2015 *)
-
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
Comments