UVA 10375 Choose and divide（组合数运算）

题意：运算两个组合数的商。

解题思路：根据唯一分解定理分解阶乘，相互约掉素因子，最后再根据素因子各个位置上的指数算出答案。

https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1316

#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<cctype>
#include<list>
#include<iostream>
#include<map>
#include<queue>
#include<set>
#include<stack>
#include<vector>

using namespace std;

#define FOR(i, s, t) for(int i = (s) ; i <= (t) ; ++i)
#define REP(i, n) for(int i = 0 ; i < (n) ; ++i)

int buf[10];
inline long long read()
{
    long long x=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}

inline void writenum(int i)
{
    int p = 0;
    if(i == 0) p++;
    else while(i)
        {
            buf[p++] = i % 10;
            i /= 10;
        }
    for(int j = p - 1 ; j >= 0 ; --j) putchar('0' + buf[j]);
}
/**************************************************************/
#define MAX_N 10005
const int INF = 0x3f3f3f3f;
int prime[MAX_N];
int isprime[MAX_N];
int e[MAX_N];
//int maxx =0;

int cnt = 0;

void make_prime()
{
    cnt = 0;
    memset(isprime, 1, sizeof(isprime));
    isprime[1] = 0;
    for(long long i =2  ; i <= MAX_N ; i++)
    {
        if(isprime[i])
        {
            prime[cnt++] = i;
            for(long long j = i * i ; j <= MAX_N ; j += i)
            {
                isprime[j] = 0;
            }
        }
    }
}

void f(int n, int d)
{
    for(int i = 0 ; i < cnt ; i++)
    {
        int c = 1;
        if(prime[i] > n)
        {
            return;
        }
        int p = prime[i];
        int tmp = p;
        while(p <= n)
        {
            e[i] += d * (n / p);
            c++;
            p = pow(tmp, c);
        }
    }
}

int main()
{
    int p,q,r,s;
    make_prime();
    while(~scanf("%d%d%d%d", &p, &q, &r, &s))
    {
        memset(e, 0, sizeof(e));
//        f(20, 1);
//        for(int i = 0 ; i < 8 ; i++)
//        {
//            cout<<e[i]<<endl;
//        }
        f(p, 1);
        f(q, -1);
        f(p - q, -1);
        f(r, -1);
        f(s, 1);
        f(r - s, 1);
        double ans = 1.0;
        for(int i = 0 ; i < cnt ; i++)
        {
            ans *= pow(prime[i], e[i]);
        }
        printf("%.5lf\n", ans);
    }


    return 0;
}