计算1+2!+3!+...+n!的和
Code
点击查看代码
#include<iostream>
#include<vector>
using namespace std;
typedef vector<int> VI; //动态数组
int n;
VI add(VI &A, VI &B) //高精度加法
{
if (A.size() < B.size()) return add(B, A); //使A长度大于B
VI C; //接收计算结果
int t=0; //保存进位
for (int i = 0; i < A.size(); i ++ ) //以长度长的作循环
{
t += A[i];
if (i < B.size()) t += B[i]; //如果B没有加完
C.push_back(t % 10); //结果mod 10
t /= 10; //进位
}
if (t) C.push_back(t); //结果长度可能长于长度最长的数
return C; //返回结果
}
VI mul(VI &A,int b){ //高精度与低精度乘法
VI C; //接收结果
int t = 0; //保存进位
for(int i = 0; i <= A.size()-1 || t;i++){ //如果没有乘完或者进位不为0
if(i <= A.size()-1)t += A[i] * b; //相乘
C.push_back(t % 10); //结果mod 10
t /= 10; //进位
}
while(C.size() > 1&&C.back() == 0)C.pop_back(); //除去前导0(可能高进度乘0)
return C; //返回结果
}
int main(){
VI A; //每次保存阶乘
VI C; //保存阶乘之和
C.push_back(0); //初始化为0
cin >> n;
for(int i = 1; i <= n; i ++){
A.erase(A.begin(),A.end()); //每次用完要清空
A.push_back(1); //初始化为1再求阶乘
for(int j = i;j >= 1; j--){
A = mul(A,j); //求阶乘
}
C = add(C,A); //阶乘之和
}
for(int i = C.size() - 1; i >= 0; i -- )cout << C[i]; //逆序输出
return 0;
}