代码:

---

/* 
   可以把n!的结果放在数组中，数组中每个元素都表示n!值的一位. 
   对整数范围内的n,求n!. 
   对于输入的n想办法昼精确地估计出n!所占的位数.就能确定数组元素的个数 
   可以将n!表示成10的次幂,即n!=10^M(10的M次方)则不小于M的最小整数就是 
   n!的位数，对该式两边取对数，有�M=log10^n!即： 
   M = log10^1+log10^2+log10^3...+log10^n 
   循环求和,就能算得M值，该M是n!的精确位数。 
   数组初始化时,令数组第一个元素(n!的第一位)为整数1,其余为0. 
   在数组中计算n!时是通过将数组中的值乘2,3,4,。。一直到乘n的方式得的 
   把数组的第一个元素看做是n!的最低位，最后一个元素看做是最高位。 
   */ 
  
#include <iostream> 
using namespace std; 
#include <cmath> 
#include <cstdlib> 
#include <iomanip> 
  
int getN();  //输入n 
int getBitNum(int n);  //求n!的位数 
char *init(int size); 
void calc(char *a,int n);  //求n! 
void display(char *a,int size); 
  
int main() 
{ 
    int n = getN(); 
    int size = getBitNum(n); 
    char *pa = init(size); 
    calc(pa,n); 
    display(pa,size); 
    delete []pa;  //note: 
} 
  
int getN() 
{ 
    int n; 
    cout << "Please input n of n! : "; 
    cin >> n; 
    while(n < 0) { 
        cout << "input error,again: "; 
        cin >> n; 
    } 
    if(n == 0) 
        exit(1); 
    return n; 
} 
  
int getBitNum(int n) 
{ 
    double sum = 1.0; 
    for(int i=1; i<=n; i++) 
        sum += log10((long double)i); //参数为double &和long double& 
                                      //不转换重载函数不能解析 
    return int(sum); //不转换会有warning 
} 
  
char * init(int size) 
{ 
    char *pa = new char[size]; 
    if(!pa) { //n太大时考虑申请内存不成功情况的执行动作 
        cout << "too large factor of " << size << endl; 
        exit(1); 
    } 
    pa[0] = 1; 
    for(int i=1; i<size; i++) 
        pa[i] = 0; 
    return pa; 
} 
  
void calc(char *a, int n) 
{ 
    double bitcount = 1; 
    int begin = 0; 
  
    for(int i=2; i<=n; ++i) { 
        long add = 0;  //源为long and=0;由于and是C++关键字，所以出错。 
                       //故所有为add的源文为and.不小心的出错点。 
        bitcount += log10((long double)i); 
        if(a[begin] == 0) 
            begin++; 
  
        for(int j=begin; j<int(bitcount); ++j) { 
            add += (i*a[j]); 
            a[j] = char(add % 10); 
            add /= 10; 
        } 
    } 
} 
  
void display(char *a, int size) 
{ 
    int bit = 0; 
    for(int i=size-1; i>=0; i--) { 
        if(bit % 50 == 0) 
            cout << endl << "The" 
                << setw(5) << (bit/50+1) << " \'s 50 bit: "; 
        cout << int (a[i]);  //思考：为什么需要进行显式的整形转换呢? 
        bit++; 
    } 
    cout << endl; 
} 

---