Euler's totient function
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. In other words, it is 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 ℤ/nℤ).[6] It is also used for defining 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] although 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.
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 is multiplicative
This means that if gcd(m, n) = 1, then φ(mn) = φ(m) φ(n). (Outline of proof: let A, B, C be the sets of nonnegative integers, which are, respectively, coprime to and less than m, n, and mn; then there is a bijection between A × B and C, by the Chinese remainder theorem.)
Value for a prime power argument
If p is prime and k ≥ 1, then
Proof: since p is a prime number the only possible values of gcd(p**k, m) are 1, p, p2, ..., p**k, and the only way for gcd(p**k, 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 p**k are p, 2p, 3p, ..., p**k − 1p = p**k, and there are p**k − 1 of them. Therefore, the other p**k − p**k − 1 numbers are all relatively prime to p**k.
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 < ... < p**r are prime numbers and each k**i ≥ 1. (The case n = 1 corresponds to the empty product.)
Repeatedly using the multiplicative property of φ and the formula for φ(p**k) 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 xk = gcd(k,n) for k ∈ {1, …, n}. Then
The real part of this formula is
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 C**d; specifically, if C**d = ⟨g⟩, then g**k is a generator for every k coprime to d. Since every element of C**n generates a cyclic subgroup, and all subgroups of C**d ≤ C**n are generated by some element of C**n, 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 (1/20, 3/20, 7/20, 9/20, 11/20, 13/20, 17/20, 19/20). By definition this is φ(20) fractions. Similarly, there are φ(10) = 4 fractions with denominator 10 (1/10, 3/10, 7/10, 9/10), φ(5) = 4 fractions with denominator 5 (1/5, 2/5, 3/5, 4/5), and so on.
In detail, we are 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(k, n) = 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
Some values of the function
Graph of the first 100 values
0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | |
---|---|---|---|---|---|---|---|---|---|---|---|---|
0 | N/A | 1 | 1 | 2 | 2 | 4 | 2 | 6 | 4 | 6 | 4 | 10 |
12 | 4 | 12 | 6 | 8 | 8 | 16 | 6 | 18 | 8 | 12 | 10 | 22 |
24 | 8 | 20 | 12 | 18 | 12 | 28 | 8 | 30 | 16 | 20 | 16 | 24 |
36 | 12 | 36 | 18 | 24 | 16 | 40 | 12 | 42 | 20 | 24 | 22 | 46 |
48 | 16 | 42 | 20 | 32 | 24 | 52 | 18 | 40 | 24 | 36 | 28 | 58 |
60 | 16 | 60 | 30 | 36 | 32 | 48 | 20 | 66 | 32 | 44 | 24 | 70 |
72 | 24 | 72 | 36 | 40 | 36 | 60 | 24 | 78 | 32 | 54 | 40 | 82 |
84 | 24 | 64 | 42 | 56 | 40 | 88 | 24 | 72 | 44 | 60 | 46 | 72 |
96 | 32 | 96 | 42 | 60 | 40 | 100 | 32 | 102 | 48 | 48 | 52 | 106 |
108 | 36 | 108 | 40 | 72 | 48 | 112 | 36 | 88 | 56 | 72 | 58 | 96 |
120 | 32 | 110 | 60 | 80 | 60 | 100 | 36 | 126 | 64 | 84 | 48 | 130 |
132 | 40 | 108 | 66 | 72 | 64 | 136 | 44 | 138 | 48 | 92 | 70 | 120 |
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]
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 a ↦ ae mod n, where e is the (public) encryption exponent, is the function b ↦ bd mod n, where d, the (private) decryption exponent, is the multiplicative inverse of e modulo φ(n). The difficulty of computing φ(n) without knowing the factorization of n is thus the difficulty of computing d: this is known as the RSA problem which can be solved by factoring n. The owner of the private key knows the factorization, since an RSA private key is constructed by choosing n as the product of two (randomly chosen) large primes p and q. Only n is publicly disclosed, and given the difficulty to factor large numbers we have the guarantee that no-one else knows the factorization.
Other formulae
- Note the special cases
- Compare this to the formula
φ(n) is even for n ≥ 3. Moreover, if n has r distinct odd prime factors, 2r | φ(n)
For any a > 1 and n > 6 such that 4 ∤ n there exists an l ≥ 2n such that l | φ(an − 1).
- whererad(n)is theradical ofn.
- (whereγis theEuler–Mascheroni constant).
- wherem > 1is a positive integer andω(m)is the number of distinct prime factors ofm.[25]
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[26] found a pair of identities connecting the totient function, the golden ratio and the Möbius function μ(n). In this section φ(n) is the totient function, and ϕ = 1 + √5/2 = 1.618… 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 e:
The proof is based on the two formulae
Generating functions
The Dirichlet series for φ(n) may be written in terms of the Riemann zeta function as:[27]
The Lambert series generating function is[28]
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’.”[29]
First[30]
but as n goes to infinity,[31] 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..., so eγ = 1.7810724... and e−γ = 0.56145948....
and
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."[36]
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[38] that the probability of two randomly chosen numbers being relatively prime is 6/π2.
Ratio of consecutive values
is dense in the positive real numbers. They also proved[39] 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.[41] 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,[42] and indeed every positive integer has a multiple which is an even nontotient.[43]
The number of totient numbers up to a given limit x is
for a constant C = 0.8178146....[44]
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 x/(log x)k for any positive k.[45]
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,[48] and it had been obtained as a consequence of Schinzel's hypothesis H.[44] Indeed, each multiplicity that occurs, does so infinitely often.[44][47]
However, no number m is known with multiplicity k = 1. Carmichael's totient function conjecture is the statement that there is no such m.[49]
Perfect totient numbers
Applications
Cyclotomy
In the last section of the Disquisitiones[50][51] Gauss proves[52] 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 n is a first power and 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[53]
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 = m**e (mod n).
It is decrypted by computing t = S**d (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.[54]
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.[41]
See also
Carmichael function
Duffin–Schaeffer conjecture
Generalizations of Fermat's little theorem
Highly composite number
Multiplicative group of integers modulo n
Ramanujan sum
Totient summatory function