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,2

Comments

See A362031 for further details.

Links

Table of n, a(n) for n=1..32.

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

Cf. A362031, A001222.

Sequence in context: A268649 A316983 A136326 * A059973 A030035 A123431

Adjacent sequences: A362030 A362031 A362032 * A362034 A362035 A362036

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