Featured image of post Kann man die ursprünglichen Zahlen zurückgewinnen, wenn nur die Summe bekannt ist? Die Mathematik der Möbius-Inversion

Kann man die ursprünglichen Zahlen zurückgewinnen, wenn nur die Summe bekannt ist? Die Mathematik der Möbius-Inversion

Wenn sich reine, primitive Größen nicht direkt beobachten lassen und nur ihre kumulative Summe über Teiler oder Teilmengen bekannt ist: Wie schält man doppelt gezählte Beiträge systematisch heraus? Eine umfassende Einführung in die Möbius-Funktion, die Dirichlet-Faltung, die Riemannsche Zeta-Funktion, klassische Anwendungen sowie Gian-Carlo Rotas große Vereinheitlichung über Halbordnungen (Posets) und Inzidenzalgebren.

1. Einleitung: Das fundamentale inverse Problem der Mathematik

In der Zahlentheorie, Kombinatorik, Algebra und statistischen Physik begegnen wir immer wieder demselben inversen Problem: Die reine, primitive Größe $f(n)$, die wir eigentlich bestimmen möchten, lässt sich kaum direkt zählen oder messen, während die kumulative Summe $g(n)$, in der alle untergeordneten Beiträge einer Hierarchie zusammenfließen, erstaunlich leicht zu berechnen ist.

Betrachten wir drei klassische Beispiele:

  1. Periodische Zeichenketten und aperiodische Lyndon-Wörter: Über einem Alphabet mit $q$ Zeichen gibt es insgesamt $g(n) = q^n$ Wörter der Länge $n$. Darunter befinden sich jedoch degenerierte Wörter, die durch Wiederholung eines kürzeren Blocks der Periodenlänge $d \mid n$ entstehen. Bezeichnet $f(d)$ die Anzahl der echt aperiodischen Wörter mit minimaler Periode $d$, so gilt die Teilersummen-Relation: $$g(n) = \sum_{d \mid n} f(d)$$
  2. Kürzbare und unkürzbare Brüche (Eulersche $\varphi$-Funktion): Kürzt man die $n$ Brüche $\frac{1}{n}, \frac{2}{n}, \dots, \frac{n}{n}$ vollständig, so wird jeder Nenner zu einem Teiler $d \mid n$. Da es zu jedem Nenner $d$ genau $\varphi(d)$ unkürzbare Brüche gibt, folgt die Gaußsche Identität: $$n = \sum_{d \mid n} \varphi(d)$$
  3. Teilmengensummen und das Prinzip von Inklusion und Exklusion (PIE): Ist für jede Teilmenge $A \subseteq S$ ein Gewicht $f(A)$ gegeben, aber nur die kumulative Summe $g(B) = \sum_{A \subseteq B} f(A)$ beobachtbar, wie gewinnt man $f(B)$ zurück?
$$f(6) = g(6) - g(3) - g(2) + g(1)$$$$f(12) = g(12) - g(6) - g(4) + g(2)$$

Überraschenderweise tauchen die Teiler $d = 3$ und $d = 1$ (für die der Quotient $12/d$ durch das Quadrat $2^2 = 4$ teilbar ist) überhaupt nicht auf — ihr Koeffizient ist $0$!

Im Jahr 1832 löste der deutsche Mathematiker und Astronom August Ferdinand Möbius dieses Problem durch die Einführung der Möbius-Funktion $\mu(n)$. Im Jahr 1964 verallgemeinerte Gian-Carlo Rota dieses Prinzip auf partiell geordnete Mengen (Posets) und deren Inzidenzalgebren und zeigte, dass die zahlentheoretische Möbius-Inversion, das Prinzip von Inklusion und Exklusion, der Hauptsatz der diskreten Analysis und die Kumulantenentwicklung der Wahrscheinlichkeitstheorie allesamt Spezialfälle der Gleichung $\mu = \zeta^{-1}$ sind.

  flowchart TD
    A["Reine primitive Größe f(n)<br/>(Schwer direkt zu beobachten oder zu zählen)"] -->|"Kumulative Summe über alle Teiler<br/>(Dirichlet-Faltung mit 1: g = f * 1)"| B["Beobachtbare Gesamtsumme g(n) = Σ_{d|n} f(d)<br/>(Leicht berechenbar, enthält Überlappungen)"]
    B -->|"Möbius-Inversionsformel<br/>(Faltung mit dem Inversen μ = 1⁻¹: f = g * μ)"| C["Exakte Rekonstruktion von f(n)<br/>f(n) = Σ_{d|n} μ(d) g(n/d)"]
    D["Möbius-Funktion μ(n) ∈ {-1, 0, 1}<br/>0 bei Quadratfaktoren; (-1)^k bei k verschiedenen Primfaktoren"] -.->|"Hebt Mehrfachzählungen<br/>exakt auf"| C

