大数的阶乘
没办法,今天是贼不爽的一天,就这个大数阶乘搞了我一天,搞得我头昏脑胀,所以今天必须得吐槽一下。不得不说,要做出这题,逻辑思维一定要缜密。废话少说,直接上题。
题目描述
输入一个正整数n,输出n!的值。
其中n!=123*…*n。
算法描述
n!可能很大,而计算机能表示的整数范围有限,需要使用高精度计算的方法。使用一个数组A来表示一个大整数a,A[0]表示a的个位,A[1]表示a的十位,依次类推。
将a乘以一个整数k变为将数组A的每一个元素都乘以k,请注意处理相应的进位。
首先将a设为1,然后乘2,乘3,当乘到n时,即得到了n!的值。
输入
输入包含一个正整数n,n< =1000。
输出
输出n!的准确值。
样例输入
10
样例输出
3628800
题目已经提示的很清楚了,要做出这题,就必须得用数组来计数相乘N。否则我也没想到其他的法子。解题思想就是你把每次计算出的结果用数组存储,a[0]存个位,a[1]存十位,a[2]存百位~~~~~~再把存取的结果也就是数组的每一位分别与N相乘,在把每一位相乘结果相加又得到一个新的数组,循环这个过程,直到乘以n;一开始想的就是如何让这些结果相加,然后发现了一个规律,b[i] = b[i] + (a[j]/(int)pow(10,cnt) ) % 10,也就是说,假如你要得到百位结果a[2],就可以用a[2]%10+a[1]/10 + a[0]/100 % 10;通过这种规律就可以写出相加的代码。
c++代码如下:
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int a[100001];
int b[100001];
int n;
int flag;
int length(int num)
{
int len = 0;
if(num == 0)
return 1;
while(num)
{
num /= 10;
len ++;
}
return len;
}
int pow(int num)
{
int sta = 1;
while(num)
{
sta = sta * 10;
num --;
}
return sta;
}
void mult(int N,int len)
{
if(N == n + 1)
{
int temp = 0;
for(int i = len - 1; i >= 0; i --)
{
if(a[i] != 0)
temp = 1;
if(temp == 1 || a[i] != 0)
cout << a[i];
}
return;
}
for(int i = 0; i < len; i ++)
{
a[i] = a[i] * N;
}
len = length(a[len - 1]) + len + 1;
int temp = 0;
for(int i = 0; i < len ; i ++)
{
int cnt = i;
int flag;
if(i - 5 >= 0)
{
cnt = 5;
flag = i - 5;
}
else flag = 0;
for(int j = flag; j <= i; j ++)
{
//b[i] = b[i] + (a[j]/(int)pow(10,cnt) + temp) % 10;
if(flag == 0)
b[i] = b[i] + (a[j]/pow(cnt)) % 10;
else
b[i] = b[i] + (a[j]/pow(cnt)) % 10;
cnt --;
}
b[i] += temp;
if(b[i] >= 10)
{
temp = b[i] / 10;
b[i] %= 10;
}else
temp = 0;
}
for(int i = 0; i < len; i ++)
{
a[i] = b[i];
}
memset(b,0,sizeof(b));
mult(N + 1,len);
return;
}
我这里有些东西做的不够精确,也只是为了过测试点,就比如结果的位数我是没有仔细考虑的,只能让位数增大,去除前面的0,得到结果,还有就是用了两个数组,但毕竟是为了刷题,也没有什么好法子了,直接AC就对了。