POJ - 2992 Divisors （求阶乘的因数，）

 转载他人的博客https://www.cnblogs.com/chenxiwenruo/p/3569491.html

先说一个n!的性质：
n!的素因子分解中的素数p的个数为
n/p+n/(p^2)+...+n/(p^k)+...

《ACM-ICPC程序设计系列 数论及应用》上的方法，200+ms：
首先先求解435以内的素因子。
然后预处理出j!中每个素因子的个数，公式如下：
num[j][i]=j/prime[i]+num[j/prime[i]][i];

设n!中素因子p的个数为：a=n/p+n/(p^2)+...+n/(p^k)+...
那么(n/p)!中素因子p的个数为：b=n/(p^2)+...+n/(p^k)+...
很显然a=b+n/p，因此可以利用上述递推公式预处理出所有的j!中每个素因子的个数。

接下来就可以预处理出C(i,j)的因子个数，然后一切就好办了。

 这个是 800 ms 过的。

没有预处理，每次读取n和k后，利用公式计算n!,k!,(n-k)!的各个因素的个数，最后再总的求

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#define mem(x,v) memset(x,v,sizeof(x))
using namespace std;
typedef long long LL;
int n,k,num = 0;
int f[500],c[500];
bool vis[500];
void init(){
	mem(vis,1);
	vis[1] = 0;
	for (int i = 2; i <= 431; i++){
		if (vis[i]){
			f[num++] = i;
			for (int j = 2; i*j <= 431; j++)
				vis[i*j] = 0;
		}
	}
	return;
}
void slove(int x,int y){
	for (int i = 0; i < num; i++){
		int cc = 0,p = f[i];
		while(x/p){
			cc += x/p;
			p *= f[i];
		}
		c[i] = c[i] + y*cc;
	}
}
int main(){
	init();
	while(scanf("%d%d",&n,&k) != EOF){
		mem(c,0);
		slove(n,1);
		slove(k,-1);
		slove(n-k,-1);
		LL ans = 1;
		for (int i = 0; i < num; i++)
		ans *= (c[i]+1);		
		printf("%I64d\n",ans);
	}


	return 0;
}

这个是 200 ms 过的‘

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#define mem(x,v) memset(x,v,sizeof(x))
using namespace std;
typedef long long LL;
int n,k,cnt = 0;
int f[500],num[500][500];
LL C[500][500];
bool vis[500];
void init(){
	mem(vis,1);
	vis[1] = 0;
	for (int i = 2; i <= 435; i++){
		if (vis[i]){
			f[cnt++] = i;
			for (int j = 2; i*j <= 435; j++)
				vis[i*j] = 0;
		}
	}
	mem(num,0);
	for (int i = 0; i < cnt; i++)
		for (int j = 1; j < 435; j++)  //num[i][j]表示i!中素因子prime[j]的个数。
			num[j][i] = j/f[i] + num[j/f[i]][i];
//预处理出C(i,j)的因子个数
	for (int i = 1; i < 450; i++)
		for (int j = 1; j < i; j++){
			C[i][j] = 1;
			for (int k = 0; k < cnt; k++){
				int d = num[i][k] -num[i-j][k] - num[j][k];
					C[i][j] *= (d+1);  //C[i][j](0<=j<=i)表示组合C(i,j)的因子个数。
			}
		}
	return;
}

int main(){
	init();
	while(scanf("%d%d",&n,&k) != EOF){
		if (n == k || k == 0) printf("1\n"); else
		printf("%I64d\n",C[n][k]);
	}


	return 0;
}