題意：

求解C（p,q）/ C（r，s），保留五位小數

題解：

劉汝佳老師的程式碼真的很好，思路很清晰，好好品味，好好學習

#include<math.h>
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;

#define MAXN 10010
#define LL long long
const LL mod=1000000007LL;
LL flag[MAXN];
int primer[MAXN],cnt;
void get_primer()
{
    flag[0]=flag[1]=1;
    int sqrtmaxn=(int)sqrt((double)MAXN);
    for(int i=2;i<MAXN;i++){
        if(primer[i]==0){
            flag[i]=flag[i-1];

            primer[cnt++]=i;
            if(i<=sqrtmaxn)
                for(int j=i*i;j<MAXN;j+=i)
                primer[j]=1;
        }
        else flag[i]=(flag[i-1]*i)%mod;
    }
}

int e[MAXN];
void add_integer(int n,int d){
    for(int i=0;i<cnt;i++){
        while(n%primer[i]==0){
            n/=primer[i];
            e[i]+=d;
        }
        if(n==1) break;
    }
}
void add_factorial(int n,int d){
    for(int i=1;i<=n;i++)
        add_integer(i,d);
}
int main()
{
    int p,q,r,s;
    get_primer();
    //freopen("in.txt","r",stdin);
    while(scanf("%d%d%d%d",&p,&q,&r,&s)!=EOF)
    {
        memset(e,0,sizeof(e));
        add_factorial(p,1);
        add_factorial(q,-1);
        add_factorial(p-q,-1);
        add_factorial(r,-1);
        add_factorial(s,1);
        add_factorial(r-s,1);
        double ans=1.0;
        for(int i=0;i<cnt;i++)
            ans*=pow(primer[i],e[i]);
        printf("%0.5lf\n",ans);
    }
    return 0;
}