P1445 [Violet]樱花

传送门

看到题目就要开始愉快地推式子

原式 $\frac{1}{x}+\frac{1}{y}=\frac{1}{n!}$

$\rightarrow \frac{x+y}{xy}=\frac{1}{n!} \rightarrow (x+y)n!=xy \rightarrow xy-(x+y)n!=0$

两边同时加上 $(n!)^2$ 得

$xy-(x+y)n!+(n!)^2=(n!)^2\rightarrow (x-n!)(y-n!)=(n!)^2$

设$a=(x-n!),b=(y-n!)$，则原式化为 $ab=(n!)^2$

如果 a 确定了，那么 b 也确定了，那么 x,y 也都确定了

所以就变成了求 a 的取值方案数

设$z=(n!)^2$

运用唯一分解定理把$z$分解成几个质数的乘积

$z=p_1^{k_1}p_2^{k_2}...p_n^{k_n}$

a可以在每个质数中取 0 到 k1 个作为质因子，那么每个质数有 $k_i+1$ 种取法

根据乘法原理，共有 $(k_1+1)(k_2+1)...(k_n+1)$ 种取法

所以答案就是这个

把阶乘分解就只要把 1 到 n 每个数都分解一遍，复杂度O(nlogn)

分解一个数的质因子都会吧..就不讲了

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<iostream>
using namespace std;
typedef long long ll;
inline int read()
{
    int x=,f=; char ch=getchar();
    while(ch<''||ch>'') { if(ch=='-') f=-; ch=getchar(); }
    while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
    return x*f;
}
const int N=1e6+,mo=1e9+;
int n,cnt[N],ans=;//cnt存 (n!)^2 的质因子数量
int pri[N],tot,p[N];//p存每个数最小的质因子
bool not_pri[N];
void pre()
{
    not_pri[]=; p[]=;
    for(int i=;i<=n;i++)
    {
        if(!not_pri[i]) { pri[++tot]=i; p[i]=i; }
        for(int j=;j<=tot;j++)
        {
            ll g=pri[j]*i; if(g>n) break;
            not_pri[g]=; p[g]=pri[j];
            if(!(i%pri[j])) break;
        }
    }
}
int main()
{
    n=read();
    pre();
    for(int i=;i<=n;i++)
    {
        int t=i;
        while(t!=) cnt[p[t]]++,t/=p[t];//分解质因子
    }
    for(int i=;i<=tot;i++) ans=1ll*ans*(cnt[pri[i]]*+)%mo;//记得cnt要乘2加1
    printf("%d",ans);
    return ;
}