2. Die Möbius-Funktion $\mu(n)$: Definition, Multiplikativität und Fundamentalidentität

2.1 Definition und Wertetabelle

$$\mu(n) = \begin{cases} 1 & \text{falls } n = 1, \\[4pt] (-1)^k & \text{falls } n = p_1 p_2 \cdots p_k \text{ ein Produkt von } k \text{ verschiedenen Primzahlen ist (quadratfrei)}, \\[4pt] 0 & \text{falls } p^2 \mid n \text{ für eine Primzahl } p \text{ gilt (nicht quadratfrei)}. \end{cases}$$
$n$PrimfaktorzerlegungQuadratfrei?Anzahl $k$ versch. PrimfaktorenMöbius-Funktion $\mu(n)$Mertens-Funktion $M(n) = \sum_{m=1}^n \mu(m)$
$1$$1$ (leeres Produkt)Ja$0$$+1$$+1$
$2$$2$Ja$1$$-1$$0$
$3$$3$Ja$1$$-1$$-1$
$4$$2^2$Nein ($2^2 \mid 4$)—$0$$-1$
$5$$5$Ja$1$$-1$$-2$
$6$$2 \cdot 3$Ja$2$$+1$$-1$
$7$$7$Ja$1$$-1$$-2$
$8$$2^3$Nein ($2^2 \mid 8$)—$0$$-2$
$9$$3^2$Nein ($3^2 \mid 9$)—$0$$-2$
$10$$2 \cdot 5$Ja$2$$+1$$-1$
$11$$11$Ja$1$$-1$$-2$
$12$$2^2 \cdot 3$Nein ($2^2 \mid 12$)—$0$$-2$

2.2 Multiplikativität und die Fundamentalidentität

  1. Multiplikativität: Für teilerfremde $a, b \in \mathbb{N}^+$ mit $\gcd(a, b) = 1$ gilt $\mu(ab) = \mu(a)\mu(b)$. Beweis: Besitzt $a$ oder $b$ einen quadratischen Primfaktor $p^2$, so auch $ab$, und beide Seiten sind $0$. Sind beide quadratfrei mit $r$ bzw. $s$ Primfaktoren, so haben sie wegen $\gcd(a, b) = 1$ keine gemeinsamen Primteiler; daher hat $ab$ genau $r+s$ verschiedene Primfaktoren und es gilt $\mu(ab) = (-1)^{r+s} = (-1)^r (-1)^s = \mu(a)\mu(b)$. $\blacksquare$

  2. $$\sum_{d \mid n} \mu(d) = \varepsilon(n) := \begin{cases} 1 & n = 1, \\ 0 & n > 1. \end{cases}$$$$\sum_{d \mid n} \mu(d) = \sum_{j = 0}^{k} \binom{k}{j} (-1)^j 1^{k - j} = (1 - 1)^k = 0. \quad \blacksquare$$

3. Der kommutative Ring arithmetischer Funktionen unter der Dirichlet-Faltung

Auf dem Raum $\mathcal{A} = \{f : \mathbb{N}^+ \to \mathbb{C}\}$ aller zahlentheoretischen Funktionen ist die Dirichlet-Faltung definiert durch:

$$(f * g)(n) := \sum_{d \mid n} f(d)\,g\!\left(\frac{n}{d}\right) = \sum_{ab = n} f(a)\,g(b)$$

$(\mathcal{A}, +, *)$ bildet einen kommutativen Ring mit dem Einselement $\varepsilon(n) = [n = 1]$ ($f * \varepsilon = f$). Mit der konstanten Funktion $\mathbf{1}(n) = 1$ schreibt sich jede Teilersumme als Faltung $g(n) = \sum_{d \mid n} f(d) = (f * \mathbf{1})(n)$, also $g = f * \mathbf{1}$. Die Fundamentalidentität $\sum_{d \mid n} \mu(d) = \varepsilon(n)$ lautet in dieser Sprache:

$$\mu * \mathbf{1} = \mathbf{1} * \mu = \varepsilon \quad \Longleftrightarrow \quad \mu = \mathbf{1}^{-1}$$

