Euler's totient function

thumb|The first thousand values of . The points on the top line represent when is a prime number, which is

In number theory, Euler's totient function counts the positive integers up to a given integer <math>n</math> that are relatively prime to <math>n</math>. It is written using the Greek letter phi as <math>\varphi(n)</math> or <math>\phi(n)</math>, and may also be called Euler's phi function. In other words, it is the number of integers <math>k</math> in the range <math>1\leq k\leq n</math> for which the greatest common divisor <math>\gcd(n,k)</math> is equal to 1. The integers <math>k</math> of this form are sometimes referred to as totatives of <math>n</math>.

For example, the totatives of <math>n=9</math> 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, since <math>\gcd(9,3)=\gcd(9,6)=3</math> and <math>\gcd(9,9)=9</math>. Therefore, <math>\varphi(9)=6</math>. As another example, <math>\varphi(1)=1</math> since for <math>n=1</math> the only integer in the range from 1 to <math>n</math> is 1 itself, and <math>\gcd(1,1)=1</math>.

Euler's totient function is a multiplicative function, meaning that if two numbers <math>m</math> and <math>n</math> are relatively prime, then <math>\varphi(mn)=\varphi(m)\varphi(n)</math>.

This function gives the order of the multiplicative group of integers modulo (the group of units of the ring <math>\Z/n\Z</math>). It is also used for defining the RSA encryption system.

History, terminology, and notation

Leonhard Euler introduced the function in 1763. 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 <math>\pi</math> to denote it: he wrote <math>\pi D</math> for "the multitude of numbers less than <math>D</math>, and which have no common divisor with it". This definition varies from the current definition for the totient function at <math>D=1</math> but is otherwise the same. The now-standard notation <math>\varphi(A)</math> comes from Gauss's 1801 treatise Disquisitiones Arithmeticae, although Gauss did not use parentheses around the argument and wrote <math>\varphi A</math>. 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, 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 <math>n</math> is defined as <math>n-\varphi(n)</math>. It counts the number of positive integers less than or equal to <math>n</math> that have at least one prime factor in common with <math>n</math>.

Computing Euler's totient function

There are several formulae for computing <math>\varphi(n)</math>.

Euler's product formula

It states

:<math>\varphi(n) =n \prod_{p\mid n} \left(1-\frac{1}{p}\right),</math>

where the product is over the distinct prime numbers dividing .

An equivalent formulation is

:<math>\varphi(n) = p_1^{k_1-1}(p_1{-}1)\,p_2^{k_2-1}(p_2{-}1)\cdots p_r^{k_r-1}(p_r{-}1),</math>

where <math>n = p_1^{k_1} p_2^{k_2} \cdots p_r^{k_r}</math> is the prime factorization of <math>n</math> (that is, <math>p_1, p_2,\ldots, p_r</math> are distinct prime numbers).

The proof of these formulae depends on two important facts.

Phi is a multiplicative function

This means that if <math>\gcd(m,n) = 1</math>, then <math>\varphi(m) \varphi(n) = \varphi(mn)</math>. Proof outline: Let <math>A,B,C</math> be the sets of positive integers which are coprime to and less than , , , respectively, so that <math>|A| = \varphi(m)</math>, etc. Then there is a bijection between <math>A\times B</math> and by the Chinese remainder theorem.

Value of phi for a prime power argument

If is prime and <math>k\geq1</math>, then

:<math>\varphi \left(p^k\right) = p^k-p^{k-1} = p^{k-1}(p-1) = p^k \left( 1 - \tfrac{1}{p} \right).</math>

Proof: Since is a prime number, the only possible values of <math>\gcd(p^{k},m)</math> are <math>1,p,p^{2},\dots,p^{k}</math>, and the only way to have <math>\gcd(p^{k},m)>1</math> is if is a multiple of , that is, <math> m \in \{ p, 2p, 3p, \ldots, p^{k-1} p= p^{k}\}</math>, and there are <math>p^{k-1}</math> such multiples not greater than <math>p^{k}</math>. Therefore, the other <math>p^{k}-p^{k-1}</math> numbers are all relatively prime to <math>p^{k}</math>.

Proof of Euler's product formula

The fundamental theorem of arithmetic states that if there is a unique expression <math>n = p_1^{k_1} p_2^{k_2} \cdots p_r^{k_r}, </math> where are prime numbers and each . (The case corresponds to the empty product.) Repeatedly using the multiplicative property of and the formula for gives

:<math>\begin{array} {rcl}

\varphi(n)&=& \varphi(p_1^{k_1})\, \varphi(p_2^{k_2})

