这题刚开始拿起来就果断用唯一分解定理,结果发现TLE。。。。。发现貌似数据超大,给做一个预处理。。。可是怎嘛预处理?不知道。后来看了别人的评论惊奇的发现前几天学过的一个定理用上了场,就是判断N!中有几个质因子:P(N!)=N/i+N/i^2+N/i^3+.....N/i^m(i为一个质因子,m是使N/i^m=0的最小值)
这样的化就可以预处理了,然后利用公式C(N,K)=N!/(K!*(N-K)!)和约数定理,就可以得到答案。。。。下面附代码:
#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cstring>
using namespace std;
int a[450],b[450],p[450][450];
int m;
int pj(int n,int x)
{
int r=0,y=x;
while (n>=y)
{
r+=n/y;
y*=x;
}
return r;
}
int main()
{
int n,k,i,j;
for (i=2;i<=431;i++)
if (!a[i])
{
for (j=i*2;j<=431;j+=i) a[j]=1;
b[++m]=i;
}
for (i=1;i<=431;i++)
for (j=1;j<=m;j++)
p[i][j]=pj(i,b[j]);
while (scanf("%d",&n)!=EOF)
{
scanf("%d",&k);
__int64 r=1;
for (i=1;i<=m;i++)
if (n>=b[i])
{
r*=(p[n][i]-p[n-k][i]-p[k][i]+1);
}
printf("%I64d\n",r);
}
return 0;
}