In number theory, an arithmetic, arithmetical, or number-theoretic function is for most authors any function ''f''(''n'') whose domain is the positive integers and whose range is a subset of the complex numbers. Hardy & Wright include in their definition the requirement that an arithmetical function "expresses some arithmetical property of ''n''". An example of an arithmetic function is the divisor function whose value at a positive integer ''n'' is equal to the number of divisors of ''n''. There is a larger class of number-theoretic functions that do not fit the above definition, for example, the prime-counting functions. This article provides links to functions of both classes. Arithmetic functions are often extremely irregular (see table), but some of them have series expansions in terms of Ramanujan's sum.

Multiplicative and additive functions

An arithmetic function ''a'' is * completely additive if ''a''(''mn'') = ''a''(''m'') + ''a''(''n'') for all natural numbers ''m'' and ''n''; * completely multiplicative if ''a''(''mn'') = ''a''(''m'')''a''(''n'') for all natural numbers ''m'' and ''n''; Two whole numbers ''m'' and ''n'' are called coprime if their greatest common divisor is 1, that is, if there is no prime number that divides both of them. Then an arithmetic function ''a'' is * additive if ''a''(''mn'') = ''a''(''m'') + ''a''(''n'') for all coprime natural numbers ''m'' and ''n''; * multiplicative if ''a''(''mn'') = ''a''(''m'')''a''(''n'') for all coprime natural numbers ''m'' and ''n''.


\sum_p f(p)   and   \prod_p f(p)   mean that the sum or product is over all prime numbers: :\sum_p f(p) = f(2) + f(3) + f(5) + \cdots :\prod_p f(p)= f(2)f(3)f(5)\cdots. Similarly,   \sum_ f(p^k)   and   \prod_ f(p^k)   mean that the sum or product is over all prime powers with strictly positive exponent (so ''k'' = 0 is not included): :\sum_ f(p^k) = \sum_p\sum_ f(p^k) = f(2) + f(3) + f(4) +f(5) +f(7)+f(8)+f(9)+\cdots \sum_ f(d)   and   \prod_ f(d)   mean that the sum or product is over all positive divisors of ''n'', including 1 and ''n''. For example, if ''n'' = 12, :\prod_ f(d) = f(1)f(2) f(3) f(4) f(6) f(12).\ The notations can be combined:   \sum_ f(p)   and   \prod_ f(p)   mean that the sum or product is over all prime divisors of ''n''. For example, if ''n'' = 18, :\sum_ f(p) = f(2) + f(3),\ and similarly   \sum_ f(p^k)   and   \prod_ f(p^k)   mean that the sum or product is over all prime powers dividing ''n''. For example, if ''n'' = 24, :\prod_ f(p^k) = f(2) f(3) f(4) f(8).\

Ω(''n''), ''ω''(''n''), ''ν''''p''(''n'') – prime power decomposition

The fundamental theorem of arithmetic states that any positive integer ''n'' can be represented uniquely as a product of powers of primes:   n = p_1^\cdots p_k^   where ''p''1 < ''p''2 < ... < ''p''''k'' are primes and the ''aj'' are positive integers. (1 is given by the empty product.) It is often convenient to write this as an infinite product over all the primes, where all but a finite number have a zero exponent. Define the ''p''-adic valuation ν''p''(''n'') to be the exponent of the highest power of the prime ''p'' that divides ''n''. That is, if ''p'' is one of the ''p''''i'' then ''ν''''p''(''n'') = ''a''''i'', otherwise it is zero. Then :n=\prod_p p^. In terms of the above the prime omega functions ω and Ω are defined by :''ω''(''n'') = ''k'', :Ω(''n'') = ''a''1 + ''a''2 + ... + ''a''''k''. To avoid repetition, whenever possible formulas for the functions listed in this article are given in terms of ''n'' and the corresponding ''p''''i'', ''a''''i'', ω, and Ω.

Multiplicative functions

σ''k''(''n''), τ(''n''), ''d''(''n'') – divisor sums

σ''k''(''n'') is the sum of the ''k''th powers of the positive divisors of ''n'', including 1 and ''n'', where ''k'' is a complex number. σ1(''n''), the sum of the (positive) divisors of ''n'', is usually denoted by σ(''n''). Since a positive number to the zero power is one, σ0(''n'') is therefore the number of (positive) divisors of ''n''; it is usually denoted by ''d''(''n'') or τ(''n'') (for the German ''Teiler'' = divisors). :\sigma_k(n) = \prod_^ \frac = \prod_^ \left(1 + p_i^k + p_i^ + \cdots + p_i^\right). Setting ''k'' = 0 in the second product gives :\tau(n) = d(n) = (1 + a_)(1+a_)\cdots(1+a_).