Die Möbius-Funktion $\mu$ ist exakt das bezüglich der Dirichlet-Faltung inverse Element der konstanten Funktion $\mathbf{1}$! Damit folgt die klassische Möbius-Inversionsformel in einer einzigen Zeile:

$$g = f * \mathbf{1} \iff g * \mu = f * (\mathbf{1} * \mu) = f * \varepsilon = f \iff f(n) = \sum_{d \mid n} \mu(d)\,g\!\left(\frac{n}{d}\right) = \sum_{d \mid n} \mu\!\left(\frac{n}{d}\right) g(d)$$

Verbindung zur analytischen Zahlentheorie: Dirichlet-Reihen und $\zeta(s)$

Ordnet man jeder Funktion $f \in \mathcal{A}$ ihre Dirichlet-Reihe $D_f(s) = \sum_{n=1}^{\infty} \frac{f(n)}{n^s}$ zu, so geht die Faltung in das gewöhnliche Produkt über: $D_{f * g}(s) = D_f(s) D_g(s)$. Da $D_{\mathbf{1}}(s) = \zeta(s) = \sum_{n=1}^{\infty} \frac{1}{n^s}$ die Riemannsche Zeta-Funktion und $D_{\varepsilon}(s) = 1$ ist, liefert $\mathbf{1} * \mu = \varepsilon$ sofort:

$$\zeta(s) \sum_{n = 1}^{\infty} \frac{\mu(n)}{n^s} = 1 \quad \Longrightarrow \quad \frac{1}{\zeta(s)} = \sum_{n = 1}^{\infty} \frac{\mu(n)}{n^s} \qquad (\operatorname{Re}(s) > 1)$$
  • Der Primzahlsatz ist äquivalent zu $\sum_{n=1}^{\infty} \frac{\mu(n)}{n} = 0$ bzw. $M(x) = \sum_{n \le x} \mu(n) = o(x)$.
  • Die Riemannsche Vermutung (Riemann Hypothesis) ist äquivalent dazu, dass $\frac{1}{\zeta(s)}$ für $\operatorname{Re}(s) > \frac{1}{2}$ holomorph ist, was wiederum exakt der Quadratwurzel-Auslöschung der Mertens-Funktion entspricht: $$\text{Riemannsche Vermutung} \quad \Longleftrightarrow \quad M(x) = \sum_{n \le x} \mu(n) = O\!\left(x^{\frac{1}{2} + \varepsilon}\right) \quad (\forall \varepsilon > 0)$$

4. Drei klassische Anwendungen

  1. Eulersche $\varphi$-Funktion: Aus $\sum_{d \mid n} \varphi(d) = n$ ($\varphi * \mathbf{1} = \mathrm{id}$) folgt durch Möbius-Inversion $\varphi = \mathrm{id} * \mu$: $$\varphi(n) = \sum_{d \mid n} \mu(d)\,\frac{n}{d} = n \prod_{p \mid n} \left(1 - \frac{1}{p}\right)$$
  2. Aperiodische Halsketten (Lyndon-Wörter) und irreduzible Polynome über $\mathbb{F}_q$: Sowohl die Anzahl $N_q(d)$ aperiodischer Halsketten der Länge $d$ über $q$ Symbolen als auch die Anzahl $I_q(d)$ normierter irreduzibler Polynome vom Grad $d$ über dem endlichen Körper $\mathbb{F}_q$ (da $x^{q^n} - x = \prod_{d \mid n} \prod_{\deg P = d} P(x)$) erfüllen $q^n = \sum_{d \mid n} d \cdot N_q(d) = \sum_{d \mid n} d \cdot I_q(d)$. Möbius-Inversion liefert die Gaußsche Formel: $$N_q(n) = I_q(n) = \frac{1}{n} \sum_{d \mid n} \mu\!\left(\frac{n}{d}\right) q^d$$
  3. Kreisteilungspolynome (Zyklotomische Polynome) $\Phi_n(x)$: Aus $x^n - 1 = \prod_{d \mid n} \Phi_d(x)$ folgt durch multiplikative Möbius-Inversion: $$\Phi_n(x) = \prod_{d \mid n} (x^d - 1)^{\mu(n/d)}$$ Für $n = 12$ mit Teilern $d \in \{1, 2, 3, 4, 6, 12\}$ und $\mu(12/d) \in \{0, +1, 0, -1, -1, +1\}$ ergibt sich sofort: $$\Phi_{12}(x) = \frac{(x^{12} - 1)(x^2 - 1)}{(x^6 - 1)(x^4 - 1)} = \frac{x^6 + 1}{x^2 + 1} = x^4 - x^2 + 1$$

