OFFSET
1,1
COMMENTS
The height of a Collatz run starting at a number m is the number of steps to reach 1, A006577(m).
There are three additional blocks with starting values less than 5000000: a(35, 40, 49) = (2054827, 596310, 4330040); a(34) = 8363804. Altogether there are 55 blocks with starting values at most 50000000, the highest of which is a(47) = 37669696 while a(45) > 50000000.
After searching up to about k = 5.3595*10^12, the largest-indexed term observed in the sequence thus far is a(1770) = 2490262807816 which begins a string of 1770 numbers whose Collatz sequence height is 331. - Kevin P. Thompson, Aug 27 2022
After searching up to about k = 3.293*10^13, the largest-indexed term observed in the sequence thus far is a(2225) = 23969528245354 which begins a string of 2225 numbers whose Collatz sequence height is 373. - Kevin P. Thompson, May 28 2023
LINKS
Kevin P. Thompson, Table of n, a(n) for n = 1..784
Kevin P. Thompson, Table of n, a(n) for n = 1..1001 with missing values
EXAMPLE
a(2) = 14 since the 2 adjacent numbers 14 and 15 are the first consecutive 2 whose height in their Collatz runs is the same prime number, in this case 17, while the heights for the Collatz runs at 13 and 16 are the nonprimes 9 and 4, respectively.
a(9) = 24408 since the 9 adjacent numbers 24408 .. 24416 are the first consecutive 9 whose height in their Collatz runs is the same prime number, in this case 157, while the heights for the Collatz runs at 24407 and 24417 are the nonprimes 64 and 69, respectively.
MATHEMATICA
collatz[n_] := If[EvenQ[n], n/2, 3n+1]
height[n_] := Length[NestWhileList[collatz, n, #!=1&]] - 1
(* b is an estimate on the size of the list being computed *)
a339773[n_, b_] := Module[{k=2, c, d, j, pList=Table[0, {b}]}, While[k<=n, c=height[k-1]; d=height[k]; j=k+1; If[!PrimeQ[c]&&PrimeQ[d], While[height[j]==d, j++]; If[!PrimeQ[height[j]]&&pList[[j-k]]==0, pList[[j-k]]=k]]; k=j]; pList]
Take[a339773[5000000, 50], 33] (* sequence data *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Hartmut F. W. Hoft, Dec 16 2020
EXTENSIONS
a(34)-a(35) from Kevin P. Thompson, Aug 27 2022
STATUS
approved