On solutions of $\sum_{i=1}^n 1/x_i = 1$ in integers of the form $2^a k^b$
Joel Louwsma
TL;DR
This work addresses representing $1$ as a sum of unit fractions with restricted denominators $x_i \in S(k)=\{2^a k^b : a\in\{0,1,2\}, b\ge0\}$ for $k$ not a power of $2$. It introduces a count-based array encoding $c_{a,b}$ and a set of reduction/expansion moves that manipulate these arrays to generate all such solutions, with reductions guaranteed to decrease the term count. The main contributions are a complete algorithmic framework to enumerate all solutions for given $k$ and $n$, and a precise criterion determining when nontrivial solutions exist, depending on $k$ modulo 4, plus explicit constructive cases. Computational results from a Python implementation illustrate counts for small $k$ and $n$, and the method yields a practical tool for exploring structured Egyptian fraction decompositions and their combinatorial structure.
Abstract
We give an algorithm that produces all solutions of the equation $\sum_{i=1}^n 1/x_i = 1$ in integers of the form $2^a k^b$, where $k$ is a fixed positive integer that is not a power of $2$, $a$ is an element of $\{0,1,2\}$ that can vary from term to term, and $b$ is a nonnegative integer that can vary from term to term. We also completely characterize the pairs $(k,n)$ for which this equation has a nontrivial solution in integers of this form.
