%I #49 May 24 2022 17:50:59
%S 3,6,12,13,24,26,48,52,53,96,104,106,113,192,208,212,213,226,227,384,
%T 416,424,426,452,453,454,768,832,848,852,853,904,906,908,909,1536,
%U 1664,1696,1704,1706,1808,1812,1813,1816,1818,3072,3328,3392,3408,3412,3413,3616
%N Numbers with 3 odd integers in their Collatz (or 3x+1) trajectory.
%C The Collatz (or 3x+1) function is f(x) = x/2 if x is even, 3x+1 if x is odd (A006370).
%C The Collatz trajectory of n is obtained by applying f repeatedly to n until 1 is reached.
%C A078719(a(n)) = 3; A006667(a(n)) = 2.
%D J. R. Goodwin, Results on the Collatz Conjecture, Sci. Ann. Comput. Sci. 13 (2003) pp. 1-16
%D J. Shallit and D. Wilson, The "3x+1" Problem and Finite Automata, Bulletin of the EATCS #46 (1992) pp. 182-185.
%H Reinhard Zumkeller and David A. Corneth, <a href="/A062053/b062053.txt">Table of n, a(n) for n = 1..16191</a> (first 250 terms from Reinhard Zumkeller, terms < 10^25)
%H J. Shallit and D. Wilson, <a href="http://www.cs.uwaterloo.ca/~shallit/Papers/wilson.ps">The "3x+1" Problem and Finite Automata</a>, Bulletin of the EATCS #46 (1992) pp. 182-185.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/CollatzProblem.html">Collatz Problem</a>
%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Collatz_conjecture">Collatz conjecture</a>
%H <a href="/index/3#3x1">Index entries for sequences related to 3x+1 (or Collatz) problem</a>
%H <a href="/index/Ar#2-automatic">Index entries for 2-automatic sequences</a>.
%F The two formulas giving this sequence are listed in Corollary 3.1 and Corollary 3.2 in J. R. Goodwin with the following caveats: the value x cannot equal zero in Corollary 3.2, one must multiply the formulas by all powers of 2 (2^1, 2^2, ...) to get the evens. - _Jeffrey R. Goodwin_, Oct 26 2011
%e The Collatz trajectory of 3 is (3,10,5,16,8,4,2,1), which contains 3 odd integers.
%t Collatz[n_?OddQ] := (3n + 1)/2; Collatz[n_?EvenQ] := n/2; oddIntCollatzCount[n_] := Length[Select[NestWhileList[Collatz, n, # != 1 &], OddQ]]; Select[Range[4000], oddIntCollatzCount[#] == 3 &] (* _Alonso del Arte_, Oct 28 2011 *)
%o (Haskell)
%o import Data.List (elemIndices)
%o a062053 n = a062053_list !! (n-1)
%o a062053_list = map (+ 1) $ elemIndices 3 a078719_list
%o -- _Reinhard Zumkeller_, Oct 08 2011
%Y Cf. A006370, A078719, A006667.
%Y Cf. A000079, A062052, A062054, A062055, A062056, A062057, A062058, A062059, A062060.
%Y Cf. A198584 (this sequence without the even numbers).
%Y See also A198587.
%Y Column k=3 of A354236.
%K nonn,easy
%O 1,1
%A _David W. Wilson_