%I #23 Jul 14 2021 01:59:52
%S 3,13,37,151,577,811,1453,1873,12697,14401,18253,27901,44101,75853,
%T 87121,93151,106033,151687,178753,188443,242173,291853,319057,333271,
%U 362953,410701,643453,666073,712891,787153,1040503,1379953,1742401
%N Primes of the form k^3 + k^2 + 1.
%C Primes in A098547. - _Michel Marcus_, Jan 21 2015
%H Harvey P. Dale, <a href="/A120479/b120479.txt">Table of n, a(n) for n = 1..1000</a>
%e 1^3 + 1^2 + 1 = 3 (prime), so 3 is in the sequence.
%e 2^3 + 2^2 + 1 = 13 (prime), so 13 is in the sequence.
%e 3^3 + 3^2 + 1 = 37 (prime), so 37 is in the sequence.
%e 4^3 + 4^2 + 1 = 81 = 3^4, so 81 is not in the sequence.
%e 5^3 + 5^2 + 1 = 151 (prime), so 151 is in the sequence.
%t Select[Table[n^3+n^2+1,{n,200}],PrimeQ] (* _Harvey P. Dale_, Oct 23 2020 *)
%o (PARI) for(n=1,10^3,if(isprime(p=n^3+n^2+1),print1(p,", "))) \\ _Derek Orr_, Jan 21 2015
%Y Cf. A119863 (corresponding k).
%Y Cf. A000040, A098547.
%K easy,nonn
%O 1,1
%A _Jonathan Vos Post_, Jul 21 2006
|