UVA10375唯一分解定理运用

本题用到了唯一分解定理，唯一分解定理就是指任意一个数可以被分解为若干素数相乘的形式。数学描述不太严谨，但我想大家都懂的。
话不多说贴代码吧。

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
using namespace std;

const int maxn = 1e4+20;
int p,q,r,s;
int prime[maxn],e[maxn],tot;
bool v[maxn];

void init(){
    memset(v,1,sizeof(v));
    int m = sqrt(maxn+0.5);
    tot = 0;
    for (int i=2; i<=m; i++) if (v[i]){
        for (int j=i*i; j<maxn; j+=i) v[j] = 0;
    }
    for (int i=2; i<maxn; i++) if (v[i]) prime[++tot] = i;
}

void add_integer(int n, int d){
   for (int i=1; i<=tot; i++) {
      while (n%prime[i] == 0) {e[i]+=d; n/=prime[i];}
      if (n==1) break;
   }
}

void add_factorial(int n, int d){
    for (int i=2; i<=n; i++) add_integer(i,d);
}

int main(){
    init();
    while (scanf("%d %d %d %d",&p,&q,&r,&s)!=EOF){
        memset(e,0,sizeof(e));
        add_factorial(p,1);
        add_factorial(q,-1);
        add_factorial(p-q,-1);
        add_factorial(r,-1);
        add_factorial(r-s,1);
        add_factorial(s,1);
        double ans = 1;
        for (int i=1; i<=tot; i++) ans *= pow(prime[i],e[i]);
        printf("%.5lf\n",ans);
    }
    return 0;
}