在刷题时偶然遇到了高精度阶乘的问题,于是立即自学了一下,在此作简要总结。
提出背景:普通的阶乘算法(不断用i*(i+1)循环)对于较大数的阶乘,由于运算结果位数过大,会造成溢出,使结果精度低。此时引出高精度计算阶乘的方法。
实现思路:用数组倒序存储每一位的数,手动调整进位,再倒序输出
程序代码:
#include <stdio.h>
#define LEN 100000
void factorial(int* p,int x);
int main(int argc,char* argv[])
{
int num[LEN]={0};
int x,i,j;
int* p=num;
scanf("%d",&x);
factorial(p,x);
for(j=LEN-1;j>=0;j--)
if(*(p+j)>0) break;
for(i=j;i>=0;i--)
printf("%d",*(p+i));
}
void factorial(int* p,int x)//计算x的阶乘,存入数组num
{
int sum,i,j,c;
*p=1;
for(i=2;i<=x;i++)
{
c=0;//c表示进位
for(j=0;j<LEN;j++)
{
sum=*(p+j)*i+c;
*(p+j)=sum%10;
c=sum/10;
}
}
}