login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

Numbers k for which A363763(k) = -1.
7

%I #26 Jul 13 2023 07:26:50

%S 46,55,62,71,80,86,107,130,172,187,195,208,222,247,259,263,268,272,

%T 280,297,314,330,358,363,370,372,379,394,400,405,429,449,450,462,489,

%U 500,529,534,587,629,641,646,652,667,668,672,704,715,733,736,749,769,775,776,778,785,793,799

%N Numbers k for which A363763(k) = -1.

%H Hugo Pfoertner, <a href="/A363762/b363762.txt">Table of n, a(n) for n = 1..2400</a>

%o (Python)

%o from itertools import count, islice

%o from sympy import factorint

%o def A363762_gen(startvalue=1): # generator of terms >= startvalue

%o for n in count(max(startvalue,1)):

%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 break

%o else:

%o yield n

%o A363762_list = list(islice(A363762_gen(),20)) # _Chai Wah Wu_, Jun 22 2023

%Y Cf. A001481, A004018, A363522, A363761.

%Y Numbers not occurring as terms of A077773.

%K nonn

%O 1,1

%A _Hugo Pfoertner_, Jun 20 2023