UVA 10375 唯一分解定理

题意：求c(m,n)/c(r,s)

思路：先对素数打表，构造阶乘的素数，设置一个e数组用来存储素数对应序号i的个数，阶乘作为分子是+1，作为分母是-1，最后计算的时候直接ans*=pow（素数，对应的个数）即可

#include<bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
#define ll long long
#define fo freopen("in.txt","r",stdin)
#define fc fclose(stdin)
#define fu0(i,n) for(i=0;i<n;i++)
#define fu1(i,n) for(i=1;i<=n;i++)
#define fd0(i,n) for(i=n-1;i>=0;i--)
#define fd1(i,n) for(i=n;i>0;i--)
#define mst(a,b) memset(a,b,sizeof(a))
#define sd(n) scanf("%d",&n)
#define sdd(n,m) scanf("%d %d",&n,&m)
#define ss(s) scanf("%s",s)
#define sddd(n,m,k) scanf("%d %d %d",&n,&m,&k)
#define pans() printf("%d\n",ans)
#define all(a) a.begin(),a.end()
#define sc(c) scanf("%c",&c)
#define we(a) while(scanf("%d",&a)!=EOF)
const int maxn=200005;
const double eps=1e-8;
int prime[maxn],check[maxn],e[maxn];
void solve(int n,int d)
{
    for(int i=1;i<=n;i++)
    {
        int t=i;
        for(int j=0;;j++)
        {
            while(t%prime[j]==0)
            {
                t/=prime[j];
                e[j]+=d;
            }
            if(t==1)
                break;
        }
    }
}
int main()
{
    int p,q,r,s,tot=0;
    memset(check,0,sizeof(check));
    for(int i=2;i<maxn;i++)
    {
        if(!check[i])
        {
            prime[tot++]=i;
        }
        for(int j=0;j<tot;j++)
        {
            if(i*prime[j]>maxn)
            {
                break;
            }
            check[i*prime[j]]=1;
            if(i%prime[j]==0)
                break;
        }
    }
    while(cin>>p>>q>>r>>s)
    {
        memset(e,0,sizeof(e));
        solve(p,1);
        solve(q,-1);
        solve(p-q,-1);
        solve(r-s,1);
        solve(s,1);
        solve(r,-1);
        double ans=1.0;
        for(int i=0;i<maxn;i++)
        {
            ans*=(pow(prime[i],e[i]));
        }
        cout<<fixed<<setprecision(5)<<ans<<endl;
    }
    return 0;
}