Farey sequence

Last updated November 12, 2024

Farey diagram to F₉.

Symmetrical pattern made by the denominators of the Farey sequence, F₉.

Symmetrical pattern made by the denominators of the Farey sequence, F₂₅.

In mathematics, the Farey sequence of order n is the sequence of completely reduced fractions, either between 0 and 1, or without this restriction,^[a] which when in lowest terms have denominators less than or equal to n, arranged in order of increasing size.

Examples
Farey sunburst
History
Properties
Sequence length and index of a fraction
Farey neighbours
Farey neighbours and continued fractions
Farey fractions and the least common multiple
Farey fractions and the greatest common divisor
Applications
Ford circles
Riemann hypothesis
Other sums involving Farey fractions
Next term
See also
Footnotes
References
Further reading
External links

With the restricted definition, each Farey sequence starts with the value 0, denoted by the fraction ⁠0/1⁠, and ends with the value 1, denoted by the fraction ⁠1/1⁠ (although some authors omit these terms).

A Farey sequence is sometimes called a Farey series, which is not strictly correct, because the terms are not summed.^[2]

Examples

The Farey sequences of orders 1 to 8 are :

F₁ = { ⁠01⁠_,⁠11⁠ }

F₂ = { ⁠01⁠_,⁠12⁠_,⁠11⁠ }

F₃ = { ⁠01⁠_,⁠13⁠_,⁠12⁠_,⁠23⁠_,⁠11⁠ }

F₄ = { ⁠01⁠_,⁠14⁠_,⁠13⁠_,⁠12⁠_,⁠23⁠_,⁠34⁠_,⁠11⁠ }

F₅ = { ⁠01⁠_,⁠15⁠_,⁠14⁠_,⁠13⁠_,⁠25⁠_,⁠12⁠_,⁠35⁠_,⁠23⁠_,⁠34⁠_,⁠45⁠_,⁠11⁠ }

F₆ = { ⁠01⁠_,⁠16⁠_,⁠15⁠_,⁠14⁠_,⁠13⁠_,⁠25⁠_,⁠12⁠_,⁠35⁠_,⁠23⁠_,⁠34⁠_,⁠45⁠_,⁠56⁠_,⁠11⁠ }

F₇ = { ⁠01⁠_,⁠17⁠_,⁠16⁠_,⁠15⁠_,⁠14⁠_,⁠27⁠_,⁠13⁠_,⁠25⁠_,⁠37⁠_,⁠12⁠_,⁠47⁠_,⁠35⁠_,⁠23⁠_,⁠57⁠_,⁠34⁠_,⁠45⁠_,⁠56⁠_,⁠67⁠_,⁠11⁠ }

F₈ = { ⁠01⁠_,⁠18⁠_,⁠17⁠_,⁠16⁠_,⁠15⁠_,⁠14⁠_,⁠27⁠_,⁠13⁠_,⁠38⁠_,⁠25⁠_,⁠37⁠_,⁠12⁠_,⁠47⁠_,⁠35⁠_,⁠58⁠_,⁠23⁠_,⁠57⁠_,⁠34⁠_,⁠45⁠_,⁠56⁠_,⁠67⁠_,⁠78⁠_,⁠11⁠ }

Centered
F₁ = { ⁠0/1⁠_,⁠1/1⁠ }
F₂ = { ⁠0/1⁠_,⁠1/2⁠_,⁠1/1⁠ }
F₃ = { ⁠0/1⁠_,⁠1/3⁠_,⁠1/2⁠_,⁠2/3⁠_,⁠1/1⁠ }
F₄ = { ⁠0/1⁠_,⁠1/4⁠_,⁠1/3⁠_,⁠1/2⁠_,⁠2/3⁠_,⁠3/4⁠_,⁠1/1⁠ }
F₅ = { ⁠0/1⁠_,⁠1/5⁠_,⁠1/4⁠_,⁠1/3⁠_,⁠2/5⁠_,⁠1/2⁠_,⁠3/5⁠_,⁠2/3⁠_,⁠3/4⁠_,⁠4/5⁠_,⁠1/1⁠ }
F₆ = { ⁠0/1⁠_,⁠1/6⁠_,⁠1/5⁠_,⁠1/4⁠_,⁠1/3⁠_,⁠2/5⁠_,⁠1/2⁠_,⁠3/5⁠_,⁠2/3⁠_,⁠3/4⁠_,⁠4/5⁠_,⁠5/6⁠_,⁠1/1⁠ }
F₇ = { ⁠0/1⁠_,⁠1/7⁠_,⁠1/6⁠_,⁠1/5⁠_,⁠1/4⁠_,⁠2/7⁠_,⁠1/3⁠_,⁠2/5⁠_,⁠3/7⁠_,⁠1/2⁠_,⁠4/7⁠_,⁠3/5⁠_,⁠2/3⁠_,⁠5/7⁠_,⁠3/4⁠_,⁠4/5⁠_,⁠5/6⁠_,⁠6/7⁠_,⁠1/1⁠ }
F₈ = { ⁠0/1⁠_,⁠1/8⁠_,⁠1/7⁠_,⁠1/6⁠_,⁠1/5⁠_,⁠1/4⁠_,⁠2/7⁠_,⁠1/3⁠_,⁠3/8⁠_,⁠2/5⁠_,⁠3/7⁠_,⁠1/2⁠_,⁠4/7⁠_,⁠3/5⁠_,⁠5/8⁠_,⁠2/3⁠_,⁠5/7⁠_,⁠3/4⁠_,⁠4/5⁠_,⁠5/6⁠_,⁠6/7⁠_,⁠7/8⁠_,⁠1/1⁠ }

