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 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Hieronymus Fischer, Dec 31 2008
EXTENSIONS
a(22)-a(25) from Robert Price, May 10 2012
STATUS
approved