Choose and divide UVA - 10375 (数论+唯一分解定理)

直击题目：https://vjudge.net/problem/19601/origin

题目大意：

有一个组合公式为C(m,n)=m!(n!(m-n)!),然后输入p,q,r,s,（<=10000）求C(p,q)/C(r,s);结果在10^8内。

题解：

把公式拆开写成分数形式然后分析。很明显直接乘会溢出，那么会想到怎么才能在算的过程中给其约分，可以想到每个数都能分解成某些素数的成绩，也就是用唯一分解定理，对其分解的素数的指数进行相加减就ok了。

#include <iostream>
#include <bits/stdc++.h>

using namespace std;

int prime[2000];
int e[2000];

int vis[11111];
int top;
void f2(int x,int d)
{
    for(int i=0; i<top; i++)
    {
        while(x%prime[i]==0)
        {
            e[i]+=d;
            x=x/prime[i];
        }
        if(x==1)return ;
    }
}
void f1(int n,int d)
{
    for(int i=2; i<=n; i++)
    {
        f2(i,d);
    }
}
int main()
{
    int p,q,r,s;
    memset(vis,0,sizeof(vis));
    for(int i=2; i<=sqrt(10005+0.5); i++)
    {
        if(!vis[i])
        {
            for(int j=i*i; j<=10005; j=j+i)
            {
                vis[j]=1;
            }
        }
    }
    top=0;
    for(int i=2; i<=10005; i++)
    {
        if(!vis[i])prime[top++]=i;
    }
    while(~scanf("%d%d%d%d",&p,&q,&r,&s))
    {
        memset(e,0,sizeof(e));
        f1(p,1);
        f1(r-s,1);
        f1(s,1);
        f1(q,-1);
        f1(p-q,-1);
        f1(r,-1);
        double sum=1.0;
        for(int i=0; i<top; i++)
        {
            sum*=pow(prime[i],e[i]);
        }
        printf("%.5lf\n",sum);
    }

    return 0;
}