Jump to content

Möbius inversion formula: Difference between revisions

From Wikipedia, the free encyclopedia
Content deleted Content added
m Contributions of Weisner, Hall, and Rota: Journal cites, Added 1 doi to a journal cite using AWB (10870)
 
(38 intermediate revisions by 31 users not shown)
Line 1: Line 1:
{{short description|Relation between pairs of arithmetic functions}}
:''Möbius transform redirects here. It should not be confused with [[Möbius transformation]].''
{{redirect-distinguish|Möbius transform|Möbius transformation}}
In [[mathematics]], the classic '''Möbius inversion formula''' was introduced into [[number theory]] during the 19th century by [[August Ferdinand Möbius]].
In [[mathematics]], the classic '''Möbius inversion formula''' is a relation between pairs of [[arithmetic function]]s, each defined from the other by sums over [[divisor]]s. It was introduced into [[number theory]] in 1832 by [[August Ferdinand Möbius]].<ref>{{Harvnb|Möbius|1832|pp=105-123}}</ref>


Other Möbius inversion formulas are obtained when different [[Locally finite poset|locally finite partially ordered set]]s replace the classic case of the natural numbers ordered by divisibility; for an account of those, see [[incidence algebra]].
A large generalization of this formula applies to summation over an arbitrary [[Locally finite poset|locally finite partially ordered set]], with Möbius' classical formula applying to the set of the natural numbers ordered by divisibility: see [[incidence algebra]].


==Statement of the formula==
==Statement of the formula==
The classic version states that if ''g'' and ''f'' are [[arithmetic function]]s satisfying
The classic version states that if {{mvar|g}} and {{mvar|f}} are [[arithmetic function]]s satisfying


: <math>g(n)=\sum_{d\,\mid \,n}f(d)\quad\text{for every integer }n\ge 1</math>
: <math>g(n)=\sum_{d \mid n}f(d)\quad\text{for every integer }n\ge 1</math>


then
then


:<math>f(n)=\sum_{d\,\mid\, n}\mu(d)g(n/d)\quad\text{for every integer }n\ge 1</math>
:<math>f(n)=\sum_{d \mid n}\mu(d)g\left(\frac{n}{d}\right)\quad\text{for every integer }n\ge 1</math>


where μ is the [[Möbius function]] and the sums extend over all positive [[divisor]]s ''d'' of ''n''. In effect, the original ''f''(''n'') can be determined given ''g''(''n'') by using the inversion formula. The two sequences are said to be '''Möbius transforms''' of each other.
where {{mvar|μ}} is the [[Möbius function]] and the sums extend over all positive [[divisor]]s {{mvar|d}} of {{mvar|n}} (indicated by <math>d \mid n</math> in the above formulae). In effect, the original {{math|''f''(''n'')}} can be determined given {{math|''g''(''n'')}} by using the inversion formula. The two sequences are said to be '''Möbius transforms''' of each other.


The formula is also correct if ''f'' and ''g'' are functions from the positive integers into some [[abelian group]] (viewed as a <math>\mathbb{Z}</math>-[[module (mathematics)|module]]).
The formula is also correct if {{mvar|f}} and {{mvar|g}} are functions from the positive integers into some [[abelian group]] (viewed as a {{math|'''Z'''}}-[[module (mathematics)|module]]).


In the language of [[Dirichlet convolution]]s, the first formula may be written as
In the language of [[Dirichlet convolution]]s, the first formula may be written as


:<math>g=f*1</math>
:<math>g=\mathit{1}*f</math>


where ''*'' denotes the Dirichlet convolution, and ''1'' is the [[constant function]] <math>1(n)=1</math>. The second formula is then written as
where {{math|∗}} denotes the Dirichlet convolution, and {{math|''1''}} is the [[constant function]] {{math|1=''1''(''n'') = 1}}. The second formula is then written as


:<math>f=\mu * g.</math>
:<math>f=\mu * g.</math>
Line 27: Line 28:
Many specific examples are given in the article on [[multiplicative function]]s.
Many specific examples are given in the article on [[multiplicative function]]s.


