%I #15 May 08 2018 02:42:07
%S 1,32,2048,8192,2199023255552,144115188075855872,576460752303423488,
%T 2305843009213693952
%N Denominators of records low values of the ratio n*3^A006667(n)/2^A006666(n).
%e For n=1 to 10 the ratios are 1, 1, 27/32, 1, 15/16, 27/32, 1701/2048, 1, 6561/8192, 15/16, so the low records are 1, 27/32, 1701/2048, 6561/8192, ...
%t q=1;Collatz[n_]:=NestWhileList[If[EvenQ[#],#/2,3 #+1]&,n,#>1&];nn=5000;t={};n=0;While[Length[t]<nn,n++;c=Collatz[n];ev=Length[Select[c,EvenQ]];od=Length[c]-ev-1;If[n*3^od/2^ev<q,Print[n," ",od," ",ev," ",n*3^od/2^ev];q=n*3^od/2^ev]]
%o (PARI) ht(n) = my(h, t); while(n>1, if(n%2, n=3*n+1; t++, n>>=1; h++)); return([h, t]);
%o lista(nn) = {m = 2; for (n=1, nn, v = ht(n); newm = n*3^v[2]/2^v[1]; if (newm < m, print1(denominator(newm), ", "); m = newm)); } \\ _Michel Marcus_, May 06 2018
%Y Cf. A127789 (for the indices where these records occur).
%Y Cf. A304119 (for the corresponding numerators).
%Y Cf. A006666, A006667.
%K nonn,frac,more
%O 1,2
%A _Michel Lagneau_, May 06 2018
|