uva 10375

这个题参照紫书。惨惨惨，由于数组越界，一直wa，为什么找了很久，因为我坚信我数组一定没有越界。

回到问题，C(m, n) = m! / ((m − n)! n!)， 数好大！ 参考紫书知道了用唯一分解定理，N=P1a1P2a2P3a3......Pnan，这里P1<P2<P3......<Pn均为质数，其中指数ai是正整数。这样的分解称为 N 的标准分解式。比如7 = 7^1, 9 = 3 ^2 , 15 = 3^1 * 5 ^1。 即任何数都可以用素数表示。那么我们可以进行拆分。把C(m, n)  / C(r, s) 表示为一个数。以m!, (m - n) !, n ! ......分开计算。因为他们是乘除关系。分别找出他们的质因数，维护每个质因数的指数。最后pow(质因数，指数)就能算出P1a1P2a2P3a3......Pnan。

当然，得小心他们是阶乘。

/*
ID: DickensTone
LANG: C++
PROB: namenum
*/
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
const int maxn = 10000 + 7;
int prime[maxn], e[maxn];
int num[maxn], cnt = 0;
void print()
{
    for(int i = 0; i < cnt; i++)
        printf("%d\n", prime[i]);
}
void add_integer(int n, int d)
{
    for(int i = 0; i < cnt; i++)
    {
        if(n == 1) return ;
        while(n % prime[i] == 0)
        {
            n = n / prime[i];
            e[i] = e[i] + d;
        }
    }
}
void add_factor(int n, int d)
{
    for(int i = 2; i <= n; i++)
        add_integer(i, d);
}
int main()
{
    //freopen("namenum.in","r",stdin);
    //freopen("namenum.out","w",stdout);
    for(int i = 1; i < maxn; i++)
        num[i] = i;
    for(int i = 2; i < maxn; i++)
    {
        if(num[i] != -1)
        {
            //printf("%d\n", cnt);
            prime[cnt++] = i;
            int t = 1;
            while(t * i < maxn)
                num[i * t++] = -1;
        }
    }
    //printf("%d\n", cnt);
    //print();
    int p, q, r, s;
    while(scanf("%d%d%d%d", &p, &q, &r, &s) == 4)
    {
        //printf("%d %d %d %d\n", p, q, r, s);
        memset(e, 0, sizeof(e));
        add_factor(p, 1);
        add_factor(s, 1);
        add_factor(r - s, 1);
        add_factor(p - q, -1);
        add_factor(q, -1);
        add_factor(r, -1);
        double ans = 1;
        for(int i = 0; i < cnt; i++)
            ans = ans * pow(prime[i], e[i]);
        printf("%.5lf\n", ans);
    }
    return 0;
}