login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A098294 a(n) = ceiling(n*log_2(3/2)). 11
0, 1, 2, 2, 3, 3, 4, 5, 5, 6, 6, 7, 8, 8, 9, 9, 10, 10, 11, 12, 12, 13, 13, 14, 15, 15, 16, 16, 17, 17, 18, 19, 19, 20, 20, 21, 22, 22, 23, 23, 24, 24, 25, 26, 26, 27, 27, 28, 29, 29, 30, 30, 31, 32, 32, 33, 33, 34, 34, 35, 36, 36, 37, 37, 38, 39, 39, 40, 40, 41, 41, 42, 43, 43, 44 (list; graph; refs; listen; history; text; internal format)
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:
A227048(n,a(n)) = A001047(n). - Reinhard Zumkeller, Jun 30 2013
Differs from A005378 at indices n = 0,17,20,22,25,27,29,30,... - M. F. Hasler, Jun 29 2014
LINKS
Sonic Studio, Pythagorean Scale.
Eric Weisstein's World of Music, Pythagorean Scale
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
Sequence in context: A214656 A194979 A366502 * A195119 A077467 A005378
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, Oct 18 2004
STATUS
approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 7 06:58 EDT 2024. Contains 375008 sequences. (Running on oeis4.)