Solutions to $\sum_{i=1}^n 1/x_i=1$ in integers $p^a\,q^b$ with $p$ and $q$ two set primes
Claire I. Levaillant
TL;DR
The paper studies representations of the unit as a sum of $n$ unit fractions with denominators $x_i$ of the form $p^a q^b$, where $p$ and $q$ are distinct primes and each appears at least once. It develops a structured, algorithmic framework based on a two-dimensional $p$–$q$ exponent table, admissibility criteria for the top row, and push/left/right moves to enumerate all such solutions of length $n$ with bounded $p$-valuation, including the special case $p=2$. Key results include a fundamental admissibility bound $q < p^{eta} n$, a constructive existence result for $p=2$ with large $n$, and corollaries that bound the admissible primes $q$ for fixed $n$ and width; finally, the algorithm is fully described to generate all solutions while controlling for duplicates. These contributions extend prior work by handling arbitrary two-prime denominators with unbounded exponents and provide a practical enumeration method for Egyptian-fraction decompositions under the stated restrictions.
Abstract
We present an algorithm for computing all the solutions in not necessarily distinct integers to the decomposition of the unit into a sum of unit fractions with denominators $p^a.q^b$ where $p$ and $q$ are two distinct primes, each appearing at least once in the solution.
