A362033 The indices where A362031(n) = 1.

1, 2, 4, 9, 17, 37, 88, 185, 387, 783, 1611, 3354, 6959, 14469, 29983, 62386, 127992, 256022, 513110, 1029534, 2069248, 4161860, 8380507, 16880924, 34031112, 68622510, 138439065, 279365954, 563874347, 1138322112, 2298288343, 4640647720

Offset: 1

Author

Scott R. Shannon, Apr 06 2023

Keywords

• nonn
• more

Comments

See A362031 for further details.

Crossrefs

Cf. A362031, A001222.

Programs

Python

# uses imports and generator in A362031
def agen():
    yield from (n for n, an in enumerate(A362031gen(), 1) if an == 1)
print(list(islice(agen(), 20))) # Michael S. Branicky, Apr 06 2023

Python

from itertools import count, islice
from sympy import primeomega
def A362033_gen(): # generator of terms
    a, b, c = {}, {}, 1
    for n in count(1):
        if c == 1: yield n
        d = b[c] = b.get(c,primeomega(c))
        c = a[d] = a.get(d,0)+1
A362033_list = list(islice(A362033_gen(),10)) # Chai Wah Wu, Apr 06 2023

Extensions

a(20)-a(28) from Michael S. Branicky, Apr 06 2023

a(29)-a(31) from Chai Wah Wu, Apr 06 2023

a(32) from Chai Wah Wu, Apr 08 2023