大数阶乘的算法

前几天朋友问我一个问题：“10000的阶乘怎么算？”当时我就有点懵，“10000”这个数字太大了，无论用什么数据类型保存结果都会溢出。这可怎么办呢？一时间束手无策。然后被一顿鄙视。后来经朋友的提醒，才恍然大悟，终于知道怎么实现了，原来是使用数组来模拟数字，这样无论结果数字有多大，只要数组的长度够长就能表示出来，用这个办法可以进行大数据的运算。看起来还是挺有用的。我把它用程序实现出来，如果有用到的地方还可以借鉴一下。（最起码还可以拿来鄙视别人^_^）
首先定义一个足够长的数组。
拿10000的阶乘为例，最后的结果长度是35660位，所以我们定义一个40000个成员的数组就可以了。
int result[40000];
其核心思想就是把计算结果每一位上的数字保存到一个数组成员中，例如：
把124保存至数组中，保存结果应该是
result[0] 4
result[1] 2
result[2] 1
这样肯定是没有问题的,一个int型数据存放一个小于10的数是绝对不会溢出。但是处理起来就稍微有点麻烦。
把整个数组看成一个数字，这个数字和一个数相乘的时候，需要每一位都和这个乘数进行相乘运算还需要把前一为的进位加上。运算方法和小学数学是一样的，乘积的个位是当前位上应该表示的数字，10位以上的需要进位。因为乘数不可能大于10000，所以乘数和一个小于10的书相乘的时候不会大于100000，再加上前一位的进位用一个int型数据来保持这个结果就没有问题。写法如下：
int 结果 = result[x] * 乘数 + 进位;
每一位的计算结果有了，把这个结果的个位数拿出来放到这个数组元素上：
result[x] = 结果%10;
接下来的工作就是计算出进位：
进位 = 结果 / 10;
这样一位一位的把整个数组计算一遍，最后可能还有进位，用同样的方法，把进位的数值拆成单个数字，放到相应的数组元素中。
最后输出一下结果，从最高位吧数字打印一遍就OK了。完整程序如下，在VC6.0下调试通过。

 

#include <iostream> using namespace std; void main() { int result[40000]; //保存结算结果的数组 int height = 1; //结果的最高位 int num; //计算阶乘的数字 cout<<"Input num:"; //控制台输入数字 cin>>num; if (num > 10000 || num < 0) { cout<<"Wrong num!"; return; } result[0] = 1; for (int i=1;i<=num;i++) { int res = 0; //进位 for (int j=0;j<height;j++) { int buf = result[j] * i + res; //计算结果 result[j] = buf % 10; //取当前位 res = buf / 10; //计算进位 } while (res) { result[height++] = res % 10; //取当前位 res /= 10; //计算进位 } } cout<<num<<"'s factorial is:"<<endl; //输出结果 for (int k=height-1;k>=0;k--) { cout<<result[k]; } cout<<endl; cout<<"length"<<height<<endl; }

完了。