\cdots\varphi(p_r^{k_r})\\[.1em]

&=& p_1^{k_1} \left(1- \frac{1}{p_1} \right) p_2^{k_2} \left(1- \frac{1}{p_2} \right) \cdots p_r^{k_r}\left(1- \frac{1}{p_r} \right)\\[.1em]

&=& p_1^{k_1} p_2^{k_2} \cdots p_r^{k_r} \left(1- \frac{1}{p_1} \right) \left(1- \frac{1}{p_2} \right) \cdots \left(1- \frac{1}{p_r} \right)\\[.1em]

&=&n \left(1- \frac{1}{p_1} \right)\left(1- \frac{1}{p_2} \right) \cdots\left(1- \frac{1}{p_r} \right).

\end{array}</math>

This gives both versions of Euler's product formula.

An alternative proof that does not require the multiplicative property instead uses the inclusion-exclusion principle applied to the set <math>\{1,2,\ldots,n\}</math>, excluding the sets of integers divisible by the prime divisors.

Example

:<math>\varphi(20)=\varphi(2^2 5)=20\,(1-\tfrac12)\,(1-\tfrac15)

=20\cdot\tfrac12\cdot\tfrac45=8.</math>

In words: the distinct prime factors of 20 are 2 and 5; half of the twenty integers from 1 to 20 are divisible by 2, leaving ten; a fifth of those are divisible by 5, leaving eight numbers coprime to 20; these are: 1, 3, 7, 9, 11, 13, 17, 19.

The alternative formula uses only integers:<math display="block">\varphi(20) = \varphi(2^2 5^1)= 2^{2-1}(2{-}1)\,5^{1-1}(5{-}1) = 2\cdot 1\cdot 1\cdot 4 = 8.</math>

Fourier transform

The totient is the discrete Fourier transform of the gcd, evaluated at 1. Let

:<math> \mathcal{F} \{ \mathbf{x} \}[m] = \sum\limits_{k=1}^n x_k \cdot e^</math>

where for . Then

