login
A153681
Greatest number m such that the fractional part of (1024/1000)^A153677(n) <= 1/m.
2
41, 60, 76, 116, 233, 463, 718, 1350, 12472, 13733, 17428, 27955, 32276, 41155, 62437, 69643, 111085, 811799, 2656810, 11462221, 56414953
OFFSET
1,1
FORMULA
a(n) = floor(1/fract((1024/1000)^A153677(n))), where fract(x) = x-floor(x).
EXAMPLE
a(2)=60 since 1/61 < fract((1024/1000)^A153677(2)) = fract((1024/1000)^68) = 0.0164... <= 1/60.
MATHEMATICA
A153677 = {1, 68, 142, 341, 395, 490, 585, 1164, 1707, 26366, 41358,
46074, 120805, 147332, 184259, 205661, 385710, 522271, 3418770,
3675376, 9424094};
Table[fp = FractionalPart[(1024/1000)^A153677[[n]]]; m = Floor[1/fp];
While[fp <= 1/m, m++]; m - 1, {n, 1, Length[A153677]}] (* Robert Price, Mar 25 2019 *)
KEYWORD
nonn,more
AUTHOR
Hieronymus Fischer, Jan 06 2009
EXTENSIONS
a(18)-a(21) from Robert Price, Mar 25 2019
STATUS
approved