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

Author

Derek Orr, Jan 26 2014

Keywords

• nonn

Comments

Primes in the sequence A236475.

Examples

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

Links

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

Crossrefs

Cf. A236475.

Programs

PARI

s=[]; forprime(p=2, 2000, if(isprime(p^3+p-1), s=concat(s, p))); s \\ Colin Barker, Jan 27 2014

Python

import sympy
from sympy import isprime
{print(n) for n in range(10**4) if isprime(n) and isprime(n**3+n-1)}