OFFSET
1,1
COMMENTS
The first 17 primes correspond to n from 0 to 16, which makes n^2 + 23n + 23 a prime-generating polynomial (see the link). This is a monic polynomial of the form n^2 + pn + p, where p is prime. Among the first 10^8 primes, only two more besides 23 give rise to prime-generating polynomials of this form. They are 8693 and 50983511 and they generate only 11 primes for n = 0 to 10.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Prime-Generating Polynomial
EXAMPLE
For n = 1, we have 1^2 + 23 * 1 + 23 = 47, which is prime, so 47 is in the sequence.
For n = 2, we have 2^2 + 23 * 2 + 23 = 4 + 46 + 23 = 73, which is prime, so 73 is in the sequence.
Contrast to n = 17, which gives us 17^2 + 23 * 17 + 23 = 289 + 391 + 23 = 703 = 19 * 37, so 703 is not in the sequence.
MAPLE
select(isprime, [seq(x^2+23*x+23, x=0..1000)]); # Robert Israel, Sep 18 2017
MATHEMATICA
Select[Range[0, 100]//#^2 + 23# + 23 &, PrimeQ]
PROG
(PARI) for(n=0, 100, isprime(n^2+23*n+23)&&print1(n^2+23*n+23 ", "))
(Magma) [a: n in [0..100] | IsPrime(a) where a is n^2+23*n+23 ]; // Vincenzo Librandi, Sep 23 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Waldemar Puszkarz, Sep 17 2017
STATUS
approved