Efficient Computation and Congruences for Colored Partition Functions

Source abstract

For a positive integer $α$, let $p_α(n)$ denote the number of $α$-colored partitions of $n$. The Rademacher-type expansion of Iskander, Jain, and Talvola is valid for every real $α>0$, but when $α>24$ it involves several polar terms and a two-parameter family of exponential sums $A_k^{(α)}(n,m)$. For integral $α$, we prove multiplicativity and prime-power reduction formulas for these sums, expressing their local factors as classical or quadratically twisted Kloosterman sums. Together with explicit truncation and precision bounds, these formulas yield an efficient algorithm for computing $p_α(n)$ exactly. Our SageMath implementation computes the 1,113,767-digit integer $p_{100}(10^{10})$ in less than one hour. As an application, we use the algorithm to certify new Ramanujan-type congruences.