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,1

Links

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

Example

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

Mathematica

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

Prog

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

Crossrefs

Cf. A049407, A126421.

Sequence in context: A228643 A304066 A296994 * A287458 A050572 A249668

Adjacent sequences: A236472 A236473 A236474 * A236476 A236477 A236478

Keyword

nonn

Author

Derek Orr, Jan 26 2014

Status

approved