大数阶乘思考

  最近一场校招比试题考的是编程实现9999！，最后面试官竟然因为我写的程序最复杂给我一次面试机会，最后发现招聘的岗位是偏向于硬件，最后就没有最后了。

面试官说这个问题应该是很简单，直接long double n=1.0;for(int i=1;i<10000;i++)n*=i;printf("%e",n);回来之后我在在vs2008运行，截图如下

在这里说一下float double的最大存储吧

float 32 6-7 -3.4*10(-38)～3.4*10(38)

double 64 15-16 -1.7*10(-308)～1.7*10(308)

long double 128 18-19 -1.2*10(-4932)～1.2*10(4932)

下面废话不多说给出博主的解决方法：

#include <iostream>
using namespace std;
struct bigNum{
double base;
unsigned long eps;
};
void main(){
	bigNum plus;plus.base=1.0;plus.eps=0;
for(int i=1;i<10000;i++){
plus.base*=i;
if(plus.base>100000) {plus.eps+=5;plus.base/=100000;}
else if(plus.base>10000) {plus.eps+=4;plus.base/=10000;}
else if(plus.base>1000) {plus.eps+=3;plus.base/=1000;}
else if(plus.base>100) {plus.eps+=2;plus.base/=100;}
     else if(plus.base>10) {plus.eps+=1;plus.base/=10;}
}
cout<<"1*2*3*4*5*....*9999="<<plus.base<<"乘以10的"<<plus.eps<<"次方"<<endl;
}

输出结果：

与真实结果还是表近：

2846259680917054518906413212119868890148051401702799230794179994274411340003764443772990786757784775815884062142317........

虽然博主的结果不是精确的，但是结果还是可以接受的。博主的程序还可以改进实现更大数的阶乘。下面是一些牛人们实现的大数阶乘，但个人觉得都比较复杂，大家觉得好理解就行不多说了。

http://www.oschina.net/search?scope=code&q=%E5%A4%A7%E6%95%B0%E9%98%B6%E4%B9%98

http://www.cnblogs.com/confach/archive/2005/07/14/192703.html

希望大神指导