|
%I #17 Jan 01 2024 09:08:23
%S 1,1,4,1,36,1,4,9,4,1,36,1,4,9,4,1,36,1,4,9,16,1,36,49,4,81,4,1,144,1,
%T 16,9,4,9,36,1,4,9,4,1,576,1,4,9,16,1,36,25,4,9,16,1,36,25,4,81,4,1,
%U 324,1,36,9,4,9,36,1,4,81,4,1,36,1,16,9,4,25,36,1,4,9,16,1,144,25,4,81
%N a(n) = smallest nonzero square x^2 such that n+x^2 is prime.
%C a(n) = 1 for n in A006093. - _Robert Israel_, Dec 31 2023
%H Robert Israel, <a href="/A059844/b059844.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) + n is the smallest prime of the form x^2 + n.
%e a(24) = 49 because 49 + 24 = 73 is prime and 1 + 24 = 25, 4 + 24 = 28, 9 + 24 = 33, 16 + 24 = 40, 25 + 24 = 49, and 36 + 24 = 60 are composite.
%p f:= proc(n) local x;
%p for x from 1 + (n mod 2) by 2 do
%p if isprime(n+x^2) then return x^2 fi;
%p od
%p end proc:
%p f(1):= 1:
%p map(f, [$1..100]); # _Robert Israel_, Dec 31 2023
%t sqs[n_]:=Module[{q=1},While[!PrimeQ[n+q],q=(Sqrt[q]+1)^2];q]; Array[ sqs,90] (* _Harvey P. Dale_, Aug 11 2017 *)
%Y Cf. A002496, A006093, A056899, A049423, A005473, A056905, A056909.
%K nonn
%O 1,3
%A _Labos Elemer_, Feb 26 2001
|
