login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A078373 n sets a record for the number of primes in {n, f(n), f(f(n)), ..., 1}, where f is the Collatz function defined by f(x) = x/2 if x is even; f(x) = 3x + 1 if x is odd. 2
2, 3, 7, 19, 27, 97, 171, 231, 487, 763, 1071, 4011, 6171, 10971, 17647, 47059, 99151, 117511, 202471, 260847, 481959, 963919, 1564063, 1805311, 1993215, 6991599, 8400511, 11200681, 36791535, 46564287, 103359483, 206718967, 359502063 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

LINKS

Jud McCranie, Table of n, a(n) for n = 1..40

J. C. Lagarias, The 3x+1 problem and its generalizations, Amer. Math. Monthly, 92 (1985), 3-23.

Carlos Rivera, Puzzle 634. Primes in Collatz trajectory

Carlos Rivera, Puzzle 922. Follow up to Puzzle 634

EXAMPLE

The sequence n, f(n), f(f(n)), ..., 1 for n = 7 is: 7, 22, 11, 34, 17, 52, 26, 13, 40, 20, 10, 5, 16, 8, 4, 2, 1, which has six prime terms, more prime terms than for any n < 7. Hence 7 sets a record and so belongs to a(n).

MATHEMATICA

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]; high = 0; a = {}; For[j = 1, j <= 10^5, j++, k = g[j]; If[k > high, high = k; a = Append[a, j]]]; a

(* Second program: *)

With[{s = Array[Count[NestWhileList[If[EvenQ[#], #/2, 3 # + 1] &, #, # > 1 &], _?PrimeQ] &, 10^5]}, Map[FirstPosition[s, #][[1]] &, Union@ FoldList[Max, s]]] (* Michael De Vlieger, Apr 21 2018 *)

PROG

(PARI) A078350(n)=my(s=isprime(n)); while(n>1, if(n%2, n=(3*n+1)/2, n/=2); s+=isprime(n)); s

r=0; for(n=2, 1e9, t=A078350(n); if(t>r, r=t; print1(n", "))) \\ Charles R Greathouse IV, Apr 28 2015

CROSSREFS

Cf. A055509, A181921.

Sequence in context: A178954 A138111 A218100 * A038878 A040112 A246373

Adjacent sequences:  A078370 A078371 A078372 * A078374 A078375 A078376

KEYWORD

nonn

AUTHOR

Joseph L. Pe, Dec 24 2002

EXTENSIONS

a(18)-a(30) from Donovan Johnson, Jul 02 2010

a(31)-a(33) from Carlos Rivera, Apr 15 2012

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy .

Last modified May 28 03:59 EDT 2018. Contains 304733 sequences. (Running on oeis4.)