The OEIS is supported by the many generous donors to the OEIS Foundation.
%I #13 Sep 08 2022 08:45:51
%S 7,73,271,367,379,397,739,823,859,1291,1297,1627,1747,2083,2203,2221,
%T 2269,2311,3067,3079,3109,3163,3463,3607,3733,3907,4093,4339,4603,
%U 4933,5077,5431,5527,5821,5869,6091,6301,6373,6421,6529,6883,7927,7951,8011,8467,8689,8737
%N Primes p such that p^6-p^3+1 is also prime.
%H Vincenzo Librandi, <a href="/A175205/b175205.txt">Table of n, a(n) for n = 1..1000</a>
%F A000040 INTERSECT A164989.
%t Select[Prime[Range[2000]], PrimeQ[#^6 - #^3 + 1]&] (* _Vincenzo Librandi_, Apr 16 2013 *)
%o (PARI) forprime(p=2,10000,if(isprime(p^6-p^3+1),print1(p,", ")))
%o (Magma) [p: p in PrimesUpTo(9000) | IsPrime(p^6-p^3+1)]; // _Vincenzo Librandi_, Apr 16 2013
%Y Cf. A060891, A175170, A000040, A164989.
%Y Subsequence of A002476.
%K nonn,easy
%O 1,1
%A _Vincenzo Librandi_, Dec 03 2010