大数阶乘

原创 2015年07月10日 20:46:14
/*【分析】
为了保存结果,先分析1000!大约等于4×102567,因此可以用一个3000个元素的数组f保存。
让f[0]保存结果的个位,f[1]是十位,f[2]是百位,…,则每次只需要模似手算即可完成n!。
在输出时需要忽略前导0。
注意,如果结果本身就是0,那么忽略所有前导0后将什么都不输出。所幸n!肯定不等于0,因本题可以忽略这个细节。
*/
#include <stdio.h>
#include <string.h>
const int maxn = 3000;
int f[maxn];
int main()
{
	  int i, j, n;
	  scanf("%d", &n);
	  memset(f, 0, sizeof(f));
	  f[0] = 1;
	  for(i = 2; i <= n; i++) 
	  { /* 乘以i */
	    int c = 0;
	    for(j = 0; j < maxn; j++) 
		{
	      int s = f[j] * i + c;
	      f[j] = s % 10;
	      c = s / 10;
	    }
	  }
	  for(j = maxn-1; j >= 0; j--) 
	   if(f[j]) 
	     break;  /* 忽略前导0 */
	  for(i = j; i >= 0; i--)
	    printf("%d", f[i]);
	  printf("\n");
	  return 0;
}



这段代码:
for(i = 2; i <= n; i++) { /* 乘以i */
  int c = 0;
  for(j = 0; j < maxn; j++) {
    int s = f[j] * i + c;
    f[j] = s % 10;
    c = s / 10; 
  }
}

实际上就是求阶乘的核心代码。由于计算机的字长有限,能表示的整数范围远远小于1000阶乘的结果。目前PC机通常是32位,只能表示最大4294967295的(无符号)整数。而15阶乘的结果就已经是:1307674368000,远远超过32位的表示范围,更不用说1000的阶乘了。因此,程序里采用整数数组的方式存储“超长整数”,以保持精度。程序里采用int f[]存储阶乘的结果。例如,15阶乘为1307674368000,采用int f[]表示为:f[0] = 0 f[1] = 0 f[2] = 0 f[3] = 8 f[4] = 6 f[5] = 3 f[5] = 4 f[5] = 7 f[5] = 6 f[5] = 7 f[5] = 0 f[5] = 3 f[5] = 1 程序里,f[j] = s % 10;c = s / 10; 就是模仿手工计算的方式,采用10进制,一位位计算,大于10就进到下一位。你自己手工演算一下就明白了。
版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

donet 求 大数阶乘

  • 2015-03-06 15:19
  • 74KB
  • 下载

大数阶乘C++代码

大数阶乘

1.防止溢出的一种方式,有些数字超出int范围,这时候用int就会溢出。这时候要用到大数的思想。 例子:#include int main() { int n;//阶乘大小 prin...

大数阶乘问题

  • 2013-10-10 12:54
  • 974B
  • 下载

C/C++编程小练习 大数阶乘

我的这个大数阶乘是完全基于模拟手动运算的,速度略慢,其中的核心函数是之前我写的一个大数乘法算法{CSDN:CODE:2031929} 最大支持4万位的number输出(10000!是35660位,也可...

C# 大数阶乘

  • 2013-04-14 20:16
  • 99KB
  • 下载

One Day One Step 之大数阶乘

     在说大数阶乘之前,我先来说一下大数乘法。一般来说,当操作数的数量级很大的时候,肯定是要用数组进行模拟的!之所以说是模拟,就是因为他模仿了我们我们平时做乘法时的步骤!不过我想了很久都...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)