阶乘计算——蓝桥杯基础练习

问题描述

	输入一个正整数n，输出n!的值。
	其中n! = 1*2*3*...*n。

输入格式

	输入包含一个正整数n，n<=1000。

输出格式

	输出n!的准确值。

样例输入

	10

样例输出

	3628800

问题分析

	当看到阶乘计算问题时，首先想到的就是使用for循环进行累乘，然而再回头看
    一下数据规模与约定，显然就算用64位的无符号长整型也存不下如此庞大的数据。这
    时我们可以运用数组，进行分块存储运算结果。
    	这里我选择方式的是数组一个空间存放一位十进制数。

代码

#include <stdio.h>

#define Max_Length 3000  //数组最大长度

int length_num(int *k)   //返回当前数组有效长度的下标
{
     int i = Max_Length - 1;
     while(k[i] == 0)
         i--;
     return i;
}

int main()
{
     int a[Max_Length] = {1,0}; //这样初始化省去了计算1! 
     int n,i,j,len;    //i，j为循环变量；len存放数组有效长度的最高下标 
     scanf("%d",&n);
     if(n > 0)      //问题中规定 n 为正整数，所以忽略0的阶乘。 
     {
          for(i = 2;i <= n;i++)
          {
               len = length_num(a);     //获得当前数组有效长度的最高下标 
               for(j = 0;j <= len; j++)   
               a[j] *= i;     //对每一个数据乘以i， 
               for(j = 0;j <=len; j++)
               {
                    if(a[j] >= 10)        //如果有数据大于10，则需要进位 
                    {
                         a[j+1] += a[j]/10%10;   //求a[j]的十位，同时在a[j]高一位处加上这个值
                         a[j+2] += a[j]/100%10;  //求a[j]的百位，同时在a[j+1]高一位处加上这个值
                         a[j+3] += a[j]/1000%10; //求a[j]的千位，同时在a[j+2]高一位处加上这个值
                         a[j] = a[j]%10;         //求a[j]的个位，并赋给a[j]
                    }
                    len = length_num(a);        //获得当前数组有效长度的最高下标
               }
          }
          len = length_num(a);  //获得当前数组有效长度的最高下标
          while(len != EOF)   //从最高位开始打印 
          printf("%d",a[len--]);
          printf("\n");
     }
     return 0;
}

代码分析

数组最大长度的确定

由系统自带的计算器算出1000!的位数是2568位，所以这里取了3000，稍微大了一些。
数组中存放的数据，从下标为0开始，依次为个位、十位、百位、千位………。

length_num函数

 在计算的过程中需要频繁运用数组实时的有效长度，所以单独定义了函数用于计算其长度。

运算过程

    第一层循环中，首先求得有效长度，再进行内嵌的第一个循环，对数组的每一个数据
乘以当前所求的阶数值，这样做的原理是乘法分配律，对前一次求得的阶乘值中每一位数
据进行拆解
例如：120 × 6 = (100 + 20 + 0) × 6，其中100所对应的是数组中的第3个元素1，即
前一次算得的阶乘值的百位；同理20对应的是数组的第2个元素2，0对应的是数组的第1个
元素0。
    如此乘完以后，有的数据会超过10，而我们数组中每一个数据仅代表十进制中的一位
数，不能超过10，所以此时需要进位，接下来执行内嵌的第二个循环。
    如果数据大于10，则依次计算十位、百位、千位、个位。计算四位是因为数组中存放
的最大值为9，同时要求的最大阶乘为1000，所以数组中会存放的最大数据为9000。个位
放在最后计算，如果首先就计算个位并且赋值给自己，就会破坏当前的数据，使得接下来
计算的十位、百位、千位都为0，首先计算个位那么则需要多定义一个变量来存放当前处理
的数据。
    当我们进位完成之后，需要重新确定数组有效长度，因为可能会出现如下这种情况：

数组下标	…	20	21	22
数据	…	72	108	180

其中22为当前最高位，21为次高位，20为次次高位。（这只是一个举例，可能与实际不符）
进位后：

数组下标	…	21	22	23	24
数据	…	5	0	10	1

此时新增了两位23和24，如果此时不重新获取数组有效长度，而如果这恰好又是最后一次
计算阶乘，那么打印出来的结果为10105...，而实际为2005...导致结果错误。
    同理在最后打印结果之前也应该重新获取当前数组长度，因为最后一次进位可能会有
新出现的更高位数据，不重新获取会漏打印，也导致结果出错。

写在最后

这种方法完全是作用于当前n所设定的范围内设计的，同时由于频繁的计算数组长度，时间效率也不高，不过时间效率的问题可以通过增加if语句来减少时间的耗费；空间上，由于变量定义的比较随意，针对此题还可以减小很多消耗。总之，还有很多不足之处。