UVA-10375 唯一分解定理

#include<iostream>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<stdio.h>
#define rep(i,j,k) for(int i=j;i<=k;i++)

using namespace std;
const int N = ;
bool is_prime[N];
int prime[N];
int cnt=;
void get_prime(){
  int m=sqrt(N);
  memset(is_prime,,sizeof(is_prime));
  rep(i,,m){
     if (!is_prime[i]){
       for (int j=i*i;j<=N;j+=i){
         is_prime[j]=;
       }
     }
  }
  rep(i,,N){
    if (is_prime[i]==){
        prime[cnt++]=i;
    }
  }
}
int e[N];
void add_integer(int num,int d){
  for (int i=;i<cnt;i++){
     while(num%prime[i]==){
        num/=prime[i];
        e[i]+=d;
     }
     if (num==)break;
  }
}
void add_factorial(int n,int d){
    for(int i=;i<=n;i++)
        add_integer(i,d);
}
int main(){
  int p,q,r,s;
  get_prime();
  while(~scanf("%d%d%d%d",&p,&q,&r,&s)){
    memset(e,,sizeof(e));
    add_factorial(p,);
        add_factorial(q,-);
        add_factorial(p-q,-);
        add_factorial(r,-);
        add_factorial(s,);
        add_factorial(r-s,);
        int maxx = max(r,p);
        double ans=;
        for (int i=;i<=maxx;i++){
            ans*=pow(prime[i],e[i]);
        }
        printf("%.5f\n",ans);
  }
  return ;
}