Sorted

 F1 = {0/1,                                                                                                          1/1}  F2 = {0/1,                                                   1/2,                                                   1/1}  F3 = {0/1,                               1/3,                1/2,                2/3,                               1/1}  F4 = {0/1,                     1/4,      1/3,                1/2,                2/3,      3/4,                     1/1}  F5 = {0/1,                1/5, 1/4,      1/3,      2/5,      1/2,      3/5,      2/3,      3/4, 4/5,                1/1}  F6 = {0/1,           1/6, 1/5, 1/4,      1/3,      2/5,      1/2,      3/5,      2/3,      3/4, 4/5, 5/6,           1/1}  F7 = {0/1,      1/7, 1/6, 1/5, 1/4, 2/7, 1/3,      2/5, 3/7, 1/2, 4/7, 3/5,      2/3, 5/7, 3/4, 4/5, 5/6, 6/7,      1/1}  F8 = {0/1, 1/8, 1/7, 1/6, 1/5, 1/4, 2/7, 1/3, 3/8, 2/5, 3/7, 1/2, 4/7, 3/5, 5/8, 2/3, 5/7, 3/4, 4/5, 5/6, 6/7, 7/8, 1/1}

Farey sunburst

Plotting the numerators versus the denominators of a Farey sequence gives a shape like the one to the right, shown for $F 6 .$

Reflecting this shape around the diagonal and main axes generates the Farey sunburst, shown below. The Farey sunburst of order $n$ connects the visible integer grid points from the origin in the square of side $2 n$ , centered at the origin. Using Pick's theorem, the area of the sunburst is $4(| F n | - 1)$ , where $| F n |$ is the number of fractions in $F n$ .

History

The history of 'Farey series' is very curious — Hardy & Wright (1979)^[3]

... once again the man whose name was given to a mathematical relation was not the original discoverer so far as the records go. — Beiler (1964)^[4]

Farey sequences are named after the British geologist John Farey, Sr., whose letter about these sequences was published in the Philosophical Magazine in 1816. Farey conjectured, without offering proof, that each new term in a Farey sequence expansion is the mediant of its neighbours. Farey's letter was read by Cauchy, who provided a proof in his Exercices de mathématique, and attributed this result to Farey. In fact, another mathematician, Charles Haros, had published similar results in 1802 which were not known either to Farey or to Cauchy.^[4] Thus it was a historical accident that linked Farey's name with these sequences. This is an example of Stigler's law of eponymy.

Properties

Sequence length and index of a fraction

The Farey sequence of order $n$ contains all of the members of the Farey sequences of lower orders. In particular $F n$ contains all of the members of $F n - 1$ and also contains an additional fraction for each number that is less than $n$ and coprime to $n$ . Thus $F 6$ consists of $F 5$ together with the fractions ⁠1/6⁠ and ⁠5/6⁠.

The middle term of a Farey sequence $F n$ is always ⁠1/2⁠, for $n > 1$ . From this, we can relate the lengths of $F n$ and $F n - 1$ using Euler's totient function $φ (n)$ :

$|F_{n}|=|F_{n-1}|+\varphi (n).$

Using the fact that $| F 1 | = 2$ , we can derive an expression for the length of $F n$ :^[5]

$|F_{n}|=1+\sum _{m=1}^{n}\varphi (m)=1+\Phi (n),$ where $Φ (n)$ is the summatory totient.

We also have : $|F_{n}|={\frac {1}{2}}\left(3+\sum _{d=1}^{n}\mu (d)\left\lfloor {\tfrac {n}{d}}\right\rfloor ^{2}\right),$ and by a Möbius inversion formula : $|F_{n}|={\frac {1}{2}}(n+3)n-\sum _{d=2}^{n}|F_{\lfloor n/d\rfloor }|,$ where $μ (d)$ is the number-theoretic Möbius function, and $\lfloor n/d\rfloor$ is the floor function.

The asymptotic behaviour of $| F n |$ is : $|F_{n}|\sim {\frac {3n^{2}}{\pi ^{2}}}.$

The number of Farey fractions with denominators equal to $k$ in $F n$ is given by $φ (k)$ when $k \leq n$ and zero otherwise. Concerning the numerators one can define the function ${\mathcal {N}}_{n}(h)$ that returns the number of Farey fractions with numerators equal to $h$ in $F n$ . This function has some interesting properties as^[6]

{\mathcal {N}}_{n}(1)=n

,

{\mathcal {N}}_{n}(p^{m})=\left\lceil (n-p^{m})\left(1-1/p\right)\right\rceil

for any prime number

p

,

{\mathcal {N}}_{n+mh}(h)={\mathcal {N}}_{n}(h)+m\varphi (h)

for any integer

m \geq 0

,

{\mathcal {N}}_{n}(4h)={\mathcal {N}}_{n}(2h)-\varphi (2h).

In particular, the property in the third line above implies ${\mathcal {N}}_{mh}(h)=(m-1)\varphi (h)$ and, further, ${\mathcal {N}}_{2h}(h)=\varphi (h).$ The latter means that, for Farey sequences of even order $n$ , the number of fractions with numerators equal to $⁠ n / 2 ⁠$ is the same as the number of fractions with denominators equal to $⁠ n / 2 ⁠$ , that is ${\mathcal {N}}_{n}(n/2)=\varphi (n/2)$ .