The theorem follows because <math>*</math> is (commutative and) associative, and <math>1 * \mu = \epsilon</math>, where <math>\epsilon</math> is the identity function for the Dirichlet convolution, taking values <math>\epsilon(1) = 1, \epsilon(n) = 0</math> for all <math>n > 1</math>. Thus <math>\mu * g = \mu * (1 * f) = (\mu * 1) * f = \epsilon * f = f</math>.
The theorem follows because {{math|∗}} is (commutative and) associative, and {{math|1=''1'' ''μ'' = ''ε''}}, where {{mvar|ε}} is the identity function for the Dirichlet convolution, taking values {{math|1=''ε''(1) = 1}}, {{math|1=''ε''(''n'') = 0}} for all {{math|''n'' > 1}}. Thus
:<math>\mu * g = \mu * (\mathit{1} * f) = (\mu * \mathit{1}) * f = \varepsilon * f = f</math>.
Replacing <math>f, g</math> by <math>\ln f, \ln g</math>, we obtain the product version of the Möbius inversion formula:

:<math>g(n) = \prod_{d|n} f(d) \iff f(n) = \prod_{d|n} g\left(\frac{n}{d}\right)^{\mu(d)}, \forall n \geq 1.</math>


==Series relations==
==Series relations==
Line 40: Line 45:
is its transform. The transforms are related by means of series: the [[Lambert series]]
is its transform. The transforms are related by means of series: the [[Lambert series]]


:<math>\sum_{n=1}^\infty a_n x^n =
:<math>\sum_{n=1}^\infty a_n x^n = \sum_{n=1}^\infty b_n \frac{x^n}{1-x^n}</math>
\sum_{n=1}^\infty b_n \frac{x^n}{1-x^n}</math>


and the [[Dirichlet series]]:
and the [[Dirichlet series]]:


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


where <math>\zeta(s)</math> is the [[Riemann zeta function]].
where {{math|''ζ''(''s'')}} is the [[Riemann zeta function]].


==Repeated transformations==
==Repeated transformations==
Given an arithmetic function, one can generate a bi-infinite sequence of other arithmetic functions by repeatedly applying the first summation.
Given an arithmetic function, one can generate a bi-infinite sequence of other arithmetic functions by repeatedly applying the first summation.


