The binomial coefficient C ( m;n ) is de ned as C ( m;n ) = m ! ( m
直接算的话结果会爆double。利用唯一分解定理,先利用每个质因数的次数进行计算,最后算出答案。
#include<cstdio>
#include<cstring>
#include<cmath>
const int maxn=10000;
int prm[10010],cnt[10010],tot;
bool have[10010];
void add(int x,int d)
{
int i,j,tem;
for (i=2;i<=x;i++)
{
tem=i;
for (j=1;tem>1;j++)
while (tem%prm[j]==0)
{
tem/=prm[j];
cnt[j]+=d;
}
}
}
int main()
{
int i,j,k,m,n,p,q,r,s,x,y,z;
double ans;
for (i=2;i<=maxn;i++)
{
if (!have[i]) prm[++tot]=i;
for (j=1;j<=tot&&i*prm[j]<=maxn;j++)
{
have[i*prm[j]]=1;
if (i%prm[j]==0) break;
}
}
while (scanf("%d%d%d%d",&p,&q,&r,&s)==4)
{
for (i=1;i<=tot;i++)
cnt[i]=0;
add(p,1);
add(q,-1);
add(p-q,-1);
add(r,-1);
add(s,1);
add(r-s,1);
ans=1;
for (i=1;i<=tot;i++)
ans*=pow(prm[i],cnt[i]);
printf("%.5f\n",ans);
}
}