poj2249

果然是推荐题啊~~~果然有收获啊~~~

这题就是让你求c(n,m)开始的时候我是这样想的不就是求(n*(n-1)*...*(k+1))/(n-k)!吗？觉得应该很是简单。。用数组进行记录n-k个数是否被除用两个循环不就好了啊记得离散的上机的那道题也是这样的。于是开心的敲完代码，提交WA~~~开始茫然，怎么会错呢...进入DISCUSS...看见DISCUSS 发现我的疑惑早就有人解答过了。。。n的范围是2^31..也就是说你的数组肯定开不了那样大。于是这条思路被舍弃。再看发现有人是用的那个m/n*c(m,n)=c(m-1,n-1)那个公式，突然想到原来做组合的不是敲过这样一个递归的吗。于是又一轮 开心的敲代码。提交，RE。。无语。。又逛进DISCUSS看到有个好心的孩子告诉大家第二种方法肯定会RE不用去试了。。。于是木有办法了。。再看了DISCUSS居然有人在用double类型。。好奇之下进去看了，只用了一个循环。最简单的思路。只是一个一个的乘起来就好了。。能够除掉的就除掉，不能的话最后再来除就好了。

当然完全可以不用double。。只是有一个小细节，觉得需要注意：就是从大到小的这样除。这样保存的数据相对来讲要小一点~~~

好了~~附代码（3份：1,WA，2.RE，3,AC）：

1.

#include<stdio.h>
#include<string.h>
main(){
long long n,k,sum,i,p,j;
int a[100000];
while(scanf("%lld%lld",&n,&k)!=-1 && n!=0){
sum=1;
memset(a,0,sizeof(a));
if(n-k>k)
k=n-k;
for(i=n;i>=k+1;i--){
p=i;
for(j=n-k;j>=1;j--){
if(a[j]==0 &&  p%j==0){
a[j]=1;
p=p/j;
}
}
sum*=p;
}
printf("%lld\n",sum);
}

}

2.

#include<stdio.h>
#include<string.h>
long long c(long long n,long long m){
if(m==0)
return 1;
else
return c(n-1,m-1)*n/m;
}
main(){
long long n,k;


while(scanf("%lld%lld",&n,&k)!=-1 && n!=0){

printf("%lld\n",c(n,k));
}
}

3.

#include<stdio.h>
#include<string.h>
main(){
long long n,k,sum,p,i,j,ans,num;
while(scanf("%lld%lld",&n,&k)!=-1 && n!=0){
p=1;
if(n-k>k)
k=n-k;
num=n-k;
for(i=n;i>=k+1;i--){
p*=i;
  if(num>=1 && p%num==0){
   p/=num;
   num--;
   }
} 
for(i=1;i<=num;i++)
p/=i;
printf("%lld\n",p);

}
}