UVa 10375 Choose and divide (组合数相除)

10375 - Choose and divide

Time limit: 3.000 seconds 

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

The binomial coefficient C(m,n) is defined as

         m!C(m,n) = --------         n!(m-n)!

Given four natural numbers p, q, r, and s, compute the the result of dividingC(p,q) by C(r,s).

The Input

Input consists of a sequence of lines. Each line contains four non-negative integer numbers giving values forp, q, r, and s, respectively, separated by a single space. All the numbers will be smaller than 10,000 withp>=q and r>=s.

The Output

For each line of input, print a single line containing a real number with 5 digits of precision in the fraction, giving the number as described above. You may assume the result is not greater than 100,000,000.

Sample Input

10 5 14 993 45 84 59145 95 143 92995 487 996 4882000 1000 1999 9999998 4999 9996 4998

Output for Sample Input

0.12587505606.460551.282230.489962.000003.99960

思路：不要被题目那个公式误导了~自己算一算几个组合数，比如C(100,4)，发现分子分母都是4个乘数，那就好办了，写一个for循环每次一除一乘就是。

完整代码：

/*0.022s*/#include<cstdio>#include<algorithm>using namespace std;int main(){int p, q, r, s,maxn,i;double ans;while (~scanf("%d%d%d%d", &p, &q, &r, &s)){q = max(q, p - q), s = max(s, r - s);maxn = max(q, s);ans = 1.0;for (i = 1; i <= maxn; ++i){if (i <= q) ans = ans / i * (p - q + i);if (i <= s) ans = ans / (r - s + i) * i;}printf("%.5f\n", ans);}return 0;}