5. Gian-Carlo Rotas Große Vereinheitlichung (1964): Halbordnungen (Posets) und Inzidenzalgebren

Sei $(P, \le)$ eine lokal endliche partiell geordnete Menge (Poset), sodass jedes Intervall $[x, y] = \{z \in P \mid x \le z \le y\}$ endlich ist. Die Inzidenzalgebra $I(P)$ ist mit der Intervall-Faltung ausgestattet:

$$(\alpha * \beta)(x, y) := \sum_{x \le z \le y} \alpha(x, z)\,\beta(z, y) \qquad (x \le y)$$

Das Einselement ist das Kronecker-Delta $\delta(x, y) = [x = y]$. Die Zeta-Funktion $\zeta(x, y) = 1$ (für $x \le y$) besitzt ein eindeutiges Inverses $\mu = \zeta^{-1}$, das rekursiv durch $\mu(x, x) = 1$ und $\mu(x, y) = -\sum_{x \le z < y} \mu(x, z)$ für $x < y$ bestimmt ist. Damit gilt der allgemeine Möbius-Inversionssatz auf Posets:

$$g(y) = \sum_{x \le y} f(x) \quad \Longleftrightarrow \quad f(y) = \sum_{x \le y} g(x)\,\mu(x, y)$$
  flowchart TD
    R["Rotas allgemeine Möbius-Inversion (1964)<br/>μ = ζ⁻¹ in der Inzidenzalgebra I(P)<br/>g = f * ζ  ⇔  f = g * μ"]
    R --> P1["1. Lineare Kette (ℕ₀, ≤)<br/>0 < 1 < 2 < 3 < ..."]
    R --> P2["2. Boolescher Teilmengenverband (𝒫(S), ⊆)<br/>Inklusion A ⊆ B"]
    R --> P3["3. Teilbarkeitsverband (ℕ⁺, |)<br/>Teilbarkeit a | b"]
    R --> P4["4. Partitionsverband (Π_n, ≤)<br/>Verfeinerung π ≤ σ"]
    P1 --> T1["μ(i, i) = 1, μ(i, i+1) = -1, sonst 0<br/>⇒ Hauptsatz der diskreten Analysis<br/>f(n) = g(n) - g(n-1) (Δ = Σ⁻¹)"]
    P2 --> T2["μ(A, B) = (-1)^|B \ A|<br/>⇒ Prinzip von Inklusion und Exklusion (PIE)<br/>& Fixpunktfreie Permutationen D_n"]
    P3 --> T3["μ(a, b) = μ_klassisch(b/a)<br/>⇒ Klassische zahlentheoretische Möbius-Inversion<br/>(Direktes Produkt von Primzahlpotenz-Ketten!)"]
    P4 --> T4["μ(0̂, 1̂) = (-1)^(n-1) (n-1)!<br/>⇒ Moment-Kumulanten-Inversion<br/>in der Wahrscheinlichkeitstheorie"]
Halbordnung (Poset) $(P, \le)$Ordnung $x \le y$Möbius-Funktion $\mu(x, y)$Resultierender Fundamentalsatz ($f = g * \mu$)
Lineare Kette $(\mathbb{N}_0, \le)$Größe $i \le j$$1$ ($j=i$), $-1$ ($j=i+1$), $0$ (sonst)Hauptsatz der diskreten Analysis $f(n) = g(n) - g(n-1)$
Boolescher Verband $(\mathcal{P}(S), \subseteq)$Teilmenge $A \subseteq B$$(-1)^{B \setminus A
Teilbarkeitsverband $(\mathbb{N}^+, \mid)$Teilbarkeit $a \mid b$$\mu_{\text{klassisch}}(b/a) \in \{-1, 0, 1\}$Klassische Möbius-Inversion $f(n) = \sum_{d \mid n} \mu(n/d) g(d)$
Partitionsverband $(\Pi_n, \le)$Verfeinerung $\pi \le \sigma$$\mu(\pi, \hat{1}) = (-1)^{\pi

Insbesondere zeigt die Produktstruktur des Teilbarkeitsverbandes als direktes Produkt linearer Primexponent-Ketten sofort, warum $\mu(n) = 0$ für $p^2 \mid n$ gilt: Der diskrete Differenzoperator $\Delta g(m) = g(m) - g(m-1)$ blickt entlang jeder Primfaktor-Achse nur genau einen Schritt zurück und weist Abständen $\ge 2$ das Gewicht $0$ zu!

comments powered by Disqus