\sigma(n) = (1 + p_1 + p_1^2 + \cdots + p_1^{a_1})(1 + p_2 + \cdots + p_2^{a_2}) \cdots (1 + p_k + \cdots + p_k^{a_k})