UVa 10375 - Choose and divide

唯一分解， 筛出素数即可。

int p, q, r, s;
int e[maxn];

void add_int(int n, int d) {
    for(int i = 0; i < 1229; i++) {
        while(!(n % primes[i])) {
            n /= primes[i];
            e[i] += d;
        }
        if(n == 1) break;
    }
}

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

int main() {
    while(~scanf("%d%d%d%d", &p, &q, &r, &s)) {
        memset(e, 0, sizeof(e));
        add_fac(p, 1);
        add_fac(q, -1);
        add_fac(p - q, -1);
        add_fac(r, -1);
        add_fac(s, 1);
        add_fac(r - s, 1);
        double ans = 1;
        for(int i = 0; i < 1229; i++)
            ans *= pow(primes[i], e[i]);
        printf("%.5f\n", ans);
    }
    return 0;
}