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;
}