1. Introduction : Le problème inverse fondamental en mathématiques

En théorie des nombres, en combinatoire, en algèbre et en physique statistique, nous rencontrons sans cesse le même problème inverse : la quantité pure et primitive $f(n)$ que nous cherchons à déterminer est très difficile à compter directement, tandis que la somme cumulée $g(n)$, qui mélange les contributions de tous les éléments inférieurs d’une hiérarchie, se calcule avec une facilité déconcertante.

Considérons trois exemples fondamentaux :

1. Mots périodiques et mots de Lyndon apériodiques : Sur un alphabet à $q$ lettres, il existe $g(n) = q^n$ mots de longueur $n$. Mais ce total inclut les mots périodiques obtenus en répétant un motif plus court de période minimale $d \mid n$. Si $f(d)$ désigne le nombre de mots véritablement apériodiques de période minimale $d$, alors : $$g(n) = \sum_{d \mid n} f(d)$$
2. Fractions réductibles et irréductibles (Indicatrice d’Euler $\varphi(n)$) : En réduisant les $n$ fractions $\frac{1}{n}, \frac{2}{n}, \dots, \frac{n}{n}$ sous forme irréductible, chaque dénominateur devient un diviseur $d \mid n$. Comme il y a exactement $\varphi(d)$ fractions irréductibles de dénominateur $d$, on obtient l’identité de Gauss : $$n = \sum_{d \mid n} \varphi(d)$$
3. Sommes sur les sous-ensembles et principe d’inclusion-exclusion (PIE) : Si chaque sous-ensemble $A \subseteq S$ possède un poids intrinsèque $f(A)$ mais que seule la somme cumulée $g(B) = \sum_{A \subseteq B} f(A)$ est observable, comment isoler $f(B)$ ?

De manière remarquable, bien que $3$ et $1$ divisent $12$, les termes $g(3)$ et $g(1)$ ont un coefficient égal à $0$ (car les quotients $12/3 = 4 = 2^2$ et $12/1 = 12 = 2^2 \cdot 3$ contiennent un facteur premier au carré) !

En 1832, le mathématicien allemand August Ferdinand Möbius a résolu ce problème en introduisant la fonction de Möbius $\mu(n)$. En 1964, Gian-Carlo Rota a généralisé cette théorie aux ensembles partiellement ordonnés (posets) et à leurs algèbres d’incidence, démontrant que l’inversion arithmétique de Möbius, le principe d’inclusion-exclusion, le théorème fondamental du calcul discret et l’inversion moments-cumulants en probabilités ne sont que des cas particuliers de l’équation universelle $\mu = \zeta^{-1}$.

  flowchart TD
    A["Grandeur primitive pure f(n)<br/>(Difficile à observer ou compter directement)"] -->|"Somme cumulée sur les diviseurs<br/>(Convolution de Dirichlet avec 1 : g = f * 1)"| B["Somme totale observable g(n) = Σ_{d|n} f(d)<br/>(Facile à calculer, contient des redondances)"]
    B -->|"Formule d'inversion de Möbius<br/>(Convolution avec l'inverse μ = 1⁻¹ : f = g * μ)"| C["Reconstruction exacte de f(n)<br/>f(n) = Σ_{d|n} μ(d) g(n/d)"]
    D["Fonction de Möbius μ(n) ∈ {-1, 0, 1}<br/>0 si p²|n ; (-1)^k pour k facteurs premiers distincts"] -.->|"Élimine exactement les doubles comptages"| C

2. La fonction de Möbius $\mu(n)$ : Définition, multiplicativité et identité fondamentale

2.1 Définition rigoureuse

