X^Y 因子数 因子和

本篇博客讲解一下 x^y 的因子数，因子和的求解问题。

首先根据算术基本定理(唯一分解定理)将 x 进行质因数分解：
x = p₁^c₁ × p₂^c₂ × ··· × p_n^c_n

则 x 的正因数可以表示为：{p₁^b₁ p₂^b₂ ··· p_n^b_n} (0 ≤ b_i ≤ c_i)

每个 b_i 有 c_i + 1 种取法，根据乘法原理可得 x 的正因数个数为(Π为连乘积符号)：
(c₁ + 1) × (c₂ + 1) × ··· × (c_n + 1)   =   Π_i=1ⁿ (c_i + 1)

x 的所有正因数之和为:
(1 + p₁ + p₁² + ··· + p₁^c₁) × (1 + p₂ + p₂² + ··· + p₂^c₂) × ··· × (1 + p_n + p_n² + ··· + p_n^c_n)   =   Π_i=1ⁿ ( ∑_j=0^c_i p_i ^j )
解释如下：
p₁^c₁ 的因数有: p₁⁰, p₁¹, p₁², …, p₁^c₁
… …
p_n^c_n 的因数有: p_n⁰, p_n¹, p_n², …, p_n^c_n
而 x 的因数就是在 p₁^c₁, p₂^c₂, …, p_n^c_n 的因数中各取一个然后相乘得来，再由乘法原理得到上式。

经过上面的论述，同理可得：
x^y = (p₁^c₁ × p₂^c₂ × p₃^c₃ × … × p_n^c_n) ^y

x^y 的正因数个数为:
(yc₁ + 1) × (yc₂ + 1) × ··· × (yc_n + 1)   =   Π_i=1ⁿ (yc_i + 1)

x^y 的所有正因数之和为:
(1 + p₁ + p₁² + ··· + p₁^yc₁) × (1 + p₂ + p₂² + ··· + p₂^yc₂) × ··· × (1 + p_n + p_n² + ··· + p_n^yc_n)   =   Π_i=1ⁿ ( ∑_j=0^yc_i p_i ^j )

至于计算，快速幂 + 等比数列求和公式 ( + 逆元 + 快速乘) 即可。

相关题目会后续更新。

如有错误或不足希望大家能够及时指出，谢谢。