OFFSET
0,3
COMMENTS
Original name was: Smallest exponent of 2 which gives a power of 2 which is equal to or bigger than (3/2)^n, n = 0,1,... .
Stacking perfect fifths (the frequency ratio of a fifth is 3/2) this sequence determines into which octave the n-th fifth falls. For example, the third fifth, (3/2)^3, falls into the second octave, which means that it lies in the interval [2^1,2^2)=[2,4). The k-th octave comprises ratios in the interval [2^(k-1),2^k), k=1,2,...
Related to the initial number of sequential even terms in an "ideal" sequence under iteration of the 3x+1 Problem on a positive odd value m, where the piecewise function f is given by f(2*m)=m, f(2*m+1)=6*m+4, to ensure f^A122437(n) (m) < m, where n > 1 is the number of odds in the sequence (including m) and floor(1+n*(log(3)/log(2))) is the number of evens. An "ideal" sequence minimizes the effects of f(2*m+1) by following a certain order of even or odd terms along with the rules of the function. A representation of such sequences in terms of parity sequences for values n >= 2 follows:
n=2, (o,e,e,o,e,e)
n=3, (o,e,e,o,e,o,e,e)
n=4, (o,e,e,e,o,e,o,e,o,e,e)
n=5, (o,e,e,e,o,e,o,e,o,e,o,e,e)
n=6, (o,e,e,e,e,o,e,o,e,o,e,o,e,o,e,e)
n=7, (o,e,e,e,e,e,o,e,o,e,o,e,o,e,o,e,o,e,e)
The pattern is clear, and the formula for the initial number of sequential even terms in each sequence is given by a(n) = floor(1+n*(log(3)/log(2)))-n for n > 1, where the sum of the number of even and odd terms is given by A122437(n) for n > 1. Of course, most values m do not have sequences following this pattern of iteration under f. Also, the reason for placing an extra even term at the end of such sequences is to mitigate to some degree the effects of the possibility that the last odd term is only "slightly" larger than m, i.e., (3*m+1)/4 < m for all m > 1. - Jeffrey R. Goodwin, Aug 25 2011
a(n) gives the position in n-th row of A227048 where (3^n - 2^n) occurs:
Differs from A005378 at indices n = 0,17,20,22,25,27,29,30,... - M. F. Hasler, Jun 29 2014
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
Sonic Studio, Pythagorean Scale.
Eric Weisstein's World of Music, Pythagorean Scale
Wikipedia, Pythagorean tuning
Mike Winkler, On the structure and the behaviour of Collatz 3n + 1 sequences, 2014.
FORMULA
2^a(n) >= (3/2)^n but 2^(a(n) - 1) < (3/2)^n, n >= 0.
a(n) = ceiling(tau*n) with tau := log(3)/log(2) - 1 = 0.584962501..., n >= 0.
a(n) = floor(1 + n * log(3)/log(2)) - n, n >= 1. - Mike Winkler, Dec 31 2010
EXAMPLE
a(0) = 0 because 2^0 = 1 = (3/2)^0 but 2^(-1) = 1/2 < 1.
a(11) = 7 because 2^7 = 128 > 86.497... = (3/2)^11 but 2^6 = 64 < (3/2)^11.
MAPLE
seq(ceil(n*log[2](3/2)), n=0..100); # Robert Israel, Jul 12 2015
MATHEMATICA
With[{c=Log2[3/2]}, Ceiling[c*Range[0, 80]]] (* Harvey P. Dale, Feb 24 2024 *)
PROG
(Haskell)
import Data.List (elemIndex); import Data.Maybe (fromJust)
a098294 0 = 0
a098294 n = fromJust (a001047 n `elemIndex` a227048_row n) + 1
-- Reinhard Zumkeller, Jun 30 2013
(Magma) [0] cat [Floor(1 + n * Log(3)/Log(2)) - n: n in [1..70]]; // Vincenzo Librandi, Jul 13 2015
(PARI) a(n)=ceil(n*log(3/2)/log(2)) \\ Charles R Greathouse IV, Jul 13 2015
(PARI) a(n) = !!n + logint(3^n, 2) - n \\ Ruud H.G. van Tol, Nov 21 2023
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, Oct 18 2004
STATUS
approved