$$\mu(n) = \begin{cases} 1 & \text{si } n = 1, \\[4pt] (-1)^k & \text{si } n = p_1 p_2 \cdots p_k \text{ est le produit de } k \text{ nombres premiers distincts (sans facteur carré)}, \\[4pt] 0 & \text{s'il existe un nombre premier } p \text{ tel que } p^2 \mid n. \end{cases}$$

2.2 Multiplicativité et identité fondamentale

1. Multiplicativité : Si $\gcd(a, b) = 1$, alors $\mu(ab) = \mu(a)\mu(b)$. Démonstration : Si $a$ ou $b$ est divisible par un carré $p^2$, alors $ab$ l’est aussi et les deux membres valent $0$. S’ils sont tous deux sans facteur carré avec $r$ et $s$ facteurs premiers, l’hypothèse $\gcd(a, b) = 1$ garantit que leurs ensembles de facteurs premiers sont disjoints, d’où $\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. L’anneau commutatif des fonctions arithmétiques sous la convolution de Dirichlet

Sur l’espace $\mathcal{A} = \{f : \mathbb{N}^+ \to \mathbb{C}\}$ des fonctions arithmétiques, la convolution de Dirichlet est définie par :

$(\mathcal{A}, +, *)$ est un anneau commutatif unitaire dont l’élément neutre multiplicatif est $\varepsilon(n) = [n = 1]$ ($f * \varepsilon = f$). En notant $\mathbf{1}(n) = 1$ la fonction constante égale à $1$, la sommation sur les diviseurs s’écrit $g(n) = \sum_{d \mid n} f(d) = (f * \mathbf{1})(n)$, soit $g = f * \mathbf{1}$. L’identité fondamentale $\sum_{d \mid n} \mu(d) = \varepsilon(n)$ signifie exactement que :

La fonction de Möbius $\mu$ est l’inverse exact de la fonction constante $\mathbf{1}$ pour la convolution de Dirichlet ! La formule d’inversion de Möbius classique se démontre alors en une seule ligne :

Séries de Dirichlet, fonction zêta de Riemann $\zeta(s)$ et hypothèse de Riemann

Le passage aux séries de Dirichlet $D_f(s) = \sum_{n=1}^{\infty} \frac{f(n)}{n^s}$ transforme la convolution en produit ordinaire : $D_{f * g}(s) = D_f(s) D_g(s)$. Comme $D_{\mathbf{1}}(s) = \zeta(s) = \sum_{n=1}^{\infty} \frac{1}{n^s}$ est la fonction zêta de Riemann et $D_{\varepsilon}(s) = 1$, la relation $\mathbf{1} * \mu = \varepsilon$ donne immédiatement :

• Le théorème des nombres premiers est équivalent à $\sum_{n=1}^{\infty} \frac{\mu(n)}{n} = 0$, soit $M(x) = \sum_{n \le x} \mu(n) = o(x)$.
• L’hypothèse de Riemann (RH) est équivalente au caractère holomorphe de $\frac{1}{\zeta(s)}$ sur $\operatorname{Re}(s) > \frac{1}{2}$, ce qui équivaut exactement à la majoration en racine carrée de la fonction de Mertens : $$\text{Hypothèse de Riemann} \quad \Longleftrightarrow \quad M(x) = \sum_{n \le x} \mu(n) = O\!\left(x^{\frac{1}{2} + \varepsilon}\right) \quad (\forall \varepsilon > 0)$$

4. Trois applications classiques

1. Indicatrice d’Euler $\varphi(n)$ : De $\sum_{d \mid n} \varphi(d) = n$ ($\varphi * \mathbf{1} = \mathrm{id}$), l’inversion de Möbius donne $\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. Colliers apériodiques (mots de Lyndon) et polynômes irréductibles sur $\mathbb{F}_q$ : Le nombre $N_q(d)$ de colliers apériodiques de longueur $d$ sur un alphabet de taille $q$ et le nombre $I_q(d)$ de polynômes irréductibles unitaires de degré $d$ sur le corps fini $\mathbb{F}_q$ (via $x^{q^n} - x = \prod_{d \mid n} \prod_{\deg P = d} P(x)$) vérifient tous deux $q^n = \sum_{d \mid n} d \cdot N_q(d) = \sum_{d \mid n} d \cdot I_q(d)$. L’inversion de Möbius fournit la formule de Gauss : $$N_q(n) = I_q(n) = \frac{1}{n} \sum_{d \mid n} \mu\!\left(\frac{n}{d}\right) q^d$$
3. Polynômes cyclotomiques $\Phi_n(x)$ : À partir de $x^n - 1 = \prod_{d \mid n} \Phi_d(x)$, l’inversion de Möbius multiplicative donne : $$\Phi_n(x) = \prod_{d \mid n} (x^d - 1)^{\mu(n/d)}$$ Pour $n = 12$, avec $d \in \{1, 2, 3, 4, 6, 12\}$ et $\mu(12/d) \in \{0, +1, 0, -1, -1, +1\}$, on obtient : $$\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. La grande unification de Gian-Carlo Rota (1964) : Posets et algèbres d’incidence

Soit $(P, \le)$ un ensemble partiellement ordonné (poset) localement fini (chaque intervalle fermé $[x, y] = \{z \in P \mid x \le z \le y\}$ est fini). L’algèbre d’incidence $I(P)$ est munie du produit de convolution :

L’élément neutre est le delta de Kronecker $\delta(x, y) = [x = y]$. La fonction zêta $\zeta(x, y) = 1$ (pour $x \le y$) admet un unique inverse $\mu = \zeta^{-1}$ défini par récurrence par $\mu(x, x) = 1$ et $\mu(x, y) = -\sum_{x \le z < y} \mu(x, z)$ pour $x < y$. On obtient ainsi le théorème d’inversion de Möbius généralisé sur les posets :

  flowchart TD
    R["Inversion de Möbius généralisée de Rota (1964)<br/>μ = ζ⁻¹ dans l'algèbre d'incidence I(P)<br/>g = f * ζ  ⇔  f = g * μ"]
    R --> P1["1. Chaîne linéaire (ℕ₀, ≤)<br/>0 < 1 < 2 < 3 < ..."]
    R --> P2["2. Treillis booléen des sous-ensembles (𝒫(S), ⊆)<br/>Inclusion A ⊆ B"]
    R --> P3["3. Treillis de divisibilité (ℕ⁺, |)<br/>Divisibilité a | b"]
    R --> P4["4. Treillis des partitions (Π_n, ≤)<br/>Raffinement π ≤ σ"]
    P1 --> T1["μ(i, i) = 1, μ(i, i+1) = -1, sinon 0<br/>⇒ Théorème fondamental du calcul discret<br/>f(n) = g(n) - g(n-1) (Δ = Σ⁻¹)"]
    P2 --> T2["μ(A, B) = (-1)^|B \ A|<br/>⇒ Principe d'inclusion-exclusion (PIE)<br/>& Dérangements D_n"]
    P3 --> T3["μ(a, b) = μ_classique(b/a)<br/>⇒ Inversion de Möbius arithmétique<br/>(Produit direct de chaînes d'exposants premiers !)"]
    P4 --> T4["μ(0̂, 1̂) = (-1)^(n-1) (n-1)!<br/>⇒ Formule d'inversion moments-cumulants<br/>en théorie des probabilités"]