洛谷P1445：没占到1444的愤怒 （数学推导+线性筛）

题目传送门：https://www.luogu.org/problem/show?pid=1445

---

题目分析：这题我先自己推导了一遍，用了一种很麻烦的推法，后来发现题解的做法比我的快好多。

我自己的推导：
求：$\frac{1}{x}+\frac{1}{y}=\frac{1}{n!}$
$\frac{xy}{x+y}=n!$
不妨令$k=gcd(x,y),x=ak,y=bk$：
原式化为：$\frac{abk}{a+b}=n!$
由于$a,b$互质，$ab$与$a+b$也互质，不妨令$k=u(a+b)$，则有：
$ab=\frac{n!}{u}$
而最开始的式子中：$x=a(a+b)u,y=b(a+b)u$
令$\frac{n!}{u}=p_1^{k_1}p_2^{k_2}……p_q^{k_q}$，则对于$a$而言，它要么选某个质因子的最高次幂，要么不选（因为$gcd(a,b)=1$）。我们不妨按顺序考虑$n!$的每个质因子，设$f[i]$表示考虑完$n!$的前$i$个质因子后的方案数。则对于第$i$个质因子，如果$u$不选完$p_i$的最高次幂，剩余的$p_i$因子就要交给$a$或$b$其中一个去选，于是方案数是$2*k_i*f[i-1]$；如果$p_i^{k_i}|u$，那么方案数就是$f[i-1]$。综上：$f[i]=(2k_i+1)f[i-1]$。

还有题解的快速简便做法，比我的高到不知哪里去了（大概就是因式分解）。话说我初中奥数白学了

虽然推导方法不同，但最后代码都是一样的。用线性筛搞出所有素数，看看它在$n!$中出现几次（比如要求2在$n!$中的出现次数，就是$\lfloor\frac{n}{2}\rfloor+\lfloor\frac{n}{4}\rfloor+\lfloor\frac{n}{8}\rfloor+\lfloor\frac{n}{16}\rfloor……$），然后一通乱搞。

---

CODE：

#include<iostream>
#include<string>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<stdio.h>
#include<algorithm>
using namespace std;

const int maxn=1000100;
const long long M=1000000007;
typedef long long LL;

bool vis[maxn];
int prime[maxn];
int cur=0;

int num[maxn];
LL f[maxn];

int n;

void Linear_shaker()
{
    for (int i=2; i<=n; i++)
    {
        if (!vis[i]) prime[++cur]=i;
        for (int j=1; j<=cur && i*prime[j]<=n; j++)
        {
            vis[ i*prime[j] ]=true;
            if ( !(i%prime[j]) ) break;
        }
    }
}

int main()
{
    freopen("1445.in","r",stdin);
    freopen("1445.out","w",stdout);

    scanf("%d",&n);
    Linear_shaker();

    for (int i=1; i<=cur; i++)
    {
        LL x=prime[i];
        while ( x<=(long long)n ) num[i]+=(n/x),x*=(long long)prime[i];
    }

    f[1]=(num[1]<<1)|1;
    for (int i=2; i<=cur; i++) f[i]=f[i-1]*(long long)((num[i]<<1)|1)%M;
    printf("%lld\n",f[cur]);

    return 0;
}