处理阶乘问题

处理阶乘有三种办法：

（1）       传统意义上的直接递归，n的规模最多到20+，太小了，在本题不适用，而且非常慢

（2）       稍快一点的算法，就是利用log()化乘为加，n的规模虽然扩展到1000+，但是由于要用三重循环，一旦n规模变得更大，耗时就会非常之严重,时间复杂度达到O(n*m*(n-m))，本题规定了n，m用unsigned int32类型，就是说n,m的规模达到了21E以上，铁定TLE的。而且就算抛开时间不算，还存在一个致命的问题，就是精度损失随着n的增加会变得非常严重。

因为n有多大，就要进行n次对数运算，n规模一旦过大，就会丢失得非常严重了。所以这种方法是绝对不可取的，因为中途的精度丢失不是简单的四舍五入可以挽回的。

（3）       拆分阶乘，逐项相除，再乘以前面所有项之积。这种方法用一个循环就OK了，时间复杂度只有O(n-m)，非常可观。

 

 下面我根据程序详细说说算法（3）：

       double cnm=1.0;

       while(b>0)

              cnm*=(double)(a- -)/(double)(b- -);

 

这是我写的函数原型，计算的是 aCb

 

这种算法巧妙地利用了分子分母的关系，而不是把公示中的3个阶乘单独处理。

例如当 a=5,b=2时