OFFSET
1,2
LINKS
Rémy Sigrist, Table of n, a(n) for n = 1..8191
Rémy Sigrist, PARI program
EXAMPLE
A368427(67) = 107, so a(107) = 67.
PROG
(PARI) See Links section.
(Python)
from itertools import islice
from functools import reduce
def uniq(r): return reduce(lambda u, e: u if e in u else u+[e], r, [])
def A368427gen(): # generator of terms of A368427
R = [["1"]]
while R:
r = R.pop(0)
yield from map(lambda b: int(b, 2), r)
if len(r) > 1: R.append(uniq([r[k]+"0" for k in range(1, len(r))]))
R.append(uniq([r[0]+"0", r[0]+"1"] + [r[k]+"1" for k in range(1, len(r))]))
def agen(): # generator of terms
adict, n = dict(), 1
for i, k in enumerate(A368427gen(), 1):
if k not in adict:
adict[k] = i
while n in adict: yield adict[n]; n += 1
print(list(islice(agen(), 66))) # Michael S. Branicky, Dec 24 2023
CROSSREFS
KEYWORD
AUTHOR
Rémy Sigrist, Dec 24 2023
STATUS
approved