A362062 First index in A362061 where A001221(A362061(a(n))) = n.

1, 3, 10, 52, 398, 4734, 78428, 1326821, 26386288, 645482585

Offset

0,2

Links

Table of n, a(n) for n=0..9.

Prog

(Python)
from sympy import primefactors
from collections import Counter
from itertools import count, islice
def agen(): # generator of terms
    an, c, d, d2, k = 1, Counter(), dict(), dict(), 0
    for n in count(1):
        if an in d: dpf = d[an]
        else: dpf, d2[dpf] = len(primefactors(an)), n
        c[dpf] += 1
        an = c[dpf]
        if k in d2: yield d2[k]; k += 1
print(list(islice(agen(), 7))) # Michael S. Branicky, Apr 06 2023
(Python)
from itertools import count
from sympy import primenu
def A362062(n):
    a, b, c = {}, {}, 1
    for m in count(1):
        d = b[c] = b.get(c, primenu(c))
        if d == n: return m
        c = a[d] = a.get(d, 0)+1 # Chai Wah Wu, Apr 07 2023

Crossrefs

Cf. A362061, A001221.

Sequence in context: A309885 A363209 A052446 * A373284 A290489 A002873

Adjacent sequences: A362059 A362060 A362061 * A362063 A362064 A362065

Keyword

nonn,more

Author

Scott R. Shannon, Apr 06 2023

Extensions

a(9) from Chai Wah Wu, Apr 07 2023

Status

approved