|
%I #12 Feb 07 2017 13:50:51
%S 2,5,11,19,29,37,47,71,89,97,107,127,139,179,281,317,331,337,379,397,
%T 461,499,509,571,617,659,677,691,709,839,877,887,907,1021,1049,1051,
%U 1087,1097,1229,1279,1297,1307,1381,1429,1447,1481,1489,1499,1511,1579
%N Primes p such that p^3 + p^2 - 1 is also prime.
%H Robert Israel, <a href="/A088796/b088796.txt">Table of n, a(n) for n = 1..10000</a>
%e a(5) = 29 because (29^2 + 29^3 - 1) = (841 + 24389 - 1) = 25229, which is prime.
%p filter:= n -> isprime(n) and isprime(n^2 + n^3 - 1):
%p select(filter, [2, seq(i,i=3..10000,2)]); # _Robert Israel_, Jan 29 2017
%t Select[Prime[Range[300]],PrimeQ[#^3+#^2-1]&] (* _Harvey P. Dale_, Feb 07 2017 *)
%o (PARI) isok(n) = isprime(n) && isprime(n^2+n^3-1) \\ _Michel Marcus_, Jul 29 2013
%K nonn
%O 1,1
%A Chuck Seggelin (barkeep(AT)plastereddragon.com), Oct 16 2003
|
