大数运算之阶乘

发布时间：2021-01-18 13:30:18 所属栏目：大数据来源：网络整理

导读：Int不能存储较大数 int 代表有符号整数，也就是说，用 int 声明的变量可以是正数，可以是负数，也可以是零，但是只能是整数。标准规定 int 的最小取值范围是 -32767 到 32767。 int 的取值范围因机器而异，但是一定要大于或者等于 -32767 到 32767。一般来说

副标题[/!--empirenews.page--]

Int不能存储较大数

int 代表有符号整数，也就是说，用 int 声明的变量可以是正数，可以是负数，也可以是零，但是只能是整数。标准规定 int 的最小取值范围是 -32767 到 32767。

int 的取值范围因机器而异，但是一定要大于或者等于 -32767 到 32767。一般来说，int 占用一个字的内存空间。因此，字长为 16 位（Bit）的旧式 IBM 兼容机使用 16 位来储存整型 int ，取值范围是-32768 到 32767 。目前的个人电脑一般都是 32 位字长的，这些电脑中，int 一般也是 32 位的，取值范围是 -2147483648 到 2147483647。对于使用 64 位 CPU 的电脑，使用更多位储存 int 也是很自然的事情，取值范围当然也会更大。

Int:

16位 -2^15? ~ 2^15-1

32位 -2^31? ~2^31-1

64位 -2^63? ~2^63-1

对于大数感觉长度还是不够，大数动辄都是几百位数，最近看到一段代码，自己研究了一段时间，感觉C语言真是强大，再加上这次国庆去北京，了解到C语言其实从来都不像大家人为的那样只是作为初级语言，入门语言供大家学习的，其实C语言的强大远远超出你我的想象。好吧，闲扯了一下。

大数阶乘：1000！

精简代码：

#include<stdio.h>
#define N 1000
long a[8916]={1,0},n,i,c,len;
int main(void) {
         n=N;
         for ( len=1; n>1; n--) {
                   for (c=0,i=0; i<len; i++ ){
                            a[i]= ( c+= a[i]*n )% 10000;
                            c/=10000;
                   }
 
                   ((a[i]=c)>0)?len++:0;
 
         }
         for(len--,printf("%d",a[len--]); len>=0; len--)printf("%04d",a[len]);
 
         return 0;
}

我的调试代码：

#include<stdio.h>
#define N 1000/* 最大到10555的阶乘 */
long a[8816] = {1,0 },len;
int main( void )
{
         n = N;
         for ( len = 1; n > 1; n-- )
         {
                   printf( "第1层for:len:%d,n:%dn",len,n );
 
                   for ( c = 0,i = 0; i <len; i++ )
                   {
                            printf( "第2层for参数:c:%d,i:%dn",i );
 
                            a[i]   = (c += a[i] * n) % 10000;
                            c       /= 10000;
                            printf( "第2层for:a[i]:%d,c:%dn",a[i],c );
                   }
                   ( (a[i] = c) > 0) ? len++: 0;
                   printf( "第2层for结束:len:%dn",len );
                   printf( "n" );
         }
 
 
/*
 * for( len--,printf("for内：%dn",a[len--]); len>=0; len--)
 * {
 * //          printf("打印：len:%d,a[len--]:%dn",a[len--]);
 *     printf("%04dn",a[len]);
 * }
 */
         for ( len--,printf( "%d",a[len--]); len >= 0; len-- )
                   printf( "%04d",a[len] );
         /* %04d  4代表着一共有几位，前面加0代表着输出时不够的位数用0填充 */
         return(0);
}
 
 
/* for (len=1;n>1; n--)     //把len的长度初始为1，因为数组中已经有一个元素了a[0]=1
 * {
 *  for(c=0,i=0; i<len;i++ ) //c初始为0，一开始还没运算哪来的进位，这是个内层循环，数组a中的len
 *                         //个元素都必须参与运算，都必须*n，这就是下面a[i]*n的来由
 *
 *    {
 *         a[i]= ( c+= a[i]*n ) % 10; c/=10;  //c是进位，不用说了，c+=a[i]*n，a每个元素与n
 *                                               //相乘必须考虑低位是否有进位，有就加起来
 *    }
 *
 * ((a[i]=c)>0)?len++:0;  //最后一个元素也有进位吗，有的话就在当前的元素的下个数组位置直接等于进位值，并
 *                         // 且数组的元素值要加1,没有的话什么都不干,光一个0什么都不是
 *
 * }
 */

精简版代码编译信息：

Compiling singlefile...

--------

- Filename:D:dev++99.cpp

- Compiler Name: TDM-GCC 4.9.2 64-bit Release

Processing C++source file...

--------

- C++ Compiler:E:Dev-CppMinGW64bing++.exe

- Command: g++.exe "D:dev++99.cpp" -o"D:dev++99.exe"?-I"E:Dev-CppMinGW64include"-I"E:Dev-CppMinGW64x86_64-w64-mingw32include"-I"E:Dev-CppMinGW64libgccx86_64-w64-mingw324.9.2include"-I"E:Dev-CppMinGW64libgccx86_64-w64-mingw324.9.2includec++"-L"E:Dev-CppMinGW64lib"-L"E:Dev-CppMinGW64x86_64-w64-mingw32lib" -static-libgcc

Compilationresults...

--------

- Errors: 0

- Warnings: 0

- OutputFilename: D:dev++99.exe

- Output Size: 163.51953125 KiB

- Compilation Time: 0.34s