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
G. C. Greubel, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (0,0,6,0,0,-1).
FORMULA
For n > 0, a(3*n) = A000129(2*n+1).
a(3*n+2) = a(3*n) + A000129(2*n+2).
a(3*n+4) = a(3*n+2) + a(3*n+3).
a(3*n) = ceiling((3+2*sqrt(2))^n*(2+sqrt(2))/4).
a(3*n+2)/a(3*n+1) -> 1/sqrt(2).
a(3*n+1)/a(3*n) -> 3-sqrt(2).
a(3*n)/a(3*n-1) -> (8+5*sqrt(2))/7.
G.f.: x*(2*x^9 - 13*x^6 - x^5 - x^4 + 4*x^3 + 5*x^2 + 3*x + 1) / (x^6 - 6*x^3 + 1). - Colin Barker, Jun 16 2013
EXAMPLE
a(4) = 10 since frac(1x) + frac(3x) + frac(5x) + frac(10x) < 1, while frac(1x) + frac(3x) + frac(5x) + frac(k*x) > 1 for all k > 5 and k < 10.
MATHEMATICA
CoefficientList[Series[(1 + 3*z + 5*z^2 + 4*z^3 - z^4 - z^5 - 13*z^6 + 2*z^9)/(1 - 6*z^3 + z^6), {z, 0, 40}], z] (* Vladimir Joseph Stephan Orlovsky, Feb 16 2012 *)
PROG
(PARI) x='x+O('x^50); Vec(x*(2*x^9 -13*x^6 -x^5 -x^4 +4*x^3 +5*x^2 +3*x +1)/(x^6-6*x^3 +1)) \\ G. C. Greubel, Sep 22 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Benoit Cloitre and Paul D. Hanna, Jan 20 2003
STATUS
approved