|
%I #11 Mar 26 2019 11:46:29
%S 41,60,76,116,233,463,718,1350,12472,13733,17428,27955,32276,41155,
%T 62437,69643,111085,811799,2656810,11462221,56414953
%N Greatest number m such that the fractional part of (1024/1000)^A153677(n) <= 1/m.
%F a(n) = floor(1/fract((1024/1000)^A153677(n))), where fract(x) = x-floor(x).
%e a(2)=60 since 1/61 < fract((1024/1000)^A153677(2)) = fract((1024/1000)^68) = 0.0164... <= 1/60.
%t A153677 = {1, 68, 142, 341, 395, 490, 585, 1164, 1707, 26366, 41358,
%t 46074, 120805, 147332, 184259, 205661, 385710, 522271, 3418770,
%t 3675376, 9424094};
%t Table[fp = FractionalPart[(1024/1000)^A153677[[n]]]; m = Floor[1/fp];
%t While[fp <= 1/m, m++]; m - 1, {n, 1, Length[A153677]}] (* _Robert Price_, Mar 25 2019 *)
%Y Cf. A081464, A153669, A153677, A154130, A153689, A153697, A153705, A153713, A153721.
%K nonn,more
%O 1,1
%A _Hieronymus Fischer_, Jan 06 2009
%E a(18)-a(21) from _Robert Price_, Mar 25 2019
|
