精度计算-大数阶乘
我们的思路是把每4位数看做数组的一个元素来存储,例如:个、十、百、千存在a[0],万、十万、百万、千万存在a[1]以此类推。
我们用10的阶乘来模拟一下求结果大于4位数阶乘的过程,9的阶乘为362880,而10的阶乘为9的阶乘乘以10,在计算完9的阶乘时a[0] = 2880,a[1]=36,因为362880*10 = (36*10+2880*10/10000)*10000+28800%10000 = 3628800。同理其他高阶的数的计算过程和这个大同小异了
下面是我的C语言实现过程:
void factorial(int n) { long a[10000]; int i,j,c,m=0; a[0]=1; for(i=1;i<=n;i++) { c=0; for(j=0;j<=m;j++) { a[j]=a[j]*i+c; c=a[j]/10000; a[j]=a[j]%10000; } printf("\n"); if(c>0) { m++; a[m]=c; } } printf("\n%ld",a[m]); for(i=m-1;i>=0;i--) printf("%4.4ld",a[i]); }
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。