%I #9 Oct 08 2019 12:31:15
%S 14,21,27,42,51,53,69,78,90,104,111,128,137,149,156,159,190,219,231,
%T 247,254,289,315,322,330,331,336,344,354,397,414,432,442,449,452,456,
%U 469,473,491,511,541,551,566,581,614,621,648,683,687,692,698,699,702,707,740,771,772,775,813,820,832
%N Numbers n for which n^9+n^8-1 is prime.
%H Robert Israel, <a href="/A125967/b125967.txt">Table of n, a(n) for n = 1..10000</a>
%p filter:= k -> isprime(k^9+k^8-1):
%p select(filter, [$1..1000]); # _Robert Israel_, Oct 07 2019
%t Do[If[PrimeQ[x^9 + x^8 - 1], Print[x]], {x, 1, 400}]
%o (PARI) is(n)=isprime(n^9+n^8-1) \\ _Charles R Greathouse IV_, May 15 2013
%Y Cf. A000040, A045546, A125881-A125885, A125965-A125973.
%K nonn
%O 1,1
%A _Artur Jasinski_, Dec 14 2006
%E More terms from _Robert Israel_, Oct 07 2019
|