
TL;DR
This paper analyzes the cycle lengths of a linear map on modular integer vectors, revealing how these periods depend on the prime factorization of the modulus and characterizing cycles for prime moduli.
Contribution
It provides a formula for the period of the map based on prime factorization and characterizes cycles when the modulus is prime.
Findings
Period depends on prime factorization of m
Cycle characterization for prime m
Main theorem relates periods to prime factors
Abstract
We study the period of the linear map as a function of and , where stands for the ring of integers modulo . Since this map is a variant of the Ducci sequence, several known results are adapted in the context of . The main theorem of this paper states that the period modulo can be deduced from the prime factorization of and the periods of its prime factors. We also characterize the tuples that belong to a cycle when is prime.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Cycles of Sums of Integers
Bruno Dular
Abstract.
We study the period of the linear map as a function of and , where stands for the ring of integers modulo . Since this map is a variant of the Ducci sequence, several known results are adapted in the context of . The main theorem of this paper states that the period modulo can be deduced from the prime factorization of and the periods of its prime factors. We also characterize the tuples that belong to a cycle when is prime.
Contents
- 1 Introduction
- 2 Basic results
- 3 Multiplicity of the period function
- 4 Order of 2 and Wieferich primes
- 5 Period modulo powers of primes
- 6 Characterization of tuples in a cycle
- 7 Open questions
1. Introduction
The aim of this paper111Published in: The Fibonacci Quarterly, Volume 58 No 2 (2020) 126. is to study a variant of the well-known Ducci game of differences. In this game, one starts with a -tuple of integers and iterates the Ducci map to generate a Ducci sequence. This process suggests the name Cycles of differences of integers [10], which inspired the name of the present paper.
If is a power of , we know that every Ducci sequence eventually vanishes, i.e., reaches the zero -tuple. Else, a Ducci sequence will either vanish or enter a periodic cycle. As several authors pointed out [10, 5], studying the latter case comes down to considering n-tuples that consist only of [math] and ’s. Hence the Ducci map can be considered to be a linear map over . This map can be generalized as stated in Definition 1.1, performing sums modulo for some positive integer .
This variant has been introduced by Wong in [14] and has been extensively studied by F. Breuer in [3], who noticed a link between Ducci sequences and cyclotomic polynomials. The results we prove here are similar, but we use an elementary approach, which allows us to solve the inseparable case, i.e., with a prime and divisible by . This is the object of section 5. There we also see that our method does not generalize to certain special cases ( or a Wieferich prime), for which F. Breuer’s method works.
Here denotes the set of positive integers (with the convention ).
Definition 1.1**.**
Let . Let
[TABLE]
A T-sequence of is a sequence of the form where and it is said to be generated by the tuple a.
The tuple and the T-sequence it generates are respectively called the basic tuple and the basic T-sequence of .
For example, the T-sequence generated by the basic tuple of starts as shown below. Note that the tuple repeats, hence the T-sequence becomes periodic at that point, with a cycle length of . With the notations introduced later, we write .
[TABLE]
Remark 1.2**.**
Let a be a tuple of where for some integers and . We consider a as an element of by identifying it with the element of .
As , note that T-sequences are Ducci sequences when . Several known results can then be generalized.
To simplify notation, the components of a tuple are indexed from [math] to . We sometimes write for . Note that addition and subtraction of the indices will always be performed modulo . Thus, should be understood as .
Since is finite, a T-sequence must be eventually periodic. The goal of this paper is to study the maximal cycle length as a function of and , which we denote by . We shall detail what we mean by the length of the period.
Definition 1.3**.**
Given , and , a positive integer is the cycle length of the T-sequence if the following conditions hold:
- (1)
There exists a positive integer such that for all . 2. (2)
Every positive integer satisfying for large enough is a multiple of .
The smallest such is called the pre-period. If is the pre-period, then the finite sequences and are called pre-cycle and cycle, respectively.
In other words, the cycle length of the T-sequence is the smallest positive integer such that there exists some satisfying for all .
We define as the subset of of all tuples that belong to a cycle. It directly follows from Remark 1.2 that whenever divides .
To simplify notation, cyclic permutations of a cycle are also called cycles and thus we refer to a cycle than the cycle.
In Section 6, we give a characterization of for some and .
Definition 1.4**.**
Given and in , we define
[TABLE]
It is called the period of T-sequences in . This defines a function called the period function.
As Ehrlich pointed out in [5], studying the cycle length of basic T-sequences suffices to determine the period of T-sequences. Actually, Ehrlich proved this result for Ducci sequences but the proof is essentially the same for T-sequences.
Proposition 1.5**.**
For all , the cycle length of the basic T-sequence of equals . Cycle lengths of other T-sequences in divide .
In Section 2 we give basic results that are useful to study more interesting properties of T-sequences. Among those we prove a generalization of the known fact that Ducci sequences of -tuples eventually vanish.
In Sections 4 and 5 we give important theorems about the multiplicity of the period function, which are summed up in the following theorem. It is the main result of this paper.
Theorem**.**
Let with the prime factorization of . If are odd and non-Wieferich222See definition 4.4., then
[TABLE]
Acknowledgements
I thank my friend Lucas Michel for telling me about the question of this paper.
2. Basic results
The first result below allows us to compute iterations of in a very simple way. It will be used extensively throughout the paper.
Proposition 2.1**.**
Let , with . For all and such that ,
[TABLE]
Proof.
We prove this by induction on . For , the result is obvious. Suppose it holds for . We show that it holds for , by using Pascal’s triangle formula and manipulating the sums as follows,
[TABLE]
which completes the proof. ∎
Definition 2.2**.**
We say that a T-sequence , where , vanishes if there exists a positive integer such that .
Recall that T-sequences are Ducci sequences if . Thus it is well known that every T-sequence of vanishes if and only if is a power of . It has first been proven by Ciamberlini and Marengoni in [4], and it has been reproven many times since then [2, 5]. Actually, this result still holds when is any power of . This has been proven by Wong in [14]. We give here a shorter proof using the notations we introduced and proposition 2.1. In [1], C. Avart shows a converse to this theorem for the base case , stating that the only tuples that vanish are the tuples obtained by concatenation of several copies of a tuple of length a power of 2.
Theorem 2.3**.**
If and for some positive integers , , then every T-sequence of vanishes, that is, . Reciprocally, if every T-sequence of vanishes, then is a power of .
Proof.
We first prove the case . Let . Since is even for and by Proposition 2.1, we have
[TABLE]
for all between [math] and .
We now proceed with a proof by induction. Assume the result holds for some integer . We prove that it holds for . Let . If considered over , the T-sequence generated by a vanishes. Let be an integer such that . Therefore for some tuple u, which we can assume to consist only of [math] and ’s. It follows from the base case and by linearity of that there exists some integer such that , which concludes the proof.
The reciprocal follows from the proof of Proposition 4.3. ∎
We denote by the left-shift map [5], defined as
[TABLE]
Thus where is the identity map. The following Lemma generalizes lemma in [5].
Lemma 2.4**.**
If is a prime and is a positive integer, then as linear maps from into itself.
Proof.
We have where is the identity map. The proof follows from Proposition 2.1 and the fact that a prime divides for . ∎
The next proposition is a generalization of Corollary 3 and Theorem 2 in [5].
Proposition 2.5**.**
Let be a prime and positive integers.
- (1)
If , then divides . 2. (2)
If , then divides .
Proof.
By Lemma 2.4, we have:
- (1)
. 2. (2)
, hence .
∎
If and are coprime, then , the order of in , always satisfies (1).
3. Multiplicity of the period function
In the following sections, we focus on the main question of this paper: can we deduce from the prime factorization of , knowing for ?
The next proposition suggests a positive answer. We will often use it without reference.
Proposition 3.1**.**
If , then for all .
Proof.
Let be large enough so that . The congruence still holds modulo since . The conclusion follows from Definition 1.3. ∎
The goal of the next few sections is to study this relation more precisely.
Theorem 3.2**.**
If with pairwise coprime, then
[TABLE]
Proof.
We prove this for two coprime integers. The generalization for pairwise coprime integers follows by induction. By Proposition 1.5, we only need to consider the basic T-sequence.
Let be two coprime integers. We assume here that is large enough for to be in the different cycles. By Definition 1.3, we have
[TABLE]
where . Since and are coprime, it directly follows333If and , then for some integers , . Hence divides , so divides by Euclid’s lemma and we get for some integer . that . Hence is a multiple of the period .
We now show that satisfies of Definition 1.3. Suppose is such that . In particular, ; hence is a multiple of . Similarly, is a multiple of . Therefore, by definition of the least common multiple, , and . ∎
With this theorem in our toolbox, we can restrict our attention to the periods modulo powers of primes. The question one may ask is whether we can deduce from . We will shortly determine that it is (almost) the case.
4. Order of 2 and Wieferich primes
Definition 4.1**.**
For a tuple , we write for the sum of components of a modulo .
Definition 4.2**.**
For odd, we use to denote the order of in ; that is, it is the smallest integer such that . Its existence follows from Euler’s theorem.
Proposition 4.3**.**
If is odd, then divides .
Proof.
Let be a positive integer and a be a tuple of . By linearity, we have . If is the cycle length of the T-sequence generated by a and is greater than the pre-period, then we must have . Hence or . Note that if is not a power of , then for all . Considering the basic T-sequence thus implies that must equal . ∎
Definition 4.4**.**
A prime is a Wieferich prime if .
Wieferich primes surprisingly occur in several number theoretical subjects [8]. It is believed that there are infinitely many such numbers. What is extraordinary about these is that we only know two of them, and , and there are no other Wieferich primes below [13]. We will see in Section 5 that Wieferich primes are of considerable interest here.
We can characterize Wieferich primes by the order of modulo .
Lemma 4.5**.**
A prime is a Wieferich prime if and only if .
Proof.
We first show that equals either or . By definition, . It also holds modulo , so is a multiple of . We have
[TABLE]
so divides . Thus it equals either or .
Suppose that is a Wieferich prime, i.e., . Then divides , since we cannot have , we have . Conversely, if , then divides , so . ∎
As the following proposition shows, the first part of the previous proof also holds for all primes and positive integers , that is, is either or . In fact, as soon as it is the latter for one , it is the latter for all subsequent .
Proposition 4.6**.**
If is a prime and , then we have:
- (1)
* is either or .* 2. (2)
If , then . 3. (3)
If is a non-Wieferich prime, then .
Proof.
The proof of is similar to the first part of the proof of Lemma 4.5 and (3) follows from (1) and (2) by induction. We show (2):
Suppose . Then where . Hence for some . By the binomial theorem and the fact that divides for , we have . Since , we have , so must equal . This concludes the proof. ∎
For the known Wieferich primes, we have . Hence by (3) of the previous proposition, it follows that for all .
5. Period modulo powers of primes
Propositions 4.3 and 4.6 suggest a similar induction relation for the period function. Indeed, if is an odd prime and , then for all positive integers . The fact that eventually grows as a geometric sequence forces to behave in the same way.
Theorem 5.1**.**
If is a prime and , then we have:
- (1)
* is either or .* 2. (2)
If and , then . 3. (3)
If for some , then for all .
Proof.
In (1) and (2), we choose sufficiently large for to be in a cycle of and , respectively.
(1)
Let . Since a is in a cycle modulo , it is also in a cycle modulo and . Hence divides .
We have , so for some tuple u that we can consider to be in . By linearity of , the tuple is in a cycle modulo . It implies that u is in a cycle in , hence for some tuple v. Then, by linearity,
[TABLE]
By iterating on a, we then obtain . Hence divides .
Therefore, is either or .
(2)
Suppose and let , and . Suppose . By (1), we know that is either or . We show that it equals the latter.
Since a is in a cycle modulo , it is also in a cycle modulo and . Then , but this congruence does not hold modulo , for we assumed that . This implies that where .
By linearity of , the tuple is in a cycle modulo , hence u is in a cycle modulo . The condition implies that divides . We then have, by the same argument as in (1),
[TABLE]
Therefore, , so .
(3)
This follows directly from (1) and (2) by induction. ∎
The following proposition exhibits the fact that the condition is needed for (2) in the previous theorem to hold.
Proposition 5.2**.**
We have and for all positive integers .
Proof.
By Theorem 1.5, computing the first iterations of the basic T-sequence gives .
We prove that for all positive integers by induction. The idea of the argument is to show that, for all :
- (1)
The pre-period of the basic T-sequence of is and it has a cycle of the form
[TABLE]
where , with . 2. (2)
The pre-period of the basic T-sequence of is and it has a cycle of the form
[TABLE]
where and (Note that and are those from (1)). In other words, it means that the cycle in starts at the second tuple of the cycle in .
For the base cases, and , the cycles of the basic T-sequences are respectively
[TABLE]
and
[TABLE]
so (1) and (2) are satisfied.
Now let and suppose (1) is satisfied. Since , we have for some . Then by linearity,
[TABLE]
hence since are smaller than . If we let and , we have and . By symmetry of this argument under cyclic permutations, we have , , and , as desired. Then and (2) holds. ∎
The goal of the rest of this section is to show the main theorem of this paper, stated in the Introduction. To do that, we need to find base cases in order to use (3) of Theorem 5.1. We first prove the case in Proposition 5.3. Then we will use combinatorial congruences to deduce the case .
Proposition 5.3**.**
Let be a non-Wieferich prime and . If is not a multiple of , then and .
Proof.
By Theorem 5.1, we know that is either or . Given the fact that and are coprime, Proposition 2.5 tells us that divides , so it cannot be a multiple of . However, by Proposition 4.3, is a multiple of , which equals for is non-Wieferich (Proposition 4.6). Therefore, divides and we must have .
Similarly, is either or . Since divides and is not divisible by , we must have . ∎
Corollary 5.4**.**
If is a non-Wieferich prime and is not a multiple of , then for all .
Proof.
It follows from Proposition 5.3 and (3) in Theorem 5.1. ∎
To generalize Proposition 5.3 for any , we first need to prove a few lemmas.
The p-adic valuation of an integer is the exponent of the largest power of that divides . It is denoted by . We write the sum of the digits of when written in base . If is a prime, Legendre’s formula [9] states that
[TABLE]
Note that and for all integers .
The following Lemma generalizes the fact that divides for all but not for and . We will use the cases and later in Proposition 5.9.
Lemma 5.5**.**
Let be a prime and with . We have
[TABLE]
Proof.
It is clear that [math] and belong to both sets. Let and its decomposition in base . We show that where . Since , we have .
By Legendre’s formula, we have and where . We also have
[TABLE]
Thus,
[TABLE]
and
[TABLE]
Therefore, we obtain , hence an integer belongs to the first set if and only if , i.e., , which happens if and only if is a multiple of . This concludes the proof. ∎
The proof of the following lemma is due to Darij Grinberg and Victor Reiner ((12.69.3) in [6]).
Lemma 5.6**.**
Let and be a prime factor of . For all and such that is an integer, we have
[TABLE]
Proof.
The argument consists in counting the -elements subsets of the set . It is clearly .
At the same time, the subsets fall into two classes:
- (1)
The subsets which are invariant under the permutation of . 2. (2)
The other ones.
Say there are and subsets in the first and second class, respectively.
If a -elements subset belongs to the first class, then the intersection must have elements, which uniquely determine all of by iterating the permutation given above. Thus the first class contains elements.
Besides, the permutation of acts on the subsets of the second class, splitting them into orbits. Its -power acts trivially on the subsets of the second class. Then, the size of each orbit divides . Suppose that the size of an orbit is a proper divisor of . Then divides , so acts trivially on , hence elements of this orbit are subsets of the first class, a contradiction. Then every orbit has size .
Since the set of all second class subsets is the union of these orbit, it has size divisible by .
Therefore, we have and the proof is complete. ∎
To prove the following lemma, we shall introduce Babbage’s theorem (Theorem 1.12 in [7]). It states that for any prime and integers , we have . Note that if , we can replace by , thus strengthening the result. This case is known as Wolstenholme’s theorem.
Lemma 5.7**.**
If is a prime, and , then
[TABLE]
Proof.
It follows directly from Babbage’s theorem by induction. ∎
Lemma 5.8**.**
If is a prime, and , then
[TABLE]
Proof.
For , the proof follows from Wolstenholme’s theorem by induction.
Now we consider the cases and . If , it is direct. If , then Lemma 5.6 with , (then ) and gives
[TABLE]
and the congruence also holds modulo since . Then the proof proceeds by induction. ∎
Now we can apply these lemmas in order to prove the following proposition, which gives us base cases to apply Theorem 5.1.
Proposition 5.9**.**
Let be a prime and with , where . Then we have
- (1)
. If , it holds only if . 2. (2)
If is a non-Wieferich prime, then . 3. (3)
If is a non-Wieferich prime, then .
Proof.
The idea of this proof is to study the behavior of a T-sequence of (where , respectively) by studying the behavior of the T-sequence generated by a subtuple of a, which is a T-sequence of smaller tuples that we better understand.
For that purpose, we introduce a family of functions , , that extract an interesting subtuple from a given tuple. For , let
[TABLE]
(1)
Here we use . For , the subtuple is in .
By Lemma 2.4, we have , hence for any tuple . Considering the basic tuple e of , it gives where e’ is the basic tuple of . Note that components of e that are not components of e’ remain zero after any number of iterations of , hence the behavior of the T-sequence is entirely determined by the behavior of . Since the cycle length of the latter is , the cycle length of is .
If , note that this argument only holds if the T-sequences do not vanish. This explains the additional condition in this case.
(2)
Here we suppose is a non-Wieferich prime. We first show that and then that .
To prove the first part, we use . By Lemmas 5.5 and 5.7, we have
[TABLE]
which implies that . Since is an odd non-Wieferich prime, divides by Proposition 5.3. Thus we obtain .
We now use . For , the subtuple is in . We also have
[TABLE]
where the second and third equalities follow from Lemmas 5.5 and 5.7, respectively. Thus, . Therefore we obtain as desired, where the last equality follows from the first part.
(3)
Suppose is a non-Wieferich prime. To begin with, we suppose . We first show that and then that .
First, we use . For , the subtuple is in . By Lemmas 5.5 and 5.8, we have
[TABLE]
hence . By Proposition 5.3 (it is why we consider non-Wieferich), divides , thus we obtain .
Now we use . For , the subtuple is in . First, using Lemmas 5.5 and 5.8, we obtain
[TABLE]
hence . Together with the first part, since divides (by Proposition 5.3), we get as desired.
To conclude the proof, we consider the case . We have to show that . We use . For , the subtuple is in . As above, using Lemmas 5.5 and 5.8, we have , hence . By Proposition 5.3, divides . Thus we obtain , which concludes. ∎
Theorem 5.10**.**
If is a non-Wieferich prime, we have for all positive integers and .
Proof.
Let . Write where . We show that (1) and (2) .
(1)
Points (1) and (2) of Proposition 5.9 yield and , respectively. Since is non-Wieferich and coprime with , we have by Proposition 5.3. The conclusion follows directly.
(2)
The argument is the same, using (3) of Proposition 5.9 instead of (2).
Therefore, we complete the proof by Theorem 5.1. ∎
We are now able to prove the main result of this paper.
Theorem 5.11**.**
Let with the prime factorization of . If are odd and non-Wieferich, then
[TABLE]
Proof.
It follows from Theorems 3.2 and 5.10. ∎
At this point, a question one may ask is whether we can generalize this result for and Wieferich primes. The proofs above are considerably dependant on Proposition 5.3, which itself is dependant on Proposition 4.6. Therefore, it would not be possible to use the same method to obtain similar results for these special primes.
However, F. Breuer shows in Theorem 8.2 of [3] that a variant of Theorem 5.10 holds for and Wieferich primes, namely that for some integer . That is, these special cases eventually behave as we would expect. Whether it is possible or not to derive such results with an elementary method, similar to the ones used here, remains an open question.
6. Characterization of tuples in a cycle
In this section we try to characterize tuples of that belong to . By theorem 2.3, we already know that if and are powers of .
The linear map is represented in the standard basis by the matrix
[TABLE]
Note that when is even and when is odd. This simple fact yields the following proposition.
Proposition 6.1**.**
Let and be two odd integers. Then .
Proof.
Since is invertible444An integer is invertible modulo if and only if and are coprime. modulo , the matrix is invertible555A matrix is invertible if and only if its determinant is invertible., hence is bijective. Consequently, we can unambiguously move backward in the eventually periodic T-sequence determined by a tuple a of , so a belongs to a cycle. ∎
For even, things are a bit more complicated and require the introduction of a few new notations. We denote the alternating sum of components of a tuple by
[TABLE]
Proposition 6.2**.**
Let be odd and be even. If and are coprime, then a tuple belongs to a cycle if and only if .
Proof.
We first show that the condition is sufficient. Let be such that . Finding a preimage x of a is equivalent to solving the system
[TABLE]
All components of x are determined by the value chosen for and the last equation is satisfied since . Thus, a has exactly preimages. Moreover, and are coprime, so is invertible, hence the equation has exactly one solution . Thus, a has exactly one preimage x with .
Therefore, the map restricted to the set is a one-to-one correspondence and we can conclude as in Proposition 6.1.
The last equation of (6.1) shows that the condition is necessary. ∎
The following result, concerning , is shown in [10] and in [12] for the Ducci map, which is the same as our map in this case.
Proposition 6.3**.**
In , we have and every has exactly two preimages. For odd , a tuple belongs to a cycle if and only if .
If is even, the two preimages of a tuple are either both in or both in . Thus we have to find a way to characterize tuples of that have preimages in .
That has actually already been done by Ludington-Young in [10, 11], where the following definition and theorem 6.4 come from. We only consider here tuples of . A tuple a is even if . Suppose where is odd, we say a tuple is r-even if
[TABLE]
for . For example, if , then a is -even if
[TABLE]
Theorem 6.4**.**
Let with odd. A tuple of belongs to a cycle if and only if it is -even.
Proposition 6.3 turns out to be the special case of this theorem. Indeed, a tuple a is [math]-even if and only if .
We now generalize this characterization to odd primes.
Definition 6.5**.**
Let with . We say a tuple a of is even if . We note
[TABLE]
for . We say a is r-even if for all .
Theorem 6.6**.**
Let be a prime and even with . A tuple of belongs to a cycle if and only if it is r-even.
Proof.
The congruences below are all modulo .
First, note that is r-even if a is r-even. Indeed, if a is r-even, then for each . By Lemma 2.4, we have that is r-even. Thus, every tuple of a cycle must be r-even by linearity of r-evenness, hence the condition is necessary.
We now show that it is sufficient. Let be r-even. By (6.1) (it is here that we use the condition that is even) and since r-evenness implies evenness, the tuple a has preimages. Let be one of these, hence all preimages are given by for . To simplify notation, we write for .
Since ,
[TABLE]
and, similarly, , ,, are all and too666The negative sign comes from the equality .. Hence,
[TABLE]
Since is invertible modulo , the equation has exactly one solution , hence is the only r-even preimage of a.
Therefore, the map restricted to the set of r-even tuples of is bijective and the proof is complete. ∎
7. Open questions
We can see the iterations of the map on the set (which is bijective when restricted to ) as the action of the group on this set, splitting it into orbits. What are the possible sizes for these orbits? The tuple has an orbit of size 1, whereas the basic tuple, after enough iterations, generates an orbit of size . What are the values between and that are the size of some orbit?
Further questions arise naturally. What is the largest pre-period that can happen? How can the results of this paper be generalized to any linear map of ? Do there exist explicit formulas to find for every prime and positive integer ?
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] C. Avart, A characterization of converging Ducci sequences over Z 2 , The Fibonacci Quarterly, 49 (2011).
- 2[2] F. Breuer, Ducci sequences in higher dimensions , INTEGERS: Electronic journal of combinatorial number theory, 7.24 (2007).
- 3[3] F. Breuer, Ducci sequences and cyclotomic fields , Journal of Difference Equations and Applications, 16.7 (2010), 847–862.
- 4[4] C. Ciamberlini and A. Marengoni, Su una interessante curiosita‘ numerica , Periodiche di Matematiche, 17 (1937, 25–30.
- 5[5] A. Ehrlich, Periods in Ducci’s n-number game of differences , The Fibonacci Quaterly, 28.4 (1990), 302–305.
- 6[6] D. Grinberg and V. Reiner, Hopf Algebras in Combinatorics , https://arxiv.org/abs/1409.8356 v 5/ , Ancillary file, 2014.
- 7[7] D. Grinberg, The Lucas and Babbage congruences , http://www.cip.ifi.lmu.de/~grinberg/ , Accessed: 2019-02-09 , 2019.
- 8[8] G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers , 5th edition, Oxford University Press, 1979, 202.
