poj1306解题报告

题目大意:求C(m)(n)，输出结果的二进制不超过32位。。

思路：先将1~100的素数打表，然后把每个分子分母分解质因数，记录下分子因数分解的情况，分子质因数表加到一起，分母的表加到一起，最后把分子的质因数表减去分母的质因数表，然后结果就好球了。。

#include<iostream>#include<cmath>using namespace std;int prime[50];void find(){prime[0]=2;int i,j=1,n=3;while(j<50){for(i=0;prime[i]*prime[i]<=n;i++)if(n%prime[i]==0)break;if(prime[i]*prime[i]>n)prime[j++]=n;n++;}}int main(){find();int i,j,m,n;while(cin>>m>>n&&(m||n)){int fenzi,fenmu,lockz[50]={0},lockm[50]={0},nc=n;if(n>m/2)n=m-n;if(m!=0){for(i=1;i<=n;i++){fenzi=m-i+1;fenmu=i;for(j=0;prime[j]<=fenzi;j++)if(fenzi%prime[j]==0){lockz[j]++;fenzi/=prime[j--];}for(j=0;prime[j]<=fenmu;j++)if(fenmu%prime[j]==0){lockm[j]++;fenmu/=prime[j--];}}}__int64 o=1;for(i=0;i<50;i++)lockz[i]-=lockm[i];for(i=0;i<50;i++)if(lockz[i])o*=int(pow(double(prime[i]),double(lockz[i])));printf("%d things taken %d at a time is %I64d exactly./n",m,nc,o);}return 0;}