别让习惯害了你，注意size_t变量作循环条件时，逆序是不可行的！

欢迎转载，但请标明作者 “九天雁翎”，当然，你给出这个帖子的链接更好。

在自己编写求自然数内n的n!时碰到的问题。

#include <iostream>
#include <iomanip>
#include <cmath>
#include <cstdlib>
using namespace std;

size_t getbit(const size_t &aN);
void init(unsigned char *p,size_t &am);
void compute(unsigned char *p,const size_t &an);
void print(unsigned char *p,const size_t &am);
int main()
{
cout <<"Input the /"n/" you want to compute:";
size_t n;
cin >> n;
if(n == 0) //检验输入
{
cerr<<"Please input a integer more than zero."<<endl;
}
size_t m = getbit(n);
unsigned char *pn = new unsigned char[m];
if(!pn) //检验空间分配
{
cerr<<"The factor is too big."<<endl;
}
init(pn,m); //初始化
compute(pn,n); //计算n!
print(pn,m); //输出
delete []pn;
return 0;
}
//-----------------------------------------------------------
size_t getbit(const size_t &aN)
{
double sum = 1.0;
for(size_t i = 1;i<=aN;++i)
{
sum+=log10(double(i));
}
return size_t(sum);
}
//------------------------------------------------------------
void init(unsigned char *p,size_t &am)
{
p[0]=1;
for(size_t i = 1;i != am;++i)
{
p[i] = 0;
}
}
//------------------------------------------------------------
void compute(unsigned char *p,const size_t &an)
{
double bitcount = 1.1;
size_t begin = 0;
for(size_t i = 2;i<=an;++i)
{
size_t and = 0;
bitcount +=log10(double(i));
if(p[begin]==0)
++begin;
for(size_t j = begin ; j < size_t(bitcount);++j)
{
and += i * p[j];
p[j] = unsigned char(and % 10);
and /= 10;
}
}
}
//----------------------------------------------------------------
void print(unsigned char *p,const size_t &am)
{
size_t bit = 0;
for(size_t i = am - 1;i >= 0;--i)
{
if(bit % 50 == 0)
cout <<endl<<"第"
<<setw(3)<<(bit/50+1)<<"个50位:";
cout << size_t(p[i]);
++bit;
}
}
首先声明的是，我是在学习用数组处理这个问题，所以没有用vector，而是动态分配了数组，求位数用的是（n!=10^M）估算算法，每一位只保存结果的一位，为了节省空间用unsigned char保存一位数，在在vc2005下编译通过了，不过运行却得不出结果，还有假如在最后把输出的for(size_t i = am;i >= 0;--i)改为
for(int i = am;i >= 0;--i)，结果会好一点，但是还是得不出正确结果，正准备找高手说明下问题。还有用size_t为什么似乎会导致这里的死循环，而int不会？
最后我发现，把print（）定义成这个样子，反而又可以。所以比较肯定的说是输出函数的问题，但是我实在不知道是什么问题。
void print(unsigned char *p,const size_t &am)
{
cout <<size_t(p[2])<<size_t(p[1])<<size_t(p[0])<<endl;
}

最后我老实编了一个正常简单的逆序循环发现也不可以，我刚开始甚至还怀疑是不是size_t就没有正常的“--”自减运算符，然后编小程序证明不是，而且在循环条件〉=3，2，1时都没有问题，就到了大于0就不行了，我这时才突然醒悟过来，其实size_t 就是unsigned int啊，一个没有符号的整数在C++里面就意味着一个永远为正的数，我的行为就像是等一个数的绝对值，一个数开2次方根，一个数的平方去小于零一样可笑！呵呵，大家可不要犯我一样的错误啊。然后我把print()改成下面这个样子后，果然一切都正常了。

void print(unsigned char *p,const size_t &am)
{
size_t bit = 0;
for(size_t i = am;i != 0;--i)
{
if(bit % 50 == 0)
cout <<endl<<"第"<<setw(3)<<(bit/50+1)<<"个50位:";
cout << size_t(p[i-1]);
++bit;
}
cout<<endl;
}