poj2992 求组合数的因子个数

对 N! 进行素因子分解：

n! = 1*2*3*4*…*n,因此它一定含有小于等于n 的所有素数我们要对这些素数分别进行处理求出含有多少项。
例如我们求n!含有多少个素因子p

设f[n][p] 表示n!含有多少素因子p；

1*2*3*4*,,,,,*n = ( p*2p*3p*4p*,,,,*(n/p*p) ) * k   其中k不含有素因子p；

（注： n/p*p  =  (int)(n/p) * p）

这里把 n 以内的所有 p 的倍数拿出来，剩下的数组成 k ，很明显 k 里面不含有素因子p 

1*2*3*4*…*n = p^(n/p) * (n/p)!  * k ;  

已知 k 里面不含有 素因子 p ，那么此时只需要求 (n/p)! 包含的素因子 p 个数，这里又回到了母问题 ，所以用递归解决。

则有递归方程：f[n][p] = n/p + f[n/p][p] 
 

poj2992：http://poj.org/problem?id=2992

ac代码：

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
#define N 432
typedef long long LL;
const int cnt = 83;
int p[400] = {2,3,5,7,11,13,17,19,23,29,//10
              31,37,41,43,47,53,59,61,67,//9
              71,73,79,83,89,97,101,103,107,//9
              109,113,127,131,137,139,149,151,//8
              157,163,167,173,179,181,191,193,
              197,199,211,223,227,229,233,239,
              241,251,257,263,269,271,277,281,
              283,293,307,311,313,317,331,337,
              347,349,353,359,367,373,379,383,
              389,397,401,409,419,421,431};
int num[N+1][400];
void fac()
{
    memset(num,0,sizeof(num));
    for(int k=2;k<N;k++)
    for(int i=2;i<=k;i++)
    {
        int t = i;
        for(int j=0;j<cnt;j++)
        {
            while(t%p[j] == 0)
            {
                t /= p[j];//除以这个素因子 
                num[k][j]++;//这个素因子的数量+1 
            }
            if(t == 1)break;
        }
    }
}
int main()
{
    fac();
    LL ret[N+1][N+1];
    for(int i=0;i<N;i++)
        for(int j=0;j<N;j++)
            ret[i][j] = -1;//初始化 
    int n,k;
    while(~scanf("%d%d",&n,&k))
    {
        if(ret[n][k]!=-1){printf("%I64d\n",ret[n][k]);continue;}
        ret[n][k] = 1LL;
        for(int i=0;i<cnt;i++)
            ret[n][k]*=(LL)(num[n][i]-num[n-k][i]-num[k][i]+1);
        printf("%I64d\n",ret[n][k]);
    }
    return 0;
}