(Relax 数论1.16)POJ 2992 Divisors(C[N][K]中含素数的个数)

求组合数C(n, k)的约数的个数 (0 ≤ k ≤ n ≤ 431)。

题目链接：http://poj.org/problem?id=2992

——>>3个公式：

  1、n!中含素数p的个数为n/p + n/p^2 + n/p^3 + ...（到0停）程序中通过cal函数实现

  2、C(n, k) = n! / (n-k)! / k!

  3、n = p1^a1*p2^a2*...*pk^ak约数的个数为（a1+1）(a2+1)...(ak+1)

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
/**
* 2 2 3 17
2 2 3 29

2 2 5 37
*/
using namespace std;

int su[505];//这里不要开太大,否则会TLE
bool u[505];
int num = 0;
int n;

void prepare(){
int i,j;
memset(u,true,sizeof(u));

for(i = 2 ; i <= 500 ;++i){
if(u[i]){
su[++num] = i;
}

for(j = 1 ; j <= 500 ; ++j ){
if(i*su[j] >500){
break;
}

u[i*su[j]] = false;

if(i % su[j] ==0){
break;
}
}
}
}

int cal(int n,int pri){//n!含素数pri的个数
return (n<pri)?0:(n/pri+cal(n/pri,pri));
}

int main(){
prepare();

int n,k;
while(scanf("%d%d",&n,&k)!=EOF){
long long sum = 1;
int i;
for(i = 1 ; i <= num ; ++i){//n! / (n-k)! / k!含素数的个数为:cal(n,su[i]) - cal(n-k,su[i]) - cal(k,su[i]) + 1
sum *= (cal(n,su[i]) - cal(n-k,su[i]) - cal(k,su[i]) + 1);
}

printf("%lld\n",sum);
}

return 0;
}