View source for Stern–Brocot tree

{{Short description|Ordered binary tree of rational numbers}}
[[Image:SternBrocotTree.svg|thumb|400px|The Stern–Brocot tree, and the Stern–Brocot sequences of order {{mvar|i}} for {{math|1=''i'' = 1, 2, 3, 4}}]]

In [[number theory]], the '''Stern–Brocot tree''' is an [[Binary tree#Types of binary trees|infinite complete binary]] [[tree (graph theory)|tree]] in which the [[vertex (graph theory)|vertices]] correspond [[bijection|one-for-one]] to the [[positive number|positive]] [[rational number]]s, whose values are ordered from the left to the right as in a [[binary search tree]].
 
The Stern–Brocot tree was introduced independently by {{harvs|first=Moritz|last=Stern|authorlink=Moritz Stern|year=1858|txt}} and {{harvs|first=Achille|last=Brocot|authorlink=Achille Brocot|year=1861|txt}}. Stern was a German number theorist; Brocot was a French [[clockmaker]] who used the Stern–Brocot tree to design systems of gears with a [[gear ratio]] close to some desired value by finding a ratio of [[smooth number]]s near that value.

The root of the Stern–Brocot tree corresponds to the number 1. The parent-child relation between numbers in the Stern–Brocot tree may be defined in terms of [[simple continued fraction]]s or [[mediant (mathematics)|mediants]], and a path in the tree from the root to any other number {{mvar|q}} provides a sequence of [[Diophantine approximation|approximation]]s to {{mvar|q}} with smaller [[denominator]]s than {{mvar|q}}. Because the tree contains each positive rational number exactly once, a [[breadth first search]] of the tree provides a method of listing all positive rationals that is closely related to [[Farey sequence]]s. The left subtree of the Stern–Brocot tree, containing the rational numbers in the range {{math|(0,1)}}, is called the '''Farey tree'''.

==Generating rule==

Each vertex in the tree can be associated with a triple of fractions consisting of three fractions in the same row as the vertex, namely the fraction immediately to the left of the vertex, the fraction at the vertex itself, and the fraction immediately to the right of the vertex. (Refer to the figure above.) The left and right fractions do not correspond to vertices in the same row as the vertex, but rather to vertices in some preceding row. Each such fraction can be understood as labeling the region of the plane bounded by two infinite paths descending from the preceding vertex labeled by the same fraction. The second element of a triple will always be the mediant of the first and third elements. For example, the root is associated with <math>\bigl(\tfrac{0}{1}, \tfrac{1}{1}, \tfrac{1}{0}\bigr)</math> and its left and right descendents are associated with <math>\bigl(\tfrac{0}{1}, \tfrac{1}{2}, \tfrac{1}{1}\bigr)</math> and <math>\bigl(\tfrac{1}{1}, \tfrac{2}{1}, \tfrac{1}{0}\bigr).</math> 

The tree is generated by the following rule: 
<math display=block>\begin{array}{ccccc}
  & & \left(\dfrac{a}{b}, \dfrac{c}{d}, \dfrac{e}{f}\right) & & \\
  & \swarrow & & \searrow & \\
  \left(\dfrac{a}{b}, \dfrac{a+c}{b+d}, \dfrac{c}{d}\right) & & & & \left(\dfrac{c}{d}, \dfrac{c+e}{d+f}, \dfrac{e}{f}\right) 
\end{array}</math>

==A tree of continued fractions==

The relationship between parents and descendants can also be understood in terms of continued fractions. Every positive rational number {{mvar|q}} may be expressed as a continued fraction of the form 
<math display=block>q = a_0 + \cfrac{1}{a_1 + \cfrac{1}{a_2 + \cfrac{1}{a_3 + \cfrac{1}{\ddots+\cfrac1{a_k}}}}}= [a_0; a_1, a_2, \ldots, a_k] </math>
where {{mvar|k}} and {{math|''a''{{sub|0}}}} are non-negative integers, and each subsequent coefficient {{math|a{{sub|i}}}} is a positive integer. This representation is not unique because
<math display=block> [a_0; a_1, a_2, \ldots, a_{k-1}, 1] = [a_0; a_1, a_2, \ldots, a_{k-1}+1],</math>
but using this equivalence to replace every continued fraction ending with a one by a shorter continued fraction shows that every rational number has a unique representation in which the last coefficient is greater than one. Then, unless {{math|1=''q'' = 1}}, the number {{mvar|q}} has a parent in the Stern–Brocot tree given by the continued fraction expression 
<math display=block>[a_0; a_1, a_2, \ldots, a_k-1].</math> 
Equivalently this parent is formed by decreasing the denominator in the innermost term of the continued fraction by&nbsp;1, and contracting with the previous term if the fraction becomes {{sfrac|1|1}}. For instance, the rational number {{sfrac|23|16}} has the continued fraction representation 
<math display=block>\frac{23}{16}=1+\cfrac1{2+\cfrac1{3+\frac12}} = [1;2,3,2],</math> 
so its parent in the Stern–Brocot tree is the number 
<math display=block>[1;2,3,1]=[1;2,4]=1+\cfrac1{2+\frac14}=\frac{13}{9}.</math>

Conversely each number {{mvar|q}} in the Stern–Brocot tree has exactly two children: if 
<math display=block> q = [a_0; a_1, a_2, \ldots, a_k] = [a_0; a_1, a_2, \ldots, a_k-1,1]</math>
then one child is the number represented by the continued fraction 
<math display=block>\displaystyle [a_0; a_1, a_2, \ldots, a_k+1]</math> 
while the other child is represented by the continued fraction 
<math display=block> [a_0; a_1, a_2, \ldots, a_k-1,2].</math> 
One of these children is less than {{mvar|q}} and this is the left child; the other is greater than {{mvar|q}} and it is the right child (in fact the former expression gives the left child if {{mvar|k}} is odd, and the right child if {{mvar|k}} is even).
For instance, the continued fraction representation of {{sfrac|13|9}} is {{nowrap|[1;2,4]}} and its two children are {{nowrap|1=[1;2,5] = {{sfrac|16|11}}}} (the right child) and {{nowrap|1=[1;2,3,2] = {{sfrac|23|16}}}} (the left child).

It is clear that for each finite continued fraction expression one can repeatedly move to its parent, and reach the root {{math|1=[1;] = {{sfrac|1|1}}}} of the tree in finitely many steps (in {{math|''a''<sub>0</sub> + ... + ''a''<sub>''k''</sub> − 1}} steps to be precise). Therefore, every positive rational number appears exactly once in this tree. Moreover all descendants of the left child of any number {{mvar|q}} are less than {{mvar|q}}, and all descendants of the right child of {{mvar|q}} are greater than {{mvar|q}}. The numbers at depth {{mvar|d}} in the tree are the numbers for which the sum of the continued fraction coefficients is {{math|''d'' + 1}}.

==Mediants and binary search==

The Stern–Brocot tree forms an infinite [[binary search tree]] with respect to the usual ordering of the rational numbers.<ref>{{citation|last1=Graham|first1=Ronald L.|authorlink=Ronald Graham|last2= Knuth|first2=Donald E.|authorlink2=Donald Knuth|last3=Patashnik|first3=Oren|authorlink3=Oren Patashnik|title=Concrete mathematics|edition=Second|publisher=Addison-Wesley|year=1994|isbn=0-201-55802-5|pages=116–118}}</ref><ref>{{citation
 | last1 = Gibbons | first1 = Jeremy
 | last2 = Lester | first2 = David
 | last3 = Bird | first3 = Richard | author3-link = Richard Bird (computer scientist)
 | doi = 10.1017/S0956796806005880
 | issue = 3
 | journal = Journal of Functional Programming
 | pages = 281–291
 | title = Functional pearl: Enumerating the rationals
 | volume = 16
 | year = 2006| s2cid = 14237968
 | doi-access = free
 }}.</ref> The set of rational numbers descending from a node {{mvar|q}} is defined by the open interval {{math|(''L<sub>q</sub>'', ''H<sub>q</sub>'')}} where {{mvar|L<sub>q</sub>}} is the ancestor of {{mvar|q}} that is smaller than {{mvar|q}} and closest to it in the tree (or {{math|1=''L<sub>q</sub>'' = 0}} if {{mvar|q}} has no smaller ancestor) while {{mvar|H<sub>q</sub>}} is the ancestor of {{mvar|q}} that is larger than {{mvar|q}} and closest to it in the tree (or {{math|1=''H<sub>q</sub>'' = +∞}} if {{mvar|q}} has no larger ancestor).

The path from the root 1 to a number {{mvar|q}} in the Stern–Brocot tree may be found by a [[binary search]] algorithm, which may be expressed in a simple way using [[Mediant (mathematics)|mediants]]. Augment the non-negative rational numbers to including a value {{sfrac|1|0}} (representing +∞) that is by definition greater than all other rationals. The [[binary search algorithm]] proceeds as follows:
*Initialize two values {{mvar|L}} and {{mvar|H}} to {{sfrac|0|1}} and {{sfrac|1|0}}, respectively.
*Until {{mvar|q}} is found, repeat the following steps:
**Let {{math|1=''L'' = {{sfrac|''a''|''b''}}}} and {{math|1=''H'' = {{sfrac|''c''|''d''}}}}; compute the mediant <math>M = L \oplus H = \frac{a+c}{b+d}.</math>
**If {{mvar|M}} is less than {{mvar|q}}, then {{mvar|q}} is in the open interval {{math|(''M'', ''H'')}}; replace {{mvar|L}} by {{mvar|M}} and continue.
**If {{mvar|M}} is greater than {{mvar|q}}, then {{mvar|q}} is in the open interval {{math|(''L'', ''M'')}}; replace {{mvar|H}} by {{mvar|M}} and continue.
**In the remaining case, {{math|1=''q'' = ''M''}}; terminate the search algorithm.
The sequence of values {{mvar|M}} computed by this search is exactly the sequence of values on the path from the root to {{mvar|q}} in the Stern–Brocot tree. Each open interval {{math|(''L'', ''H'')}} occurring at some step in the search is the interval {{math|(''L<sub>M</sub>'', ''H<sub>M</sub>''}}) representing the descendants of the mediant {{mvar|M}}. The parent of {{mvar|q}} in the Stern–Brocot tree is the last mediant found that is not equal to {{mvar|q}}.

This binary search procedure can be used to convert [[floating-point]] numbers into rational numbers. By stopping once the desired precision is reached, floating-point numbers can be approximated to arbitrary precision.<ref>Sedgewick and Wayne, ''Introduction to Programming in Java''. A Java implementation of this algorithm can be found [http://www.cs.princeton.edu/introcs/92symbolic/RationalApprox.java.html here].</ref> If a real number {{mvar|x}} is approximated by any rational number {{math|{{sfrac|''a''|''b''}}}} that is not in the sequence of mediants found by the algorithm above, then the sequence of mediants contains a closer approximation to {{mvar|x}} that has a denominator at most equal to {{mvar|b}}; in that sense, these mediants form the [[best rational approximation]]s to {{mvar|x}}.

The Stern–Brocot tree may itself be defined directly in terms of mediants: the left child of any number {{mvar|q}} is the mediant of {{mvar|q}} with its closest smaller ancestor, and the right child of {{mvar|q}} is the mediant of {{mvar|q}} with its closest larger ancestor. In this formula, {{mvar|q}} and its ancestor must both be taken in lowest terms, and if there is no smaller or larger ancestor then {{sfrac|0|1}} or {{sfrac|1|0}} should be used respectively. Again, using {{sfrac|7|5}} as an example, its closest smaller ancestor is {{sfrac|4|3}}, so its left child is {{nowrap|1={{sfrac|4 + 7|3 + 5}} = {{sfrac|11|8}}}}, and its closest larger ancestor is {{sfrac|3|2}}, so its right child is {{nowrap|1={{sfrac|7 + 3|5 + 2}} = {{sfrac|10|7}}}}.

==Relation to Farey sequences==
The [[Farey sequence]] of order {{mvar|n}} is the sorted sequence of fractions in the closed interval {{math|[0,1]}} that have denominator less than or equal to {{mvar|n}}. As in the binary search technique for generating the Stern–Brocot tree, the Farey sequences can be constructed using mediants: the Farey sequence of order {{math|''n'' + 1}} is formed from the Farey sequence of order {{mvar|n}} by computing the mediant of each two consecutive values in the Farey sequence of order {{mvar|n}}, keeping the subset of mediants that have denominator exactly equal to {{math|''n'' + 1}}, and placing these mediants between the two values from which they were computed.

A similar process of mediant insertion, starting with a different pair of interval endpoints <math>\bigl[ \tfrac{0}{1}, \tfrac{1}{0} \bigr],</math> may also be seen to describe the construction of the vertices at each level of the Stern–Brocot tree. The '''Stern–Brocot sequence''' of order 0 is the sequence <math>\bigl[ \tfrac{0}{1}, \tfrac{1}{0} \bigr],</math> and the Stern–Brocot sequence of order {{mvar|i}} is the sequence formed by inserting a mediant between each consecutive pair of values in the Stern–Brocot sequence of order {{math|''i'' &minus; 1}}. The Stern–Brocot sequence of order {{mvar|i}} consists of all values at the first {{mvar|i}} levels of the Stern–Brocot tree, together with the boundary values {{sfrac|0|1}} and {{sfrac|1|0}}, in numerical order.

Thus the Stern–Brocot sequences differ from the Farey sequences in two ways: they eventually include all positive rationals, not just the rationals within the interval {{math|[0,1]}}, and at the {{mvar|n}}th step all mediants are included, not only the ones with denominator equal to {{mvar|n}}. The Farey sequence of order {{mvar|n}} may be found by an inorder traversal of the left subtree of the Stern–Brocot tree, backtracking whenever a number with denominator greater than {{mvar|n}} is reached.

==Additional properties==

If <math>\tfrac{p_1}{q_1}, \tfrac{p_2}{q_2}, \dots, \tfrac{p_n}{q_n}</math> are all the rationals at the same depth in the Stern–Brocot tree, then

<math display=block> \sum_{k=1}^n \frac{1}{p_kq_k} = 1.</math>

Moreover, if <math>\tfrac{p}{q} < \tfrac{p'}{q'}</math> are two consecutive fractions at or above a certain level in the tree (in the sense that any fraction between them, must be in a lower level of the tree), then<ref>Bogomolny credits this property to Pierre Lamothe, a Canadian music theorist.</ref>

<math display=block> p'q - pq' = 1.</math>

Along with the definitions in terms of continued fractions and mediants described above, the Stern–Brocot tree may also be defined as a [[Cartesian tree]] for the rational numbers, prioritized by their denominators. In other words, it is the unique binary search tree of the rational numbers in which the parent of any vertex {{mvar|q}} has a smaller denominator than {{mvar|q}} (or if {{mvar|q}} and its parent are both integers, in which the parent is smaller than {{mvar|q}}). It follows from the theory of Cartesian trees that the [[lowest common ancestor]] of any two numbers {{mvar|q}} and {{mvar|r}} in the Stern–Brocot tree is the rational number in the closed interval {{math|[''q'', ''r'']}} that has the smallest denominator among all numbers in this interval.

Permuting the vertices on each level of the Stern–Brocot tree by a [[bit-reversal permutation]] produces a different tree, the [[Calkin–Wilf tree]], in which the children of each number {{math|{{sfrac|''a''|''b''}}}} are the two numbers {{math|{{sfrac|''a''|''a'' + ''b''}}}} and {{math|{{sfrac|''a'' + ''b''|''b''}}}}. Like the Stern–Brocot tree, the Calkin–Wilf tree contains each positive rational number exactly once, but it is not a binary search tree.

==See also==
* [[Minkowski's question-mark function]], whose definition for rational arguments is closely related to the Stern–Brocot tree
* [[Calkin–Wilf tree]]
* [https://github.com/jconorgrogan/Stern-Brocot-Tree-exploration/tree/main Code for generating Stern-Brocot trees]

==Notes==
{{reflist}}

==References==
*{{citation
 | last = Brocot | first = Achille
 | journal = Revue Chronométrique
 | pages = 186–194
 | title = Calcul des rouages par approximation, nouvelle méthode
 | volume = 3
 | year = 1861}}.
*Brocot, Achille (1862), "Calcul des rouages par approximation, nouvelle méthode", https://gallica.bnf.fr/ark:/12148/bpt6k1661912?rk=21459;2
*{{citation
 | last = Stern | first = Moritz A. | author-link = Moritz Abraham Stern
 | journal = [[Crelle's Journal|Journal für die reine und angewandte Mathematik]]
 | pages = 193–220
 | title = Ueber eine zahlentheoretische Funktion
 | url = http://www.digizeitschriften.de/resolveppn/GDZPPN002150301
 | volume = 55
 | year = 1858}}.
* {{citation | last1=Berstel | first1=Jean | last2=Lauve | first2=Aaron | last3=Reutenauer | first3=Christophe | last4=Saliola | first4=Franco V. | title=Combinatorics on words. Christoffel words and repetitions in words | series=CRM Monograph Series | volume=27 | location=Providence, RI | publisher=[[American Mathematical Society]] | year=2009 | isbn=978-0-8218-4480-9 | zbl=1161.68043 }}

==External links==
*{{citation|first=Dhroova|last=Aiylam|arxiv=1301.6807|title=Modified Stern-Brocot Sequences|year=2013 |bibcode=2013arXiv1301.6807A}}
*{{citation|first=David|last=Austin|url=http://www.ams.org/featurecolumn/archive/stern-brocot.html|title=Trees, Teeth, and Time: The mathematics of clock making|publisher=Feature Column from the AMS}}
*{{citation | first=Alexander |last=Bogomolny | authorlink = Alexander Bogomolny | title=Stern Brocot-Tree | url=http://www.cut-the-knot.org/blue/Stern.shtml | accessdate=2008-09-03 | publisher=[[cut-the-knot]]}}.
*{{citation | first = N. J. A. | last = Sloane | authorlink = Neil Sloane | url = http://oeis.org/stern_brocot.html | title = The Stern–Brocot or Farey Tree | publisher = [[On-line Encyclopedia of Integer Sequences]]}}.
*{{citation | first = Norman | last = Wildberger | url = https://www.youtube.com/watch?v=CiO8iAYC6xI | title = MF96: Fractions and the Stern-Brocot tree}}.
*{{MathWorld|urlname=Stern-BrocotTree|title=Stern–Brocot Tree|mode=cs2}}
*{{PlanetMath|urlname=SternBrocotTree|title=Stern–Brocot tree}}
*[https://www.youtube.com/watch?v=DpwUVExX27E Infinite fractions], [[Numberphile]]
*[https://www.youtube.com/watch?v=j0o-pMIR8uk Amazing Graphs III], [[Numberphile]]
* {{OEIS el|1=A002487|2=Stern's diatomic series (or Stern-Brocot sequence)}}

{{DEFAULTSORT:Stern-Brocot tree}}
[[Category:Continued fractions]]
[[Category:Trees (data structures)]]