φ(''n'') – Euler totient function

φ(''n''), the Euler totient function, is the number of positive integers not greater than ''n'' that are coprime to ''n''. :\varphi(n) = n \prod_ \left(1-\frac\right) =n \left(\frac\right)\left(\frac\right) \cdots \left(\frac\right) .

J''k''(''n'') – Jordan totient function

J''k''(''n''), the Jordan totient function, is the number of ''k''-tuples of positive integers all less than or equal to ''n'' that form a coprime (''k'' + 1)-tuple together with ''n''. It is a generalization of Euler's totient, . :J_k(n) = n^k \prod_ \left(1-\frac\right) =n^k \left(\frac\right)\left(\frac\right) \cdots \left(\frac\right) .

μ(''n'') – Möbius function

μ(''n''), the Möbius function, is important because of the Möbius inversion formula. See Dirichlet convolution, below. :\mu(n)=\begin (-1)^=(-1)^ &\text\; \omega(n) = \Omega(n)\\ 0&\text\;\omega(n) \ne \Omega(n).\end This implies that μ(1) = 1. (Because Ω(1) = ω(1) = 0.)

τ(''n'') – Ramanujan tau function

τ(''n''), the Ramanujan tau function, is defined by its generating function identity: :\sum_\tau(n)q^n=q\prod_(1-q^n)^. Although it is hard to say exactly what "arithmetical property of ''n''" it "expresses", (''τ''(''n'') is (2π)−12 times the ''n''th Fourier coefficient in the q-expansion of the modular discriminant function) it is included among the arithmetical functions because it is multiplicative and it occurs in identities involving certain σ''k''(''n'') and ''r''''k''(''n'') functions (because these are also coefficients in the expansion of modular forms).

''c''''q''(''n'') – Ramanujan's sum

''c''''q''(''n''), Ramanujan's sum, is the sum of the ''n''th powers of the primitive ''q''th roots of unity: :c_q(n)= \sum_ e^ . Even though it is defined as a sum of complex numbers (irrational for most values of ''q''), it is an integer. For a fixed value of ''n'' it is multiplicative in ''q'': :If ''q'' and ''r'' are coprime, then c_q(n)c_r(n)=c_(n).

''ψ''(''n'') - Dedekind psi function

The Dedekind psi function, used in the theory of modular functions, is defined by the formula : \psi(n) = n \prod_\left(1+\frac\right).

Completely multiplicative functions

λ(''n'') – Liouville function

''λ''(''n''), the Liouville function, is defined by :\lambda (n) = (-1)^.

''χ''(''n'') – characters

All Dirichlet characters ''χ''(''n'') are completely multiplicative. Two characters have special notations: The principal character (mod ''n'') is denoted by ''χ''0(''a'') (or ''χ''1(''a'')). It is defined as : \chi_0(a) = \begin 1 & \text \gcd(a,n) = 1, \\ 0 & \text \gcd(a,n) \ne 1. \end The quadratic character (mod ''n'') is denoted by the Jacobi symbol for odd ''n'' (it is not defined for even ''n''.): :\Bigg(\frac\Bigg) = \left(\frac\right)^\left(\frac\right)^\cdots \left(\frac\right)^. In this formula (\tfrac) is the Legendre symbol, defined for all integers ''a'' and all odd primes ''p'' by : \left(\frac\right) = \begin \;\;\,0 & \text a \equiv 0 \pmod p, \\+1 & \texta \not\equiv 0\pmod p \textx, \;a\equiv x^2\pmod p \\-1 & \text x. \end Following the normal convention for the empty product, \left(\frac\right) = 1.

Additive functions

''ω''(''n'') – distinct prime divisors

ω(''n''), defined above as the number of distinct primes dividing ''n'', is additive (see Prime omega function).

Completely additive functions

Ω(''n'') – prime divisors

Ω(''n''), defined above as the number of prime factors of ''n'' counted with multiplicities, is completely additive (see Prime omega function).

''ν''''p''(''n'') – ''p''-adic valuation of an integer ''n''

