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 :
- 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)\)
- 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)\)
- 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}\]
| \(n\) | Factorisation première | Sans facteur carré ? | Nombre \(k\) de facteurs premiers | Fonction de Möbius \(\mu(n)\) | Fonction de Mertens \(M(n) = \sum_{m=1}^n \mu(m)\) |
|---|---|---|---|---|---|
| \(1\) | \(1\) (produit vide) | Oui | \(0\) | \(+1\) | \(+1\) |
| \(2\) | \(2\) | Oui | \(1\) | \(-1\) | \(0\) |
| \(3\) | \(3\) | Oui | \(1\) | \(-1\) | \(-1\) |
| \(4\) | \(2^2\) | Non (\(2^2 \mid 4\)) | — | \(0\) | \(-1\) |
| \(5\) | \(5\) | Oui | \(1\) | \(-1\) | \(-2\) |
| \(6\) | \(2 \cdot 3\) | Oui | \(2\) | \(+1\) | \(-1\) |
| \(7\) | \(7\) | Oui | \(1\) | \(-1\) | \(-2\) |
| \(8\) | \(2^3\) | Non (\(2^2 \mid 8\)) | — | \(0\) | \(-2\) |
| \(9\) | \(3^2\) | Non (\(3^2 \mid 9\)) | — | \(0\) | \(-2\) |
| \(10\) | \(2 \cdot 5\) | Oui | \(2\) | \(+1\) | \(-1\) |
| \(11\) | \(11\) | Oui | \(1\) | \(-1\) | \(-2\) |
| \(12\) | \(2^2 \cdot 3\) | Non (\(2^2 \mid 12\)) | — | \(0\) | \(-2\) |
2.2 Multiplicativité et identité fondamentale
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\)
- \[\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 :
\[(f * g)(n) := \sum_{d \mid n} f(d)\,g\!\left(\frac{n}{d}\right) = \sum_{ab = n} f(a)\,g(b)\]\((\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 :
\[\mu * \mathbf{1} = \mathbf{1} * \mu = \varepsilon \quad \Longleftrightarrow \quad \mu = \mathbf{1}^{-1}\]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 :
\[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)\]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 :
\[\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)\]- 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
- 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)\)
- 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\)
- 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 :
\[(\alpha * \beta)(x, y) := \sum_{x \le z \le y} \alpha(x, z)\,\beta(z, y) \qquad (x \le y)\]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 :
\[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["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"]
| Poset \((P, \le)\) | Relation d’ordre \(x \le y\) | Fonction de Möbius \(\mu(x, y)\) | Théorème fondamental issu de \(f = g * \mu\) |
|---|---|---|---|
| Chaîne linéaire \((\mathbb{N}_0, \le)\) | Ordre naturel \(i \le j\) | \(1\) (\(j=i\)), \(-1\) (\(j=i+1\)), \(0\) (sinon) | Calcul discret fondamental \(f(n) = g(n) - g(n-1)\) |
| Treillis booléen \((\mathcal{P}(S), \subseteq)\) | Inclusion \(A \subseteq B\) | $(-1)^{ | B \setminus A |
| Treillis de divisibilité \((\mathbb{N}^+, \mid)\) | Divisibilité \(a \mid b\) | \(\mu_{\text{classique}}(b/a) \in \{-1, 0, 1\}\) | Inversion de Möbius classique \(f(n) = \sum_{d \mid n} \mu(n/d) g(d)\) |
| Treillis des partitions \((\Pi_n, \le)\) | Raffinement \(\pi \le \sigma\) | $\mu(\pi, \hat{1}) = (-1)^{ | \pi |