A236476 Primes p such that p^3 + p - 1 is prime.

3, 7, 43, 73, 103, 109, 127, 139, 151, 199, 223, 241, 283, 313, 367, 379, 421, 541, 631, 661, 733, 739, 751, 769, 829, 991, 1117, 1129, 1201, 1231, 1249, 1297, 1303, 1429, 1471, 1663, 1669, 1693, 1699, 1741, 1789, 1867, 1933

Offset

1,1

Comments

Primes in the sequence A236475.

Links

Daniel Starodubtsev, Table of n, a(n) for n = 1..10000

Example

241 is prime and 241^3 + 241 - 1 = 13997761 is prime.

Prog

(Python)
import sympy
from sympy import isprime
{print(n) for n in range(10**4) if isprime(n) and isprime(n**3+n-1)}
(PARI)
s=[]; forprime(p=2, 2000, if(isprime(p^3+p-1), s=concat(s, p))); s \\ Colin Barker, Jan 27 2014

Crossrefs

Cf. A236475.

Sequence in context: A257366 A141304 A213893 * A282178 A101208 A050633

Adjacent sequences: A236473 A236474 A236475 * A236477 A236478 A236479

Keyword

nonn

Author

Derek Orr, Jan 26 2014

Status

approved