For a fixed prime ''p'', ''ν''''p''(''n''), defined above as the exponent of the largest power of ''p'' dividing ''n'', is completely additive.

Neither multiplicative nor additive

(''x''), Π(''x''), ''θ''(''x''), ''ψ''(''x'') – prime count functions

These important functions (which are not arithmetic functions) are defined for non-negative real arguments, and are used in the various statements and proofs of the prime number theorem. They are summation functions (see the main section just below) of arithmetic functions which are neither multiplicative nor additive. (''x''), the prime counting function, is the number of primes not exceeding ''x''. It is the summation function of the characteristic function of the prime numbers. :\pi(x) = \sum_1 A related function counts prime powers with weight 1 for primes, 1/2 for their squares, 1/3 for cubes, ... It is the summation function of the arithmetic function which takes the value 1/''k'' on integers which are the k-th power of some prime number, and the value 0 on other integers. :\Pi(x) = \sum_\frac. ''θ''(''x'') and ''ψ''(''x''), the Chebyshev functions, are defined as sums of the natural logarithms of the primes not exceeding ''x''. :\vartheta(x)=\sum_ \log p, : \psi(x) = \sum_ \log p. The Chebyshev function ''ψ''(''x'') is the summation function of the von Mangoldt function just below.

Λ(''n'') – von Mangoldt function

Λ(''n''), the von Mangoldt function, is 0 unless the argument ''n'' is a prime power , in which case it is the natural log of the prime ''p'': :\Lambda(n) = \begin\log p &\text n = 2,3,4,5,7,8,9,11,13,16,\ldots=p^k \text\\ 0&\text n=1,6,10,12,14,15,18,20,21,\dots \;\;\;\;\text. \end

''p''(''n'') – partition function

''p''(''n''), the partition function, is the number of ways of representing ''n'' as a sum of positive integers, where two representations with the same summands in a different order are not counted as being different: : p(n) = |\left\|.

λ(''n'') – Carmichael function

''λ''(''n''), the Carmichael function, is the smallest positive number such that a^\equiv 1 \pmod   for all ''a'' coprime to ''n''. Equivalently, it is the least common multiple of the orders of the elements of the multiplicative group of integers modulo ''n''. For powers of odd primes and for 2 and 4, ''λ''(''n'') is equal to the Euler totient function of ''n''; for powers of 2 greater than 4 it is equal to one half of the Euler totient function of ''n'': :\lambda(n) = \begin \;\;\phi(n) &\textn = 2,3,4,5,7,9,11,13,17,19,23,25,27,\dots\\ \tfrac12\phi(n)&\textn=8,16,32,64,\dots \end and for general ''n'' it is the least common multiple of λ of each of the prime power factors of ''n'': :\lambda(p_1^p_2^ \dots p_^) = \operatornamelambda(p_1^),\;\lambda(p_2^),\dots,\lambda(p_^)

''h''(''n'') – Class number

''h''(''n''), the class number function, is the order of the ideal class group of an algebraic extension of the rationals with discriminant ''n''. The notation is ambiguous, as there are in general many extensions with the same discriminant. See quadratic field and cyclotomic field for classical examples.

''r''''k''(''n'') – Sum of ''k'' squares

