A088796 Primes p such that p^3 + p^2 - 1 is also prime.

2, 5, 11, 19, 29, 37, 47, 71, 89, 97, 107, 127, 139, 179, 281, 317, 331, 337, 379, 397, 461, 499, 509, 571, 617, 659, 677, 691, 709, 839, 877, 887, 907, 1021, 1049, 1051, 1087, 1097, 1229, 1279, 1297, 1307, 1381, 1429, 1447, 1481, 1489, 1499, 1511, 1579

Offset

1,1

Links

Robert Israel, Table of n, a(n) for n = 1..10000

Example

a(5) = 29 because (29^2 + 29^3 - 1) = (841 + 24389 - 1) = 25229, which is prime.

Maple

filter:= n -> isprime(n) and isprime(n^2 + n^3 - 1):
select(filter, [2, seq(i, i=3..10000, 2)]); # Robert Israel, Jan 29 2017

Mathematica

Select[Prime[Range[300]], PrimeQ[#^3+#^2-1]&] (* Harvey P. Dale, Feb 07 2017 *)

Prog

(PARI) isok(n) = isprime(n) && isprime(n^2+n^3-1) \\ Michel Marcus, Jul 29 2013

Crossrefs

Sequence in context: A031870 A040147 A152471 * A190309 A161550 A245207

Adjacent sequences: A088793 A088794 A088795 * A088797 A088798 A088799

Keyword

nonn

Author

Chuck Seggelin (barkeep(AT)plastereddragon.com), Oct 16 2003

Status

approved