This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A097546 #19 Aug 29 2024 21:01:25
%S A097546 0,1,1,1,1,1,2,3,4,5,6,7,8,15,22,29,36,43,50,57,64,71,78,85,92,99,106,
%T A097546 113,219,332,445,558,671,784,897,1010,1123,1236,1349,1462,1575,1688,
%U A097546 1801,1914,2027,2140,2253,2366,2479,2592,2705,2818,2931,3044,3157,3270
%N A097546 Denominators of "Farey fraction" approximations to Pi.
%C A097546 Given a real number x >= 1 (here x = Pi), start with 1/0 and 0/1 and construct the sequence of fractions f_n = r_n/s_n such that:
%C A097546 f_{n+1} = (r_k + r_n)/(s_k + s_n) where k is the greatest integer < n such that f_k <= x <= f_n. Sequence gives values s_n.
%C A097546 Write a 0 if f_n <= x and a 1 if f_n > x. This gives (for x = Pi) the sequence 1, 0, 0, 0, 1, 1, 1, 1, 0 (7 times), 1 (15 times), 0, 1, ... Ignore the initial string 1, 0, 0, 0, which is always the same. Look at the run lengths of the remaining sequence, which are in this case L_1 = 4, L_2 = 7, L_3 = 15, L_4 = 1, L_5 = 292, etc. (A001203). Christoffel showed that x has the continued fraction representation (L_1 - 1) + 1/(L_2 + 1/(L_3 + 1/(L_4 + ...))).
%D A097546 C. Brezinski, History of Continued Fractions and Padé Approximants, Springer-Verlag, 1991; pp. 151-152.
%D A097546 E. B. Christoffel, Observatio arithmetica, Ann. Math. Pura Appl., (II) 6 (1875), 148-153.
%H A097546 Dave Rusin, <a href="http://www.math.niu.edu/~rusin/known-math/99/farey">Farey fractions on sci.math</a> [Broken link]
%H A097546 Dave Rusin, <a href="/A002965/a002965.txt">Farey fractions on sci.math</a> [Cached copy]
%e A097546 The fractions are 1/0, 0/1, 1/1, 2/1, 3/1, 4/1, 7/2, 10/3, 13/4, 16/5, 19/6, 22/7, 25/8, 47/15, ...
%t A097546 f[x_, n_] := (m = Floor[x]; f0 = {m, m+1/2, m+1};
%t A097546 r = ({a___, b_, c_, d___} /; b < x < c) :> {b, (Numerator[b] + Numerator[c]) / (Denominator[b] + Denominator[c]), c}; Join[{m, m+1}, NestList[# /. r &, f0, n-3][[All, 2]]]); Join[{1, 0, 1, 2}, f[Pi, 52]] // Denominator (* _Jean-François Alcover_, May 18 2011 *)
%Y A097546 Cf. A097545.
%K A097546 nonn,frac,nice,easy
%O A097546 0,7
%A A097546 _N. J. A. Sloane_, Aug 28 2004
%E A097546 Corrected and extended by _Joshua Zucker_, May 08 2006