''r''''k''(''n'') is the number of ways ''n'' can be represented as the sum of ''k'' squares, where representations that differ only in the order of the summands or in the signs of the square roots are counted as different. :r_k(n) = |\|

''D''(''n'') – Arithmetic derivative

Using the Heaviside notation for the derivative, ''D''(''n'') is a function such that : D(n) = 1 if ''n'' prime, and : D(mn) = m D(n) + D(m) n (Product rule)

Summation functions

Given an arithmetic function ''a''(''n''), its summation function ''A''(''x'') is defined by : A(x) := \sum_ a(n) . ''A'' can be regarded as a function of a real variable. Given a positive integer ''m'', ''A'' is constant along open intervals ''m'' < ''x'' < ''m'' + 1, and has a jump discontinuity at each integer for which ''a''(''m'') ≠ 0. Since such functions are often represented by series and integrals, to achieve pointwise convergence it is usual to define the value at the discontinuities as the average of the values to the left and right: : A_0(m) := \frac12\left(\sum_ a(n) +\sum_ a(n)\right) = A(m) - \frac12 a(m) . Individual values of arithmetic functions may fluctuate wildly – as in most of the above examples. Summation functions "smooth out" these fluctuations. In some cases it may be possible to find asymptotic behaviour for the summation function for large ''x''. A classical example of this phenomenon is given by the divisor summatory function, the summation function of ''d''(''n''), the number of divisors of ''n'': :\liminf_d(n) = 2 :\limsup_\frac = \log 2 :\lim_\frac = 1. An average order of an arithmetic function is some simpler or better-understood function which has the same summation function asymptotically, and hence takes the same values "on average". We say that ''g'' is an ''average order'' of ''f'' if : \sum_ f(n) \sim \sum_ g(n) as ''x'' tends to infinity. The example above shows that ''d''(''n'') has the average order log(''n'').

Dirichlet convolution

Given an arithmetic function ''a''(''n''), let ''F''''a''(''s''), for complex ''s'', be the function defined by the corresponding Dirichlet series (where it converges): : F_a(s) := \sum_^\infty \frac . ''F''''a''(''s'') is called a generating function of ''a''(''n''). The simplest such series, corresponding to the constant function ''a''(''n'') = 1 for all ''n'', is ''ς''(''s'') the Riemann zeta function. The generating function of the Möbius function is the inverse of the zeta function: : \zeta(s)\,\sum_^\infty\frac=1, \;\;\mathfrak \,s >0. Consider two arithmetic functions ''a'' and ''b'' and their respective generating functions ''F''''a''(''s'') and ''F''''b''(''s''). The product ''F''''a''(''s'')''F''''b''(''s'') can be computed as follows: : F_a(s)F_b(s) = \left( \sum_^\frac \right)\left( \sum_^\frac \right) . It is a straightforward exercise to show that if ''c''(''n'') is defined by : c(n) := \sum_ a(i)b(j) = \sum_a(i)b\left(\frac\right) , then :F_c(s) =F_a(s) F_b(s). This function ''c'' is called the Dirichlet convolution of ''a'' and ''b'', and is denoted by a*b. A particularly important case is convolution with the constant function ''a''(''n'') = 1 for all ''n'', corresponding to multiplying the generating function by the zeta function: : g(n) = \sum_f(d). Multiplying by the inverse of the zeta function gives the Möbius inversion formula: : f(n) = \sum_\mu\left(\frac\right)g(d). If ''f'' is multiplicative, then so is ''g''. If ''f'' is completely multiplicative, then ''g'' is multiplicative, but may or may not be completely multiplicative.

Relations among the functions

There are a great many formulas connecting arithmetical functions with each other and with the functions of analysis, especially powers, roots, and the exponential and log functions. The page divisor sum identities contains many more generalized and related examples of identities involving arithmetic functions. Here are a few examples:

Dirichlet convolutions

: \sum_\mu(\delta)= \sum_\lambda\left(\frac\right)|\mu(\delta)|= \begin 1 & \text n=1\\ 0 & \text n\ne1 \end     where ''λ'' is the Liouville function. : \sum_\varphi(\delta) = n.       ::\varphi(n) =\sum_\mu\left(\frac\right)\delta =n\sum_\frac.       Möbius inversion : \sum_ J_k(d) = n^k.       :: J_k(n) =\sum_\mu\left(\frac\right)\delta^k =n^k\sum_\frac.       Möbius inversion : \sum_\delta^sJ_r(\delta)J_s\left(\frac\right) = J_(n)       : \sum_\varphi(\delta)d\left(\frac\right)= \sigma(n).       : \sum_|\mu(\delta)|= 2^.       ::|\mu(n)|=\sum_\mu\left(\frac\right)2^.       Möbius inversion : \sum_2^= d(n^2).       ::2^=\sum_\mu\left(\frac\right)d(\delta^2).       Möbius inversion : \sum_d(\delta^2)= d^2(n).       ::d(n^2)=\sum_\mu\left(\frac\right)d^2(\delta).       Möbius inversion : \sum_d\left(\frac\right)2^= d^2(n).       : \sum_\lambda(\delta)=\begin &1\text n \text\\ &0\text n \text \end     where λ is the Liouville function. : \sum_\Lambda(\delta) = \log n.       ::\Lambda(n)=\sum_\mu\left(\frac\right)\log(\delta).       Möbius inversion

Sums of squares

For all k \ge 4,\;\;\; r_k(n) > 0.     (Lagrange's four-square theorem). : : r_2(n) = 4\sum_\left(\frac\right), where the Kronecker symbol has the values : \left(\frac\right) = \begin +1&\textn\equiv 1 \pmod 4 \\ -1&\textn\equiv 3 \pmod 4\\ \;\;\;0&\textn\text.\\ \end There is a formula for r3 in the section on class numbers below. : r_4(n) = 8 \sum_d = 8 (2+(-1)^n)\sum_d = \begin 8\sigma(n)&\text n \text\\ 24\sigma\left(\frac\right)&\text n \text \end,     where ''ν'' = ''ν''2(''n'').     : r_6(n) = 16 \sum_ \chi\left(\frac\right)d^2 - 4\sum_ \chi(d)d^2, where \chi(n) = \left(\frac\right). Hardy & Wright, § 20.13 Define the function σ''k''*(''n'') as :\sigma_k^*(n) = (-1)^\sum_(-1)^d d^k= \begin \sum_ d^k=\sigma_k(n)&\text n \text\\ \sum_d^k -\sum_d^k&\text n \text. \end That is, if ''n'' is odd, σ''k''*(''n'') is the sum of the ''k''th powers of the divisors of ''n'', that is, σ''k''(''n''), and if ''n'' is even it is the sum of the ''k''th powers of the even divisors of ''n'' minus the sum of the ''k''th powers of the odd divisors of ''n''. : r_8(n) = 16\sigma_3^*(n).     Adopt the convention that Ramanujan's ''τ''(''x'') = 0 if ''x'' is not an integer. : r_(n) = \frac\sigma_^*(n) + \frac\left\    

Divisor sum convolutions

Here "convolution" does not mean "Dirichlet convolution" but instead refers to the formula for the coefficients of the product of two power series: : \left(\sum_^\infty a_n x^n\right)\left(\sum_^\infty b_n x^n\right) = \sum_^\infty \sum_^\infty a_i b_j x^ = \sum_^\infty \left(\sum_^n a_i b_\right) x^n = \sum_^\infty c_n x^n . The sequence c_n = \sum_^n a_i b_ is called the convolution or the Cauchy product of the sequences ''a''''n'' and ''b''''n''.
These formulas may be proved analytically (see Eisenstein series) or by elementary methods. : \sigma_3(n) = \frac\left\.    Ramanujan, ''On Certain Arithmetical Functions'', Table IV; ''Papers'', p. 146 : \sigma_5(n) = \frac\left\.    Koblitz, ex. III.2.8 : \begin \sigma_7(n) &=\frac\left\\\ &=\sigma_3(n) + 120\sum_\sigma_3(k)\sigma_3(n-k). \end     : \begin \sigma_9(n) &= \frac\left\\\ &= \frac\left\. \end     : \tau(n) = \frac\sigma_(n) + \frac\sigma_(n) - \frac\sum_\sigma_5(k)\sigma_5(n-k),     where ''τ''(''n'') is Ramanujan's function.     Since ''σ''''k''(''n'') (for natural number ''k'') and ''τ''(''n'') are integers, the above formulas can be used to prove congruences for the functions. See Ramanujan tau function for some examples. Extend the domain of the partition function by setting ''p''(0) = 1. : p(n)=\frac\sum_\sigma(k)p(n-k).       This recurrence can be used to compute ''p''(''n'').

