质因数分解。
#include <iostream>
#include <string.h>
#include <stdio.h>
using namespace std;
int isPrime[432],pChart[432];
long long cntC[432][432];
int cot;
void getPrime()
{
memset(pChart,0,sizeof(pChart));
int i,j;
for(i = 2;i < 432;i ++)
{
for(j = 2;j * i < 432;j ++)
pChart[i*j] = 1;
}
j = 0;
for(i = 2;i < 432;i ++) if(pChart[i] == 0) isPrime[j ++] = i;
cot = j;
}
int cal(int n,int p)
{
if(n < p) return 0;
else return n / p + cal(n / p, p);
}
int main()
{
int n,k,i,j;
getPrime();
//for(i = 0;i < cot;i ++) cout<<isPrime[i]<<endl;
/*for(i = 1;i < 432;i ++)
{
for(j = 1;j < i;j ++)
{
cntC[i][j] = 1;
for(k = 0;k < cot&&isPrime[k] <= i;k ++)
cntC[i][j] *= (cal(i,isPrime[k]) - cal(i-j,isPrime[k]) -cal(j,isPrime[k]) +1);
}
}*/
while(scanf("%d%d", &n, &k) != EOF)
{
__int64 cnt = 1;
//if(k > n / 2) k = n - k;
for(i = 0; isPrime[i] <= n && i < cot; i ++) {
int t = cal(n, isPrime[i]) - cal(n - k, isPrime[i]) - cal(k, isPrime[i]);
cnt *= (t + 1);
}
printf("%I64d\n", cnt);
}
return 0;
}