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

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

Offset: 1

Author

Hugo Pfoertner, Jun 20 2023

Keywords

• nonn

Links

• Hugo Pfoertner, Table of n, a(n) for n = 1..2400

Crossrefs

Cf. A001481, A004018, A363522, A363761.

Numbers not occurring as terms of A077773.

Programs

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
A363762_list = list(islice(A363762_gen(),20)) # Chai Wah Wu, Jun 22 2023