Class number related

Peter Gustav Lejeune Dirichlet discovered formulas that relate the class number ''h'' of quadratic number fields to the Jacobi symbol. An integer ''D'' is called a fundamental discriminant if it is the discriminant of a quadratic number field. This is equivalent to ''D'' ≠ 1 and either a) ''D'' is squarefree and ''D'' ≡ 1 (mod 4) or b) ''D'' ≡ 0 (mod 4), ''D''/4 is squarefree, and ''D''/4 ≡ 2 or 3 (mod 4). Extend the Jacobi symbol to accept even numbers in the "denominator" by defining the Kronecker symbol: : \left(\frac\right) = \begin \;\;\,0&\text a \text \\(-1)^&\texta \text \end Then if ''D'' < −4 is a fundamental discriminant : \begin h(D) & = \frac \sum_^r\left(\frac\right)\\ & = \frac \sum_^\left(\frac\right). \end There is also a formula relating ''r''3 and ''h''. Again, let ''D'' be a fundamental discriminant, ''D'' < −4. Then : r_3(|D|) = 12\left(1-\left(\frac\right)\right)h(D).

Prime-count related

Let H_n = 1 + \frac12 + \frac13 + \cdots +\frac   be the ''n''th harmonic number.   Then : \sigma(n) \le H_n + e^\log H_n   is true for every natural number ''n'' if and only if the Riemann hypothesis is true.     The Riemann hypothesis is also equivalent to the statement that, for all ''n'' > 5040, :\sigma(n) < e^\gamma n \log \log n     (where γ is the Euler–Mascheroni constant).     This is Robin's theorem. : \sum_\nu_p(n) = \Omega(n). : \psi(x)=\sum_\Lambda(n).     : \Pi(x)= \sum_\frac.     : e^=\prod_p.     : e^= \operatorname,2,\dots,\lfloor x\rfloor    

