Calkin–Wilf tree

In number theory, the Calkin–Wilf tree is a tree in which the vertices correspond 1-for-1 to the positive rational numbers. The tree is rooted at the number 1, and any rational number expressed in simplest terms as the fraction a/b has as its two children the numbers a/(a + b) and (a + b)/b. Every positive rational number appears exactly once in the tree.

The sequence of rational numbers in a breadth-first traversal of the Calkin–Wilf tree is known as the Calkin–Wilf sequence. Its sequence of numerators (or, offset by one, denominators) is Stern's diatomic series, and can be computed by the fusc function.

The Calkin–Wilf tree is named after Neil Calkin and Herbert Wilf, who considered it in their 2000 paper. The tree was introduced earlier by Jean Berstel and Aldo de Luca^[1] as Raney tree, since they drew some ideas from a paper by George N. Raney.^[2] Stern's diatomic series was formulated much earlier by Moritz Abraham Stern, a 19th-century German mathematician who also invented the closely related Stern–Brocot tree.

Definition and structure

The Calkin–Wilf tree may be defined as a directed graph in which each positive rational number a/b occurs as a vertex and has one outgoing edge to another vertex, its parent. We assume that a/b is in simplest terms; that is, the greatest common divisor of a and b is 1. If a/b < 1, the parent of a/b is a/(b − a); if a/b is greater than one, the parent of a/b is (a − b)/b. Thus, in either case, the parent is a fraction with a smaller sum of numerator and denominator, so repeated reduction of this type must eventually reach the number 1. As a graph with one outgoing edge per vertex and one root reachable by all other vertices, the Calkin–Wilf tree must indeed be a tree.

The children of any vertex in the Calkin–Wilf tree may be computed by inverting the formula for the parents of a vertex. Each vertex a/b has one child whose value is less than 1, a/(a + b), because this is the only value less than 1 whose parent formula leads back to a/b. Similarly, each vertex a/b has one child whose value is greater than 1, (a + b)/b.^[3]

Although it is a binary tree (each vertex has two children), the Calkin–Wilf tree is not a binary search tree: its inorder does not coincide with the sorted order of its vertices. However, it is closely related to a different binary search tree on the same set of vertices, the Stern–Brocot tree: the vertices at each level of the two trees coincide, and are related to each other by a bit-reversal permutation.^[4]

Breadth first traversal

The Calkin–Wilf sequence is the sequence of rational numbers generated by a breadth-first traversal of the Calkin–Wilf tree,

1/1, 1/2, 2/1, 1/3, 3/2, 2/3, 3/1, 1/4, 4/3, 3/5, 5/2, 2/5, 5/3, 3/4, ….

Because the Calkin–Wilf tree contains every positive rational number exactly once, so does this sequence.^[5] The denominator of each fraction equals the numerator of the next fraction in the sequence. The Calkin–Wilf sequence can also be generated directly by the formula

q_{i+1}={\frac {1}{2\lfloor q_{i}\rfloor -q_{i}+1}}

where $q_{i}$ denotes the ith number in the sequence, starting from $q_{1}=1$ , and $\lfloor q_{i}\rfloor$ represents the integral part.^[6]

It's also possible to calculate $q_{i}$ directly. If you count successively equal digits of the binary code of i you get the continued fraction of $q_{i}$ . You must start counting 1-digits at the end of the binary code. Then you continue counting 0-digits backwards, then you continue counting 1-digits backward and so on.

Example:
i = 1081 = $10000111001_{2}$ : The continued fraction is [1;2,3,4,1] hence $q_{1081}=53/37$ .
i = 1990 = $11111000110_{2}$ : The continued fraction is [0;1,2,3,5] hence $q_{1990}=37/53$ .

If you calculate the continued fraction of any $q_{i}$ you can determine i.
Example:
$q_{i}=3/4$ : The continued fraction is [0;1,3] hence i = $1110_{2}$ = 14.
$q_{i}=4/3$ : The continued fraction is [1;3]. But the length of the continued fraction must be a odd number. So you must use [1;2,1] instead which equals [1;3]. Hence i = $1001_{2}$ = 9.

Stern's diatomic sequence

Stern's diatomic sequence is the integer sequence

0, 1, 1, 2, 1, 3, 2, 3, 1, 4, 3, 5, 2, 5, 3, 4, … (sequence A002487 in the OEIS).

Using zero-based numbering, the nth value in the sequence is the value fusc(n) of the fusc function,^[7] defined by the recurrence relations fusc(2n) = fusc(n) and fusc(2n + 1) = fusc(n) + fusc(n + 1), with the base cases fusc(0) = 0 and fusc(1) = 1. The nth rational number in a breadth-first traversal of the Calkin–Wilf tree is the number fusc(n) / fusc(n + 1).^[8] Thus, the diatomic sequence forms both the sequence of numerators and the sequence of denominators of the numbers in the Calkin–Wilf sequence.

The function fusc(n + 1) is the number of odd binomial coefficients of the form $\scriptstyle {n-r \choose r},\ 0\leq 2r<n,$ ^[9] and also counts the number of ways of writing n as a sum of powers of two in which each power occurs at most twice. This can be seen from the recurrence defining fusc: the expressions as a sum of powers of two for an even number 2n either have no 1's in them (in which case they are formed by doubling each term an expression for n) or two 1's (in which case the rest of the expression is formed by doubling each term in an expression for n − 1), so the number of representations is the sum of the number of representations for n and for n − 1, matching the recurrence. Similarly, each representation for an odd number 2n + 1 is formed by doubling a representation for n and adding 1, again matching the recurrence.^[10] For instance,

6 = 4 + 2 = 4 + 1 + 1 = 2 + 2 + 1 + 1

has three representations as a sum of powers of two with at most two copies of each power, so fusc(6 + 1) = 3.

