uva 10375 Choose and divide

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
using namespace std;
const int maxn = 10000+5;
int e[maxn],es;
int f[maxn];

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

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

int main()
{
    int i,p,q,r,s,j,vis[maxn];
    double ans;
    memset(vis,0,sizeof(vis));
    es=0;
    for(i=2; i<maxn; i++)
    {
        if(!vis[i])
        {
            e[es++]=i;
            //printf("%d\n",i);
            for(j=i*i; j<maxn; j+=i) vis[j]=1;
        }
    }
    while(~scanf("%d%d%d%d",&p,&q,&r,&s))
    {
        memset(f,0,sizeof(f));
        add(p,1);
        add(s,1);
        add(r-s,1);
        add(p-q,-1);
        add(q,-1);
        add(r,-1);
        ans=1;

        for(i=0; i<es; i++)
        {
            ans*=pow(e[i],f[i]);
        }

        printf("%.5f\n",ans);
    }
    return 0;
}