poj2992：N!素因子个数，因子个数

题目大意：
求组合数Cnk的因子个数；
思路：
要求因子个数，由于因子可由素因子组合而成，所以先求得素因子个数
N!某素因子p的指数公式：
P(N!)=[N/p^1]+[N/p^2]+[N/p^3]+…+[N/p^n];
其中最后一个式子中n是使得p^n小于等于Ｎ的最大数；
因子个数公式：
cnt=(p1+1)（p2+1)(p3+1)….*(pn+1)；

此题要先打表，否则会超时;

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
const int maxn=440;
int prime[maxn+1];
long long p[maxn+1][maxn+1];
void getprime()
{
    memset(prime,0,sizeof prime);
    for(int i=2; i<=maxn; i++)
    {
        if(!prime[i]) prime[++prime[0]]=i;
        for(int j=1; j<=prime[0]&&i*prime[j]<=maxn; j++)
        {
            prime[i*prime[j]]=1;
            if(i%prime[j]==0) break;
        }
    }
}
long long cal(int n,int x)
{
   long long ans=0;
   while(n)
   {
       ans+=n/x;
       n/=x;
  }
  return ans;
}
int main()
{
    getprime();
    memset(p,0,sizeof p);
     for(int i=1;i<maxn;i++)
        for(int j=1;j<=prime[0];j++)
        p[i][prime[j]]=cal(i,prime[j]);
    int n,k;
    while(~scanf("%d%d",&n,&k))
    {
        long long ans=1;
        for(int i=prime[1],x=1;i<=n;i=prime[++x])
            ans*=(p[n][i]-p[n-k][i]-p[k][i]+1);
        cout<<ans<<endl;
    }
    return 0;
}