A153681 - OEIS

A153681

Greatest number m such that the fractional part of (1024/1000)^A153677(n) <= 1/m.

41, 60, 76, 116, 233, 463, 718, 1350, 12472, 13733, 17428, 27955, 32276, 41155, 62437, 69643, 111085, 811799, 2656810, 11462221, 56414953

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,1

LINKS

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 *)

CROSSREFS

KEYWORD

nonn,more

AUTHOR

EXTENSIONS

a(18)-a(21) from Robert Price, Mar 25 2019

STATUS

approved