A143514 Array D of denominators of Best Remaining Lower Approximates of x=(1+sqrt(5))/2, by antidiagonals.
1, 2, 3, 5, 4, 6, 13, 7, 9, 8, 34, 10, 12, 11, 16, 89, 18, 15, 14, 19, 21, 233, 26, 23, 17, 22, 24, 29, 610, 47, 31, 20, 25, 27, 32, 37, 1597, 68, 39, 28, 33, 30, 35, 40, 42, 4181, 123, 60, 36, 41, 38, 43, 48, 45, 50, 10946, 178, 81, 44, 49, 46, 51, 56, 53, 58
Offset: 1
Examples
Northwest corner of D: 1 2 5 13 3 4 7 10 6 9 12 15 8 11 14 17 Northwest corner of R: 1/1 3/2 8/5 21/13 4/3 6/4 11/7 16/10 9/6 14/9 19/12 24/15 12/8 17/11 22/14 27/17
Links
- Clark Kimberling, Best lower and upper approximates to irrational numbers, Elemente der Mathematik, 52 (1997) 122-126.
Programs
-
Mathematica
r = N[(1 + Sqrt[5])/2, 100]; Table[d = 1; t[k] = {}; Do[a = FractionalPart[n*r]; If[a < d && ! MemberQ[Apply[Union, Map[t[#] &, Range[k - 1]]], n], d = a; AppendTo[t[k], n]], {n, 10000}]; t[k], {k, 12}]; Column[Table[t[k], {k, 1, 12}]] (* Peter J. C. Moses, Feb 18 2021 *)
Formula
For any positive irrational number x, define an array D by successive rows as follows: D(n,k) = least positive integer q not already in D such that there exists an integer p such that 0 < x - p/q < x - c/d for every positive rational number c/d that has 0 < d < q. Thus p/q is the "best remaining lower approximate" of x when all better lower approximates are unavailable. For each q, define N(n,k)=p and R(n,k)=p/q. Then R is the "array of best remaining lower approximates of x," D is the corresponding array of denominators and N, of numerators.
Comments