OFFSET
1,2
COMMENTS
Numbers k such that fract((4/3)^k) > 1-(1/k), where fract(x) = x - floor(x).
The next term is greater than 3*10^8.
EXAMPLE
a(4) = 39, since fract((4/3)^39) = 0.9991861450... > 0.974358... = 1 - (1/39), but fract((4/3)^k) <= 1 - (1/k) for 8 < k < 39.
MATHEMATICA
Select[Range[5000], N[FractionalPart[(4/3)^#], 100] > 1 - (1/#) &] (* G. C. Greubel, Sep 02 2016 *)
PROG
(PARI) isok(n) = frac((4/3)^n) > 1 - 1/n; \\ Michel Marcus, Sep 03 2016
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Hieronymus Fischer, Jan 11 2009
EXTENSIONS
a(11)-a(20) from Robert Gerbicz, Nov 21 2010
STATUS
approved