|
|
A363762
|
|
Numbers k for which A363763(k) = -1.
|
|
7
|
|
|
46, 55, 62, 71, 80, 86, 107, 130, 172, 187, 195, 208, 222, 247, 259, 263, 268, 272, 280, 297, 314, 330, 358, 363, 370, 372, 379, 394, 400, 405, 429, 449, 450, 462, 489, 500, 529, 534, 587, 629, 641, 646, 652, 667, 668, 672, 704, 715, 733, 736, 749, 769, 775, 776, 778, 785, 793, 799
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
PROG
|
(Python)
from itertools import count, islice
from sympy import factorint
def A363762_gen(startvalue=1): # generator of terms >= startvalue
for n in count(max(startvalue, 1)):
for k in range(n>>1, ((n+1)**2<<1)+1):
c = 0
for m in range(k**2+1, (k+1)**2):
if all(p==2 or p&3==1 or e&1^1 for p, e in factorint(m).items()):
c += 1
if c>n:
break
if c==n:
break
else:
yield n
|
|
CROSSREFS
|
Numbers not occurring as terms of A077773.
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|