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

Offset: 1

Author

Hieronymus Fischer, Jan 06 2009

Keywords

• nonn
• more

Examples

			a(2)=60 since 1/61 < fract((1024/1000)^A153677(2)) = fract((1024/1000)^68) = 0.0164... <= 1/60.
		

Crossrefs

Cf. A081464, A153669, A153677, A154130, A153689, A153697, A153705, A153713, A153721.

Programs

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

Formula

a(n) = floor(1/fract((1024/1000)^A153677(n))), where fract(x) = x-floor(x).

Extensions

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