OFFSET
1,2
COMMENTS
The n-th greedy frac multiple of x is the smallest integer that does not cause sum(k=1..n,frac(a(k)*x)) to exceed unity; an infinite number of terms appear as the denominators of the convergents to the continued fraction of x.
EXAMPLE
a(4) = 14 since frac(1x) + frac(3x) + frac(7x) + frac(14x) < 1, while frac(1x) + frac(3x) + frac(7x) + frac(k*x) > 1 for all k>7 and k<14.
MAPLE
Digits := 100: a := []: s := 0: x := exp(1.0): for n from 1 to 1000000 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
Benoit Cloitre and Paul D. Hanna, Jan 21 2003
EXTENSIONS
Two more terms from Mark Hudson (mrmarkhudson(AT)hotmail.com), Jan 30 2003
STATUS
approved