|
%I #22 Sep 08 2022 08:45:47
%S 6,7,9,42,49,54,66,73,81,88,96,100,106,118,120,121,126,175,190,196,
%T 207,208,234,235,240,271,301,312,322,342,343,367,378,379,381,384,397,
%U 399,415,418,423,442,444,480,505,513,520,556,574,609,628,664,666,700,702
%N Numbers n such that n^6 - n^3 + 1 is prime.
%H G. C. Greubel, <a href="/A164989/b164989.txt">Table of n, a(n) for n = 1..10000</a>
%t Select[Range[10^3], PrimeQ[#^6 - #^3 + 1] &] (* _Michael De Vlieger_, Dec 02 2017 *)
%o (PARI) for(k=1,300,if(isprime(k^6-k^3+1),print1(k, ", ")))
%o (Magma) [n: n in [0..1000] |IsPrime(n^3*(n^3-1)+1)] // _Vincenzo Librandi_, Dec 03 2010
%Y Cf. A153439, A175170.
%K easy,nonn
%O 1,1
%A _Serge Batalov_, Sep 03 2009
%E More terms from _Vincenzo Librandi_, Mar 25 2010
|
