Uva-10375 Choose and divide

题意：已知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…)表示2¹*5²=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;
}