樱花

问题 Q: 樱花
时间限制: 1 Sec 内存限制: 128 MB
提交: 19 解决: 4

提交

题目描述
求不定方程：

1/x+1/y=1/n!

的正整数解 (x,y)(x,y)(x,y) 的数目。

输入
一个整数 n。

对于 30% 的数据，n≤100；
对于全部数据，1≤n≤106。

输出
一个整数，表示有多少对 (x,y) 满足题意。答案对 109+7 取模。
样例输入
2
样例输出
3
提示
[提交][状态]

分析：转化为求 (n!)^2的质因数个数

对于n!，可以发现小于等于n的质数一定都是n!的质因数，大于n的质数一定不是。那么打出质数表。对于质数ei，在n!的分解式中的指数ki=n/ei+n/ei^2 +……+n/ei^x， 其中n>=ei^x。

#include<bits/stdc++.h>
#define maxn 1000111
#define mod 1000000007
#define ll long long
using namespace std;
bool prime[maxn];
ll Prime[maxn];
ll num, ans = 1, n, cnt[maxn];
void make_prime(){
   memset(prime, true, sizeof(prime));
   prime[0] = prime[1] = false;
   for(int i = 2 ; i < maxn ; ++ i){
        if(prime[i]){
            Prime[num++] = i;
        }
        for(int j = 0 ; j < num && Prime[j] * i < maxn ; ++ j){
            prime[i * Prime[j]] = false;
            if(!(i%Prime[j])){
                break;
            }
        }
   }
   return;
}
int main()
{
    make_prime();
    cin>>n;
    for(int i = 0 ; i < num ; ++ i){
        for(ll j = Prime[i] ; j <= n ; j *= Prime[i]){
            cnt[i] += n / j;
            cnt[i] %= mod;
        }
    }
    for(int i = 0 ; i < num ; ++ i){
        ans = (((cnt[i] * 2 + 1) % mod) * (ans % mod)) % mod;
    }
    cout<<ans % mod;
    return 0;
}