UVA 10375 选择与除法（唯一分解定理）

该方法使得最后统一算出结果，不用担心中间过程溢出。

#include<bits/stdc++.h>

using namespace std;

vector<int> primes;
int vis[10005];
int e[10000];

void add_integer(int n,int d);
void add_factorial(int n,int d){
    for(int i=1;i<=n;i++)
    add_integer(i,d);
}
void add_integer(int n,int d){           //d为正表示相乘，为负表示除 
    for(int i=0;i<primes.size();i++){
        while((n%primes[i])==0){
            n/=primes[i];
            e[i]+=d;
        }        
        if(n==1) break;
        
    }
    
}

int main(){
    int p,q,r,s;
    int c=0;
    memset(vis,0,sizeof(vis));

    
    for(int i=2;i<=10000;i++)
    if(!vis[i]){
    primes.push_back(i);
    for(int j=i*i;j<=10000;j+=i){
        vis[j]=1;
        
    }
}
    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;
    for(int i=0;i<primes.size();i++)
    ans*=pow(primes[i],e[i]);
        
        printf("%.5lf\n",ans);    
    }
    
}