A236475 Numbers k such that k^3 + k - 1 is prime.

3, 4, 7, 10, 15, 16, 18, 21, 25, 27, 33, 36, 39, 43, 46, 51, 52, 55, 63, 73, 78, 81, 87, 93, 94, 96, 100, 103, 105, 109, 112, 115, 117, 120, 124, 127, 129, 135, 139, 145, 150, 151, 165, 166, 171, 178, 189, 192, 198, 199

Offset: 1

Author

Derek Orr, Jan 26 2014

Keywords

• nonn

Examples

			46^3 + 46 - 1 = 97381 is prime. So 46 is a member of this sequence.
		

Links

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

Crossrefs

Cf. A049407, A126421.

Programs

Mathematica

Select[Range[200],PrimeQ[#^3+#-1]&] (* Harvey P. Dale, Sep 02 2022 *)

PARI

s=[]; for(n=1, 500, if(isprime(n^3+n-1), s=concat(s, n))); s \\ Colin Barker, Jan 27 2014

Python

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