login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A100662
Primes of the form n^3 + (n+1)^2.
6
5, 17, 43, 89, 593, 829, 2393, 2969, 3631, 5237, 11177, 12743, 14449, 16301, 27961, 33857, 40529, 44171, 56393, 60919, 75937, 93241, 127601, 198593, 219721, 266369, 278981, 499439, 578843, 621521, 689393, 737281, 787337, 839609, 950993, 980299
OFFSET
1,1
LINKS
FORMULA
A100705 INTERSECT A000040.
EXAMPLE
a(1)=5 because 5=1^3+2^2;
a(2)=17 because 17=2^3+3^2;
a(3)=43 because 43=3^3+4^2.
MATHEMATICA
Select[Table[n^3+(n+1)^2, {n, 200}], PrimeQ] (* Vladimir Joseph Stephan Orlovsky, Feb 01 2012 *)
PROG
(Magma) [ a: n in [0..100] | IsPrime(a) where a is n^3 + (n+1)^2 ]; // Vincenzo Librandi, Jul 18 2012
(PARI) list(lim)=my(v=List(), t); for(n=1, lim, t=n^3+(n+1)^2; if(t>lim, break); if(isprime(t), listput(v, t))); Vec(v) \\ Charles R Greathouse IV, Dec 23 2016
CROSSREFS
Sequence in context: A146720 A146640 A100705 * A168617 A046916 A089527
KEYWORD
nonn,easy
AUTHOR
Giovanni Teofilatto, Jan 03 2005
EXTENSIONS
More terms from Mark Hudson, Jan 04 2005
STATUS
approved