BZOJ2721 Violet5樱花（数论）

　　有(x+y)n!=xy。套路地提出x和y的gcd，设为d，令ad=x，bd=y。则有(a+b)n!=abd。此时d已是和a、b无关的量。由a与b互质，得a+b与ab互质，于是将a+b除过来得n!=abd/(a+b)。d/(a+b)可取的值不受a、b限制，那么只要满足ab|n!(a⊥b)就可以了。

　　将n!分解质因数，答案就很容易统计了。枚举质数数一下在n!中有几个即可。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
#define N 1000010
#define P 1000000007
int n,prime[N],cnt=,ans=;
bool flag[N];
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj2721.in","r",stdin);
    freopen("bzoj2721.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read();
    flag[]=;
    for (int i=;i<=n;i++)
    {
        if (!flag[i]) prime[++cnt]=i;
        for (int j=;j<=cnt&&prime[j]*i<=n;j++)
        {
            flag[prime[j]*i]=;
            if (i%prime[j]==) break;
        }
    }
    for (int i=;i<=cnt;i++)
    {
        int w=;
        for (int j=n;j;j/=prime[i]) w+=j/prime[i];
        ans=1ll*ans*(w<<|)%P;
    }
    cout<<ans;
    return ;
}