数论 因子和问题

因子和公式 推导

我们设 $f(x)$ 为 x 的因数和， $x$ 为正数
通过唯一分解定理，可以得到：
$x=p_1^{e_1}\times p_2^{e_2}\times ...\times p_k^{e_k}$
设 $x$ 的最小质因子为 $p$, 因子 $p$ 的个数为 $e$，那么 $x=p^e\ast x\prime$

易得：
当 $x$ 的因子中 $p$ 的个数等于 0 的时候， 因子和为 $f(x\prime )$。推广至一般情况，当 $x$ 的因子中 $e$ 的时候，因子和为 $p^e\ast f(x\prime )$，所以 $x$ 的所以因子和可以表示成：
$f(x)=(p^0+p^1+...+p^e)\ast f(x\prime )=\frac{p^{e+1}-1}{p-1}\ast f(x\prime )$

$f(x\prime )$ 可以通过相同方法递归计算。最后可以表示成一系列等比数列和的乘积:
$f(x)=\coprod _{i=1}^k\frac{p_i^{e_i+1}-1}{p_i-1}$

个人板子

#define int long long
int n;
int p[N], num[N], idx;
void div(int x)
{
    for (int i = 2; i * i <= x; i++)
    {
        if (x % i == 0)
        {
            p[++idx] = i;
            while (x % i == 0)
                num[idx]++, x /= i;
        }
    }
    if (x > 1)
        p[++idx] = x, num[idx] = 1;
}
int ksm(int a, int b)
{
    int ans = 1;
    while (b)
    {
        if (b & 1)
            ans = ans * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return ans;
}
void solve()
{
    idx = 0;
    cin >> n;
    if (n == 0)
    {
        cout << 0 << '\n';
        return;
    }
    div(n);
    int ans = 1;
    for (int i = 1; i <= idx; i++)
    {
        int a = p[i], b = num[i] + 1;
        if ((a - 1) % mod)
        {
            int tt = (ksm(a, b) - 1 + mod) % mod;
            tt = tt * ksm(a - 1, mod - 2) % mod;
            ans = ans * tt % mod;
        }
        else
            ans = ans * b % mod;
    }
    cout << ans << '\n';
}
signed main()
{
    buff;
    solve();
}