Rigorous methods for computational number theory
Koen de Boer, Alice Pellet-Mary, Benjamin Wesolowski
TL;DR
The paper presents the first rigorous algorithm for computing class groups and unit groups of arbitrary number fields in probabilistic subexponential time under ERH, by introducing a general, provable ideal-sampling strategy that handles smooth and near-prime families. It blends Arakelov ray class group theory, lattice basis reduction (BKZ), and discrete Gaussian sampling to convert heuristic sampling steps into rigorously analyzed processes, and then applies it to obtain a full computation of class groups and unit groups, with complexity bounds tied to L-functions and the Dedekind zeta residue. A two-part structure is developed: Part I proves a general sampling theorem ensuring densities drive success rates, while Part II provides the provable lattice-based machinery (log-S-unit lattices, BKZ variants, and post-processing) that yields the final rigorous algorithm. The work thus transforms long-standing heuristic subexponential class-group computations into a framework with explicit probabilistic guarantees, contingent on ERH, and with broad implications for related problems in computational number theory like principal ideals and discrete logarithms.
Abstract
We present the first algorithm for computing class groups and unit groups of arbitrary number fields that provably runs in probabilistic subexponential time, assuming the Extended Riemann Hypothesis (ERH). Previous subexponential algorithms were either restricted to imaginary quadratic fields, or relied on several heuristic assumptions that have long resisted rigorous analysis. The heart of our method is a new general strategy to provably solve a recurring computational problem in number theory (assuming ERH): given an ideal class $[\mathfrak{a}]$ of a number field $K$, sample an ideal $\mathfrak b \in [\mathfrak{a}]$ belonging to a particular family of ideals (e.g., the family of smooth ideals, or near-prime ideals). More precisely, let $\mathcal{S}$ be an arbitrary family of ideals, and $\mathcal{S}_B$ the family of $B$-smooth ideals. We describe an efficient algorithm that samples ideals $\mathfrak b \in [\mathfrak{a}]$ such that $\mathfrak b \in \mathcal{S} \cdot\mathcal{S}_B$ with probability proportional to the density of $\mathcal{S}$ within the set of all ideals. The case where $\mathcal{S}$ is the set of prime ideals yields the family $\mathcal{S}\cdot\mathcal{S}_B$ of near-prime ideals, of particular interest in that it constitutes a dense family of efficiently factorable ideals. The case of smooth ideals $\mathcal{S} = \mathcal{S}_B$ regularly comes up in index-calculus algorithms (notably to compute class groups and unit groups), where it has long constituted a theoretical obstacle overcome only by heuristic arguments.