The index $I_{n}(a_{k,n})=k$ of a fraction $a_{k,n}$ in the Farey sequence $F_{n}=\{a_{k,n}:k=0,1,\ldots ,m_{n}\}$ is simply the position that $a_{k,n}$ occupies in the sequence. This is of special relevance as it is used in an alternative formulation of the Riemann hypothesis, see below. Various useful properties follow: ${\begin{aligned}I_{n}(0/1)&=0,\\[6pt]I_{n}(1/n)&=1,\\[2pt]I_{n}(1/2)&={\frac {|F_{n}|-1}{2}},\\[2pt]I_{n}(1/1)&=|F_{n}|-1,\\[2pt]I_{n}(h/k)&=|F_{n}|-1-I_{n}\left({\frac {k-h}{k}}\right).\end{aligned}}$

The index of $⁠ 1 / k ⁠$ where $⁠ n / i +1 ⁠ < k \leq ⁠ n / i ⁠$ and $n$ is the least common multiple of the first $i$ numbers, $n = lcm([2, i])$ , is given by:^[7] $I_{n}(1/k)=1+n\sum _{j=1}^{i}{\frac {\varphi (j)}{j}}-k\Phi (i).$

Farey neighbours

Fractions which are neighbouring terms in any Farey sequence are known as a Farey pair and have the following properties.

If $⁠ a / b ⁠$ and $⁠ c / d ⁠$ are neighbours in a Farey sequence, with $⁠ a / b ⁠ < ⁠ c / d ⁠$ , then their difference $⁠ c / d ⁠ - ⁠ a / b ⁠$ is equal to $⁠ 1 / bd ⁠$ . Since ${\frac {c}{d}}-{\frac {a}{b}}={\frac {bc-ad}{bd}},$

this is equivalent to saying that $bc-ad=1.$

Thus ⁠1/3⁠ and ⁠2/5⁠ are neighbours in $F 5$ , and their difference is ⁠1/15⁠.

The converse is also true. If $bc-ad=1$

for positive integers $a, b, c, d$ with $a < b$ and $c < d$ , then $⁠ a / b ⁠$ and $⁠ c / d ⁠$ will be neighbours in the Farey sequence of order $max(b,d)$ .

If $⁠ p / q ⁠$ has neighbours $⁠ a / b ⁠$ and $⁠ c / d ⁠$ in some Farey sequence, with $⁠ a / b ⁠ < ⁠ p / q ⁠ < ⁠ c / d ⁠$ , then $⁠ p / q ⁠$ is the mediant of $⁠ a / b ⁠$ and $⁠ c / d ⁠$ – in other words, ${\frac {p}{q}}={\frac {a+c}{b+d}}.$

This follows easily from the previous property, since if ${\begin{aligned}&&bp-aq&=qc-pd=1,\\[4pt]\implies &&bp+pd&=qc+aq,\\[4pt]\implies &&p(b+d)&=q(a+c),\\\implies &&{\frac {p}{q}}&={\frac {a+c}{b+d}}.\end{aligned}}$

It follows that if $⁠ a / b ⁠$ and $⁠ c / d ⁠$ are neighbours in a Farey sequence then the first term that appears between them as the order of the Farey sequence is incremented is ${\frac {a+c}{b+d}},$

which first appears in the Farey sequence of order $b + d$ .

Thus the first term to appear between ⁠1/3⁠ and ⁠2/5⁠ is ⁠3/8⁠, which appears in $F 8$ .

The total number of Farey neighbour pairs in $F n$ is $2 | F n | - 3$ .

The Stern–Brocot tree is a data structure showing how the sequence is built up from 0 (= ⁠0/1⁠) and 1 (= ⁠1/1⁠), by taking successive mediants.

Equivalent-area interpretation

Every consecutive pair of Farey rationals have an equivalent area of 1.^[8] See this by interpreting consecutive rationals $r_{1}={\frac {p}{q}}\qquad r_{2}={\frac {p'}{q'}}$ as vectors $(p, q)$ in the xy-plane. The area is given by $A\left({\frac {p}{q}},{\frac {p'}{q'}}\right)=qp'-q'p.$ As any added fraction in between two previous consecutive Farey sequence fractions is calculated as the mediant (⊕), then ${\begin{aligned}A(r_{1},r_{1}\oplus r_{2})&=A(r_{1},r_{1})+A(r_{1},r_{2})\\&=A(r_{1},r_{2})\\&=1\end{aligned}}$ (since $r 1 = ⁠ 1 / 0 ⁠$ and $r 2 = ⁠ 0 / 1 ⁠$ , its area must be 1).

Farey neighbours and continued fractions

Fractions that appear as neighbours in a Farey sequence have closely related continued fraction expansions. Every fraction has two continued fraction expansions — in one the final term is 1; in the other the final term is greater by 1. If $⁠ p / q ⁠$ , which first appears in Farey sequence $F q$ , has the continued fraction expansions ${\begin{aligned}&[0;\ a_{1},\ a_{2},\ \ldots ,\ a_{n-1},\ a_{n},\ 1]\\{}&[0;\ a_{1},\ a_{2},\ \ldots ,\ a_{n-1},\ a_{n}+1]\end{aligned}}$

then the nearest neighbour of $⁠ p / q ⁠$ in $F q$ (which will be its neighbour with the larger denominator) has a continued fraction expansion $[0;\ a_{1},\ a_{2},\ \ldots ,\ a_{n}]$

and its other neighbour has a continued fraction expansion $[0;\ a_{1},\ a_{2},\ \ldots ,\ a_{n-1}]$

For example, ⁠3/8⁠ has the two continued fraction expansions $[0; 2, 1, 1, 1]$ and $[0; 2, 1, 2]$ , and its neighbours in $F 8$ are ⁠2/5⁠, which can be expanded as $[0; 2, 1, 1]$ ; and ⁠1/3⁠, which can be expanded as $[0; 2, 1]$ .

Farey fractions and the least common multiple

