login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A219696 Numbers k such that the trajectory of 3k + 1 under the '3x + 1' map reaches k. 2

%I #44 Sep 05 2023 14:45:20

%S 1,2,4,8,10,14,16,20,22,26,40,44,52,106,184,206,244,274,322,526,650,

%T 668,790,866,976,1154,1300,1438,1732,1780,1822,2308,2734,3238,7288

%N Numbers k such that the trajectory of 3k + 1 under the '3x + 1' map reaches k.

%C This sequence seems complete; there are no other terms <= 10^9. - _T. D. Noe_, Dec 03 2012

%C If the 3x+1 step is replaced with (3x+1)/2, the sequence becomes {1, 2, 4, 8, 10, 14, 20, 22, 26, 40, 44, 206, 244, 650, 668, 866, 1154, 1822, 2308, ...}. - _Robert G. Wilson v_, Jan 13 2015

%C From _Andrew Slattery_, Aug 03 2023: (Start)

%C For most terms k, the trajectory of 3k + 1 reaches 310 or the trajectory of 310 reaches k.

%C For the rest of the terms k, the trajectory of 3k + 1 reaches 22 or the trajectory of 22 reaches k.

%C With the exception of k = 1, k is reached after S steps,

%C where S = c*8 + d*13 + e*44 + f*75, with c, d, e and f in {0, 1, 2}; in particular, S is in {8, 13, 8+8, 8+13, 13+13, 44, 75, 44+44, 75+13+13, 75+44, 75+75}. (End)

%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>

%e For k = 4, the Collatz trajectory of 3k + 1 is (13, 40, 20, 10, 5, 16, 8, 4, 2, 1), which includes 4; thus, 4 is in the sequence.

%e For k = 5, the Collatz trajectory of 3k + 1 is (16, 8, 4, 2, 1), which does not include 5; thus, 5 is not in the sequence.

%t Collatz[n_] := NestWhileList[If[EvenQ[#], #/2, 3 # + 1] &, n, # > 1 &]; Select[Range[10000], MemberQ[Collatz[3 # + 1], #] &] (* _T. D. Noe_, Dec 03 2012 *)

%o (Haskell)

%o a219696 n = a219696_list !! (n-1)

%o a219696_list = filter (\x -> collatz'' x == x) [1..] where

%o collatz'' x = until (`elem` [1, x]) a006370 (3 * x + 1)

%o -- _Reinhard Zumkeller_, Aug 11 2014

%o (Python)

%o def ok(n):

%o if n==1: return [1]

%o N=3*n + 1

%o l=[N, ]

%o while True:

%o if N%2==1: N = 3*N + 1

%o else: N/=2

%o l+=[N, ]

%o if N<2: break

%o if n in l: return 1

%o return 0 # _Indranil Ghosh_, Apr 22 2017

%o (PARI) a006370(n) = if(n%2==0, n/2, 3*n+1)

%o is(n) = my(x=3*n+1); while(1, x=a006370(x); if(x==n, return(1), if(x==1, return(0)))) \\ _Felix Fröhlich_, Jun 10 2021

%Y Cf. A014682, A070991, A006370, A070165.

%K nonn,nice,more

%O 1,2

%A _Robert C. Lyons_, Nov 25 2012

%E Initial 1 from Clark R. Lyons, Dec 02 2012

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 20:05 EDT 2024. Contains 371254 sequences. (Running on oeis4.)