cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

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

This page as a plain text file.
%I A236475 #15 Sep 02 2022 17:41:33
%S A236475 3,4,7,10,15,16,18,21,25,27,33,36,39,43,46,51,52,55,63,73,78,81,87,93,
%T A236475 94,96,100,103,105,109,112,115,117,120,124,127,129,135,139,145,150,
%U A236475 151,165,166,171,178,189,192,198,199
%N A236475 Numbers k such that k^3 + k - 1 is prime.
%H A236475 Daniel Starodubtsev, <a href="/A236475/b236475.txt">Table of n, a(n) for n = 1..10000</a>
%e A236475 46^3 + 46 - 1 = 97381 is prime. So 46 is a member of this sequence.
%t A236475 Select[Range[200],PrimeQ[#^3+#-1]&] (* _Harvey P. Dale_, Sep 02 2022 *)
%o A236475 (Python)
%o A236475 import sympy
%o A236475 from sympy import isprime
%o A236475 {print(n) for n in range(10**3) if isprime(n**3+n-1)}
%o A236475 (PARI)
%o A236475 s=[]; for(n=1, 500, if(isprime(n^3+n-1), s=concat(s, n))); s \\ _Colin Barker_, Jan 27 2014
%Y A236475 Cf. A049407, A126421.
%K A236475 nonn
%O A236475 1,1
%A A236475 _Derek Orr_, Jan 26 2014