题意:已知C(m,n)=m!/(n!(m-n)!),求C(p,q)/C(r,s)。(p≥q,r≥s,p,q,r,s≤10000)
首先,求出10000以内的所有素数primes,然后用数组e表示当时结果的唯一分解式中的各个素数的指数。例如,e={1,0,2,0,0,0…)表示21*52=50。 (虽然紫书上有详解了,但是讲的真是太帅了!抄一下膜拜之!)
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int Max_n=1e4;
int p,q,r,s;
int ans=0;
int pri[Max_n];
bool vis[Max_n];
int e[Max_n];
void getPrimes(){
memset(vis,0,sizeof(vis));
for(int i=2;i<=Max_n;i++){
if(!vis[i]){
pri[ans++]=i;
for(int j=i*i;j<=Max_n;j+=i)
vis[j]=1;
}
}
}
void factorial(int n,int d){ //求n!(d=1表示在分子上 d=-1在分母)
for(int i=2;i<=n;i++){
int t=i;
for(int j=0;j<ans;j++){
while(t%pri[j]==0){
t/=pri[j];
e[j]+=d;
}
if(t==1)break;
}
}
}
int main()
{
getPrimes();
while(~scanf("%d%d%d%d",&p,&q,&r,&s)){
memset(e,0,sizeof(e));
factorial(p,1);
factorial(s,1);
factorial(r-s,1);
factorial(q,-1);
factorial(r,-1);
factorial(p-q,-1);
double sum=1;
for(int i=0;i<ans;i++){
sum*=pow(pri[i],e[i]);
}≤≤
printf("%.5f\n",sum);
}
return 0;
}