The lcm can be expressed as the products of Farey fractions as ${\text{lcm}}[1,2,...,N]=e^{\psi (N)}={\frac {1}{2}}\left(\prod _{r\in F_{N},0<r\leq 1/2}2\sin(\pi r)\right)^{2}$

where $ψ (N)$ is the second Chebyshev function.^[9]^[10]

Farey fractions and the greatest common divisor

Since the Euler's totient function is directly connected to the gcd so is the number of elements in $F n$ , $|F_{n}|=1+\sum _{m=1}^{n}\varphi (m)=1+\sum \limits _{m=1}^{n}\sum \limits _{k=1}^{m}\gcd(k,m)\cos {2\pi {\frac {k}{m}}}.$

For any 3 Farey fractions $⁠ a / b ⁠, ⁠ c / d ⁠, ⁠ e / f ⁠$ the following identity between the gcd's of the 2x2 matrix determinants in absolute value holds:^[11] $\gcd \left({\begin{Vmatrix}a&c\\b&d\end{Vmatrix}},{\begin{Vmatrix}a&e\\b&f\end{Vmatrix}}\right)=\gcd \left({\begin{Vmatrix}a&c\\b&d\end{Vmatrix}},{\begin{Vmatrix}c&e\\d&f\end{Vmatrix}}\right)=\gcd \left({\begin{Vmatrix}a&e\\b&f\end{Vmatrix}},{\begin{Vmatrix}c&e\\d&f\end{Vmatrix}}\right)$

^[7]

Applications

Farey sequences are very useful to find rational approximations of irrational numbers.^[12] For example, the construction by Eliahou^[13] of a lower bound on the length of non-trivial cycles in the 3x+1 process uses Farey sequences to calculate a continued fraction expansion of the number $log 2 (3)$ .

In physical systems with resonance phenomena, Farey sequences provide a very elegant and efficient method to compute resonance locations in 1D^[14] and 2D.^[15]

Farey sequences are prominent in studies of any-angle path planning on square-celled grids, for example in characterizing their computational complexity^[16] or optimality.^[17] The connection can be considered in terms of $r$ -constrained paths, namely paths made up of line segments that each traverse at most $r$ rows and at most $r$ columns of cells. Let $Q$ be the set of vectors $(q, p)$ such that $1\leq q\leq r$ , $0\leq p\leq q$ , and $p$ , $q$ are coprime. Let $Q*$ be the result of reflecting $Q$ in the line $y = x$ . Let $S=\{(\pm x,\pm y):(x,y)\in Q\cup Q*\}$ . Then any $r$ -constrained path can be described as a sequence of vectors from $S$ . There is a bijection between $Q$ and the Farey sequence of order $r$ given by $(q, p)$ mapping to ${\tfrac {p}{q}}$ .

Ford circles

There is a connection between Farey sequence and Ford circles.

For every fraction $⁠ p / q ⁠$ (in its lowest terms) there is a Ford circle $C [p / q]$ , which is the circle with radius ${\tfrac {1}{2q^{2}}}$ and centre at ${\bigl (}{\tfrac {p}{q}},{\tfrac {1}{2q^{2}}}{\bigr )}.$ Two Ford circles for different fractions are either disjoint or they are tangent to one another—two Ford circles never intersect. If $0 < ⁠ p / q ⁠ < 1$ then the Ford circles that are tangent to $C [p / q]$ are precisely the Ford circles for fractions that are neighbours of $⁠ p / q ⁠$ in some Farey sequence.

Thus $C [2/5]$ is tangent to $C [1/2]$ , $C [1/3]$ , $C [3/7]$ , $C [3/8]$ , etc.

Ford circles appear also in the Apollonian gasket $(0,0,1,1)$ . The picture below illustrates this together with Farey resonance lines.^[18]

Riemann hypothesis

Farey sequences are used in two equivalent formulations of the Riemann hypothesis. Suppose the terms of $F n$ are $\{a_{k,n}:k=0,1,\ldots ,m_{n}\}.$ Define $d_{k,n}=a_{k,n}-{\tfrac {k}{m_{n}}},$ in other words $d_{k,n}$ is the difference between the $k$ th term of the $n$ th Farey sequence, and the $k$ th member of a set of the same number of points, distributed evenly on the unit interval. In 1924 Jérôme Franel ^[19] proved that the statement

$\sum _{k=1}^{m_{n}}d_{k,n}^{2}=O(n^{r})\quad \forall r>-1$