Relation to Stern–Brocot tree

The Calkin–Wilf tree resembles the Stern–Brocot tree in that both are binary trees with each positive rational number appearing exactly once. Additionally, the top levels of the two trees appear very similar, and in both trees the same numbers appear at the same levels. One tree can be obtained from the other by performing a bit-reversal permutation on the numbers at each level of the trees.^[4] Alternatively, the number at a given node of the Calkin–Wilf tree can be converted into the number at the same position in the Stern–Brocot tree, and vice versa, by a process involving the reversal of the continued fraction representations of these numbers.^[11] However, in other ways they have different properties: for instance, the Stern–Brocot tree is a binary search tree: the left-to-right traversal order of the tree is the same as the numerical order of the numbers in it. This property is not true of the Calkin–Wilf tree.

Notes

^ Berstel & de Luca (1997), Section 6.
^ Raney (1973).
^ The description here is dual to the original definition by Calkin and Wilf, which begins by defining the child relationship and derives the parent relationship as part of a proof that every rational appears once in the tree. As defined here, every rational appears once by definition, and instead the fact that the resulting structure is a tree requires a proof.
^ ^a ^b Gibbons, Lester & Bird (2006).
^ Calkin & Wilf (2000): "a list of all positive rational numbers, each appearing once and only once, can be made by writing down 1/1, then the fractions on the level just below the top of the tree, reading from left to right, then the fractions on the next level down, reading from left to right, etc." Gibbons, Lester & Bird (2006) discuss efficient functional programming techniques for performing this breadth first traversal.
^ Aigner & Ziegler (2004) credit this formula to Moshe Newman.
^ The fusc name was given to it in 1976 by Edsger W. Dijkstra; see EWD570 and EWD578.
^ Calkin & Wilf (2000), Theorem 1.
^ Carlitz (1964).
^ The OEIS entry credits this fact to Carlitz (1964) and to uncited work of Lind. However, Carlitz' paper describes a more restricted class of sums of powers of two, counted by fusc(n) instead of by fusc(n + 1).
^ Bates, Bunder & Tognetti (2010).

References

Aigner, Martin; Ziegler, Günter M. (2004), Proofs from THE BOOK (3rd ed.), Berlin; New York: Springer, pp. 94–97, ISBN 978-3-540-40460-6
Bates, Bruce; Bunder, Martin; Tognetti, Keith (2010), "Linking the Calkin-Wilf and Stern-Brocot trees", European Journal of Combinatorics, 31 (7): 1637–1661, doi:10.1016/j.ejc.2010.04.002, MR 2673006
Berstel, Jean; de Luca, Aldo (1997), "Sturmian words, Lyndon words and trees", Theoretical Computer Science, 178: 171–203, doi:10.1016/S0304-3975(96)00101-6
Calkin, Neil; Wilf, Herbert (2000), "Recounting the rationals" (PDF), American Mathematical Monthly, 107 (4), Mathematical Association of America: 360–363, doi:10.2307/2589182, JSTOR 2589182.
Carlitz, L. (1964), "A problem in partitions related to the Stirling numbers", Bulletin of the American Mathematical Society, 70 (2): 275–278, doi:10.1090/S0002-9904-1964-11118-6, MR 0157907.
Dijkstra, Edsger W. (1982), Selected Writings on Computing: A Personal Perspective, Springer-Verlag, ISBN 0-387-90652-5. EWD 570: An exercise for Dr.R.M.Burstall, pp. 215–216, and EWD 578: More about the function "fusc" (A sequel to EWD570), pp. 230–232, reprints of notes originally written in 1976.
Gibbons, Jeremy; Lester, David; Bird, Richard (2006), "Functional pearl: Enumerating the rationals", Journal of Functional Programming, 16 (3): 281–291, doi:10.1017/S0956796806005880.
Raney, George N. (1973), "On continued fractions and finite automata", Mathematische Annalen, 206: 265–283, doi:10.1007/BF01355980
Stern, Moritz A. (1858), "Ueber eine zahlentheoretische Funktion", Journal für die reine und angewandte Mathematik, 55: 193–220.

External links

Weisstein, Eric W. "Calkin–Wilf Tree". MathWorld.
Weisstein, Eric W. "Stern's Diatomic Series". MathWorld.
Bogomolny, Alexander, Fractions on a Binary Tree II, Cut-the-knot

[1] Berstel & de Luca (1997), Section 6.

[2] Raney (1973).

[3] The description here is dual to the original definition by Calkin and Wilf, which begins by defining the child relationship and derives the parent relationship as part of a proof that every rational appears once in the tree. As defined here, every rational appears once by definition, and instead the fact that the resulting structure is a tree requires a proof.

[FOOTNOTEGibbonsLesterBird2006-4] Gibbons, Lester & Bird (2006).

[5] Calkin & Wilf (2000): "a list of all positive rational numbers, each appearing once and only once, can be made by writing down 1/1, then the fractions on the level just below the top of the tree, reading from left to right, then the fractions on the next level down, reading from left to right, etc." Gibbons, Lester & Bird (2006) discuss efficient functional programming techniques for performing this breadth first traversal.

[6] Aigner & Ziegler (2004) credit this formula to Moshe Newman.

[7] The fusc name was given to it in 1976 by Edsger W. Dijkstra; see EWD570 and EWD578.

[8] Calkin & Wilf (2000), Theorem 1.

[9] Carlitz (1964).

[10] The OEIS entry credits this fact to Carlitz (1964) and to uncited work of Lind. However, Carlitz' paper describes a more restricted class of sums of powers of two, counted by fusc(n) instead of by fusc(n + 1).

[FOOTNOTEBatesBunderTognetti2010-11] Bates, Bunder & Tognetti (2010).

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]