x = p1^a1 * p2^a2 * p3^a3…;
s = (1 + p1^1 +…+ p1^a1) * (1 + p2^1 +…+ p2^a2) * … * (1 + pn^1 +…+ pn^an);