題目大意：給出  p ,q, r, s這四個數，C(m, n) = m! / (m − n)! n!   ，讓你求解   C(p, q)by C(r, s)  ，即兩個階乘相除。

思路：(   p!*s!*(r-s)!  ) /(  q!*(p-q)!*r!  )   

篩法求素數，唯一分解定理，用函式實現，從而求其各種階乘，程式碼如下

#include<iostream>//唯一分解定理
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int maxn = 1000005;
int prime[maxn],flag[maxn],e[maxn],num=0;
//篩法求素數
void Init()
{
    for(int i=2;i<=maxn;i++)
    {
        if(flag[i])
            continue;
        prime[num++]=i;
        for(int j=2;i*j<=maxn;j++)
            flag[i]=1;
    }
}
//計算因子的總個數
void add_integer(int n,int d)
{
    for(int i=0;i<num;i++)
    {
        while(n%prime[i]==0)
        {
            n/=prime[i];
            e[i]+=d;//因為其指數為d，所以每個因子的指數為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()
{
    Init();
    int p,q,r,s;
    while(~scanf("%d%d%d%d",&p,&q,&r,&s))
    {
        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<num;i++)
            ans*=pow(prime[i],e[i]);
        printf("%.5f\n",ans);
    }
    return 0;
}