poj：2992 因子数量

题意：

求 组合数c(n,k)的因子数量

由算术基本定理很容易求得，不过第一次却T了，加了好多预处理，o1查询，才过

#include <iostream>
#include <stdio.h>
#include<string.h>
#include<algorithm>
#include<string>
#include<ctype.h>
using namespace std;
#define I64d lld
int prime[];
int isnotprime[];
int fac[][];
int nfac[][];
long long ans[][];
int np;
void setprime()
{
    np=;
    for(int i=;i<=;i++)
    {
        if(!isnotprime[i])
        {
            prime[np++]=i;
        }
        for(int j=;j<np&&i*prime[j]<=;j++)
        {
            isnotprime[i*prime[j]]=;
            if(i%prime[j]==)
                break;
        }
    }
    return;
}
void setfac()
{
    memset(fac,,sizeof(fac));
    for(int i=;i<=;i++)
    {
        int p=i;
        for(int j=;j<np;j++)
        {
            while(p%prime[j]==)
            {
                fac[i][j]++;
                p/=prime[j];
            }
        }
    }
    return;
}
void setans()
{
    memset(nfac,,sizeof(nfac));
    for(int i=;i<=;i++)
    {
        for(int j=;j<np;j++)
            nfac[i][j]=nfac[i-][j]+fac[i][j];
    }
    for(int n=;n<=;n++)
    {
        for(int k=;*k<=n;k++)
        {
            long long res=;
            for(int i=;i<np;i++)
            {
                res*=nfac[n][i]-nfac[k][i]-nfac[n-k][i]+;
            }
            ans[n][k]=ans[n][n-k]=res;
        }
    }
    return;
}
int main()
{
    setprime();
    setfac();
    setans();
    int n,k;
    while(scanf("%d%d",&n,&k)!=EOF)
    {
        printf("%I64d\n",ans[n][k]);
    }
    return ;
}