51nod1189 阶乘分数(数学,好题)

题目链接

1189 阶乘分数

题目来源：  Spoj

基准时间限制：1 秒 空间限制：131072 KB 分值: 80  难度：5级算法题

 收藏

 关注

1/N! = 1/X + 1/Y （0<x<=y），给出N，求满足条件的整数解的数量。例如：N = 2，1/2 = 1/3 + 1/6，1/2 = 1/4 + 1/4。由于数量可能很大，输出Mod 10^9 + 7。

Input

输入一个数N（1 <= N <= 1000000)。

Output

输出解的数量Mod 10^9 + 7。

Input示例

2

Output示例

2

李陶冶  (题目提供者)

题解：

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#include<queue>
#include<stack>
using namespace std;
#define rep(i,a,n) for (int i=a;i<n;i++)
#define per(i,a,n) for (int i=n-1;i>=a;i--)
#define pb push_back
#define fi first
#define se second
typedef vector<int> VI;
typedef long long ll;
typedef pair<int,int> PII;
const int inf=0x3fffffff;
const ll mod=1000000007;
const int maxn=1e6+100;
const int maxp=1e3+100;
bool is[maxn];
int a[maxn],cnt[maxn];
int tol=0;
void init()
{
    memset(is,false,sizeof(is));
    for(int i=2;i<maxp;i++)  //注意此处处理到1e3就可以了,处理到1e6会T
    {
        if(is[i]) continue;
        a[++tol]=i;
        for(int j=2*i;j<maxp;j+=i)
            is[j]=1;
    }
}
int main()
{
    ll ans=0;
    init();
    int n;
    scanf("%d",&n);
    if(n==1)
    {
        puts("1");
        return 0;
    }
    rep(i,2,n+1)
    {
        int t=i;
        int j=1;
        while(t!=1)
        {
            while(t%a[j]==0&&t!=1)
                t/=a[j],cnt[a[j]]++;
            j++;
            if(j>tol||a[j]*a[j]>t) break;
        }
        if(t!=1) cnt[t]++;   //注意此处,可能含有大于1e3的因子
    }
    ll res=1;
    rep(i,2,n+1) res=(res*(2*cnt[i]+1))%mod;
    ll tt=500000004;
    ans=(res+1)*tt%mod;
    printf("%lld\n",ans);
    return 0;
}