Английская Википедия:Farey sequence
In mathematics, the Farey sequence of order n is the sequence of completely reduced fractions, either between 0 and 1, or without this restriction,Шаблон:Efn which when in lowest terms have denominators less than or equal to n, arranged in order of increasing size.
With the restricted definition, each Farey sequence starts with the value 0, denoted by the fraction Шаблон:Sfrac, and ends with the value 1, denoted by the fraction Шаблон:Sfrac (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.[1]
Examples
The Farey sequences of orders 1 to 8 are :
- F1 = { Шаблон:Sfrac, Шаблон:Sfrac }
- F2 = { Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac }
- F3 = { Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac }
- F4 = { Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac }
- F5 = { Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac }
- F6 = { Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac }
- F7 = { Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac }
- F8 = { Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac, Шаблон:Sfrac }
| 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 Шаблон:Mvar6.
Reflecting this shape around the diagonal and main axes generates the Farey sunburst, shown below. The Farey sunburst of order Шаблон:Mvar connects the visible integer grid points from the origin in the square of side 2Шаблон:Mvar, centered at the origin. Using Pick's theorem, the area of the sunburst is 4(|Шаблон:Mvarn|−1), where |Шаблон:Mvarn| is the [[#Sequence_length_and_index_of_a_fraction|number of fractions in Шаблон:Mvarn]].
History
- The history of 'Farey series' is very curious — Hardy & Wright (1979)[2]
- ... 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)[3]
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.[3] 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 Fn contains all of the members of Fn−1 and also contains an additional fraction for each number that is less than n and coprime to n. Thus F6 consists of F5 together with the fractions Шаблон:Sfrac and Шаблон:Sfrac.
The middle term of a Farey sequence Fn is always Шаблон:Sfrac, for n > 1. From this, we can relate the lengths of Fn and Fn−1 using Euler's totient function <math>\varphi(n)</math> :
- <math>|F_n| = |F_{n-1}| + \varphi(n).</math>
Using the fact that |F1| = 2, we can derive an expression for the length of Fn:[4]
- <math>|F_n| = 1 + \sum_{m=1}^n \varphi(m) = 1 + \Phi(n),</math>
where <math>\Phi(n)</math> is the summatory totient.
We also have :
- <math>|F_n| = \frac{1}{2}\left(3+\sum_{d=1}^n\mu(d)\left\lfloor\tfrac{n}{d}\right\rfloor^2\right),</math>
and by a Möbius inversion formula :
- <math>|F_n| = \frac{1}{2}(n+3)n-\sum_{d=2}^n|F_{\lfloor n/d\rfloor}|,</math>
where µ(d) is the number-theoretic Möbius function, and <math>\lfloor \tfrac{n}{d} \rfloor </math> is the floor function.
The asymptotic behaviour of |Fn| is :
- <math>|F_n| \sim \frac {3n^2}{\pi^2}.</math>
The index <math>I_n(a_{k,n})=k</math> of a fraction <math>a_{k,n}</math> in the Farey sequence <math>F_n=\{a_{k,n} : k = 0, 1, \ldots, m_n\}</math> is simply the position that <math>a_{k,n}</math> 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:
- <math>I_n(0/1) = 0,</math>
- <math>I_n(1/n) = 1,</math>
- <math>I_n(1/2) = (|F_n|-1)/2,</math>
- <math>I_n(1/1) = |F_n|-1 ,</math>
- <math>I_n(h/k) = |F_n|-1-I_n((k-h)/k).</math>
The index of <math>1/k</math> where <math>n/(i+1) < k \leq n/i </math> and <math>n</math> is the least common multiple of the first <math>i</math> numbers, <math>n={\rm lcm}([2,i]) </math>, is given by:[5]
- <math>I_n(1/k) = 1 + n \sum_{j=1}^{i} \frac{\varphi(j)}{j} - k\Phi(i).</math>
Farey neighbours
Fractions which are neighbouring terms in any Farey sequence are known as a Farey pair and have the following properties.
If Шаблон:Sfrac and Шаблон:Sfrac are neighbours in a Farey sequence, with Шаблон:Sfrac < Шаблон:Sfrac, then their difference Шаблон:Sfrac − Шаблон:Sfrac is equal to Шаблон:Sfrac. Since
- <math>\frac{c}{d} - \frac{a}{b} = \frac{bc - ad}{bd},</math>
this is equivalent to saying that
- <math>bc - ad = 1</math>.
Thus Шаблон:Sfrac and Шаблон:Sfrac are neighbours in F5, and their difference is Шаблон:Sfrac.
The converse is also true. If
- <math>bc - ad = 1</math>
for positive integers a,b,c and d with a < b and c < d then Шаблон:Sfrac and Шаблон:Sfrac will be neighbours in the Farey sequence of order max(b,d).
If Шаблон:Sfrac has neighbours Шаблон:Sfrac and Шаблон:Sfrac in some Farey sequence, with
- <math>\frac{a}{b} < \frac{p}{q} < \frac{c}{d} </math>
then Шаблон:Sfrac is the mediant of Шаблон:Sfrac and Шаблон:Sfrac – in other words,
- <math>\frac{p}{q} = \frac{a + c}{b + d}.</math>
This follows easily from the previous property, since if Шаблон:Nowrap, then Шаблон:Nowrap, Шаблон:Nowrap, Шаблон:Nowrap.
It follows that if Шаблон:Sfrac and Шаблон:Sfrac are neighbours in a Farey sequence then the first term that appears between them as the order of the Farey sequence is incremented is
- <math>\frac{a+c}{b+d},</math>
which first appears in the Farey sequence of order Шаблон:Nowrap.
Thus the first term to appear between Шаблон:Sfrac and Шаблон:Sfrac is Шаблон:Sfrac, which appears in F8.
The total number of Farey neighbour pairs in Fn is 2|Fn| − 3.
The Stern–Brocot tree is a data structure showing how the sequence is built up from 0 (= Шаблон:Sfrac) and 1 (= Шаблон:Sfrac), by taking successive mediants.
Equivalent-area interpretation
Every consecutive pair of Farey rationals have an equivalent area of 1.[6] See this by interpreting consecutive rationals r1 = p/q and r2 = p′/q′ as vectors (p, q) in the x–y plane. The area of A(p/q, p′/q′) is given by qp′ − q′p. As any added fraction in between two previous consecutive Farey sequence fractions is calculated as the mediant (⊕), then A(r1, r1 ⊕ r2) = A(r1, r1) + A(r1, r2) = A(r1, r2) = 1 (since r1 = 1/0 and r2 = 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 Шаблон:Sfrac, which first appears in Farey sequence Fq, has continued fraction expansions
- [0; a1, a2, ..., an − 1, an, 1]
- [0; a1, a2, ..., an − 1, an + 1]
then the nearest neighbour of Шаблон:Sfrac in Fq (which will be its neighbour with the larger denominator) has a continued fraction expansion
- [0; a1, a2, ..., an]
and its other neighbour has a continued fraction expansion
- [0; a1, a2, ..., an − 1]
For example, Шаблон:Sfrac has the two continued fraction expansions Шаблон:Nowrap and Шаблон:Nowrap, and its neighbours in F8 are Шаблон:Sfrac, which can be expanded as Шаблон:Nowrap; and Шаблон:Sfrac, which can be expanded as Шаблон:Nowrap.
Farey fractions and the least common multiple
The lcm can be expressed as the products of Farey fractions as
- <math> \text{lcm}[1,2,...,N] = e^{\psi(N)}=\frac{1}{2} \left( \prod_{r \in F_N, 0<r \le 1/2} 2 \sin(\pi r) \right)^2 </math>
where <math>\psi(N)</math> is the second Chebyshev function.[7][8]
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 Fn,
- <math>|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}} .</math>
For any 3 Farey fractions Шаблон:Sfrac, Шаблон:Sfrac and Шаблон:Sfrac the following identity between the gcd's of the 2x2 matrix determinants in absolute value holds:[9]
- <math>\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) </math>
Applications
Farey sequences are very useful to find rational approximations of irrational numbers.[10] For example, the construction by Eliahou[11] 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 log2(3).
In physical systems with resonance phenomena, Farey sequences provide a very elegant and efficient method to compute resonance locations in 1D[12] and 2D.[13]
Farey sequences are prominent in studies of any-angle path planning on square-celled grids, for example in characterizing their computational complexity[14] or optimality.[15] The connection can be considered in terms of r-constrained paths, namely paths made up of line segments that each traverse at most <math>r</math> rows and at most <math>r</math> columns of cells. Let <math>Q</math> be the set of vectors <math>(q,p)</math> such that <math>1 \leq q \leq r</math>, <math>0 \leq p \leq q</math>, and <math>p</math>, <math>q</math> are coprime. Let <math>Q*</math> be the result of reflecting <math>Q</math> in the line <math>y = x</math>. Let <math>S = \{ (\pm x, \pm y) : (x, y) \in Q \cup Q* \}</math>. Then any r-constrained path can be described as a sequence of vectors from <math>S</math>. There is a bijection between <math>Q</math> and the Farey sequence of order <math>r</math> given by <math>(q,p)</math> mapping to <math>\tfrac{p}{q}</math>.
Ford circles
There is a connection between Farey sequence and Ford circles.
For every fraction Шаблон:Sfrac (in its lowest terms) there is a Ford circle C[[[:Шаблон:Sfrac]]], which is the circle with radius 1/(2q2) and centre at (Шаблон:Sfrac, Шаблон:Sfrac). Two Ford circles for different fractions are either disjoint or they are tangent to one another—two Ford circles never intersect. If 0 < Шаблон:Sfrac < 1 then the Ford circles that are tangent to C[[[:Шаблон:Sfrac]]] are precisely the Ford circles for fractions that are neighbours of Шаблон:Sfrac in some Farey sequence.
Thus C[[[:Шаблон:Sfrac]]] is tangent to C[[[:Шаблон:Sfrac]]], C[[[:Шаблон:Sfrac]]], C[[[:Шаблон:Sfrac]]], C[[[:Шаблон:Sfrac]]], etc.
Ford circles appear also in the Apollonian gasket (0,0,1,1). The picture below illustrates this together with Farey resonance lines.[16]
Riemann hypothesis
Farey sequences are used in two equivalent formulations of the Riemann hypothesis. Suppose the terms of <math>F_n</math> are <math>\{a_{k,n} : k = 0, 1, \ldots, m_n\}</math>. Define <math>d_{k,n} = a_{k,n} - k/m_n</math>, in other words <math>d_{k,n}</math> is the difference between the kth term of the nth Farey sequence, and the kth member of a set of the same number of points, distributed evenly on the unit interval. In 1924 Jérôme Franel[17] proved that the statement
- <math>\sum_{k=1}^{m_n} d_{k,n}^2 = O (n^r)\quad\forall r>-1</math>
is equivalent to the Riemann hypothesis, and then Edmund Landau[18] remarked (just after Franel's paper) that the statement
- <math>\sum_{k=1}^{m_n} |d_{k,n}| = O (n^r)\quad\forall r>1/2</math>
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:
- <math>\sum_{r\in F_n} r = \frac{1}{2} |F_n| .</math>
The sum of the denominators in the Farey sequence is twice the sum of the numerators and relates to Euler's totient function:
- <math>\sum_{a/b \in F_n} b = 2 \sum_{a/b \in F_n} a = 1 + \sum_{i=1}^{n} i\varphi(i) , </math>
which was conjectured by Harold L. Aaron in 1962 and demonstrated by Jean A. Blake in 1966.[19] A one line proof of the Harold L. Aaron conjecture is as follows. The sum of the numerators is <math>{\displaystyle 1+ \sum_{ 2 \le b \le n} \sum_{(a,b)=1} a = 1+\sum_{ 2 \le b \le n} b\frac{\varphi(b)}{2}}</math>. The sum of denominators is <math>{\displaystyle 2+ \sum_{ 2 \le b \le n} \sum_{(a,b)=1} b = 2+\sum_{ 2 \le b \le n} b\varphi(b) }</math>. The quotient of the first sum by the second sum is <math>\frac{1}{2}</math>.
Let bj be the ordered denominators of Fn, then:[20]
- <math>\sum_{j=0}^{|F_n|-1} \frac{b_j}{b_{j+1}} = \frac{3|F_n|-4}{2} </math>
and
- <math>\sum_{j=0}^{|F_n|-1} \frac{1}{b_{j+1}b_{j}} = 1 .</math>
Let aj/bj the jth Farey fraction in Fn, then
- <math>\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 , </math>
which is demonstrated in.[21] Also according to this reference the term inside the sum can be expressed in many different ways:
- <math> 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, </math>
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
- <math>\sum_{j=1}^{\lfloor|F_n|/2\rfloor} (a_{j-1} b_{j+1} - a_{j+1} b_{j-1}) = 3(|F_n|-1)/2 - n- \lceil n/2 \rceil , </math>
The Mertens function can be expressed as a sum over Farey fractions as
- <math>M(n)= -1+ \sum_{a\in \mathcal{F}_n} e^{2\pi i a}</math> where <math> \mathcal{F}_n</math> is the Farey sequence of order n.
This formula is used in the proof of the Franel–Landau theorem.[22]
Next term
A surprisingly simple algorithm exists to generate the terms of Fn 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 Шаблон:Sfrac and Шаблон:Sfrac are the two given entries, and Шаблон:Sfrac is the unknown next entry, then Шаблон:Sfrac = Шаблон:Sfrac. Since Шаблон:Sfrac 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
- <math> \frac{p(k)}{q(k)}- \frac{c}{d} = \frac{cb - da}{d(kd - b)}</math>
so the larger k gets, the closer Шаблон:Sfrac gets to Шаблон:Sfrac.
To give the next term in the sequence k must be as large as possible, subject to kd − b ≤ n (as we are only considering numbers with denominators not greater than n), so k is the greatest integer ≤ Шаблон:Sfrac. Putting this value of k back into the equations for p and q gives
- <math> p = \left\lfloor\frac{n+b}{d}\right\rfloor c - a</math>
- <math> q = \left\lfloor\frac{n+b}{d}\right\rfloor d - b</math>
This is implemented in Python as follows:
def farey_sequence(n: int, descending: bool = False) -> None:
"""Print the n'th Farey sequence. Allow for either ascending or descending."""
a, b, c, d = 0, 1, 1, n
if descending:
a, c = 1, n - 1
print(f"{a}/{b}")
while 0 <= c <= n:
k = (n + b) // d
a, b, c, d = c, d, k * c - a, k * d - b
print(f"{a}/{b}")
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.[23]
See also
Footnotes
References
Further reading
- Шаблон:Citation
- Шаблон:Cite book — in particular, see §4.5 (pp. 115–123), Bonus Problem 4.61 (pp. 150, 523–524), §4.9 (pp. 133–139), §9.3, Problem 9.3.6 (pp. 462–463).
- Шаблон:Cite web — reviews the isomorphisms of the Stern-Brocot Tree.
- Шаблон:Cite web — reviews connections between Farey Fractions and Fractals.
- Шаблон:Cite journal Шаблон:Cite web Шаблон:Cite web
- Шаблон:Cite book Errata + CodeШаблон:Cbignore
External links
- Шаблон:Cite web Online copy of book
- Шаблон:Cite web
- Шаблон:Cite web
- Шаблон:Cite web
- Шаблон:Springer
- Шаблон:MathWorld
- Шаблон:OEIS el
- Шаблон:OEIS el
- Шаблон:OEIS el
- Archived at GhostarchiveШаблон:Cbignore and the Wayback MachineШаблон:Cbignore: Шаблон:Cite AV mediaШаблон:Cbignore
- ↑ Шаблон:Cite book
- ↑ Шаблон:Cite book
- ↑ 3,0 3,1 Шаблон:Cite book Cited in Шаблон:Cite web
- ↑ Шаблон:Cite OEIS
- ↑ 5,0 5,1 Шаблон:Cite journal
- ↑ Шаблон:Cite web
- ↑ Шаблон:Cite arXiv
- ↑ Шаблон:Cite arXiv
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite web
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite arXiv
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite journal
- ↑ Шаблон:Cite book
- ↑ Шаблон:Cite magazine
