Proofs involving the totient function

Proofs involving the totient function

This page provides proofs for identities involving the totient function varphi(k) and the Möbius function mu(k).

um of integers relatively prime to and less than or equal to "n"

Claim::sum_{1le kle n atop {gcd(k,n)=1 k = frac{1}{2} , varphi(n) , nquad ext{ for integers }n ge 2.Since gcd(n,n) e 1and gcd(k,n)=gcd(n-k,n) for all integers k, a change of index yields:sum_{1le kle n atop {gcd(k,n)=1 k ,= sum_{1le kle n atop {gcd(k,n)=1 (n-k).Therefore:2sum_{1le kle n atop {gcd(k,n)=1 k ,=sum_{1le kle n atop {gcd(k,n)=1 (k+n-k) = varphi(n),n.

Proofs of totient identities involving the floor function

The proof of the identity

:sum_{k=1}^nfrac{varphi(k)}{k} = sum_{k=1}^nfrac{mu(k)}{k}leftlfloorfrac{n}{k} ight floor

is by mathematical induction on n. The base case is n=1 and we see that the claim holds:

:varphi(1)/1 = 1 = frac{mu(1)}{1} leftlfloor 1 ight floor.

For the induction step we need to prove that

:frac{varphi(n+1)}{n+1} =

sum_{k=1}^nfrac{mu(k)}{k}left(leftlfloorfrac{n+1}{k} ight floor - leftlfloorfrac{n}{k} ight floor ight)+ frac{mu(n+1)}{n+1}.

The key observation is that

:leftlfloorfrac{n+1}{k} ight floor - leftlfloorfrac{n}{k} ight floor; = ; egin{cases} 1, & mbox{if }k|(n+1) \ 0, & mbox{otherwise, }end{cases}

so that the sum is:sum_{k|n+1,; k

Now the fact that:sum_{k|n+1} frac{mu(k)}{k} = frac{varphi(n+1)}{n+1}

is a basic totient identity. To see that it holds, let p_1^{v_1} p_2^{v_2} ldots p_q^{v_q}be the prime factorization of n+1. Then

:frac{varphi(n+1)}{n+1} = prod_{l=1}^q left( 1 - frac{1}{p_l} ight)= sum_{k|n+1} frac{mu(k)}{k}

by definition of mu(k). This concludes the proof.cn|date=December 2007

An alternate proof proceeds by substitutingfrac{varphi(k)}{k} = sum_{d|k}frac{mu(d)}{d}directly into the left side of the identity, givingsum_{k=1}^n sum_{d|k}frac{mu(d)}{d}.

Now we ask how often the term egin{matrix}frac{mu(d)}{d}end{matrix} occurs in the double sum. The answer is that it occurs for every multiple k of d, but there are precisely egin{matrix}leftlfloorfrac{n}{d} ight floorend{matrix} such multiples, which means that the sum is

:sum_{d=1}^nfrac{mu(d)}{d}leftlfloorfrac{n}{d} ight floor

as claimed.cn|date=December 2007

----

The trick where zero values ofcn|date=December 2007 egin{matrix}leftlfloorfrac{n+1}{k} ight floor - leftlfloorfrac{n}{k} ight floorend{matrix} are filtered out may also be used to prove the identity

:sum_{k=1}^nvarphi(k) = frac{1}{2}left(1+ sum_{k=1}^n mu(k)leftlfloorfrac{n}{k} ight floor^2 ight).

The base case is n=1 and we have

:varphi(1) = 1 = frac{1}{2} left(1+ mu(1) leftlfloorfrac{1}{1} ight floor^2 ight)

and it holds. The induction step requires us to show that

:varphi(n+1) = frac{1}{2} sum_{k=1}^n mu(k)left(leftlfloorfrac{n+1}{k} ight floor^2 - leftlfloorfrac{n}{k} ight floor^2 ight); + ; frac{1}{2} ; mu(n+1) ; leftlfloorfrac{n+1}{n+1} ight floor^2.

Next observe thatcn|date=December 2007

:leftlfloorfrac{n+1}{k} ight floor^2 - leftlfloorfrac{n}{k} ight floor^2; = ; egin{cases} 2frac{n+1}{k} - 1, & mbox{if }k|(n+1) \ 0, & mbox{otherwise.}end{cases}

This gives the following for the sum:frac{1}{2} sum_{k|n+1, ; k

Treating the two inner terms separately, we get

:(n+1) sum_{k|n+1} frac{mu(k)}{k} - frac{1}{2} sum_{k|n+1} mu(k).

The first of these two is precisely varphi(n+1) as we saw earlier, and the second is zero, by a basic property of the Möbius function (using the same factorization of n+1 as above, we have egin{matrix} sum_{k|n+1} mu(k) = prod_{l=1}^q (1-1) = 0 end{matrix}.) This concludes the proof.

This result may also be proved by inclusion-exclusion. Rewrite the identity as:-1 + 2sum_{k=1}^nvarphi(k) = sum_{k=1}^n mu(k)leftlfloorfrac{n}{k} ight floor^2.

Now we see that the left side counts the number of lattice points ("a", "b") in [1,"n"] × [1,"n"] where "a" and "b" are relatively prime to each other. Using the sets A_p where "p" is a prime less than or equal to "n" to denote the set of points where both coordinates are divisible by "p" we have

: left| igcup_p A_p ight| = sum_p left| A_p ight| ;- ; sum_{p

This formula counts the number of pairs where "a" and "b" are not relatively prime to each other.The cardinalities are as follows:

:left| A_p ight| = leftlfloor frac{n}{p} ight floor^2, ;left| A_p cap A_q ight| = leftlfloor frac{n}{pq} ight floor^2, ;left| A_p cap A_q cap A_r ight| = leftlfloor frac{n}{pqr} ight floor^2, ; ldots

and the signs are egin{matrix}-mu(pqrcdots)end{matrix}, hence the number of points with relatively prime coordinates is

mu(1), n^2; + ; sum_p mu(p) leftlfloor frac{n}{p} ight floor^2; + ; sum_{p

but this is precisely sum_{k=1}^n mu(k) leftlfloor frac{n}{k} ight floor^2 and we have the claim.

Average order of the totient

We will use the last formula of the preceding section to prove the following result::frac{1}{n^2} sum_{k=1}^n varphi(k)= frac{3}{pi^2} + mathcal{O}left(frac{log n }{n} ight)

Using x-1 < lfloor x floor le x we have the upper bound:frac{1}{2 n^2} left(1+ sum_{k=1}^n mu(k)frac{n^2}{k^2} ight) =frac{1}{2 n^2} + frac{1}{2} sum_{k=1}^n frac{mu(k)}{k^2}

and the lower bound:frac{1}{2 n^2} left(1+ sum_{k=1}^n mu(k)left(frac{n^2}{k^2} - 2frac{n}{k} + 1 ight) ight)which is:frac{1}{2 n^2} + frac{1}{2} sum_{k=1}^n frac{mu(k)}{k^2}- frac{1}{n} sum_{k=1}^n frac{mu(k)}{k}+ frac{1}{2 n^2} sum_{k=1}^n mu(k)

Working with the last two terms and using the asymptotic expansion of the nth
harmonic number we have:- frac{1}{n} sum_{k=1}^n frac{mu(k)}{k} >- frac{1}{n} sum_{k=1}^n frac{1}{k} = - frac{1}{n} H_n> -frac{1}{n} (log n +1)and:frac{1}{2 n^2} sum_{k=1}^n mu(k) > - frac{1}{2 n}.

Now we check the order of the terms in the upper and lower bound. The termegin{matrix}sum_{k=1}^n frac{mu(k)}{k^2}end{matrix} ismathcal{O}(1) by comparison with zeta(2), wherezeta(s) is the Riemann zeta function. The next largest term is the logarithmic term from the lower bound.

So far we have shown that:frac{1}{n^2} sum_{k=1}^n varphi(k)= frac{1}{2} sum_{k=1}^n frac{mu(k)}{k^2} + mathcal{O}left(frac{log n }{n} ight)

It remains to evaluate egin{matrix}sum_{k=1}^n frac{mu(k)}{k^2}end{matrix}asymptotically, which we have seen converges. The Euler product for the Riemann zeta function is:zeta(s) = prod_p left(1 - frac{1}{p^s} ight)^{-1}mbox{ for } Re(s)>1.

Now it follows immediately from the definition of the Möbius function that:frac{1}{zeta(s)} = prod_p left(1 - frac{1}{p^s} ight)= sum_{n ge 1} frac{mu(n)}{n^s}.

This means that:frac{1}{2} sum_{k=1}^n frac{mu(k)}{k^2} = frac{1}{2} frac{1}{zeta(2)} + mathcal{O}left(frac{1}{n} ight)where the integralegin{matrix} int_{n+1}^infty frac{1}{t^2} dtend{matrix}was used to estimateegin{matrix} sum_{k>n} frac{mu(k)}{k^2}end{matrix}.But egin{matrix}frac{1}{2} frac{1}{zeta(2)} = frac{3}{pi^2}end{matrix}and we have established the claim.

Average order of &phi;("n")/"n"

The material of the preceding section, together with the identity

:sum_{k=1}^nfrac{varphi(k)}{k} = sum_{k=1}^nfrac{mu(k)}{k}leftlfloorfrac{n}{k} ight floor

also yields a proof that

:frac{1}{n} sum_{k=1}^n frac{varphi(k)}{k} = frac{6}{pi^2} + mathcal{O}left(frac{log n }{n} ight).

Reasoning as before, we have the upper bound:frac{1}{n} sum_{k=1}^nfrac{mu(k)}{k}frac{n}{k} = sum_{k=1}^n frac{mu(k)}{k^2}and the lower bound:-frac{1}{n} sum_{k=1}^nfrac{mu(k)}{k} +sum_{k=1}^n frac{mu(k)}{k^2}.

Now apply the estimates from the preceding section to obtain the result.

Inequalities

We first show that

:lim inf frac{varphi (n)}{n}=0 mbox{ and } lim sup frac{varphi (n)}{n}=1.

The latter holds because when "n" is a power of a prime "p", we have

:frac{varphi (n)}{n} = 1 - frac{1}{p},

which gets arbitrarily close to 1 for "p" large enough (and we can take "p" as large as we please since there are infinitely many primes).

To see the former, let "n""k" be the product of the first "k" primes, call them p_1, p_2, ..., p_k. Let

:r_k = frac{varphi (n_k)}{n_k} = prod_{i=1}^k left( 1 - frac{1}{p_i} ight).

Then

:frac{1}{r_k} = prod_{i=1}^k left( 1 - frac{1}{p_i} ight)^{-1} >sum_{m=1}^{p_k} frac{1}{m} = H_{p_k},

a harmonic number. Hence, by the well-known bound H_n > log n, we have

:frac{1}{r_k} > log p_k.

Since the logarithm is unbounded, taking "k" arbitrarily large ensures that "r""k" achieves values arbitrarily close to zero.

We use the same factorization of "n" as in the first section to prove that

:frac {6 n^2} {pi^2} < varphi(n) sigma(n) < n^2.

Note that

: varphi(n) sigma(n) = n prod_{l=1}^q left(1 - frac{1}{p_l} ight) prod_{l=1}^q frac{p_l^{v_l+1}-1}{p_l-1} =nprod_{l=1}^q frac{p_l-1}{p_l} ; frac{p_l^{v_l+1}-1}{p_l-1}

which is

:n prod_{l=1}^q left(p_l^{v_l}-frac{1}{p_l} ight) =n^2 prod_{l=1}^q left(1 - frac{1}{p_l^{v_l+1 ight).

The upper bound follows immediately since

:prod_{l=1}^q left(1 - frac{1}{p_l^{v_l+1 ight) < 1.

We come arbitrarily close to this bound when "n" is prime. For the lower bound, note that

: prod_{l=1}^q left(1 - frac{1}{p_l^{v_l+1 ight) geprod_{l=1}^q left(1 - frac{1}{p_l^2} ight) >prod_p left(1 - frac{1}{p^2} ight),

where the product is over all primes. We have already seen this product, as in

: prod_p left(1 - frac{1}{p^s} ight) =sum_{nge 1} frac{mu(n)}{n^s} = frac{1}{zeta(s)}

so that

:prod_p left(1 - frac{1}{p^2} ight) = frac{1}{zeta(2)}= frac{6}{pi^2}

and we have the claim. The values of "n" that come closest to this bound are products of the first "k" primes.

External links

* Chris K. Caldwell, " [http://primes.utm.edu/notes/relprime.html What is the probability that gcd(n,m)=1?] "


Wikimedia Foundation. 2010.

Игры ⚽ Поможем написать курсовую

Look at other dictionaries:

  • Divisor function — σ0(n) up to n = 250 Sigma function σ …   Wikipedia

  • Arithmetic function — In number theory, an arithmetic (or arithmetical) function is a real or complex valued function ƒ(n) defined on the set of natural numbers (i.e. positive integers) that expresses some arithmetical property of n. [1] An example of an arithmetic… …   Wikipedia

  • List of mathematics articles (P) — NOTOC P P = NP problem P adic analysis P adic number P adic order P compact group P group P² irreducible P Laplacian P matrix P rep P value P vector P y method Pacific Journal of Mathematics Package merge algorithm Packed storage matrix Packing… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”