%I #115 Jul 18 2024 06:32:10
%S 0,1,1,3,1,3,3,7,1,3,3,7,3,7,7,15,1,3,3,7,3,7,7,15,3,7,7,15,7,15,15,
%T 31,1,3,3,7,3,7,7,15,3,7,7,15,7,15,15,31,3,7,7,15,7,15,15,31,7,15,15,
%U 31,15,31,31,63,1,3,3,7,3,7,7,15,3,7,7,15,7,15,15,31,3,7,7,15,7,15,15,31
%N a(n) = 2^A000120(n) - 1.
%C Essentially the same sequence as A001316, which has much more information, and also A159913. - _N. J. A. Sloane_, Jun 05 2009
%C Smallest number with same number of 1's in its binary expansion as n.
%C Fixed point of the morphism 0 -> 01, 1 -> 13, 3 -> 37, ... = k -> k, 2k+1, ... starting from a(0) = 0; 1 -> 01 -> 0113 -> 01131337 -> 011313371337377(15) -> ..., . - _Robert G. Wilson v_, Jan 24 2006
%C From _Gary W. Adamson_, Jun 04 2009: (Start)
%C As an infinite string, 2^n terms per row starting with "1": (1; 1,3; 1,3,3,7; 1,3,3,7,3,7,7,15; 1,3,3,7,3,7,7,15,3,7,7,15,7,15,15,31;...)
%C Row sums of that triangle = A027649: (1, 4, 14, 46, 454, ...); where the next row sum = current term of A027649 + next term in finite difference row of A027649, i.e., (1, 3, 10, 32, 100, 308, ...) = A053581. (End)
%C From _Omar E. Pol_, Jan 24 2016: (Start)
%C Partial sums give A267700.
%C a(n) is also the number of cells turned ON at n-th generation of the cellular automaton of A267700 in a 90-degree sector on the square grid.
%C a(n) is also the number of Y-toothpicks added at n-th generation of the structure of A267700 in a 120-degree sector on the triangular grid. (End)
%C Row sums of A090971. - _Nikolaos Pantelidis_, Nov 23 2022
%H Seiichi Manyama, <a href="/A038573/b038573.txt">Table of n, a(n) for n = 0..8191</a> (terms 0..1023 from T. D. Noe)
%H David Applegate, <a href="/A139250/a139250.anim.html">The movie version</a>
%H Michael Gilleland, <a href="/selfsimilar.html">Some Self-Similar Integer Sequences</a>
%H N. J. A. Sloane, <a href="/wiki/Catalog_of_Toothpick_and_CA_Sequences_in_OEIS">Catalog of Toothpick and Cellular Automata Sequences in the OEIS</a>
%H <a href="/index/Ce#cell">Index entries for sequences related to cellular automata</a>
%H <a href="/index/To#toothpick">Index entries for sequences related to toothpick sequences</a>
%H <a href="/index/Fi#FIXEDPOINTS">Index entries for sequences that are fixed points of mappings</a>
%F a(2n) = a(n), a(2n+1) = 2*a(n)+1, a(0) = 0. a(n) = A001316(n)-1 = 2^A000120(n) - 1. - _Daniele Parisse_
%F a(n) = number of positive integers k < n such that n XOR k = n-k (cf. A115378). - _Paul D. Hanna_, Jan 21 2006
%F a(n) = f(n, 1) with f(x, y) = if x = 0 then y - 1 else f(floor(x/2), y*(1 + x mod 2)). - _Reinhard Zumkeller_, Nov 21 2009
%F a(n) = (n mod 2 + 1) * a(floor(n/2)) + n mod 2. - _Reinhard Zumkeller_, Oct 10 2012
%F a(n) = Sum_{i=1..n} C(n,i) mod 2. - _Wesley Ivan Hurt_, Nov 17 2017
%F G.f.: -1/(1 - x) + Product_{k>=0} (1 + 2*x^(2^k)). - _Ilya Gutkovskiy_, Aug 20 2019
%F G.f. A(x) = x + x^2*A(x) + (1 + 2*x)*(1 - x^2)*A(x^2). - _Michael Somos_, Jul 24 2023
%e 9 = 1001 -> 0011 -> 3, so a(9)=3.
%e From _Gary W. Adamson_, Jun 04 2009: (Start)
%e Triangle read by rows:
%e 1;
%e 1, 3;
%e 1, 3, 3, 7;
%e 1, 3, 3, 7, 3, 7, 7, 15;
%e 1, 3, 3, 7, 3, 7, 7, 15, 3, 7, 7, 15, 7, 15, 15, 31;
%e ...
%e Row sums: (1, 4, 14, 46, ...) = A027649 = last row terms + new set of terms such that row 3 = (1, 3, 3, 7,) + (3, 7, 7, 15) = 14 + 32 = A027649(2) + A053581(3). (End)
%e The rows of this triangle converge to A159913. - _N. J. A. Sloane_, Jun 05 2009
%e G.f. = x + x^2 + 3*x^3 + x^4 + 3*x^5 + 3*x^6 + 7*x^7 + x^8 + 3*x^9 + 3*x^10 + 7*x^11 + ... - _Michael Somos_, Jul 24 2023
%p seq(2^convert(convert(n,base,2),`+`)-1, n=0..100); # _Robert Israel_, Jan 24 2016
%t Array[ 2^Count[ IntegerDigits[ #, 2 ], 1 ]-1&, 100 ]
%t Nest[ Flatten[ # /. a_Integer -> {a, 2a + 1}] &, {0}, 7] (* _Robert G. Wilson v_, Jan 24 2006 *)
%o (PARI) {a(n) = 2^subst(Pol(binary(n)), x, 1) - 1};
%o (PARI) a(n) = 2^hammingweight(n)-1; \\ _Michel Marcus_, Jan 24 2016
%o (Haskell)
%o a038573 0 = 0
%o a038573 n = (m + 1) * (a038573 n') + m where (n', m) = divMod n 2
%o -- _Reinhard Zumkeller_, Oct 10 2012, Feb 07 2011
%o (Python 3.10+)
%o def A038573(n): return (1<<n.bit_count())-1 # _Chai Wah Wu_, Nov 15 2022
%Y Cf. A007313, A115378, A073138.
%Y This is Guy Steele's sequence GS(3, 6) (see A135416).
%Y Cf. also A000079, A001316, A027649, A053581, A159913, A267700.
%Y Write n in b-ary, sort digits into increasing order: this sequence (b=2), A038574 (b=3), A319652 (b=4), A319653 (b=5), A319654 (b=6), A319655 (b=7), A319656 (b=8), A319657 (b=9), A004185 (b=10).
%Y Column k=0 of A340666.
%K nonn,easy,nice
%O 0,4
%A _Marc LeBrun_
%E More terms from _Erich Friedman_
%E New definition from _N. J. A. Sloane_, Mar 01 2008