Menon's identity

In 1965 P Kesava Menon proved : \sum_ \gcd(k-1,n) =\varphi(n)d(n). This has been generalized by a number of mathematicians. For example, B. Sury : \sum_ \gcd(k_1-1,k_2,\dots,k_s,n) =\varphi(n)\sigma_(n). N. Rao : \sum_ \gcd(k_1-a_1,k_2-a_2,\dots,k_s-a_s,n)^s =J_s(n)d(n), where ''a''1, ''a''2, ..., ''a''''s'' are integers, gcd(''a''1, ''a''2, ..., ''a''''s'', ''n'') = 1. László Fejes Tóth : \sum_ \gcd(k^2-1,m_1)\gcd(k^2-1,m_2) =\varphi(n)\sum_ \varphi(\gcd(d_1, d_2))2^, where ''m''1 and ''m''2 are odd, ''m'' = lcm(''m''1, ''m''2). In fact, if ''f'' is any arithmetical function : \sum_ f(\gcd(k-1,n)) =\varphi(n)\sum_\frac, where * stands for Dirichlet convolution.


Let ''m'' and ''n'' be distinct, odd, and positive. Then the Jacobi symbol satisfies the law of quadratic reciprocity: : \left(\frac\right) \left(\frac\right) = (-1)^.     Let ''D''(''n'') be the arithmetic derivative. Then the logarithmic derivative : \frac = \sum_ \frac Let ''λ''(''n'') be Liouville's function. Then :|\lambda(n)|\mu(n) =\lambda(n)|\mu(n)| = \mu(n),     and :\lambda(n)\mu(n) = |\mu(n)| =\mu^2(n).     Let ''λ''(''n'') be Carmichael's function. Then :\lambda(n)\mid \phi(n).     Further, :\lambda(n)= \phi(n) \textn=\begin1,2, 4;\\ 3,5,7,9,11, \ldots \text p^k \textp\text;\\ 6,10,14,18,\ldots \text 2p^k\textp\text. \end See Multiplicative group of integers modulo n and Primitive root modulo n.   : 2^ \le d(n) \le 2^.     : \frac<\frac < 1.     : \begin c_q(n) &=\frac\phi(q)\\ &=\sum_\mu\left(\frac\right)\delta. \end         Note that  \phi(q) = \sum_\mu\left(\frac\right)\delta.     :c_q(1) = \mu(q). :c_q(q) = \phi(q). : \sum_d^(\delta) = \left(\sum_d(\delta)\right)^2.       Compare this with 13 + 23 + 33 + ... + ''n''3 = (1 + 2 + 3 + ... + ''n'')2 :d(uv) = \sum_\mu(\delta)d\left(\frac\right)d\left(\frac\right).     :\sigma_k(u)\sigma_k(v) = \sum_\delta^k\sigma_k\left(\frac\right).     :\tau(u)\tau(v) = \sum_\delta^\tau\left(\frac\right),     where ''τ''(''n'') is Ramanujan's function.    Apostol, ''Modular Functions ...'', ch. 6 eq. 3

First 100 values of some arithmetic functions



* * * * * * * * * * * * * * * * * *

Further reading


External links

* * Matthew Holden, Michael Orrison, Michael Varbl
Yet another Generalization of Euler's Totient Function
* Huard, Ou, Spearman, and Williams
Elementary Evaluation of Certain Convolution Sums Involving Divisor Functions
* Dineva, Rosica
The Euler Totient, the Möbius, and the Divisor Functions
* László Tóth
Menon's Identity and arithmetical sums representing functions of several variables
{{Authority control Category:Functions and mappings