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!)
A078440 Numbers n with property that n is not a power of 2 and the finite sequence n, f(n), f(f(n)), ...., 1 in the Collatz (or 3x + 1) problem contains exactly one prime. (The earliest "1" is meant.) 2

%I #21 Mar 19 2019 13:25:57

%S 21,42,84,85,168,170,336,340,341,453,672,680,682,906,909,1344,1360,

%T 1364,1365,1812,1813,1818,2688,2720,2728,2730,3624,3626,3636,5376,

%U 5440,5456,5460,5461,7248,7252,7272,7281,9669

%N Numbers n with property that n is not a power of 2 and the finite sequence n, f(n), f(f(n)), ...., 1 in the Collatz (or 3x + 1) problem contains exactly one prime. (The earliest "1" is meant.)

%C f(n) = n/2 if n is even, = 3n + 1 if n is odd. Powers 2^n trivially have exactly one prime in n, f(n), f(f(n)), ..., 2, 1, namely 2 and so are excluded from the sequence.

%C A055509(a(n)) = 0; A078350(a(n)) <= 1.

%H Reinhard Zumkeller, <a href="/A078440/b078440.txt">Table of n, a(n) for n = 1..1000</a>

%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 n, f(n), f(f(n)), .... for n = 21 is: 21, 64, 32, 16, 8, 4, 2, 1, which has exactly one prime, that is, 2. Hence 21 belongs to the sequence.

%t f[n_] := n/2 /; Mod[n, 2] == 0 f[n_] := 3 n + 1 /; Mod[n, 2] == 1 g[n_] := Module[{i, p}, i = n; p = 0; While[i > 1, If[PrimeQ[i], p = p + 1]; i = f[i]]; p]; Select[Range[10^4], g[ # ] == 1 && ! IntegerQ[Log[2, # ]] &]

%t pQ[n_]:=Count[NestWhileList[If[EvenQ[#],#/2,3#+1]&,n, #>1&], _?PrimeQ] == 1; With[ {nn=10000},Complement[Select[Range[nn],pQ],2^Range[Floor[ Log[ 2,nn]]]]] (* _Harvey P. Dale_, Oct 19 2011 *)

%o (Haskell)

%o a078440 n = a078440_list !! (n-1)

%o a078440_list = filter notbp a196871_list where

%o notbp x = m > 0 && x > 1 || m == 0 && notbp x' where

%o (x',m) = divMod x 2

%o -- _Reinhard Zumkeller_, Oct 08 2011

%Y A006370; subsequence of A196871 (with binary powers).

%K nonn

%O 1,1

%A _Joseph L. Pe_, Dec 31 2002

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 09:04 EDT 2024. Contains 371240 sequences. (Running on oeis4.)