做一下N阶乘的小结:
#include<iostream>
#define N 500000
using namespace std;
//使用long long时候N=21溢出
long long fac(long long n)
{
if(1==n)return 1;
else return n*fac(n-1);
}
//当阶乘数N>21的时候,但是这种方法受到VS2008编译器栈空间大小的限制,
//VS2008栈空间大小为1M(默认,可以修改),定义数组保存阶乘结果,
//sizeof(int)=4,2^21/4=2^19;阶乘结果最大的位数为2^19;
//改进的办法:申请动态数组,将结果保存在数组中,数组可以是整型数组和字符数组
void fac2(int f[],int n)
{
int i,j,up,sum;
for(i=2;i<n;++i)
{
up=0;
for(j=0;j<N;++j)
{
sum=f[j]*i+up;
f[j]=sum%10;
up=sum/10;
}
}
}
//在这N可以设定大于2^19的数,上限为内存的大小.
void fac3(char f[],int n)
{
int i,j,up,sum;
for(i=2;i<n;++i)
{
up=0;
for(j=0;j<N;++j)
{
f[j]=f[j]-'0';
sum=f[j]*i+up;
f[j]=(sum%10)+'0';
up=sum/10;
}
}
}
int main()
{
int n;
int j=0;
/*定义整型数组来保存结果
int *f;
f=(int *)malloc(sizeof(int)*N);
for(j=1;j<N;++j)f[j]=0;
f[0]=1;
*/
//定义字符串来保存结果,字符串一个字符占1字节,一个整型占4字节
char *f;
f=(char *)malloc(sizeof(char)*N);
for(j=0;j<N;++j)f[j]='0';
f[0]='1';
//注意不能这样初始化空间
/*
while(f)
{*f=0;f++; }
*/
cin>>n;
fac3(f,n);
for(j=N-1;f[j]=='0';--j);
for(;j>=0;--j)cout<<f[j]-'0';
cout<<endl;
cout<<sizeof(int)<<endl;
cout<<sizeof(char)<<endl;
return 0;
}
有没有人没有听说long long型,在VS2008中sizeof(long long)=8,占64位,本例测试了当N=50000,n!=20000!的时候。