|
%I #18 Sep 08 2022 08:45:04
%S 2,3,7,13,67,79,139,151,163,193,337,349,379,457,541,613,643,727,769,
%T 919,991,1021,1093,1117,1201,1231,1381,1423,1549,1567,1597,1621,1693,
%U 1747,1789,1801,1933,1987,2011,2017,2113,2137,2143,2239,2281,2557,2647
%N Primes p such that p^2 - p + 1 is prime.
%C The primes p^2 - p + 1 are in A074268.
%H T. D. Noe, <a href="/A065508/b065508.txt">Table of n, a(n) for n = 1..1000</a>
%t Select[Prime[Range[500]],PrimeQ[#^2-#+1]&] (* _Harvey P. Dale_, Oct 06 2015 *)
%o (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
%o (Magma) [n: n in [0..10000]| IsPrime(n) and IsPrime(n^2 - n + 1)] // _Vincenzo Librandi_, Aug 07 2010
%Y Cf. A053182, A053184, A091567, A074268.
%K easy,nonn
%O 1,1
%A _Vladeta Jovovic_, Nov 26 2001
|
