login
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.
8

%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