:<math>\varphi (n) = \mathcal{F} \{ \mathbf{x} \}[1] = \sum\limits_{k=1}^n \gcd(k,n) e^{-2\pi i\frac{k}{n.</math>

The real part of this formula is

:<math>\varphi (n)=\sum\limits_{k=1}^n \gcd(k,n) \cos {\tfrac{2\pi k}{n

.</math>

For example, using <math>\cos\tfrac{\pi}5 = \tfrac{\sqrt 5+1}4

</math> and <math>\cos\tfrac{2\pi}5 = \tfrac{\sqrt 5-1}4

</math>:<math display="block">\begin{array}{rcl}

\varphi(10) &=& \gcd(1,10)\cos\tfrac{2\pi}{10} + \gcd(2,10)\cos\tfrac{4\pi}{10}

+ \gcd(3,10)\cos\tfrac{6\pi}{10}+\cdots+\gcd(10,10)\cos\tfrac{20\pi}{10}\\

&=& 1\cdot(\tfrac{\sqrt5+1}4) + 2\cdot(\tfrac{\sqrt5-1}4) +

1\cdot(-\tfrac{\sqrt5-1}4) + 2\cdot(-\tfrac{\sqrt5+1}4) + 5\cdot (-1) \\

&& +\ 2\cdot(-\tfrac{\sqrt5+1}4) + 1\cdot(-\tfrac{\sqrt5-1}4) + 2\cdot(\tfrac{\sqrt5-1}4) +

1\cdot(\tfrac{\sqrt5+1}4) + 10 \cdot (1) \\

&=& 4 .

\end{array}

</math>Unlike the Euler product and the divisor sum formula, this one does not require knowing the factors of . However, it does involve the calculation of the greatest common divisor of and every positive integer less than , which suffices to provide the factorization anyway.

Divisor sum

The property established by Gauss, that

:<math>\sum_{d\mid n}\varphi(d)=n,</math>

where the sum is over all positive divisors of , can be proven in several ways. (See Arithmetical function for notational conventions.)

One proof is to note that is also equal to the number of possible generators of the cyclic group ; specifically, if with , then is a generator for every coprime to . Since every element of generates a cyclic subgroup, and each subgroup is generated by precisely elements of , the formula follows. Equivalently, the formula can be derived by the same argument applied to the multiplicative group of the th roots of unity and the primitive th roots of unity.

The formula can also be derived from elementary arithmetic. For example, let and consider the positive fractions up to 1 with denominator 20:

:<math>

\tfrac{ 1}{20},\,\tfrac{ 2}{20},\,\tfrac{ 3}{20},\,\tfrac{ 4}{20},\,

\tfrac{ 5}{20},\,\tfrac{ 6}{20},\,\tfrac{ 7}{20},\,\tfrac{ 8}{20},\,

\tfrac{ 9}{20},\,\tfrac{10}{20},\,\tfrac{11}{20},\,\tfrac{12}{20},\,

\tfrac{13}{20},\,\tfrac{14}{20},\,\tfrac{15}{20},\,\tfrac{16}{20},\,

\tfrac{17}{20},\,\tfrac{18}{20},\,\tfrac{19}{20},\,\tfrac{20}{20}.

</math>

Put them into lowest terms:

:<math>

\tfrac{ 1}{20},\,\tfrac{ 1}{10},\,\tfrac{ 3}{20},\,\tfrac{ 1}{ 5},\,

\tfrac{ 1}{ 4},\,\tfrac{ 3}{10},\,\tfrac{ 7}{20},\,\tfrac{ 2}{ 5},\,

\tfrac{ 9}{20},\,\tfrac{ 1}{ 2},\,\tfrac{11}{20},\,\tfrac{ 3}{ 5},\,

\tfrac{13}{20},\,\tfrac{ 7}{10},\,\tfrac{ 3}{ 4},\,\tfrac{ 4}{ 5},\,

\tfrac{17}{20},\,\tfrac{ 9}{10},\,\tfrac{19}{20},\,\tfrac{1}{1}

</math>

These twenty fractions are all the positive ≤ 1 whose denominators are the divisors . The fractions with 20 as denominator are those with numerators relatively prime to 20, namely , , , , , , , ; by definition this is fractions. Similarly, there are fractions with denominator 10, and fractions with denominator 5, etc. Thus the set of twenty fractions is split into subsets of size for each dividing 20. A similar argument applies for any n.

Möbius inversion applied to the divisor sum formula gives

:<math> \varphi(n) = \sum_{d\mid n} \mu\left( d \right) \cdot \frac{n}{d} = n\sum_{d\mid n} \frac{\mu (d)}{d},</math>

where is the Möbius function, the multiplicative function defined by <math>\mu(p) = -1</math> and <math> \mu(p^k) = 0</math> for each prime and . This formula may also be derived from the product formula by multiplying out <math display="inline"> \prod_{p\mid n} (1 - \frac{1}{p}) </math> to get <math display="inline"> \sum_{d \mid n} \frac{\mu (d)}{d}. </math>

An example:<math display="block">

\begin{align}

\varphi(20) &= \mu(1)\cdot 20 + \mu(2)\cdot 10 +\mu(4)\cdot 5 +\mu(5)\cdot 4 + \mu(10)\cdot 2+\mu(20)\cdot 1\\[.5em]

&= 1\cdot 20 - 1\cdot 10 + 0\cdot 5 - 1\cdot 4 + 1\cdot 2 + 0\cdot 1 = 8.

\end{align}

</math>

Some values

The first 100 values are shown in the table and graph below:

thumb|Graph of the first 100 values

:{| class="wikitable" style="text-align: right"

|+ for

! +

! 1 || 2 || 3 || 4 || 5 || 6 || 7 || 8 || 9 || 10

! 0

| 1 || 1 || 2 || 2 || 4 || 2 || 6 || 4 || 6 || 4

! 10

| 10 || 4 || 12 || 6 || 8 || 8 || 16 || 6 || 18 || 8

! 20

| 12 || 10 || 22 || 8 || 20 || 12 || 18 || 12 || 28 || 8

! 30

| 30 || 16 || 20 || 16 || 24 || 12 || 36 || 18 || 24 || 16

! 40

| 40 || 12 || 42 || 20 || 24 || 22 || 46 || 16 || 42 || 20

! 50

| 32 || 24 || 52 || 18 || 40 || 24 || 36 || 28 || 58 || 16

! 60

| 60 || 30 || 36 || 32 || 48 || 20 || 66 || 32 || 44 || 24

! 70

| 70 || 24 || 72 || 36 || 40 || 36 || 60 || 24 || 78 || 32

! 80

| 54 || 40 || 82 || 24 || 64 || 42 || 56 || 40 || 88 || 24

! 90

| 72 || 44 || 60 || 46 || 72 || 32 || 96 || 42 || 60 || 40

In the graph at right the top line is an upper bound valid for all other than one, and attained if and only if is a prime number. A simple lower bound is <math>\varphi(n) \ge \sqrt{n/2} </math>, which is rather loose: in fact, the lower limit of the graph is proportional to .

<math>\sum_{1\le k\le n-1 \atop gcd(k,n)=1}\!\!k = \tfrac12 n\varphi(n) \quad \text{for }n>1</math>
<math>\sum_{k=1}^n\varphi(k) = \tfrac12 \left(1+ \sum_{k=1}^n \mu(k)\left\lfloor\frac{n}{k}\right\rfloor^2\right)

=\frac3{\pi^2}n^2+O\left(n(\log n)^\frac23(\log\log n)^\frac43\right)</math> ( cited in)

<math>\sum_{k=1}^n\varphi(k) =\frac3{\pi^2}n^2+O\left(n(\log n)^\frac23(\log\log n)^\frac13\right) </math> [Liu (2016)]
<math>\sum_{k=1}^n\frac{\varphi(k)}{k} = \sum_{k=1}^n\frac{\mu(k)}{k}\left\lfloor\frac{n}{k}\right\rfloor=\frac6{\pi^2}n+O\left((\log n)^\frac23(\log\log n)^\frac43\right)</math> 
<math>\sum_{k=1}^n\frac{1}{\varphi(k)} = \frac{315\,\zeta(3)}{2\pi^4}\left(\log n+\gamma-\sum_{p\text{ prime\frac{\log p}{p^2-p+1}\right)+O\left(\frac{(\log n)^\frac23}n\right)</math>  has important consequences. For instance it rules out uniform distribution of the values of <math>\varphi(n)</math> in the arithmetic progressions modulo <math>q</math> for any integer <math>q>1</math>.

For every fixed positive integer <math>q</math>, the relation <math>q|\varphi(n)</math> holds for almost all <math>n</math>, meaning for all but <math>o(x)</math> values of <math>n\le x</math> as <math>x\rightarrow\infty</math>.

This is an elementary consequence of the fact that the sum of the reciprocals of the primes congruent to 1 modulo <math>q</math> diverges, which itself is a corollary of the proof of Dirichlet's theorem on arithmetic progressions.

Generating functions

The Dirichlet series for may be written in terms of the Riemann zeta function as:

:<math>\sum_{n=1}^\infty \frac{\varphi(n)}{n^s}=\frac{\zeta(s-1)}{\zeta(s)}</math>

where the left-hand side converges for <math>\Re (s)>2</math>.

The Lambert series generating function is

:<math>\sum_{n=1}^{\infty} \frac{\varphi(n) q^n}{1-q^n}= \frac{q}{(1-q)^2}</math>

which converges for .

Both of these are proved by elementary series manipulations and the formulae for .

Growth rate

In the words of Hardy & Wright, the order of is "always 'nearly '."

First

:<math>\lim\sup \frac{\varphi(n)}{n}= 1,</math>

but as n goes to infinity, for all

:<math>\frac{\varphi(n)}{n^{1-\delta\rightarrow\infty.</math>

These two formulae can be proved by using little more than the formulae for and the divisor sum function .

In fact, during the proof of the second formula, the inequality

:<math>\frac {6}{\pi^2} < \frac{\varphi(n) \sigma(n)}{n^2} < 1,</math>

true for , is proved.

We also have

:<math>\lim\inf\frac{\varphi(n)}{n}\log\log n = e^{-\gamma}.</math>

Here is Euler's constant, , so and .

Proving this does not quite require the prime number theorem. Since goes to infinity, this formula shows that

:<math>\lim\inf\frac{\varphi(n)}{n}= 0.</math>

In fact, more is true.

:<math>\varphi(n) > \frac {n} {e^\gamma\; \log \log n + \frac {3} {\log \log n \quad\text{for } n>2</math>

and

:<math>\varphi(n) < \frac {n} {e^{ \gamma}\log \log n} \quad\text{for infinitely many } n.</math>

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."

:<math>\varphi(1)+\varphi(2)+\cdots+\varphi(n) = \frac{3n^2}{\pi^2}+O\left(n(\log n)^\frac23(\log\log n)^\frac43\right) \quad\text{as }n\rightarrow\infty,</math>

due to Arnold Walfisz, its proof exploiting estimates on exponential sums due to I. M. Vinogradov and N. M. Korobov.

By a combination of van der Corput's and Vinogradov's methods, H.-Q. Liu (On Euler's function.Proc. Roy. Soc. Edinburgh Sect. A 146 (2016), no. 4, 769–775)

improved the error term to

:<math>

O\left(n(\log n)^\frac23(\log\log n)^\frac13\right)

</math>

(this is currently the best known estimate of this type). The "Big " stands for a quantity that is bounded by a constant times the function of inside the parentheses (which is small compared to ).

This result can be used to prove that the probability of two randomly chosen numbers being relatively prime is .

Ratio of consecutive values

In 1950 Somayajulu proved

:<math>\begin{align}

\lim\inf \frac{\varphi(n+1)}{\varphi(n)}&= 0 \quad\text{and} \\[5px]

\lim\sup \frac{\varphi(n+1)}{\varphi(n)}&= \infty.

\end{align}</math>

In 1954 Schinzel and Sierpiński strengthened this, proving 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, and indeed every positive integer has a multiple which is an even nontotient.

The first few totient numbers are <math>1, 2, 4, 6, 8, 10, 12, 16, 18, 20</math>, see sequence .

The number of totient numbers up to a given limit is

:<math>\frac{x}{\log x}e^{ \big(C+o(1)\big)(\log\log\log x)^2 } </math>

for a constant .

If counted accordingly to multiplicity, the number of totient numbers up to a given limit is

:<math>\Big\vert\{ n : \varphi(n) \le x \}\Big\vert = \frac{\zeta(2)\zeta(3)}{\zeta(6)} \cdot x + R(x)</math>

where the error term is of order at most for any positive .

It is known that the multiplicity of exceeds infinitely often for any .

Ford's theorem

proved that for every integer there is a totient number of multiplicity : that is, for which the equation has exactly solutions; this result had previously been conjectured by Wacław Sierpiński, and it had been obtained as a consequence of Schinzel's hypothesis H.

Perfect totient numbers

A perfect totient number is an integer that is equal to the sum of its iterated totients. That is, we apply the totient function to a number n, apply it again to the resulting totient, and so on, until the number 1 is reached, and add together the resulting sequence of numbers; if the sum equals n, then n is a perfect totient number.

Applications

Cyclotomy

In the last section of the Disquisitiones Gauss proves that a regular -gon can be constructed with straightedge and compass if is a power of 2. If is a power of an odd prime number the formula for the totient says its totient can be a power of two only if is a first power and 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 -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 are

:2, 3, 4, 5, 6, 8, 10, 12, 15, 16, 17, 20, 24, 30, 32, 34, 40,... .

Prime number theorem for arithmetic progressions

The RSA cryptosystem

Setting up an RSA system involves choosing large prime numbers and , computing and , and finding two numbers and such that . The numbers and (the "encryption key") are released to the public, and (the "decryption key") is kept private.

A message, represented by an integer , where , is encrypted by computing .

It is decrypted by computing . Euler's Theorem can be used to show that if , then .

The security of an RSA system would be compromised if the number could be efficiently factored or if could be efficiently computed without factoring .

Unsolved problems

Lehmer's conjecture

If is prime, then . In 1932 D. H. Lehmer asked if there are any composite numbers such that divides . None are known.

In 1933 he proved that if any such exists, it must be odd, square-free, and divisible by at least seven primes (i.e. ). In 1980 Cohen and Hagis proved that and that . Further, Hagis showed that if 3 divides then and .

Carmichael's conjecture

This states that there is no number <math>n</math> with the property that for all other numbers <math>m</math>, <math>m\neq n</math>, <math>\varphi(m) \neq \varphi(n)</math>. See Ford's theorem above.

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.

Notes

References

The Disquisitiones Arithmeticae has been translated from Latin into English and German. The German edition includes all of Gauss's 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.

. See paragraph 24.3.2.
Dickson, Leonard Eugene, "History Of The Theory Of Numbers", vol 1, chapter 5 "Euler's Function, Generalizations; Farey Series", Chelsea Publishing 1952
.
.
.

==External links==

Euler's Phi Function and the Chinese Remainder Theorem — proof that is multiplicative
Euler's totient function calculator in JavaScript — up to 20 digits
Dineva, Rosica, The Euler Totient, the Möbius, and the Divisor Functions
Plytage, Loomis, Polhill Summing Up The Euler Phi Function

History, terminology, and notation

Computing Euler's totient function

Euler's product formula

Phi is a multiplicative function

Value of phi for a prime power argument

Proof of Euler's product formula

Example

Fourier transform

Divisor sum

Some values

Generating functions

Growth rate

Ratio of consecutive values

Ford's theorem

Perfect totient numbers

Applications

Cyclotomy

Prime number theorem for arithmetic progressions

The RSA cryptosystem

Unsolved problems

Lehmer's conjecture

Carmichael's conjecture

See also

Notes

References