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”).
%I #45 Sep 08 2022 08:46:19
%S 23,47,73,101,131,163,197,233,271,311,353,397,443,491,541,593,647,761,
%T 821,883,947,1013,1151,1223,1297,1373,1451,1531,1613,1697,1783,1871,
%U 2053,2243,2341,2441,2543,2647,2753,2861,2971,3083,3313,3673,3797,3923,4051,4447
%N Primes of the form n^2 + 23*n + 23.
%C 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.
%H Robert Israel, <a href="/A292509/b292509.txt">Table of n, a(n) for n = 1..10000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/Prime-GeneratingPolynomial.html">Prime-Generating Polynomial</a>
%e For n = 1, we have 1^2 + 23 * 1 + 23 = 47, which is prime, so 47 is in the sequence.
%e For n = 2, we have 2^2 + 23 * 2 + 23 = 4 + 46 + 23 = 73, which is prime, so 73 is in the sequence.
%e 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.
%p select(isprime, [seq(x^2+23*x+23, x=0..1000)]); # _Robert Israel_, Sep 18 2017
%t Select[Range[0, 100]//#^2 + 23# + 23 &, PrimeQ]
%o (PARI) for(n=0, 100, isprime(n^2+23*n+23)&&print1(n^2+23*n+23 ","))
%o (Magma) [a: n in [0..100] | IsPrime(a) where a is n^2+23*n+23 ]; // _Vincenzo Librandi_, Sep 23 2017
%Y Cf. A005846, A300473.
%K nonn
%O 1,1
%A _Waldemar Puszkarz_, Sep 17 2017