问题描述
输入一个正整数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};
int n,i,j,len;
scanf("%d",&n);
if(n > 0)
{
for(i = 2;i <= n;i++)
{
len = length_num(a);
for(j = 0;j <= len; j++)
a[j] *= i;
for(j = 0;j <=len; j++)
{
if(a[j] >= 10)
{
a[j+1] += a[j]/10%10;
a[j+2] += a[j]/100%10;
a[j+3] += a[j]/1000%10;
a[j] = a[j]%10;
}
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,首先计算个位那么则需要多定义一个变量来存放当前处理
的数据。
当我们进位完成之后,需要重新确定数组有效长度,因为可能会出现如下这种情况:
其中22为当前最高位,21为次高位,20为次次高位。(这只是一个举例,可能与实际不符)
进位后:
此时新增了两位23和24,如果此时不重新获取数组有效长度,而如果这恰好又是最后一次
计算阶乘,那么打印出来的结果为10105...,而实际为2005...导致结果错误。
同理在最后打印结果之前也应该重新获取当前数组长度,因为最后一次进位可能会有
新出现的更高位数据,不重新获取会漏打印,也导致结果出错。
写在最后
这种方法完全是作用于当前n所设定的范围内设计的,同时由于频繁的计算数组长度,时间效率也不高,不过时间效率的问题可以通过增加if语句来减少时间的耗费;空间上,由于变量定义的比较随意,针对此题还可以减小很多消耗。总之,还有很多不足之处。