is equivalent to the Riemann hypothesis, and then Edmund Landau ^[20] remarked (just after Franel's paper) that the statement $\sum _{k=1}^{m_{n}}|d_{k,n}|=O(n^{r})\quad \forall r>{\frac {1}{2}}$

is also equivalent to the Riemann hypothesis.

Other sums involving Farey fractions

The sum of all Farey fractions of order $n$ is half the number of elements: $\sum _{r\in F_{n}}r={\frac {1}{2}}|F_{n}|.$

The sum of the denominators in the Farey sequence is twice the sum of the numerators and relates to Euler's totient function:

$\sum _{a/b\in F_{n}}b=2\sum _{a/b\in F_{n}}a=1+\sum _{i=1}^{n}i\varphi (i),$

which was conjectured by Harold L. Aaron in 1962 and demonstrated by Jean A. Blake in 1966.^[21] A one line proof of the Harold L. Aaron conjecture is as follows. The sum of the numerators is $1+\sum _{2\leq b\leq n}\ \sum _{(a,b)=1}a=1+\sum _{2\leq b\leq n}b{\frac {\varphi (b)}{2}}.$ The sum of denominators is $2+\sum _{2\leq b\leq n}\ \sum _{(a,b)=1}b=2+\sum _{2\leq b\leq n}b\varphi (b).$ The quotient of the first sum by the second sum is ⁠1/2⁠.

Let $b j$ be the ordered denominators of $F n$ , then:^[22]

$\sum _{j=0}^{|F_{n}|-1}{\frac {b_{j}}{b_{j+1}}}={\frac {3|F_{n}|-4}{2}}$ and $\sum _{j=0}^{|F_{n}|-1}{\frac {1}{b_{j+1}b_{j}}}=1.$

Let ${\tfrac {a_{j}}{b_{j}}}$ the $j$ th Farey fraction in $F n$ , then $\sum _{j=1}^{|F_{n}|-1}(a_{j-1}b_{j+1}-a_{j+1}b_{j-1})=\sum _{j=1}^{|F_{n}|-1}{\begin{Vmatrix}a_{j-1}&a_{j+1}\\b_{j-1}&b_{j+1}\end{Vmatrix}}=3(|F_{n}|-1)-2n-1,$

which is demonstrated in.^[23] Also according to this reference the term inside the sum can be expressed in many different ways: $a_{j-1}b_{j+1}-a_{j+1}b_{j-1}={\frac {b_{j-1}+b_{j+1}}{b_{j}}}={\frac {a_{j-1}+a_{j+1}}{a_{j}}}=\left\lfloor {\frac {n+b_{j-1}}{b_{j}}}\right\rfloor ,$

obtaining thus many different sums over the Farey elements with same result. Using the symmetry around 1/2 the former sum can be limited to half of the sequence as

$\sum _{j=1}^{\left\lfloor {\frac {|F_{n}|}{2}}\right\rfloor }(a_{j-1}b_{j+1}-a_{j+1}b_{j-1})={\frac {3(|F_{n}|-1)}{2}}-n-\left\lceil {\frac {n}{2}}\right\rceil ,$

The Mertens function can be expressed as a sum over Farey fractions as $M(n)=-1+\sum _{a\in {\mathcal {F}}_{n}}e^{2\pi ia}$ where ${\mathcal {F}}_{n}$ is the Farey sequence of order $n$ .

This formula is used in the proof of the Franel–Landau theorem.^[24]

Next term

A surprisingly simple algorithm exists to generate the terms of F_n in either traditional order (ascending) or non-traditional order (descending). The algorithm computes each successive entry in terms of the previous two entries using the mediant property given above. If ⁠a/b⁠ and ⁠c/d⁠ are the two given entries, and ⁠p/q⁠ is the unknown next entry, then $⁠ c / d ⁠ = ⁠ a + p / b + q ⁠$ . Since ⁠c/d⁠ is in lowest terms, there must be an integer k such that $kc = a + p$ and $kd = b + q$ , giving $p = kc - a$ and $q = kd - b$ . If we consider p and q to be functions of k, then

{\frac {p(k)}{q(k)}}-{\frac {c}{d}}={\frac {cb-da}{d(kd-b)}}

so the larger k gets, the closer ⁠p/q⁠ gets to ⁠c/d⁠.

To give the next term in the sequence k must be as large as possible, subject to $kd - b \leq n$ (as we are only considering numbers with denominators not greater than n), so k is the greatest integer ≤ ⁠n + b/d⁠. Putting this value of k back into the equations for p and q gives

p=\left\lfloor {\frac {n+b}{d}}\right\rfloor c-a

q=\left\lfloor {\frac {n+b}{d}}\right\rfloor d-b

This is implemented in Python as follows:

fromfractionsimportFractionfromcollections.abcimportGeneratordeffarey_sequence(n:int,descending:bool=False)->Generator[Fraction]:"""    Print the n'th Farey sequence. Allow for either ascending or descending.    >>> print(*farey_sequence(5), sep=' ')    0 1/5 1/4 1/3 2/5 1/2 3/5 2/3 3/4 4/5 1    """a,b,c,d=0,1,1,nifdescending:a,c=1,n-1yieldFraction(a,b)while0<=c<=n:k=(n+b)//da,b,c,d=c,d,k*c-a,k*d-byieldFraction(a,b)if__name__=="__main__":importdoctestdoctest.testmod()

Brute-force searches for solutions to Diophantine equations in rationals can often take advantage of the Farey series (to search only reduced forms). While this code uses the first two terms of the sequence to initialize a, b, c, and d, one could substitute any pair of adjacent terms in order to exclude those less than (or greater than) a particular threshold.^[25]

Footnotes

↑ “The sequence of all reduced fractions with denominators not exceeding n, listed in order of their size, is called the Farey sequence of order n.” With the comment: “This definition of the Farey sequences seems to be the most convenient. However, some authors prefer to restrict the fractions to the interval from 0 to 1.” — Niven & Zuckerman (1972)^[1]

Related Research Articles

In mathematics, the Fibonacci sequence is a sequence in which each number is the sum of the two preceding ones. Numbers that are part of the Fibonacci sequence are known as Fibonacci numbers, commonly denoted $F n$ . Many writers begin the sequence with 0 and 1, although some authors start it from 1 and 1 and some from 1 and 2. Starting from 0 and 1, the sequence begins

In mathematics, two quantities are in the golden ratio if their ratio is the same as the ratio of their sum to the larger of the two quantities. Expressed algebraically, for quantities $and with, is in a golden ratio to if$

In number theory, a multiplicative function is an arithmetic function f(n) of a positive integer n with the property that f(1) = 1 and $whenever a and b are coprime.$

<span class="mw-page-title-main">Euler's totient function</span> Number of integers coprime to and less than n

In number theory, Euler's totient function counts the positive integers up to a given integer $n$ that are relatively prime to $n$ . It is written using the Greek letter phi as $or, and may also be called Euler's phi function . In other words, it is the number of integers k in the range 1 \leq k \leq n for which the greatest common divisor gcd(n, k) is equal to 1. The integers k of this form are sometimes referred to as totatives of n .$

<span class="mw-page-title-main">Fourier series</span> Decomposition of periodic functions into sums of simpler sinusoidal forms

A Fourier series is an expansion of a periodic function into a sum of trigonometric functions. The Fourier series is an example of a trigonometric series, but not all trigonometric series are Fourier series. By expressing a function as a sum of sines and cosines, many problems involving the function become easier to analyze because trigonometric functions are well understood. For example, Fourier series were first used by Joseph Fourier to find solutions to the heat equation. This application is possible because the derivatives of trigonometric functions fall into simple patterns. Fourier series cannot be used to approximate arbitrary functions, because most functions have infinitely many terms in their Fourier series, and the series do not always converge. Well-behaved functions, for example smooth functions, have Fourier series that converge to the original function. The coefficients of the Fourier series are determined by integrals of the function multiplied by trigonometric functions, described in Common forms of the Fourier series below.

Noether's theorem states that every continuous symmetry of the action of a physical system with conservative forces has a corresponding conservation law. This is the first of two theorems published by mathematician Emmy Noether in 1918. The action of a physical system is the integral over time of a Lagrangian function, from which the system's behavior can be determined by the principle of least action. This theorem only applies to continuous and smooth symmetries of physical space.

In mathematics, a generating function is a representation of an infinite sequence of numbers as the coefficients of a formal power series. Generating functions are often expressed in closed form, by some expression involving operations on the formal series.

In linear algebra, the Cayley–Hamilton theorem states that every square matrix over a commutative ring satisfies its own characteristic equation.

<span class="mw-page-title-main">Hamiltonian mechanics</span> Formulation of classical mechanics using momenta

In physics, Hamiltonian mechanics is a reformulation of Lagrangian mechanics that emerged in 1833. Introduced by Sir William Rowan Hamilton, Hamiltonian mechanics replaces (generalized) velocities $used in Lagrangian mechanics with (generalized) momenta . Both theories provide interpretations of classical mechanics and describe the same physical phenomena.$

The Lucas sequence is an integer sequence named after the mathematician François Édouard Anatole Lucas (1842–1891), who studied both that sequence and the closely related Fibonacci sequence. Individual numbers in the Lucas sequence are known as Lucas numbers. Lucas numbers and Fibonacci numbers form complementary instances of Lucas sequences.

In mathematics, theta functions are special functions of several complex variables. They show up in many topics, including Abelian varieties, moduli spaces, quadratic forms, and solitons. As Grassmann algebras, they appear in quantum field theory.

In mathematics, the Gauss–Kuzmin–Wirsing operator is the transfer operator of the Gauss map that takes a positive number to the fractional part of its reciprocal. It is named after Carl Gauss, Rodion Kuzmin, and Eduard Wirsing. It occurs in the study of continued fractions; it is also related to the Riemann zeta function.

In mathematics, the mediant of two fractions, generally made up of four positive integers

In mathematics, the resultant of two polynomials is a polynomial expression of their coefficients that is equal to zero if and only if the polynomials have a common root, or, equivalently, a common factor. In some older texts, the resultant is also called the eliminant.

In algebra, the Bring radical or ultraradical of a real number a is the unique real root of the polynomial

In mathematics, in the area of wavelet analysis, a refinable function is a function which fulfils some kind of self-similarity. A function $is called refinable with respect to the mask if$

$<span class="mw-page-title-main">Rogers–Ramanujan continued fraction</span> Continued fraction closely related to the Rogers–Ramanujan identities$

The Rogers–Ramanujan continued fraction is a continued fraction discovered by Rogers (1894) and independently by Srinivasa Ramanujan, and closely related to the Rogers–Ramanujan identities. It can be evaluated explicitly for a broad class of values of its argument.

In mathematics, Ramanujan's master theorem, named after Srinivasa Ramanujan, is a technique that provides an analytic expression for the Mellin transform of an analytic function.

In mathematics, infinite compositions of analytic functions (ICAF) offer alternative formulations of analytic continued fractions, series, products and other infinite expansions, and the theory evolving from such compositions may shed light on the convergence/divergence of these expansions. Some functions can actually be expanded directly as infinite compositions. In addition, it is possible to use ICAF to evaluate solutions of fixed point equations involving infinite expansions. Complex dynamics offers another venue for iteration of systems of functions rather than a single function. For infinite compositions of a single function see Iterated function. For compositions of a finite number of functions, useful in fractal theory, see Iterated function system.

In modern valence bond (VB) theory calculations, Chirgwin–Coulson weights are the relative weights of a set of possible VB structures of a molecule. Related methods of finding the relative weights of valence bond structures are the Löwdin and the inverse weights.

References

↑ Niven, Ivan M.; Zuckerman, Herbert S. (1972). An Introduction to the Theory of Numbers (Third ed.). John Wiley and Sons. Definition 6.1.
↑ Guthery, Scott B. (2011). "1. The Mediant". A Motif of Mathematics: History and Application of the Mediant and the Farey Sequence. Boston: Docent Press. p. 7. ISBN 978-1-4538-1057-6. OCLC 1031694495 . Retrieved 28 September 2020.
↑ Hardy, G.H.; Wright, E.M. (1979). An Introduction to the Theory of Numbers (Fifth ed.). Oxford University Press. Chapter III. ISBN 0-19-853171-0.
1 2 Beiler, Albert H. (1964). Recreations in the Theory of Numbers (Second ed.). Dover. Chapter XVI. ISBN 0-486-21096-0. Cited in "Farey Series, A Story". Cut-the-Knot.
↑ Sloane, N. J. A. (ed.). "SequenceA005728". The On-Line Encyclopedia of Integer Sequences . OEIS Foundation.
↑ Tomas Garcia, Rogelio (July 2024). "Farey Fractions with Equal Numerators and the Rank of Unit Fractions" (PDF). Integers. 24.
1 2 Tomas, Rogelio (January 2022). "Partial Franel sums" (PDF). Journal of Integer Sequences. 25 (1).
↑ Austin, David (December 2008). "Trees, Teeth, and Time: The mathematics of clock making". American Mathematical Society . Rhode Island. Archived from the original on 4 February 2020. Retrieved 28 September 2020.
↑ Martin, Greg (2009). "A product of Gamma function values at fractions with the same denominator". arXiv: 0907.4384 [math.CA].
↑ Wehmeier, Stefan (2009). "The LCM(1,2,...,n) as a product of sine values sampled over the points in Farey sequences". arXiv: 0909.1838 [math.CA].
↑ Tomas Garcia, Rogelio (August 2020). "Equalities between greatest common divisors involving three coprime pairs" (PDF). Notes on Number Theory and Discrete Mathematics. 26 (3): 5–7. doi: 10.7546/nntdm.2020.26.3.5-7 . S2CID 225280271.
↑ "Farey Approximation". NRICH.maths.org. Archived from the original on 19 November 2018. Retrieved 18 November 2018.
↑ Eliahou, Shalom (August 1993). "The 3x+1 problem: new lower bounds on nontrivial cycle lengths". Discrete Mathematics. 118 (1–3): 45–56. doi: 10.1016/0012-365X(93)90052-U .
↑ Zhenhua Li, A.; Harter, W.G. (2015). "Quantum Revivals of Morse Oscillators and Farey–Ford Geometry". Chem. Phys. Lett. 633: 208–213. arXiv: 1308.4470 . Bibcode:2015CPL...633..208L. doi:10.1016/j.cplett.2015.05.035. S2CID 66213897.
↑ Tomas, R. (2014). "From Farey sequences to resonance diagrams" (PDF). Physical Review Special Topics - Accelerators and Beams. 17 (1): 014001. Bibcode:2014PhRvS..17a4001T. doi: 10.1103/PhysRevSTAB.17.014001 .
↑ Harabor, Daniel Damir; Grastien, Alban; Öz, Dindar; Aksakalli, Vural (26 May 2016). "Optimal Any-Angle Pathfinding In Practice". Journal of Artificial Intelligence Research. 56: 89–118. doi: 10.1613/jair.5007 .
↑ Hew, Patrick Chisan (19 August 2017). "The Length of Shortest Vertex Paths in Binary Occupancy Grids Compared to Shortest r-Constrained Ones". Journal of Artificial Intelligence Research. 59: 543–563. doi: 10.1613/jair.5442 .
↑ Tomas, Rogelio (2020). "Imperfections and corrections". arXiv: 2006.10661 [physics.acc-ph].
↑ Franel, Jérôme (1924). "Les suites de Farey et le problème des nombres premiers". Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen. Mathematisch-Physikalische Klasse (in French): 198–201.
↑ Landau, Edmund (1924). "Bemerkungen zu der vorstehenden Abhandlung von Herrn Franel". Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen. Mathematisch-Physikalische Klasse (in German): 202–206.
↑ Blake, Jean A. (1966). "Some Characteristic Properties of the Farey Series". The American Mathematical Monthly. 73 (1): 50–52. doi:10.2307/2313922. JSTOR 2313922.
↑ Kurt Girstmair; Girstmair, Kurt (2010). "Farey Sums and Dedekind Sums". The American Mathematical Monthly. 117 (1): 72–78. doi:10.4169/000298910X475005. JSTOR 10.4169/000298910X475005. S2CID 31933470.
↑ Hall, R. R.; Shiu, P. (2003). "The Index of a Farey Sequence". Michigan Math. J. 51 (1): 209–223. doi: 10.1307/mmj/1049832901 .
↑ Edwards, Harold M. (1974). "12.2 Miscellany. The Riemann Hypothesis and Farey Series" . In Smith, Paul A.; Ellenberg, Samuel (eds.). Riemann's Zeta Function. Pure and Applied Mathematics. New York: Academic Press. pp. 263–267. ISBN 978-0-08-087373-2. OCLC 316553016 . Retrieved 30 September 2020.
↑ Routledge, Norman (March 2008). "Computing Farey series". The Mathematical Gazette . Vol. 92, no. 523. pp. 55–62.

External links

Hatcher, Allen. "Topology of Numbers" (PDF). Online copy of book
Bogomolny, Alexander. "Farey series". Cut-the-Knot.
Bogomolny, Alexander. "Stern-Brocot Tree". Cut-the-Knot.
Pennestri, Ettore. "A Brocot table of base 120".
"Farey series", Encyclopedia of Mathematics , EMS Press, 2001 [1994]
Weisstein, Eric W. "Stern-Brocot Tree". MathWorld .
OEIS sequenceA005728(Number of fractions in Farey series of order n)
OEIS sequenceA006842(Numerators of Farey series of order n)
OEIS sequenceA006843(Denominators of Farey series of order n)
Archived at Ghostarchive and the Wayback Machine : Bonahon, Francis. Funny Fractions and Ford Circles (video). Brady Haran . Retrieved 9 June 2015– via YouTube.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[2] “The sequence of all reduced fractions with denominators not exceeding n, listed in order of their size, is called the Farey sequence of order n.” With the comment: “This definition of the Farey sequences seems to be the most convenient. However, some authors prefer to restrict the fractions to the interval from 0 to 1.” — Niven & Zuckerman (1972)^[1]

[1] Niven, Ivan M.; Zuckerman, Herbert S. (1972). An Introduction to the Theory of Numbers (Third ed.). John Wiley and Sons. Definition 6.1.

[3] Guthery, Scott B. (2011). "1. The Mediant". A Motif of Mathematics: History and Application of the Mediant and the Farey Sequence. Boston: Docent Press. p. 7. ISBN 978-1-4538-1057-6. OCLC 1031694495 . Retrieved 28 September 2020.

[4] Hardy, G.H.; Wright, E.M. (1979). An Introduction to the Theory of Numbers (Fifth ed.). Oxford University Press. Chapter III. ISBN 0-19-853171-0.

[Beiler-5] 1 2 Beiler, Albert H. (1964). Recreations in the Theory of Numbers (Second ed.). Dover. Chapter XVI. ISBN 0-486-21096-0. Cited in "Farey Series, A Story". Cut-the-Knot.

[6] Sloane, N. J. A. (ed.). "SequenceA005728". The On-Line Encyclopedia of Integer Sequences . OEIS Foundation.

[Tomas2024-7] Tomas Garcia, Rogelio (July 2024). "Farey Fractions with Equal Numerators and the Rank of Unit Fractions" (PDF). Integers. 24.

[Tomas2018-8] 1 2 Tomas, Rogelio (January 2022). "Partial Franel sums" (PDF). Journal of Integer Sequences. 25 (1).

[Austin2008-9] Austin, David (December 2008). "Trees, Teeth, and Time: The mathematics of clock making". American Mathematical Society . Rhode Island. Archived from the original on 4 February 2020. Retrieved 28 September 2020.

[10] Martin, Greg (2009). "A product of Gamma function values at fractions with the same denominator". arXiv: 0907.4384 [math.CA].

[11] Wehmeier, Stefan (2009). "The LCM(1,2,...,n) as a product of sine values sampled over the points in Farey sequences". arXiv: 0909.1838 [math.CA].

[TomasGarcia2020-12] Tomas Garcia, Rogelio (August 2020). "Equalities between greatest common divisors involving three coprime pairs" (PDF). Notes on Number Theory and Discrete Mathematics. 26 (3): 5–7. doi: 10.7546/nntdm.2020.26.3.5-7 . S2CID 225280271.

[13] "Farey Approximation". NRICH.maths.org. Archived from the original on 19 November 2018. Retrieved 18 November 2018.

[14] Eliahou, Shalom (August 1993). "The 3x+1 problem: new lower bounds on nontrivial cycle lengths". Discrete Mathematics. 118 (1–3): 45–56. doi: 10.1016/0012-365X(93)90052-U .

[15] Zhenhua Li, A.; Harter, W.G. (2015). "Quantum Revivals of Morse Oscillators and Farey–Ford Geometry". Chem. Phys. Lett. 633: 208–213. arXiv: 1308.4470 . Bibcode:2015CPL...633..208L. doi:10.1016/j.cplett.2015.05.035. S2CID 66213897.

[16] Tomas, R. (2014). "From Farey sequences to resonance diagrams" (PDF). Physical Review Special Topics - Accelerators and Beams. 17 (1): 014001. Bibcode:2014PhRvS..17a4001T. doi: 10.1103/PhysRevSTAB.17.014001 .

[17] Harabor, Daniel Damir; Grastien, Alban; Öz, Dindar; Aksakalli, Vural (26 May 2016). "Optimal Any-Angle Pathfinding In Practice". Journal of Artificial Intelligence Research. 56: 89–118. doi: 10.1613/jair.5007 .

[18] Hew, Patrick Chisan (19 August 2017). "The Length of Shortest Vertex Paths in Binary Occupancy Grids Compared to Shortest r-Constrained Ones". Journal of Artificial Intelligence Research. 59: 543–563. doi: 10.1613/jair.5442 .

[Tomas2020-19] Tomas, Rogelio (2020). "Imperfections and corrections". arXiv: 2006.10661 [physics.acc-ph].

[20] Franel, Jérôme (1924). "Les suites de Farey et le problème des nombres premiers". Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen. Mathematisch-Physikalische Klasse (in French): 198–201.

[21] Landau, Edmund (1924). "Bemerkungen zu der vorstehenden Abhandlung von Herrn Franel". Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen. Mathematisch-Physikalische Klasse (in German): 202–206.

[22] Blake, Jean A. (1966). "Some Characteristic Properties of the Farey Series". The American Mathematical Monthly. 73 (1): 50–52. doi:10.2307/2313922. JSTOR 2313922.

[23] Kurt Girstmair; Girstmair, Kurt (2010). "Farey Sums and Dedekind Sums". The American Mathematical Monthly. 117 (1): 72–78. doi:10.4169/000298910X475005. JSTOR 10.4169/000298910X475005. S2CID 31933470.

[24] Hall, R. R.; Shiu, P. (2003). "The Index of a Farey Sequence". Michigan Math. J. 51 (1): 209–223. doi: 10.1307/mmj/1049832901 .

[25] Edwards, Harold M. (1974). "12.2 Miscellany. The Riemann Hypothesis and Farey Series" . In Smith, Paul A.; Ellenberg, Samuel (eds.). Riemann's Zeta Function. Pure and Applied Mathematics. New York: Academic Press. pp. 263–267. ISBN 978-0-08-087373-2. OCLC 316553016 . Retrieved 30 September 2020.

[26] Routledge, Norman (March 2008). "Computing Farey series". The Mathematical Gazette . Vol. 92, no. 523. pp. 55–62.

[a]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]

[16]

[17]

[18]

[19]

[20]

[21]

[22]

[23]

[24]

[25]

[1]