A174207 Depleted stream of the natural numbers written in base 2: delete even occurrences of digit 0 and odd occurrences of digit 1.
0, 1, 1, 0, 1, 0, 1, 1, 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 0, 0, 1, 0, 0
Offset: 0
Programs
-
Maple
nb := 100 ; L := [0] ; for n from 1 to nb do nn := ListTools[Reverse](convert(n,base,2)) ; L := [op(L),op(nn)] ; end do; Lout := [] ; ie :=1 ; io :=1 ; for i from 1 to nops(L) do if op(i,L) =0 then if ie>0 then if type(ie,'odd') then printf("%d,",0) ; end if; end if; ie := ie+1 ; else if io>0 then if type(io,'even') then printf("%d,",1) ; end if; end if; io := io+1 ; end if; end do: # R. J. Mathar, Dec 06 2011
-
Mathematica
nMax = 60; bits = Join @@ Table[IntegerDigits[n, 2], {n, 0, nMax}]; pos0 = Position[bits, 0] // Flatten; even0 = Partition[pos0, 2][[All, 2]]; bits2 = ReplacePart[bits, Alternatives @@ even0 -> "~"]; pos1 = Position[bits2, 1] // Flatten; odd1 = Partition[pos1, 2][[All, 1]]; DeleteCases[ ReplacePart[ bits2, Alternatives @@ odd1 -> "~"], "~"] (* Jean-François Alcover, Nov 18 2016 *)
Comments