Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A080158
Greedy frac multiples of Catalan's constant, G: a(1)=1, sum(n>0,frac(a(n)*x))=1 at x=G, where "frac(y)" denotes the fractional part of y.
0
1, 11, 107, 10579, 21158, 53014, 106028, 625708, 721157, 1442314, 2163471, 2884628, 3605785, 4326942
OFFSET
1,2
COMMENTS
For definition of how the "Greedy Frac" sequence is defined, see other sequences in index.
EXAMPLE
a(3) = 107 since frac(1x) + frac(11x) + frac(107x) < 1, while frac(1x) + frac(11x) + frac(k*x) > 1 for all k>11 and k<107.
MAPLE
Digits := 1000: a := []: s := 0: x := evalf(Catalan): for n from 1 to 5000000 do: temp := evalf(s+frac(n*x)): if (temp<1.0) then a := [op(a), n]: print(n): s := s+evalf(frac(n*x)): fi: od: a;
CROSSREFS
KEYWORD
more,nonn
AUTHOR
Mark Hudson (mrmarkhudson(AT)hotmail.com), Jan 31 2003
STATUS
approved