通过分解质因数来求约数个数+阶乘分解

---

1)如果直接通过这个式子很明显是无法求得，我们应该要给这个式子化简：
1.$\frac{x+y}{x\ast y}=\frac{1}{n!}$
2.$(x+y)\ast n!=x\ast y$
3.$x=\frac{y\ast n!}{(y-n!)}$
4.$x=\frac{(y-n!+n!)\ast n!}{(y-n!)}$
5.$x=n!+\frac{n!\ast n!}{(y-n!)}$
6.通过式子5我们知道$n!$一定是整数，如果要$x$是整数那么后面那部分也要是整数说明$(y-n!)|n!\ast n!$
7.那么我们就可以运用质因数分解来求约数个数那个公式$sum=(1+c_1)\ast (1+c_2)+...(1+c_n)$
8.其中$n!=p_1^{c_1}\ast p_2^{c_2}\ast ....\ast p_n^{c_n}$那么$n!\ast n!=p_1^{2c_1}\ast p_2^{2c_2}\ast ....\ast p_n^{2c_n}$
9.那么我们就可以联系上次快速阶乘分解得方法
不懂戳这里

---

#include <bits/stdc++.h>

using namespace std;
const int N = 1e6 + 10, mod = 1e9 + 7;
int prime[N], cnt;
bool vis[N];
int n;
void init(int n)
{
    for(int i = 2; i <= n; ++ i)
    {
        if(!vis[i]) prime[cnt ++] = i;
        for(int j = 0; j < cnt && i * prime[j] <= n;  ++ j)
        {
            vis[i * prime[j]] = true;
            if(i % prime[j] == 0) break;
        }
    }
}



int main()
{
    cin >> n;
    init(n);
    int res = 1;
    for(int i = 0; i < cnt; ++ i)
    {
        int p = prime[i];
        int s = 0;
        for(int j = n; j; j /= p) s += j / p;
        res = 1ll * res * (2 * s + 1) % mod;
    }
    cout <<res;
    return 0;
}