AtCoder 2286

题意：求N！的约数的个数

---

以下部分内容转自：http://blog.csdn.net/ji414341055/article/details/5759157 
先说基本定理：

若正整数n可分解为p1^a1*p1^a2*...*pk^ak

其中pi为两两不同的素数,ai为对应指数，则n的约数个数为(1+a1)*(1+a2)*....*(1+ak)
    如180=2*2*3*3*5=2^2*3^2*5
    180的约数个数为(1+2)*(1+2)*(1+1)=18个。

    若求A/B的约数个数，A可分解为p1^a1*p2^a2*...*pk^ak，B可分解为q1^b1*q1^b2*...*qk^bk,

    则A/B的约数个数 为(a1-b1+1)*(a2-b2+1)*(a3-b3+1)...*(ak-bk+1).

然后说N的阶乘：

例如:20!
1.先求出20以内的素数,(2,3,5,7,11,13,17,19)
2.再求各个素数的阶数
    e(2)=[20/2]+[20/4]+[20/8]+[20/16]=18;
    e(3)=[20/3]+[20/9]=8;
    e(5)=[20/5]=4;
    ...
    e(19)=[20/19]=1;
    所以
    20!=2^18*3^8*5^4*...*19^1

解释： 
2、4、6、8、10、12、14、16、18、20能被2整除 
4、8、12、16、20能被4整除（即被2除一次后还能被2整除） 
8、16能被8整除（即被2除两次后还能被2整除） 
16能被16整除（即被2除三次后还能被2整除） 
这样就得到了2的阶。其它可以依次递推。

所以在求N的阶乘质数因数个数时，从最小的质数开始

#include<bits/stdc++.h>
#define LL long long 
using namespace std;

const int MOD = 1e9+7;
const int MAX = 1e6+10;
const int INF = 0x3fffffff;
bool a[MAX];

int C(int a,int b){
    if(a<b) return 0;
    else return a/b+C(a/b,b);
}//求n！的情况下， b的阶数
void init(){
    memset(a,true,sizeof(a));
    a[0]=a[1]=false;
    for(int i=2;i<=MAX;i++){
        if(a[i]){
            for(int j=i+i;j<MAX;j+=i){
                a[j]=false;
            }
        }
    }
}//把质数打表
int main(){
    int n;
    init();
    while(scanf("%d",&n)!=EOF){
        LL sum=1;
        int i=2;
        while(i<=n){
            if(a[i]){/若a[i]是质数
                sum=sum*(LL)(C(n,i)+1)%MOD;//相乘。
            }
            i++;
        } 
        printf("%lld\n",sum);
    }

    return 0;
}