%I #73 Aug 06 2023 04:08:11
%S 0,1,2,4,5,7,8,10,13,12,15,17,19,23,21,24,25,28,32,31,34,37,39,44,41,
%T 43,45,50,51,48,57,55,56,59,64,63,68,69,74,77,78,75,72,80,88,84,-1,94,
%U 89,96,93,99,97,102,108,-1,106,111,110,113,117,120,-1,123,133,127,130,137,142,138,139,-1,135
%N a(n) is the least k such that there are exactly n distinct numbers j that can be expressed as the sum of two squares with k^2 < j < (k+1)^2, or -1 if such a k does not exist.
%C Index of first occurrence of n in A077773 if there is any, otherwise -1. - _Rainer Rosenthal_, Jul 07 2023
%H Hugo Pfoertner, <a href="/A363763/a363763.gp.txt">PARI program for calculating a single term</a>, Jul 2023.
%H Rainer Rosenthal, <a href="/A363763/a363763.png">Illustrating a(5) = 7</a>.
%H Rainer Rosenthal, <a href="/A363763/a363763_1.png">First terms of A363763 illustrated</a>.
%F If a(n) != -1, then a(n) >= n/2. - _Chai Wah Wu_, Jun 22 2023
%F a(n) < (n+1)^2/2. - _Jon E. Schoenfield_ and _Chai Wah Wu_, Jun 24-26 2023
%e From _Rainer Rosenthal_, Jul 09 2023: (Start)
%e a(5) = 7, since A077773(7) = 5 and A077773(n) != 5 for n < 7.
%e a(46) = -1, since a(46) < ((46+1)^2)/2 < 1105 and A077773(k) != 46 for all k < 1105.
%e See illustrations in the links section. (End)
%o (PARI) \\ a4018(n) after _Michael Somos_
%o a4018(n) = if( n<1, n==0, 4 * sumdiv( n, d, (d%4==1) - (d%4==3)));
%o a363763 (upto) = {for (n=0, upto, my(kfound=-1); for (k=0, (n+1)^2\2+1, my(kp=k^2+1, km=(k+1)^2-1, m=0); for (j=kp, km, if (a4018(j), m++); if (m>n, break)); if (m==n, kfound=k; break)); print1 (kfound,", ");)};
%o a363763(75)
%o (Python)
%o from sympy import factorint
%o def A363763(n):
%o for k in range(n>>1,((n+1)**2<<1)+1):
%o c = 0
%o for m in range(k**2+1,(k+1)**2):
%o if all(p==2 or p&3==1 or e&1^1 for p, e in factorint(m).items()):
%o c += 1
%o if c>n:
%o break
%o if c==n:
%o return k
%o return -1 # _Chai Wah Wu_, Jun 20-26 2023
%Y Cf. A000415, A001481, A004018, A077773, A364341.
%Y A363762 gives the positions of terms = -1.
%Y Identical with A363761 up to a(11459) = 33864, but increasingly different afterwards, i.e., a(11460) = 34451, whereas A363761(11460) = -1.
%K sign
%O 0,3
%A _Hugo Pfoertner_, Jun 20 2023