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”).

A153667
Greatest number m such that the fractional part of (3/2)^A153663(n) >= 1-(1/m).
8
2, 2, 2, 2, 3, 14, 31, 33, 69, 137, 222, 318, 901, 1772, 2747, 12347, 16540, 18198, 135794, 222246, 570361, 2134829, 6901329, 75503109, 814558605
OFFSET
1,1
FORMULA
a(n) = floor(1/(1-fract((3/2)^A153663(n)))), where fract(x) = x-floor(x).
EXAMPLE
a(5)=3, since 1-(1/4)=0.75>fract((3/2)^A153663(5))=fract((3/2)^12)=0.746...>=1-(1/3).
MATHEMATICA
A153663 = {1, 5, 8, 10, 12, 14, 46, 58, 105, 157, 163, 455, 1060, 1256, 2677, 8093, 28277, 33327, 49304, 158643, 164000, 835999, 2242294, 25380333, 92600006};
Table[fp = FractionalPart[(3/2)^A153663[[n]]]; m = Floor[1/(1-fp)];
While[fp >= 1 - (1/m), m++]; m - 1, {n, 1, Length[A153663]}] (* Robert Price, Mar 26 2019 *)
KEYWORD
nonn,more
AUTHOR
Hieronymus Fischer, Dec 31 2008
EXTENSIONS
a(22)-a(25) from Robert Price, May 10 2012
STATUS
approved