大数阶乘
思路:
数组模拟阶乘的过程,处理好进位的问题与优化的问题即可。
代码一如下:
1 //大数阶乘 2 #include<iostream> 3 #include<cstring> 4 #define N 17000 //此处只能计算到5000,如有需要,可以修改N 5 using namespace std; 6 void jiecheng(int n); 7 int ans[N]; 8 int main() 9 { 10 int n; 11 while(cin>>n) 12 { 13 jiecheng(n); 14 } 15 return 0; 16 } 17 18 void jiecheng(int n) 19 { 20 int i,j,len,t,temp; 21 memset(ans,0,sizeof(ans));//初始化 22 ans[0]=1; 23 len=1; 24 for(i=2;i<=n;i++) 25 { 26 t=0;//进位的标志 27 for(j=0;j<=len;j++)//逐个相乘 28 { 29 temp=ans[j]*i+t; 30 ans[j]=temp%10; 31 t=temp/10; 32 if(j==len&&t!=0)//渐进式增加,可以减少循环的次数 33 len++; 34 } 35 }
int k=len;
while(ans[k]==0&&k>=0)--k; 36 for(i=k;i>=0;i--)//输出 37 cout<<ans[i]; 38 cout<<endl; 39 }
代码二 (位压缩--5位)
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #define N 10000 5 6 using namespace std; 7 8 int a[N]; 9 10 int main() 11 { 12 int n,t,i,j,k; 13 while(cin>>n) 14 { 15 memset(a,0,sizeof(a)); 16 a[0]=1;k=0; 17 for(i=2;i<=n;i++) 18 { 19 t=0; 20 for(j=0;j<=k;j++) 21 { 22 a[j]=a[j]*i+t; 23 t=a[j]/100000; 24 a[j]=a[j]%100000; 25 } 26 if(t>0) 27 { 28 k++;a[k]=t; 29 } 30 } 31 printf("%d",a[k]); 32 for(i=k-1;i>=0;i--) 33 { 34 printf("%05d",a[i]); 35 } 36 printf("\n"); 37 } 38 }