poj 2992(因子分解。。。。)

一直在想用线性筛素数的过程中，来处理组合的计算。。。。想了一下午，无果。。。。。感觉应该是可以的。。。。组后还是用了很土的方法过的。。。200ms++。。

这里用了C[ n ] [ k ] = n! / ( k!  *  ( n - k ) ! ) 然后把n！=p1^a1+p2^a2。。。。中a1,a2。。。预处理出来。。。同理把k！ 和（n-k）！ 也预处理出来 然后指数加减得出

C[ n ] [ k ]=( 总的a1+1 )*（总的a2+1）*。。。。。就可以了。。。

#include <iostream>
#include<cmath>
using namespace std;
#define LL long long
LL c[432][432]={0};
int isprim[432]={0},pri[400],p=0;
int count1(int x,int y){
    int sum=0;
    while(x){
        sum+=(x/y); x/=y;
    }
    return sum;
}
void cal(){
    int i,j;
    for(i=2;i<432;i++){
        if(!isprim[i]) pri[p++]=i;
        for(j=0;j<p && i*pri[j]<432;j++){
            isprim[i*pri[j]]=1;
            if(i%pri[j]==0)
                break;
        }
    }
    for(i=2;i<432;i++)
        for(j=0;j<p && pri[j]<=i;j++){
            c[i][pri[j]]=count1(i,pri[j]);
        }
}
int main(int argc, char** argv) {
    int n,k,i;
    long long ans;
    cal();
    while(scanf("%d%d",&n,&k)!=-1){
        ans=1LL;
        for(i=0;i<p && pri[i]<=n;i++)
            ans*=(long long)(c[n][pri[i]]-c[n-k][pri[i]]-c[k][pri[i]]+1);
        printf("%lld\n",ans);
    }
    return 0;
}