cp's OEIS Frontend

A362033 The indices where A362031(n) = 1.

%I A362033 #27 Apr 09 2023 02:05:43
%S A362033 1,2,4,9,17,37,88,185,387,783,1611,3354,6959,14469,29983,62386,127992,
%T A362033 256022,513110,1029534,2069248,4161860,8380507,16880924,34031112,
%U A362033 68622510,138439065,279365954,563874347,1138322112,2298288343,4640647720
%N A362033 The indices where A362031(n) = 1.
%C A362033 See A362031 for further details.
%o A362033 (Python) # uses imports and generator in A362031
%o A362033 def agen():
%o A362033     yield from (n for n, an in enumerate(A362031gen(), 1) if an == 1)
%o A362033 print(list(islice(agen(), 20))) # _Michael S. Branicky_, Apr 06 2023
%o A362033 (Python)
%o A362033 from itertools import count, islice
%o A362033 from sympy import primeomega
%o A362033 def A362033_gen(): # generator of terms
%o A362033     a, b, c = {}, {}, 1
%o A362033     for n in count(1):
%o A362033         if c == 1: yield n
%o A362033         d = b[c] = b.get(c,primeomega(c))
%o A362033         c = a[d] = a.get(d,0)+1
%o A362033 A362033_list = list(islice(A362033_gen(),10)) # _Chai Wah Wu_, Apr 06 2023
%Y A362033 Cf. A362031, A001222.
%K A362033 nonn,more
%O A362033 1,2
%A A362033 _Scott R. Shannon_, Apr 06 2023
%E A362033 a(20)-a(28) from _Michael S. Branicky_, Apr 06 2023
%E A362033 a(29)-a(31) from _Chai Wah Wu_, Apr 06 2023
%E A362033 a(32) from _Chai Wah Wu_, Apr 08 2023