BZOJ 2721 [Violet 5]樱花

推式子+质因数分解

神犇题解：http://blog.csdn.net/popoqqq/article/details/45095755

我太弱了，推到$x=\frac{y * n!}{y-n!}$的时候没有想到把分母设成t，导致看不出什么关系。由于换元之后有$x=\frac{m^2}{t}+n!$，要求x和t是整数，那么关系就很明了了。。。这一题告诉我们以后推式子要熟练换元，特别是分母位置，常常就和因数挂钩了。。。

#include<cstdio>
#define N 1000005
#define ll long long
#define MOD 1000000007
using namespace std;
bool notprime[N];
int prime[N], pcnt;
int main()
{
    int n;
    scanf("%d",&n);
    for(int i = 2; i <= n; i++)
    {
        if(!notprime[i])
            prime[++pcnt]=i;
        for(int j = 1; j <= pcnt && 1ll*prime[j]*i <= n; j++)
        {
            notprime[prime[j]*i]=1;
            if(i%prime[j]==0)break;
        }
    }
    ll ans = 1;
    for(int i = 1; i <= pcnt; i++)
    {
        ll x=prime[i], cnt = 0;
        while(x<=n)
        {
            cnt+=n/x;
            x*=prime[i];
        }
        ans=ans*(cnt*2+1)%MOD;
    }
    printf("%lld\n",ans);
}