Choose and divide

题意

其实就是两个组合数相除C(a, b) 除C(c, d).,最后化简一下发现其实就是
a! x (c-d)! xd! 除以 c! x b! (a-b) !,虽然题目中说了最后结果不会太大int范围，但是我们如果直接循环来求阶乘，那么过程中会出现比较大的数，一种思路是一边乘一边除。

我一开始错误的思路

我认为是三个数的阶乘除以三个数的阶乘，分子的话每个数所在的数组加1，比如说分子上有6！，那么a[1]–a[6]都+1，还有一个8！，那么数组a1到a8再每一个都加1，对于分母上的数，每个数所对应的数组序列都减一，分子上的3个数和分母上的3个数如果都算完了的话，再开始对数组a操作，正数相乘，比如说a[3]=2;那么ans=ans33,如果说a[4]=-1;那么ans=ans/4;
大体上就是这样来求，题目中的样例都能过，但是最后wa，后来想了想，可能过程中还是有特别大的数，比如某一个样例中有a[9000]-a[9100]都是1，那么这个过程相乘时肯定越界int

唯一分解定理

概念：任意一个大于0的正整数都能被表示成若干个素数的乘积且表示方法是唯一的；整理可以将相同素数的合并；可以得到公式：n = P1^a1 * P2^a2 * …………* (P1 < P2 < ……Pn);
比如说 18，可以分解成233，36分解成2的2次方乘以3*3，所以在我们求 的时候只需要统计有多少个素数，然后再求幂，因为是相乘，再一个求素数时用的是埃氏筛法。

/*
* Copyright(c)
* All rights reserved.
* Author : Michael
* Date : 2020-02-21-19.02.31
*/
#include<iostream>
#include<algorithm>//快排sort()
#include<cstdio>//能不用cin就不用
#include<cstring>
#include<cmath>
#include<map>
#include<vector>
#include<queue>
#include<set>
#define ll long long
using namespace std;
bool p[10003];//判定素数的数组
int pr[10003];
int e[10002];
int cnt;
void primes(int n)//true表示不是素数
{
    cnt=0;
    memset(p,false,sizeof(p));
    for(int i=2;i<=n;i++)
    {
        if(p[i]) continue;//如果不是素数
        pr[cnt++]=i;
        for(int j=i;j<=n/i;j++)//之所以从x^2开始是因为前面的其实已经排除过了
        {
            p[i*j]=true;//不再是素数
        }
    }
    //cout<<pr[0]<<endl;
}
void solve(int n,int d)//唯一分解定理
{
    for(int i=0;i<=cnt;i++)
    {
        while(n%pr[i]==0)
        {
            n=n/pr[i];
            e[i]+=d;
        }
        if(n==1)
            break;
    }

}
void Factorial(int n,int d)
{
    int i;
    for( i=1;i<=n;i++)
    {
        solve(i,d);
    }
}
int main(){
    int a,b,c,d;
    double ans=1.0;
    primes(10000);//算一下10000内的素数
    while(scanf("%d%d%d%d",&a,&b,&c,&d)!=EOF)
    {
        ans=1.0;
        memset(e,0,sizeof(e));
        Factorial(a,1);
        Factorial(c-d,1);
        Factorial(d,1);
        Factorial(c,-1);
        Factorial(b,-1);
        Factorial(a-b,-1);
        for(int i=0;i<=cnt;i++)
        {
            ans*=pow(pr[i],e[i]);
        }
        printf("%.5f\n",ans);
    }
    return 0;
}