|
|
A079938
|
|
Greedy frac multiples of Pi: a(1)=1, Sum_{n>=1} frac(a(n)*x) = 1 at x = Pi.
|
|
4
|
|
|
|
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.
|
|
LINKS
|
|
|
EXAMPLE
|
a(4) = 8 since frac(1x*) + frac(2*x) + frac(3*x) + frac(8*x) < 1, while frac(1*x) + frac(2*x) + frac(3*x) + frac(k*x) > 1 for all k > 3 and k < 8.
|
|
MAPLE
|
Digits := 100: a := []: s := 0: x := Pi: for n from 1 to 10000000 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
|
|
|
EXTENSIONS
|
One more term from Mark Hudson, Jan 30 2003
|
|
STATUS
|
approved
|
|
|
|