|
%I #18 Feb 26 2023 17:41:46
%S 1,2,3,4,5,7,8,9,10,11,13,15,16,17,18,19,20,21,22,23,25,26,27,29,31,
%T 32,33,34,35
%N Records in A261461.
%C a(n) = A261461(A261467(n)).
%o (Haskell)
%o a261466 n = a261466_list !! (n-1)
%o (a261466_list, a261467_list) = unzip $ (0, 1) : f 0 1 where
%o f i x | y > x = (y, i) : f (i + 1) y
%o | otherwise = f (i + 1) x
%o where y = a261461 i
%o (Python)
%o from itertools import count, islice
%o def f(n):
%o b, k = bin(n)[2:], 1
%o return next(k for k in count(1) if bin(k)[2:] not in b)
%o def agen(record=-1):
%o yield from (((record:=f(k)), k)[0] for k in count(0) if f(k) > record)
%o print(list(islice(agen(), 18))) # _Michael S. Branicky_, Feb 26 2023
%Y Cf. A261461, A261467.
%K nonn,more,base
%O 1,2
%A _Reinhard Zumkeller_, Aug 30 2015
%E Initial terms corrected and more terms from _Rémy Sigrist_, Mar 10 2018
%E a(26)-a(29) from _Rémy Sigrist_, Feb 26 2023
|
