1.题目描述
用高精度计算出 S = 1! + 2! + 3! + ....+ n!(n ≤ 50)。
其中
!
表示阶乘,定义为n! = n × (n−1) × (n−2) × ⋯ ×1。例如,5! = 5 × 4 × 3 × 2 × 1 = 120。
2.输入格式
一个正整数。
3.输出格式
一个正整数x,表示计算结果。
4.输入输出样例
输入:3
输出:9
输入:22
输出:1177652997443428940313
说明/提示【数据范围】:对于100% 的数据,1 ≤ n ≤ 50。
5.关键点
1.数组a[]作为每次接受新(i)的值的数组
2.数组c[]赋值给b[],b[]作为上一个(i-1)阶乘的值,在新的一轮中再与a[]做高精乘。
(做到了保留(i-1)!的结果,省事儿了)
3.将赋值后的c[]用memset函数清空数组,接收下一次a[]*b[]的结果
3.f[]作为最后的结果接受数组,与每个阶乘的结果b[]做高精加,最后输出。
6.代码实现
#include<bits/stdc++.h>
using namespace std;
const int maxx = 100;
int main()
{
int a[maxx]= {},b[maxx]= {},c[maxx]= {},f[maxx]= {};
int len_a,len_b = 1,len_c = 1,t,m=1;
int n,i,j,k;
cin >> n;
b[0] = 1;//b[0]要赋初值
for(i = 1; i <= n; i++)
{
len_a = 0;
int p = i;
while(p > 0)
{
a[len_a++] = p % 10;//存储i到数组a[]
p = p / 10;
}
int jw;
for(k = 0; k < len_a; k++)
{
jw=0;
for(j = 0; j < len_b; j++)
{
c[k+j] = a[k]*b[j] + jw + c[k+j];
jw=c[k+j]/10;
c[j+k] = c[k+j]%10;
}
c[k+len_b] = jw;
}
len_c = len_a+len_b;
//错在了len_c的长度啊啊啊啊啊啊啊啊啊,检查好长时间呜呜呜呜
for(int test =len_c-1; test >=0; test--)
{
if(0 == c[i] && test > 1)
test--;
else
break;
}
t = len_b;
//c[]求出的两数之积赋给b[]继续运算。
for(k = len_c - 1; k >=0; k--)
{
b[k] = c[k];
}
len_b = len_c;
//高精加,m是和的长度
m = max(m,len_c)+1;//+1以防进位
len_c = len_a+t;//以防万一,取最长清零
memset(c,0,sizeof(c));
//加法运算
for(j = 0; j < m; j++)
{
f[j] = f[j]+ b[j];
if(f[j] > 9)
{
f[j+1] += f[j]/10;
f[j] = f[j]%10;
}
}
}
//删除前导0
for(i = m - 1; i >=0; i--)
{
if(0 == f[i] && m > 1)
m--;
else
break;
}
//倒序输出
for(i = m-1; i >=0; i--)
cout << f[i];
return 0;
}