|
%I #21 Sep 08 2022 08:46:13
%S 33,34,37,42,49,50,51,53,56,58,60,65,67,68,69,71,72,75,78,82,83,84,85,
%T 86,88,91,94,95,97,100,101,102,105,106,107,110,111,113,114,116,117,
%U 118,119,122,123,124,128,129,132,133,134,135,136,139,141,143,148,151,152,153
%N Numbers n>0 for which n+(17-n)^2 is not prime.
%C Motivated by the fact that n+(17-n)^2 = 1+16^2, 2+15^2, ..., 16+1^2, 17+0^2, 18+1^2, 19+2^2, ..., 32+15^2 are all prime. This has an explanation through Heegener numbers, similar to Euler's prime-generating polynomial, cf. A002837 and related crossrefs.
%H Robert Israel, <a href="/A260678/b260678.txt">Table of n, a(n) for n = 1..10000</a>
%p remove(t -> isprime(t+(17-t)^2), [$1..200]); # _Robert Israel_, May 02 2017
%t Select[Range[200], !PrimeQ[# + (17 - #)^2] &] (* _Vincenzo Librandi_, Nov 16 2015 *)
%o (PARI) for(n=1,999,isprime(n+(17-n)^2)||print1(n","))
%o (Magma) [n: n in [1..180] | not IsPrime(n+(17-n)^2)]; // _Vincenzo Librandi_, Nov 16 2015
%Y Cf. A260679 (n+(17-n)^2), A007635 (primes in that sequence = primes of the form n^2+n+17).
%Y Cf. A002837 (n^2-n+41 is prime), A005846 (primes of form n^2+n+41), A007634 (n^2+n+41 is composite), A097823 (n^2+n+41 is not squarefree).
%K nonn,easy
%O 1,1
%A _M. F. Hasler_, Nov 15 2015
|
