login
A362033
The indices where A362031(n) = 1.
2
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,2
COMMENTS
See A362031 for further details.
PROG
(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
CROSSREFS
Sequence in context: A268649 A316983 A136326 * A059973 A030035 A123431
KEYWORD
nonn,more
AUTHOR
Scott R. Shannon, 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
STATUS
approved