转载他人的博客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;
}