Euler's totient function
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 φ(n) or ϕ(n), and may also be called Euler's phi function. It can be defined more formally as the number of integers k in the range 1 ≤ k ≤ n for which the greatest common divisor gcd(n, k) is equal to 1.[2][3] The integers k of this form are sometimes referred to as totatives of n.
For example, the totatives of n = 9 are the six numbers 1, 2, 4, 5, 7 and 8. They are all relatively prime to 9, but the other three numbers in this range, 3, 6, and 9 are not, because gcd(9, 3) = gcd(9, 6) = 3 and gcd(9, 9) = 9. Therefore, φ(9) = 6. As another example, φ(1) = 1 since for n = 1 the only integer in the range from 1 to n is 1 itself, and gcd(1, 1) = 1.
Euler's totient function is a multiplicative function, meaning that if two numbers m and n are relatively prime, then φ(mn) = φ(m) φ(n).[4][5] This function gives the order of the multiplicative group of integers modulo n (the group of units of the ring Z/nZ).[6] It also plays a key role in the definition of the RSA encryption system.
History, terminology, and notation
Leonhard Euler introduced the function in 1763.[7][8][9] However, he did not at that time choose any specific symbol to denote it. In a 1784 publication, Euler studied the function further, choosing the Greek letter π to denote it: he wrote πD for "the multitude of numbers less than D, and which have no common divisor with it".[10] This definition varies from the current definition for the totient function at D = 1 but is otherwise the same. The now-standard notation[8][11] φ(A) comes from Gauss's 1801 treatise Disquisitiones Arithmeticae.[12] (However, Gauss didn't use parentheses around the argument and wrote φA.) Thus, it is often called Euler's phi function or simply the phi function.
In 1879, J. J. Sylvester coined the term totient for this function,[13][14] so it is also referred to as Euler's totient function, the Euler totient, or Euler's totient. Jordan's totient is a generalization of Euler's.
The cototient of n is defined as n − φ(n). It counts the number of positive integers less than or equal to n that have at least one prime factor in common with n.
Computing Euler's totient function
There are several formulas for computing φ(n).
Euler's product formula
It states
where the product is over the distinct prime numbers dividing n. (The notation is described in the article Arithmetical function.)
The proof of Euler's product formula depends on two important facts.
The function φ(n) is multiplicative
This means that if gcd(m, n) = 1, then φ(mn) = φ(m) φ(n). (Sketch of proof: let A, B, C be the sets of residue classes modulo-and-coprime to m, n, mn respectively; then there is a bijection between A × B and C, by the Chinese remainder theorem.)
φ(pk) = pk − pk−1 = pk−1(p − 1)
That is, if p is prime and k ≥ 1 then
Proof: since p is a prime number the only possible values of gcd(pk, m) are 1, p, p2, ..., pk, and the only way for gcd(pk, m) to not equal 1 is for m to be a multiple of p. The multiples of p that are less than or equal to pk are p, 2p, 3p, ..., pk − 1p = pk, and there are pk − 1 of them. Therefore, the other pk − pk − 1 numbers are all relatively prime to pk.
Proof of Euler's product formula
The fundamental theorem of arithmetic states that if n > 1 there is a unique expression for n,
where p1 < p2 < ... < pr are prime numbers and each ki ≥ 1. (The case n = 1 corresponds to the empty product.)
Repeatedly using the multiplicative property of φ and the formula for φ(pk) gives
This is Euler's product formula.
Example
In words, this says that the distinct prime factors of 36 are 2 and 3; half of the thirty-six integers from 1 to 36 are divisible by 2, leaving eighteen; a third of those are divisible by 3, leaving twelve numbers that are coprime to 36. And indeed there are twelve positive integers that are coprime with 36 and lower than 36: 1, 5, 7, 11, 13, 17, 19, 23, 25, 29, 31, and 35.
Fourier transform
The totient is the discrete Fourier transform of the gcd, evaluated at 1:[15] Let
where for Then
The real part of this formula is
Note that unlike the other two formulae (the Euler product and the divisor sum) this one does not require knowing the factors of n. However, it does involve the calculation of the greatest common divisor of n and every positive integer less than n, which suffices to provide the factorization anyway.
Divisor sum
The property established by Gauss,[16] that
where the sum is over all positive divisors d of n, can be proven in several ways. (see Arithmetical function for notational conventions.)
One way is to note that φ(d) is also equal to the number of possible generators of the cyclic group Cd; specifically, if Cd = <g>, then gk is a generator for every k coprime to d. Since every element of Cn generates a cyclic subgroup, and all subgroups of Cd ≤ Cn are generated by some element of Cn, the formula follows.[17] In the article Root of unity Euler's formula is derived by using this argument in the special case of the multiplicative group of the nth roots of unity.
This formula can also be derived in a more concrete manner.[18] Let n = 20 and consider the fractions between 0 and 1 with denominator 20:
Put them into lowest terms:
First note that all the divisors of 20 are denominators. And second, note that there are 20 fractions. Which fractions have 20 as denominator? The ones whose numerators are relatively prime to 20 By definition this is φ(20) fractions. Similarly, there are φ(10) = 4 fractions with denominator 10 φ(5) = 4 fractions with denominator 5 and so on.
In detail, we're considering the fractions of the form k/n where k is an integer from 1 to n inclusive. Upon reducing these to lowest terms, each fraction will have as its denominator some divisor of n. We can group the fractions together by denominator, and we must show that for a given divisor d of n, the number of such fractions with denominator d is φ(d).
Note that to reduce k/n to lowest terms, we divide the numerator and denominator by gcd(k, n). The reduced fractions with denominator d are therefore precisely the ones originally of the form k/n in which gcd(k, n)=n/d. The question therefore becomes: how many k are there less than or equal to n which verify gcd(k, n)=n/d? Any such k must clearly be a multiple of n/d, but it must also be coprime to d (if it had any common divisor s with d, then sn/d would be a larger common divisor of n and k). Conversely, any multiple k of n/d which is coprime to d will satisfy gcd(n, k)=n/d. We can generate φ(d) such numbers by taking the numbers less than d coprime to d and multiplying each one by n/d (these products will of course each be smaller than n, as required). This in fact generates all such numbers, as if k is a multiple of n/d coprime to d (and less than n), then k/(n/d) will still be coprime to d, and must also be smaller than d, else k would be larger than n. Thus there are precisely φ(d) values of k less than or equal to n such that gcd(k, n)=n/d, which was to be demonstrated.
Möbius inversion gives
where μ is the Möbius function.
This formula may also be derived from the product formula by multiplying out to get
Riemann zeta function limit
For the Euler totient function can be calculated as a limit involving the Riemann zeta function:
where
is the Riemann zeta function, is the Möbius function, is e (mathematical constant), and is a divisor.
Some values of the function
The first 99 values (sequence A000010 in the OEIS) are shown in the table and graph below:[19]
+0 | +1 | +2 | +3 | +4 | +5 | +6 | +7 | +8 | +9 | |
---|---|---|---|---|---|---|---|---|---|---|
0+ | 1 | 1 | 2 | 2 | 4 | 2 | 6 | 4 | 6 | |
10+ | 4 | 10 | 4 | 12 | 6 | 8 | 8 | 16 | 6 | 18 |
20+ | 8 | 12 | 10 | 22 | 8 | 20 | 12 | 18 | 12 | 28 |
30+ | 8 | 30 | 16 | 20 | 16 | 24 | 12 | 36 | 18 | 24 |
40+ | 16 | 40 | 12 | 42 | 20 | 24 | 22 | 46 | 16 | 42 |
50+ | 20 | 32 | 24 | 52 | 18 | 40 | 24 | 36 | 28 | 58 |
60+ | 16 | 60 | 30 | 36 | 32 | 48 | 20 | 66 | 32 | 44 |
70+ | 24 | 70 | 24 | 72 | 36 | 40 | 36 | 60 | 24 | 78 |
80+ | 32 | 54 | 40 | 82 | 24 | 64 | 42 | 56 | 40 | 88 |
90+ | 24 | 72 | 44 | 60 | 46 | 72 | 32 | 96 | 42 | 60 |
The top line in the graph, y = n − 1, is a true upper bound. It is attained whenever n is prime. There is no lower bound that is a straight line of positive slope; no matter how gentle the slope of a line is, there will eventually be points of the plot below the line. More precisely, the lower limit of the graph is proportional to n/log log n rather than being linear.[20]
Euler's theorem
This states that if a and n are relatively prime then
The special case where n is prime is known as Fermat's little theorem.
This follows from Lagrange's theorem and the fact that φ(n) is the order of the multiplicative group of integers modulo n.
The RSA cryptosystem is based on this theorem: it implies that the inverse of the function , where is the (public) encryption exponent, is the function , where , the (private) decryption exponent, is the multiplicative inverse of modulo . The difficulty of computing without knowing the factorization of is thus the difficulty of computing : this is known as the RSA problem which can be solved by factoring . The owner of the private key knows the factorization, since an RSA private key is constructed by choosing as the product of two (randomly chosen) large primes and . Only is publicly disclosed, and given the difficulty to factor large numbers we have the guarantee that no-one else knows the factorization.
Other formulae
- implies
- (a, n > 1)
- where d = gcd(m, n). Note the special cases
and
- Compare this to the formula (See lcm).
- is even for Moreover, if n has r distinct odd prime factors,
- For any a > 1 and n > 6 such that there exists an such that .
- where is the radical of n.
- [21]
- ([22] cited in [23])
(here γ is the Euler–Mascheroni constant).
where m > 1 is a positive integer and ω(m) is the number of distinct prime factors of m. (a, b) is a standard abbreviation for gcd(a, b).[26]
Menon's identity
In 1965 P. Kesava Menon proved
where d(n) = σ0(n) is the number of divisors of n.
Formulae involving the golden ratio
Schneider[27] found a pair of identities connecting the totient function, the golden ratio and the Möbius function . In this section is the totient function, and is the golden ratio.
They are:
and
Subtracting them gives
Applying the exponential function to both sides of the preceding identity yields an infinite product formula for Euler's number e
The proof is based on the formulae
- and valid for 0 < x < 1.
Generating functions
The Dirichlet series for φ(n) may be written in terms of the Riemann zeta function as:[28]
The Lambert series generating function is[29]
which converges for |q| < 1.
Both of these are proved by elementary series manipulations and the formulae for φ(n).
Growth rate
In the words of Hardy & Wright, the order of φ(n) is "always ‘nearly n’."[30]
First[31]
but as n goes to infinity,[32] for all δ > 0
These two formulae can be proved by using little more than the formulae for φ(n) and the divisor sum function σ(n).
In fact, during the proof of the second formula, the inequality
true for n > 1, is proved.
We also have[20]
Here γ is Euler's constant, γ = 0.577215665..., eγ = 1.7810724..., e−γ = 0.56145948... .
Proving this doesn't quite require the prime number theorem.[33][34] Since log log (n) goes to infinity, this formula shows that
In fact, more is true.[35][36][37]
- for n > 2, and
- for infinitely many n.
The second inequality was shown by Jean-Louis Nicolas. Ribenboim says "The method of proof is interesting, in that the inequality is shown first under the assumption that the Riemann hypothesis is true, secondly under the contrary assumption."[37]
For the average order, we have[22][38]
due to Arnold Walfisz, its proof exploiting estimates on exponential sums due to I. M. Vinogradov and N.M. Korobov (this is currently the best known estimate of this type). The "Big O" stands for a quantity that is bounded by a constant times the function of "n" inside the parentheses (which is small compared to n2).
This result can be used to prove[39] that the probability of two randomly chosen numbers being relatively prime is
Ratio of consecutive values
In 1950 Somayajulu proved[40][41]
- and
In 1954 Schinzel and Sierpiński strengthened this, proving[40][41] that the set
is dense in the positive real numbers. They also proved[40] that the set
is dense in the interval (0, 1).
Totient numbers
A totient number is a value of Euler's totient function: that is, an m for which there is at least one n for which φ(n) = m. The valency or multiplicity of a totient number m is the number of solutions to this equation.[42] A nontotient is a natural number which is not a totient number. Every odd integer exceeding 1 is trivially a nontotient. There are also infinitely many even nontotients,[43] and indeed every positive integer has a multiple which is an even nontotient.[44]
The number of totient numbers up to a given limit x is
for a constant C = 0.8178146... .[45]
If counted accordingly to multiplicity, the number of totient numbers up to a given limit x is
where the error term R is of order at most for any positive k.[46]
It is known that the multiplicity of m exceeds mδ infinitely often for any δ < 0.55655.[47][48]
Ford's theorem
Ford (1999) proved that for every integer k ≥ 2 there is a totient number m of multiplicity k: that is, for which the equation φ(n) = m has exactly k solutions; this result had previously been conjectured by Wacław Sierpiński,[49] and it had been obtained as a consequence of Schinzel's hypothesis H.[45] Indeed, each multiplicity that occurs, does so infinitely often.[45][48]
However, no number m is known with multiplicity k = 1. Carmichael's totient function conjecture is the statement that there is no such m.[50]
Applications
Cyclotomy
In the last section of the Disquisitiones[51][52] Gauss proves[53] that a regular n-gon can be constructed with straightedge and compass if φ(n) is a power of 2. If n is a power of an odd prime number the formula for the totient says its totient can be a power of two only if a) n is a first power and b) n − 1 is a power of 2. The primes that are one more than a power of 2 are called Fermat primes, and only five are known: 3, 5, 17, 257, and 65537. Fermat and Gauss knew of these. Nobody has been able to prove whether there are any more.
Thus, a regular n-gon has a straightedge-and-compass construction if n is a product of distinct Fermat primes and any power of 2. The first few such n are[54] 2, 3, 4, 5, 6, 8, 10, 12, 15, 16, 17, 20, 24, 30, 32, 34, 40, ... . (sequence A003401 in the OEIS)
The RSA cryptosystem
Setting up an RSA system involves choosing large prime numbers p and q, computing n = pq and k = φ(n), and finding two numbers e and d such that ed ≡ 1 (mod k). The numbers n and e (the "encryption key") are released to the public, and d (the "decryption key") is kept private.
A message, represented by an integer m, where 0 < m < n, is encrypted by computing S = me (mod n).
It is decrypted by computing t = Sd (mod n). Euler's Theorem can be used to show that if 0 < t < n, then t = m.
The security of an RSA system would be compromised if the number n could be factored or if φ(n) could be computed without factoring n.
Unsolved problems
Lehmer's conjecture
If p is prime, then φ(p) = p − 1. In 1932 D. H. Lehmer asked if there are any composite numbers n such that φ(n) | n − 1. None are known.[55]
In 1933 he proved that if any such n exists, it must be odd, square-free, and divisible by at least seven primes (i.e. ω(n) ≥ 7). In 1980 Cohen and Hagis proved that n > 1020 and that ω(n) ≥ 14.[56] Further, Hagis showed that if 3 divides n then n > 101937042 and ω(n) ≥ 298848.[57][58]
Carmichael's conjecture
This states that there is no number n with the property that for all other numbers m, m ≠ n, φ(m) ≠ φ(n). See Ford's theorem above.
As stated in the main article, if there is a single counterexample to this conjecture, there must be infinitely many counterexamples, and the smallest one has at least ten billion digits in base 10.[42]
See also
- Carmichael function
- Duffin–Schaeffer conjecture
- Generalizations of Fermat's little theorem
- Highly composite number
- Multiplicative group of integers modulo n
- Ramanujan sum
Notes
- ↑ "Euler's totient function". Khan Academy. Retrieved 2016-02-26.
- ↑ Long (1972, p. 85)
- ↑ Pettofrezzo & Byrkit (1970, p. 72)
- ↑ Long (1972, p. 162)
- ↑ Pettofrezzo & Byrkit (1970, p. 80)
- ↑ See Euler's theorem.
- ↑ L. Euler "Theoremata arithmetica nova methodo demonstrata" (An arithmetic theorem proved by a new method), Novi commentarii academiae scientiarum imperialis Petropolitanae (New Memoirs of the Saint-Petersburg Imperial Academy of Sciences), 8 (1763), 74-104. (The work was presented at the Saint-Petersburg Academy on October 15, 1759. A work with the same title was presented at the Berlin Academy on June 8, 1758). Available on-line in: Ferdinand Rudio, ed., Leonhardi Euleri Commentationes Arithmeticae, volume 1, in: Leonhardi Euleri Opera Omnia, series 1, volume 2 (Leipzig, Germany, B. G. Teubner, 1915), pages 531-555. On page 531, Euler defines n as the number of integers that are smaller than N and relatively prime to N (… aequalis sit multitudini numerorum ipso N minorum, qui simul ad eum sint primi, …), which is the phi function, φ(N).
- 1 2 Sandifer, p. 203
- ↑ Graham et al. p. 133 note 111
- ↑ L. Euler, Speculationes circa quasdam insignes proprietates numerorum, Acta Academiae Scientarum Imperialis Petropolitinae, vol. 4, (1784), pp. 18-30, or Opera Omnia, Series 1, volume 4, pp. 105-115. (The work was presented at the Saint-Petersburg Academy on October 9, 1775).
- ↑ Both and are seen in the literature. These are two forms of the lower-case Greek letter phi.
- ↑ Gauss, Disquisitiones Arithmeticae article 38
- ↑ J. J. Sylvester (1879) "On certain ternary cubic-form equations", American Journal of Mathematics, 2 : 357-393; Sylvester coins the term "totient" on page 361.
- ↑ "totient". Oxford English Dictionary (2nd ed.). Oxford University Press. 1989.
- ↑ Schramm (2008)
- ↑ Gauss, DA, art 39
- ↑ Gauss, DA art. 39, arts. 52-54
- ↑ Graham et al. pp. 134-135
- ↑ The cell for n = 0 in the upper-left corner of the table is empty, as the function φ(n) is commonly defined only for positive integers, so it is not defined for n = 0.
- 1 2 Hardy & Wright 1979, thm. 328
- ↑ Dineva (in external refs), prop. 1
- 1 2 3 Walfisz, Arnold (1963). Weylsche Exponentialsummen in der neueren Zahlentheorie. Mathematische Forschungsberichte (in German). 16. Berlin: VEB Deutscher Verlag der Wissenschaften. Zbl 0146.06003.
- ↑ Lomadse, G., "The scientific work of Arnold Walfisz" (PDF), Acta Arithmetica, 10 (3): 227–237
- ↑ R. Sitaramachandrarao. On an error term of Landau II, Rocky Mountain J. Math. 15 (1985), 579-588
- ↑ Also R. Sitaramachandrarao (loc. cit.)
- ↑ Bordellès in the external links
- ↑ All formulae in the section are from Schneider (in the external links)
- ↑ Hardy & Wright 1979, thm. 288
- ↑ Hardy & Wright 1979, thm. 309
- ↑ Hardy & Wright 1979, intro to § 18.4
- ↑ Hardy & Wright 1979, thm. 326
- ↑ Hardy & Wright 1979, thm. 327
- ↑ In fact Chebychev's theorem (Hardy & Wright 1979, thm.7) and Mertens' third theorem is all that's needed
- ↑ Hardy & Wright 1979, thm. 436
- ↑ Theorem 15 of Rosser, J. Barkley; Schoenfeld, Lowell. Approximate formulas for some functions of prime numbers. Illinois J. Math. 6 (1962), no. 1, 64–94. http://projecteuclid.org/euclid.ijm/1255631807.
- ↑ Bach & Shallit, thm. 8.8.7
- 1 2 Ribenboim, The Book of Prime Number Records, Section 4.I.C
- ↑ Sándor, Mitrinović & Crstici (2006) pp.24-25
- ↑ Hardy & Wright 1979, thm. 332
- 1 2 3 Ribenboim, p.38
- 1 2 Sándor, Mitrinović & Crstici (2006) p.16
- 1 2 Guy (2004) p.144
- ↑ Sándor & Crstici (2004) p.230
- ↑ Zhang, Mingzhi (1993). "On nontotients". Journal of Number Theory. 43 (2): 168–172. doi:10.1006/jnth.1993.1014. ISSN 0022-314X. Zbl 0772.11001.
- 1 2 3 Ford, Kevin (1998). "The distribution of totients". Ramanujan J. 2 (1-2): 67–151. doi:10.1007/978-1-4757-4507-8_8. ISSN 1382-4090. Zbl 0914.11053.
- ↑ Sándor et al (2006) p.22
- ↑ Sándor et al (2006) p.21
- 1 2 Guy (2004) p.145
- ↑ Sándor & Crstici (2004) p.229
- ↑ Sándor & Crstici (2004) p.228
- ↑ Gauss, DA. The 7th § is arts. 336-366
- ↑ Gauss proved if n satisfies certain conditions then the n-gon can be constructed. In 1837 Pierre Wantzel proved the converse, if the n-gon is constructible, then n must satisfy Gauss's conditions
- ↑ Gauss, DA, art 366
- ↑ Gauss, DA, art. 366. This list is the last sentence in the Disquisitiones
- ↑ Ribenboim, pp. 36-37.
- ↑ Cohen, Graeme L.; Hagis, Peter, jun. (1980). "On the number of prime factors of n if φ(n) divides n−1". Nieuw Arch. Wiskd., III. Ser. 28: 177–185. ISSN 0028-9825. Zbl 0436.10002.
- ↑ Hagis, Peter, jun. (1988). "On the equation M⋅φ(n)=n−1". Nieuw Arch. Wiskd., IV. Ser. 6 (3): 255–261. ISSN 0028-9825. Zbl 0668.10006.
- ↑ Guy (2004) p.142
References
The Disquisitiones Arithmeticae has been translated from Latin into English and German. The German edition includes all of Gauss' papers on number theory: all the proofs of quadratic reciprocity, the determination of the sign of the Gauss sum, the investigations into biquadratic reciprocity, and unpublished notes.
References to the Disquisitiones are of the form Gauss, DA, art. nnn.
- Abramowitz, M.; Stegun, I. A. (1964), Handbook of Mathematical Functions, New York: Dover Publications, ISBN 0-486-61272-4. See paragraph 24.3.2.
- Bach, Eric; Shallit, Jeffrey (1996), Algorithmic Number Theory (Vol I: Efficient Algorithms), MIT Press Series in the Foundations of Computing, Cambridge, MA: The MIT Press, ISBN 0-262-02405-5, Zbl 0873.11070
- Ford, Kevin (1999), "The number of solutions of φ(x) = m", Annals of Mathematics, 150 (1): 283–311, doi:10.2307/121103, ISSN 0003-486X, JSTOR 121103, MR 1715326, Zbl 0978.11053.
- Gauss, Carl Friedrich; Clarke, Arthur A. (translator into English) (1986), Disquisitiones Arithmeticae (Second, corrected edition), New York: Springer, ISBN 0-387-96254-9
- Gauss, Carl Friedrich; Maser, H. (translator into German) (1965), Untersuchungen uber hohere Arithmetik (Disquisitiones Arithmeticae & other papers on number theory) (Second edition), New York: Chelsea, ISBN 0-8284-0191-8
- Graham, Ronald; Knuth, Donald; Patashnik, Oren (1994), Concrete Mathematics: a foundation for computer science (2nd ed.), Reading, MA: Addison-Wesley, ISBN 0-201-55802-5, Zbl 0836.00001
- Guy, Richard K. (2004), Unsolved Problems in Number Theory, Problem Books in Mathematics (3rd ed.), New York, NY: Springer-Verlag, ISBN 0-387-20860-7, Zbl 1058.11001
- Hardy, G. H.; Wright, E. M. (1979), An Introduction to the Theory of Numbers (Fifth ed.), Oxford: Oxford University Press, ISBN 978-0-19-853171-5
- Long, Calvin T. (1972), Elementary Introduction to Number Theory (2nd ed.), Lexington: D. C. Heath and Company, LCCN 77-171950
- Pettofrezzo, Anthony J.; Byrkit, Donald R. (1970), Elements of Number Theory, Englewood Cliffs: Prentice Hall, LCCN 77-81766
- Ribenboim, Paulo (1996), The New Book of Prime Number Records (3rd ed.), New York: Springer, ISBN 0-387-94457-5, Zbl 0856.11001
- Sandifer, Charles (2007), The early mathematics of Leonhard Euler, MAA, ISBN 0-88385-559-3
- Sándor, József; Mitrinović, Dragoslav S.; Crstici, Borislav, eds. (2006), Handbook of number theory I, Dordrecht: Springer-Verlag, pp. 9–36, ISBN 1-4020-4215-9, Zbl 1151.11300
- Sándor, Jozsef; Crstici, Borislav (2004). Handbook of number theory II. Dordrecht: Kluwer Academic. pp. 179–327. ISBN 1-4020-2546-7. Zbl 1079.11001.
- Schramm, Wolfgang (2008), "The Fourier transform of functions of the greatest common divisor", Electronic Journal of Combinatorial Number Theory, A50 (8(1)).
External links
- Hazewinkel, Michiel, ed. (2001), "Totient function", Encyclopedia of Mathematics, Springer, ISBN 978-1-55608-010-4
- Kirby Urner, Computing totient function in Python and scheme, (2003)
- Euler's Phi Function and the Chinese Remainder Theorem — proof that Φ(n) is multiplicative
- Euler's totient function calculator in JavaScript — up to 20 digits
- Bordellès, Olivier, Numbers prime to q in
- Dineva, Rosica, The Euler Totient, the Möbius, and the Divisor Functions
- Miyata, Daisuke & Yamashita, Michinori, Derived logarithmic function of Euler's function
- Plytage, Loomis, Polhill Summing Up The Euler Phi Function
- Schneider, Robert P. A Golden Product Identity for e.