[go: up one dir, main page]

login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A330591
Number of Collatz steps to reach 1 starting from 6^n + 1.
0
16, 21, 26, 101, 83, 83, 145, 145, 220, 158, 145, 207, 114, 114, 450, 114, 357, 357, 282, 419, 419, 494, 494, 494, 494, 494, 494, 494, 543, 494, 543, 799, 799, 543, 543, 799, 543, 543, 799, 799, 791, 791, 791, 791, 861, 861, 861, 861, 998, 998, 998, 861, 861, 861
OFFSET
1,1
COMMENTS
The Collatz transform maps any positive integer k to k/2 if k is even or 3*k+1 if k is odd. There is a famous unsolved problem which says that, starting with any positive integer k, repeated application of the Collatz transform will eventually reach 1 or equivalently enter the cycle (4,2,1).
This sequence is related to A179118 and A212653, which look at the stopping times of numbers of the form 2^n+1 and 3^n+1 respectively. We note that there exist several sequences of arithmetic progressions with common difference 1 in the former and with common difference -1 in the latter. This sequence looks at stopping times of numbers of the form 2^n*3^n+1 where we see that there exist arithmetic progressions with common difference 1+(-1)=0. This is an interesting result that requires further investigation.
FORMULA
a(n) = A006577(6^n+1).
EXAMPLE
a(2)=21 because the Collatz trajectory of 6^2 + 1 = 37 is
37 -> 112 -> 56 -> 28 -> 14 -> 7 -> 22 ->
11 -> 34 -> 17 -> 52 -> 26 -> 13 -> 40 ->
20 -> 10 -> 5 -> 16 -> 8 -> 4 -> 2 -> 1, which is 21 steps.
MATHEMATICA
n=200;
For [j=1, j<n, j++,
i=6^j+1;
count=0;
While[i!=1,
count=count+1; i=If[Mod[i, 2]==0, i/2, 3*i+1]]; Print[count]]
PROG
(Python) from decimal import *
n=1000
for j in range(2, n):
i=6**j+1
count=0
while(i!=1):
if(i%2==0):
i=i//2
else:
i=3*i+1
count=count+1
print(count)
(PARI) nbsteps(n) = if(n<0, 0, my(s=n, c=0); while(s>1, s=if(s%2, 3*s+1, s/2); c++); c); \\ A006577
a(n) = nbsteps(6^n+1); \\ Michel Marcus, Dec 21 2019
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Aranya Kumar Bal, Dec 18 2019
STATUS
approved