For example, if one starts with [[Euler's totient function]] <math>\varphi</math>, and repeatedly applies the transformation process, one obtains:
For example, if one starts with [[Euler's totient function]] {{mvar|φ}}, and repeatedly applies the transformation process, one obtains:


#<math>\varphi</math> the totient function
#{{mvar|φ}} the totient function
#<math>\varphi*1=\operatorname{Id}</math> where <math>\operatorname{Id}(n)=n</math> is the [[identity function]]
#{{math|1=''φ'' ∗ ''1'' = ''I''}}, where {{math|1=''I''(''n'') = ''n''}} is the [[identity function]]
#<math>\operatorname{Id} *1 =\sigma_1 =\sigma</math>, the [[divisor function]]
#{{math|1=''I'' ∗ ''1'' = ''σ''<sub>1</sub> = ''σ''}}, the [[divisor function]]


If the starting function is the Möbius function itself, the list of functions is:
If the starting function is the Möbius function itself, the list of functions is:
#<math>\mu</math>, the Möbius function
#{{mvar|μ}}, the Möbius function
#<math>\mu*1 = \varepsilon</math> where <math>\varepsilon(n) = \begin{cases} 1, & \mbox{if }n=1 \\ 0, & \mbox{if }n>1 \end{cases} </math> is the [[unit function]]
#{{math|1=''μ'' ∗ ''1'' = ''ε''}} where <math display="block">\varepsilon(n) = \begin{cases} 1, & \text{if }n=1 \\ 0, & \text{if }n>1 \end{cases} </math> is the [[unit function]]
#<math>\varepsilon*1 = 1 </math>, the [[constant function]]
#{{math|1=''ε'' ∗ ''1'' = ''1''}}, the [[constant function]]
#<math>1*1=\sigma_0=\operatorname{d}=\tau</math>, where <math>\operatorname{d}=\tau</math> is the number of divisors of ''n'', (see [[divisor function]]).
#{{math|1=''1'' ∗ ''1'' = ''σ''<sub>0</sub> = d = ''τ''}}, where {{math|1=d = ''τ''}} is the number of divisors of {{mvar|n}}, (see [[divisor function]]).


Both of these lists of functions extend infinitely in both directions. The Möbius inversion formula enables these lists to be traversed backwards.
Both of these lists of functions extend infinitely in both directions. The Möbius inversion formula enables these lists to be traversed backwards.


As an example the sequence starting in <math>\varphi</math> is:
As an example the sequence starting with {{mvar|φ}} is:


<math>
:<math>f_n =
f_n =
\begin{cases}
\begin{cases}
\underbrace{\mu * \ldots * \mu}_{-n \text{ factors}} * \varphi & \text{if } n < 0 \\
\underbrace{\mu * \ldots * \mu}_{-n \text{ factors}} * \varphi & \text{if } n < 0 \\[8px]
\varphi & \text{if } n = 0 \\
\varphi & \text{if } n = 0 \\[8px]
\varphi * \underbrace{1* \ldots * 1}_{n \text{ factors}} & \text{if } n > 0
\varphi * \underbrace{\mathit{1}* \ldots * \mathit{1}}_{n \text{ factors}} & \text{if } n > 0
\end{cases}
\end{cases}
</math>
</math>
Line 81: Line 83:


==Generalizations==
==Generalizations==
A related inversion formula more useful in [[combinatorics]] is as follows: suppose {{math|''F''(''x'')}} and {{math|''G''(''x'')}} are [[complex number|complex]]-valued [[function (mathematics)|function]]s defined on the [[interval (mathematics)|interval]] {{closed-open|1, }} such that
{{See also|Incidence algebra}}
A related inversion formula more useful in [[combinatorics]] is as follows: suppose ''F''(''x'') and ''G''(''x'') are [[complex number|complex]]-valued [[function (mathematics)|function]]s defined on the [[interval (mathematics)|interval]] <nowiki>[1,∞)</nowiki> such that


:<math>G(x) = \sum_{1 \le n \le x}F(x/n)\quad\mbox{ for all }x\ge 1</math>
:<math>G(x) = \sum_{1 \le n \le x}F\left(\frac{x}{n}\right)\quad\mbox{ for all }x\ge 1</math>


then
then


:<math>F(x) = \sum_{1 \le n \le x}\mu(n)G(x/n)\quad\mbox{ for all }x\ge 1.</math>
:<math>F(x) = \sum_{1 \le n \le x}\mu(n)G\left(\frac{x}{n}\right)\quad\mbox{ for all }x\ge 1.</math>


Here the sums extend over all positive integers ''n'' which are less than or equal to ''x''.
Here the sums extend over all positive integers {{mvar|n}} which are less than or equal to {{mvar|x}}.


This in turn is a special case of a more general form. If <math>\alpha(n)</math> is an [[arithmetic function]] possessing a [[Dirichlet inverse]] <math>\alpha^{-1}(n)</math>, then if one defines
This in turn is a special case of a more general form. If {{math|''α''(''n'')}} is an [[arithmetic function]] possessing a [[Dirichlet inverse]] {{math|''α''<sup>−1</sup>(''n'')}}, then if one defines


:<math>G(x) = \sum_{1 \le n \le x}\alpha (n) F(x/n)\quad\mbox{ for all }x\ge 1</math>
:<math>G(x) = \sum_{1 \le n \le x}\alpha (n) F\left(\frac{x}{n}\right)\quad\mbox{ for all }x\ge 1</math>


then
then


:<math>F(x) = \sum_{1 \le n \le x}\alpha^{-1}(n)G(x/n)\quad\mbox{ for all }x\ge 1.</math>
:<math>F(x) = \sum_{1 \le n \le x}\alpha^{-1}(n)G\left(\frac{x}{n}\right)\quad\mbox{ for all }x\ge 1.</math>


The previous formula arises in the special case of the constant function <math>\alpha(n)=1</math>, whose [[Dirichlet inverse]] is <math>\alpha^{-1}(n)=\mu(n)</math>.
The previous formula arises in the special case of the constant function {{math|1=''α''(''n'') = 1}}, whose [[Dirichlet inverse]] is {{math|1=''α''<sup>−1</sup>(''n'') = ''μ''(''n'')}}.


A particular application of the first of these extensions arises if we have (complex-valued) functions ''f''(''n'') and ''g''(''n'') defined on the positive integers, with
A particular application of the first of these extensions arises if we have (complex-valued) functions {{math|''f''(''n'')}} and {{math|''g''(''n'')}} defined on the positive integers, with


:<math>g(n) = \sum_{1 \le m \le n}f\left(\left\lfloor \frac{n}{m}\right\rfloor\right)\quad\mbox{ for all } n\ge 1.</math>
:<math>g(n) = \sum_{1 \le m \le n}f\left(\left\lfloor \frac{n}{m}\right\rfloor\right)\quad\mbox{ for all } n\ge 1.</math>


By defining <math>F(x) = f(\lfloor x\rfloor)</math> and <math>G(x) = g(\lfloor x\rfloor)</math>, we deduce that
By defining {{math|1=''F''(''x'') = ''f''(⌊''x''⌋)}} and {{math|1=''G''(''x'') = ''g''(⌊''x''⌋)}}, we deduce that


:<math>f(n) = \sum_{1 \le m \le n}\mu(m)g\left(\left\lfloor \frac{n}{m}\right\rfloor\right)\quad\mbox{ for all } n\ge 1.</math>
:<math>f(n) = \sum_{1 \le m \le n}\mu(m)g\left(\left\lfloor \frac{n}{m}\right\rfloor\right)\quad\mbox{ for all } n\ge 1.</math>


A simple example of the use of this formula is counting the number of [[reduced fraction]]s 0 < ''a''/''b'' < 1, where ''a'' and ''b'' are coprime and ''b''≤''n''. If we let ''f''(''n'') be this number, then ''g''(''n'') is the total number of fractions 0 < ''a''/''b'' < 1 with ''b''≤''n'', where ''a'' and ''b'' are not necessarily coprime. (This is because every fraction ''a''/''b'' with gcd(''a'',''b'') = ''d'' and ''b''≤''n'' can be reduced to the fraction (''a''/''d'')/(''b''/''d'') with ''b''/''d'' ≤ ''n''/''d'', and vice versa.) Here it is straightforward to determine ''g''(''n'') = ''n''(''n''-1)/2, but ''f''(''n'') is harder to compute.
A simple example of the use of this formula is counting the number of [[reduced fraction]]s {{math|0 < {{sfrac|''a''|''b''}} < 1}}, where {{mvar|a}} and {{mvar|b}} are coprime and {{math|''b'' ''n''}}. If we let {{math|''f''(''n'')}} be this number, then {{math|''g''(''n'')}} is the total number of fractions {{math|0 < {{sfrac|''a''|''b''}} < 1}} with {{math|''b'' ''n''}}, where {{mvar|a}} and {{mvar|b}} are not necessarily coprime. (This is because every fraction {{math|{{sfrac|''a''|''b''}}}} with {{math|1=gcd(''a'',''b'') = ''d''}} and {{math|''b'' ''n''}} can be reduced to the fraction {{math|{{sfrac|''a''/''d''|''b''/''d''}}}} with {{math|{{sfrac|''b''|''d''}}{{sfrac|''n''|''d''}}}}, and vice versa.) Here it is straightforward to determine {{math|1=''g''(''n'') = {{sfrac|''n''(''n''1)|2}}}}, but {{math|''f''(''n'')}} is harder to compute.


Another inversion formula is (where we assume that the series involved are absolutely convergent):
Another inversion formula is (where we assume that the series involved are absolutely convergent):
Line 117: Line 118:
f(x) = \sum_{m=1}^\infty \mu(m)\frac{g(mx)}{m^s}\quad\mbox{ for all } x\ge 1.</math>
f(x) = \sum_{m=1}^\infty \mu(m)\frac{g(mx)}{m^s}\quad\mbox{ for all } x\ge 1.</math>


As above, this generalises to the case where <math>\alpha(n)</math> is an arithmetic function possessing a Dirichlet inverse <math>\alpha^{-1}(n)</math>:
As above, this generalises to the case where {{math|''α''(''n'')}} is an arithmetic function possessing a Dirichlet inverse {{math|''α''<sup>−1</sup>(''n'')}}:


:<math>g(x) = \sum_{m=1}^\infty \alpha(m)\frac{f(mx)}{m^s}\quad\mbox{ for all } x\ge 1\quad\Longleftrightarrow\quad
:<math>g(x) = \sum_{m=1}^\infty \alpha(m)\frac{f(mx)}{m^s}\quad\mbox{ for all } x\ge 1\quad\Longleftrightarrow\quad
f(x) = \sum_{m=1}^\infty \alpha^{-1}(m)\frac{g(mx)}{m^s}\quad\mbox{ for all } x\ge 1.</math>
f(x) = \sum_{m=1}^\infty \alpha^{-1}(m)\frac{g(mx)}{m^s}\quad\mbox{ for all } x\ge 1.</math>

For example, there is a well known proof relating the [[Riemann zeta function]] to the [[prime zeta function]] that uses the series-based form of
Möbius inversion in the previous equation when <math>s = 1</math>. Namely, by the [[Euler product]] representation of <math>\zeta(s)</math> for
<math>\Re(s) > 1</math>

:<math>\log\zeta(s) = -\sum_{p\mathrm{\ prime}} \log\left(1-\frac{1}{p^s}\right) = \sum_{k \geq 1} \frac{P(ks)}{k} \iff P(s) = \sum_{k \geq 1} \frac{\mu(k)}{k} \log\zeta(ks), \Re(s) > 1.</math>

These identities for alternate forms of Möbius inversion are found in.<ref>NIST Handbook of Mathematical Functions, Section 27.5.</ref>
A more general theory of Möbius inversion formulas partially cited in the next section on incidence algebras is constructed by Rota in.<ref>[On the foundations of combinatorial theory, I. Theory of Möbius Functions|https://link.springer.com/content/pdf/10.1007/BF00531932.pdf]</ref>


==Multiplicative notation==
==Multiplicative notation==
As Möbius inversion applies to any abelian group, it makes no difference whether the group operation is written as addition or as multiplication. This gives rise to the following notational variant of the inversion formula:
As Möbius inversion applies to any abelian group, it makes no difference whether the group operation is written as addition or as multiplication. This gives rise to the following notational variant of the inversion formula:


:<math>\mbox{if } F(n) = \prod_{d|n} f(d),\mbox{ then } f(n) = \prod_{d|n} F\left(\frac{n}{d}\right)^{\mu(d)}.</math>
: <math>
\mbox{If } F(n) = \prod_{d|n} f(d),\mbox{ then } f(n) = \prod_{d|n} F(n/d)^{\mu(d)}. \,
</math>


==Proofs of generalizations==
==Proofs of generalizations==


The first generalization can be proved as follows. We use Iverson's convention that [condition] is the indicator function of the condition, being 1 if the condition is true and 0 if false. We use the result that <math>\sum_{d|n}\mu(d)=i(n)</math>, that is, 1*μ=''i''.
The first generalization can be proved as follows. We use [[Iverson's convention]] that [condition] is the indicator function of the condition, being 1 if the condition is true and 0 if false. We use the result that
:<math>\sum_{d|n}\mu(d)=\varepsilon (n),</math>
that is, <math> 1 * \mu = \varepsilon</math>, where <math>\varepsilon</math> is the [[unit function]].


We have the following:
We have the following:
:<math>\begin{align}

<math>\begin{align}
\sum_{1\le n\le x}\mu(n)g\left(\frac{x}{n}\right)
\sum_{1\le n\le x}\mu(n)g\left(\frac{x}{n}\right)
&= \sum_{1\le n\le x} \mu(n) \sum_{1\le m\le x/n} f\left(\frac{x}{mn}\right)\\
&= \sum_{1\le n\le x} \mu(n) \sum_{1\le m\le \frac{x}{n}} f\left(\frac{x}{mn}\right)\\
&= \sum_{1\le n\le x} \mu(n) \sum_{1\le m\le x/n} \sum_{1\le r\le x} [r=mn] f\left(\frac{x}{r}\right)\\
&= \sum_{1\le n\le x} \mu(n) \sum_{1\le m\le \frac{x}{n}} \sum_{1\le r\le x} [r=mn] f\left(\frac{x}{r}\right)\\
&= \sum_{1\le r\le x} f\left(\frac{x}{r}\right) \sum_{1\le n\le x} \mu(n) \sum_{1\le m\le x/n} [m=r/n] \qquad\text{rearranging the summation order}\\
&= \sum_{1\le r\le x} f\left(\frac{x}{r}\right) \sum_{1\le n\le x} \mu(n) \sum_{1\le m\le \frac{x}{n}} \left[m=\frac{r}{n}\right] \qquad\text{rearranging the summation order}\\
&= \sum_{1\le r\le x} f\left(\frac{x}{r}\right) \sum_{n|r} \mu(n) \\
&= \sum_{1\le r\le x} f\left(\frac{x}{r}\right) \sum_{n|r} \mu(n) \\
&= \sum_{1\le r\le x} f\left(\frac{x}{r}\right) i(r) \\
&= \sum_{1\le r\le x} f\left(\frac{x}{r}\right) \varepsilon (r) \\
&= f(x) \qquad\text{since }i(r)=0\text{ except when }r=1
&= f(x) \qquad\text{since } \varepsilon (r)=0\text{ except when }r=1
\end{align}</math>
\end{align}</math>


The proof in the more general case where α(''n'') replaces 1 is essentially identical, as is the second generalisation.
The proof in the more general case where {{math|''α''(''n'')}} replaces 1 is essentially identical, as is the second generalisation.

==On posets==
{{See also|Incidence algebra}}
For a [[Partially ordered set|poset]] {{mvar|P}}, a set endowed with a partial order relation <math>\leq</math>, define the Möbius function <math>\mu</math> of {{mvar|P}} recursively by

:<math>\mu(s,s) = 1 \text{ for } s \in P, \qquad \mu(s,u) = - \sum_{s \leq t < u} \mu(s,t), \quad \text{ for } s < u \text{ in } P.</math>

(Here one assumes the summations are finite.) Then for <math>f,g: P \to K</math>, where {{mvar|K}} is a commutative ring, we have

:<math>g(t) = \sum_{s \leq t} f(s) \qquad \text{ for all } t \in P</math>

if and only if

:<math>f(t) = \sum_{s \leq t} g(s)\mu(s,t) \qquad \text{ for all }t \in P.</math>

(See Stanley's ''Enumerative Combinatorics'', Vol 1, Section 3.7.) The classical arithmetic Mobius function is the special case of the poset ''P'' of positive integers ordered by [[Divisor|divisibility]]: that is, for positive integers ''s, t,'' we define the partial order <math>s \preccurlyeq t </math> to mean that ''s'' is a divisor of ''t''.


==Contributions of Weisner, Hall, and Rota==
==Contributions of Weisner, Hall, and Rota==
{{Quotation|
{{Quotation|
The statement of the general Möbius inversion formula was first given independently by [[Louis Weisner|Weisner]] (1935) and [[Philip Hall]] (1936); both authors were motivated by group theory problems. Neither author seems to have been aware of the combinatorial implications of his work and neither developed the theory of Möbius functions. In a fundamental paper on Möbius functions, [[Gian-Carlo Rota|Rota]] showed the importance of this theory in combinatorial mathematics and gave a deep treatment of it. He noted the relation between such topics as inclusion-exclusion, classical number theoretic Möbius inversion, coloring problems and flows in networks. Since then, under the strong influence of Rota, the theory of Möbius inversion and related topics has become an active area of combinatorics.<ref>{{cite journal|author=Bender, Edward A.|author2=Goldman, J. R.|title=On the applications of Mö inversion in combinatorial analysis|journal=Amer. Math. Monthly|volume=82|year=1975|pages=789–803|url=http://www.maa.org/programs/maa-awards/writing-awards/on-the-applications-of-m-bius-inversion-in-combinatorial-analysis|doi=10.2307/2319793}}</ref>
The statement of the general Möbius inversion formula [for partially ordered sets] was first given independently by [[Louis Weisner|Weisner]] (1935) and [[Philip Hall]] (1936); both authors were motivated by group theory problems. Neither author seems to have been aware of the combinatorial implications of his work and neither developed the theory of Möbius functions. In a fundamental paper on Möbius functions, [[Gian-Carlo Rota|Rota]] showed the importance of this theory in combinatorial mathematics and gave a deep treatment of it. He noted the relation between such topics as inclusion-exclusion, classical number theoretic Möbius inversion, coloring problems and flows in networks. Since then, under the strong influence of Rota, the theory of Möbius inversion and related topics has become an active area of combinatorics.<ref>{{Harvnb|Bender|Goldman|1975|pp=789–803}}</ref>
|sign=|source=}}
}}


==See also==
==See also==
* [[Farey sequence]]
*[[Farey sequence]]
*[[Inclusion–exclusion principle]]

==Notes==
{{reflist|2}}


==References==
==References==
* {{Apostol IANT}}
* {{Apostol IANT}}
* {{citation|last1=Bender|first1=Edward A.|last2=Goldman|first2=J.&nbsp;R.|title=On the applications of Möbius inversion in combinatorial analysis|journal=Amer. Math. Monthly|volume=82|year=1975|issue=8|pages=789–803|url=http://www.maa.org/programs/maa-awards/writing-awards/on-the-applications-of-m-bius-inversion-in-combinatorial-analysis|doi=10.2307/2319793|jstor=2319793}}
* {{citation|first1=K.|last1=Ireland|first2=M.|last2=Rosen|title=A Classical Introduction to Modern Number Theory|date=2010|series=Graduate Texts in Mathematics (Book 84)|edition=2nd|publisher=Springer-Verlag|isbn=978-1-4419-3094-1}}
* {{SpringerEOM|id=M/m130180 |title=Möbius inversion |first=Joseph P.S. |last=Kung}}
* {{SpringerEOM|id=M/m130180 |title=Möbius inversion |first=Joseph P.S. |last=Kung}}
*{{Citation |last=Möbius |first=A. F. |author-link=August Ferdinand Möbius |year=1832 |title=Über eine besondere Art von Umkehrung der Reihen. |journal=[[Crelle's Journal|Journal für die reine und angewandte Mathematik]] |volume=9 |pages=105–123 |url=https://www.digizeitschriften.de/en/dms/img/?PID=GDZPPN002138654 }}
*K. Ireland, M. Rosen. ''A Classical Introduction to Modern Number Theory'', (1990) Springer-Verlag.
*{{Citation |last=Stanley |first=Richard P.|year=1997 |url=http://www-math.mit.edu/~rstan/ec/ |title=Enumerative Combinatorics |volume=1 |publisher=Cambridge University Press |isbn=0-521-55309-1}}

*{{Citation |last=Stanley |first=Richard P.|year=1999 |url=http://www-math.mit.edu/~rstan/ec/ |title=Enumerative Combinatorics |volume=2 |publisher=Cambridge University Press |isbn=0-521-56069-1}}
{{reflist}}


==External links==
==External links==
{{ProofWiki|id=Möbius_Inversion_Formula|title=Möbius Inversion Formula}}
*{{MathWorld|MoebiusTransform|Möbius Transform}}
*{{MathWorld|MoebiusTransform|Möbius Transform}}



Latest revision as of 21:51, 1 December 2024

In mathematics, the classic Möbius inversion formula is a relation between pairs of arithmetic functions, each defined from the other by sums over divisors. It was introduced into number theory in 1832 by August Ferdinand Möbius.[1]

A large generalization of this formula applies to summation over an arbitrary locally finite partially ordered set, with Möbius' classical formula applying to the set of the natural numbers ordered by divisibility: see incidence algebra.

Statement of the formula

[edit]

The classic version states that if g and f are arithmetic functions satisfying

then

where μ is the Möbius function and the sums extend over all positive divisors d of n (indicated by in the above formulae). In effect, the original f(n) can be determined given g(n) by using the inversion formula. The two sequences are said to be Möbius transforms of each other.

The formula is also correct if f and g are functions from the positive integers into some abelian group (viewed as a Z-module).

In the language of Dirichlet convolutions, the first formula may be written as

where denotes the Dirichlet convolution, and 1 is the constant function 1(n) = 1. The second formula is then written as

Many specific examples are given in the article on multiplicative functions.

The theorem follows because is (commutative and) associative, and 1μ = ε, where ε is the identity function for the Dirichlet convolution, taking values ε(1) = 1, ε(n) = 0 for all n > 1. Thus

.

Replacing by , we obtain the product version of the Möbius inversion formula:

Series relations

[edit]

Let

so that

is its transform. The transforms are related by means of series: the Lambert series

and the Dirichlet series:

where ζ(s) is the Riemann zeta function.

Repeated transformations

[edit]

Given an arithmetic function, one can generate a bi-infinite sequence of other arithmetic functions by repeatedly applying the first summation.

For example, if one starts with Euler's totient function φ, and repeatedly applies the transformation process, one obtains:

  1. φ the totient function
  2. φ1 = I, where I(n) = n is the identity function
  3. I1 = σ1 = σ, the divisor function

If the starting function is the Möbius function itself, the list of functions is:

  1. μ, the Möbius function
  2. μ1 = ε where is the unit function
  3. ε1 = 1, the constant function
  4. 11 = σ0 = d = τ, where d = τ is the number of divisors of n, (see divisor function).

Both of these lists of functions extend infinitely in both directions. The Möbius inversion formula enables these lists to be traversed backwards.

As an example the sequence starting with φ is:

The generated sequences can perhaps be more easily understood by considering the corresponding Dirichlet series: each repeated application of the transform corresponds to multiplication by the Riemann zeta function.

Generalizations

[edit]

A related inversion formula more useful in combinatorics is as follows: suppose F(x) and G(x) are complex-valued functions defined on the interval [1, ∞) such that

then

Here the sums extend over all positive integers n which are less than or equal to x.

This in turn is a special case of a more general form. If α(n) is an arithmetic function possessing a Dirichlet inverse α−1(n), then if one defines

then

The previous formula arises in the special case of the constant function α(n) = 1, whose Dirichlet inverse is α−1(n) = μ(n).

A particular application of the first of these extensions arises if we have (complex-valued) functions f(n) and g(n) defined on the positive integers, with

By defining F(x) = f(⌊x⌋) and G(x) = g(⌊x⌋), we deduce that

A simple example of the use of this formula is counting the number of reduced fractions 0 < a/b < 1, where a and b are coprime and bn. If we let f(n) be this number, then g(n) is the total number of fractions 0 < a/b < 1 with bn, where a and b are not necessarily coprime. (This is because every fraction a/b with gcd(a,b) = d and bn can be reduced to the fraction a/d/b/d with b/dn/d, and vice versa.) Here it is straightforward to determine g(n) = n(n − 1)/2, but f(n) is harder to compute.

Another inversion formula is (where we assume that the series involved are absolutely convergent):

As above, this generalises to the case where α(n) is an arithmetic function possessing a Dirichlet inverse α−1(n):

For example, there is a well known proof relating the Riemann zeta function to the prime zeta function that uses the series-based form of Möbius inversion in the previous equation when . Namely, by the Euler product representation of for

These identities for alternate forms of Möbius inversion are found in.[2] A more general theory of Möbius inversion formulas partially cited in the next section on incidence algebras is constructed by Rota in.[3]

Multiplicative notation

[edit]

As Möbius inversion applies to any abelian group, it makes no difference whether the group operation is written as addition or as multiplication. This gives rise to the following notational variant of the inversion formula:

Proofs of generalizations

[edit]

The first generalization can be proved as follows. We use Iverson's convention that [condition] is the indicator function of the condition, being 1 if the condition is true and 0 if false. We use the result that

that is, , where is the unit function.

We have the following:

The proof in the more general case where α(n) replaces 1 is essentially identical, as is the second generalisation.

On posets

[edit]

For a poset P, a set endowed with a partial order relation , define the Möbius function of P recursively by

(Here one assumes the summations are finite.) Then for , where K is a commutative ring, we have

if and only if

(See Stanley's Enumerative Combinatorics, Vol 1, Section 3.7.) The classical arithmetic Mobius function is the special case of the poset P of positive integers ordered by divisibility: that is, for positive integers s, t, we define the partial order to mean that s is a divisor of t.

Contributions of Weisner, Hall, and Rota

[edit]

The statement of the general Möbius inversion formula [for partially ordered sets] was first given independently by Weisner (1935) and Philip Hall (1936); both authors were motivated by group theory problems. Neither author seems to have been aware of the combinatorial implications of his work and neither developed the theory of Möbius functions. In a fundamental paper on Möbius functions, Rota showed the importance of this theory in combinatorial mathematics and gave a deep treatment of it. He noted the relation between such topics as inclusion-exclusion, classical number theoretic Möbius inversion, coloring problems and flows in networks. Since then, under the strong influence of Rota, the theory of Möbius inversion and related topics has become an active area of combinatorics.[4]

See also

[edit]

Notes

[edit]
  1. ^ Möbius 1832, pp. 105–123
  2. ^ NIST Handbook of Mathematical Functions, Section 27.5.
  3. ^ [On the foundations of combinatorial theory, I. Theory of Möbius Functions|https://link.springer.com/content/pdf/10.1007/BF00531932.pdf]
  4. ^ Bender & Goldman 1975, pp. 789–803

References

[edit]
  • Apostol, Tom M. (1976), Introduction to analytic number theory, Undergraduate Texts in Mathematics, New York-Heidelberg: Springer-Verlag, ISBN 978-0-387-90163-3, MR 0434929, Zbl 0335.10001
  • Bender, Edward A.; Goldman, J. R. (1975), "On the applications of Möbius inversion in combinatorial analysis", Amer. Math. Monthly, 82 (8): 789–803, doi:10.2307/2319793, JSTOR 2319793
  • Ireland, K.; Rosen, M. (2010), A Classical Introduction to Modern Number Theory, Graduate Texts in Mathematics (Book 84) (2nd ed.), Springer-Verlag, ISBN 978-1-4419-3094-1
  • Kung, Joseph P.S. (2001) [1994], "Möbius inversion", Encyclopedia of Mathematics, EMS Press
  • Möbius, A. F. (1832), "Über eine besondere Art von Umkehrung der Reihen.", Journal für die reine und angewandte Mathematik, 9: 105–123
  • Stanley, Richard P. (1997), Enumerative Combinatorics, vol. 1, Cambridge University Press, ISBN 0-521-55309-1
  • Stanley, Richard P. (1999), Enumerative Combinatorics, vol. 2, Cambridge University Press, ISBN 0-521-56069-1
[edit]