Choose and divide UVA - 10375

求两个组合数的商，数据大，肯定不能直接算。
根据唯一分解定理，每个数可以分解成若干个素数的幂的乘积的形式。可以把分解成的素数的幂记下来，最后知道每个素数的出现次数，依次求幂即可。
主要的收获：
1.学会用唯一分解定理
2.复习了欧拉筛
3.pow()的函数原型：double pow(double,double) 因此用来算整数的幂时，会有精度损失。
可以用重载来解决。
4.组合数的定义。
代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <cmath>
using namespace std;
const int N=1e4;
int p,q,r,s;
bool vis[N+5];
vector<int>prime;
int e[1235];
void f_prime()//打出一个素数表
{
    memset(vis,false,sizeof(vis));
    prime.clear();
    for(int i=2;i<=N;i++)
    {
        if(!vis[i])
            prime.push_back(i);
        for(int j=0;j<prime.size()&&prime[j]*i<=N;j++)
        {
            vis[i*prime[j]]=true;
            if(i%prime[j]==0)
                break;
        }
    }
}
void integer(int m,int d)
{
    for(int i=0;i<prime.size();i++)
    {
        while(m%prime[i]==0)
        {
            m/=prime[i];
            e[i]+=d;
        }
        if(m==1)
            break;
    }
}
void factorial(int n,int d)
{
    for(int i=1;i<=n;i++)
        integer(i,d);
}
int main()
{
    f_prime();
    while(scanf("%d%d%d%d",&p,&q,&r,&s)!=EOF)
    {
        memset(e,0,sizeof(e));
        factorial(p,1);
        factorial(s,1);
        factorial(p-q,-1);
        factorial(q,-1);
        factorial(r,-1);
        factorial(r-s,1);
        double ans=1.0;
        for(int i=0;i<prime.size();i++)
            ans*=pow(prime[i]*1.0,e[i]*1.0);//pow本身是一个返回double的函数
        printf("%.5lf\n",ans);
    }
    return 0;
}