P1009 [NOIP1998 普及组] 阶乘之和

题目描述

用高精度计算出 S = 1! + 2! + 3! +…+ n!=1!+2!+3!+⋯+n!
（n ≤50n≤50）。

其中“!”表示阶乘

输入格式

一个正整数 n。

输出格式

一个正整数 SS，表示计算结果。

分析一下本题，看到n的范围就知道不能使用正常的计算方式，因为会超范围。所以就要使用高精度计算，主要思路就是使用字符串来进行计算，计算逻辑就模拟手算时的逻辑。
以下是代码，可以作为大数运算的模板

#include<bits/stdc++.h>
#define mem(arr,p) memset(arr,p,sizeof(arr))
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
bool space;
inline void print(){if(space)cout<<" ";space=true;}
const int maxn =1e3;
const int INF=0x3f3f3f3f;
const double PI=acos(-1);

string toString(int n){
    stack<char> tmp;
    while (n){
        tmp.push(n%10+'0');
        n/=10;
    }
    string ans="";
    while (!tmp.empty()){
        ans.append(1,tmp.top());
        tmp.pop();
    }
    return ans;
}
string Add(string l,string r){
    if(l.size()<r.size()) swap(l,r);
    if(r.size()==0) return l;
    int i=0,flag=0;
    string ans="";
    while (i<r.size()){
        int sum=(l[i]-'0')+(r[i]-'0')+flag;
        if(sum>=10){
            flag=1;
            sum-=10;
        }
        else flag=0;
        ans.append(1,sum+'0');
        i++;
    }
    while (i<l.size()){
        int sum=(l[i]-'0')+flag;
        if(sum>=10){
            flag=1;
            sum-=10;
        }
        else flag=0;
        ans.append(1,sum+'0');
        i++;
    }
    if(flag) ans.append(1,'1');
    return ans;
}
string Mul(string l,string r){
    reverse(r.begin(),r.end());
    if(l.size()>r.size()) swap(l,r);
    string ans="";
    for(int i=0;i<l.size();i++){
        string tmp="";
        tmp.append(i,'0');
        int flag=0;
        for(int j=0;j<r.size();j++){
            int sum=(l[i]-'0')*(r[j]-'0')+flag;
            if(sum>=10){
                flag=sum/10;
                sum=sum%10;
            }
            else flag=0;
            tmp.append(1,sum+'0');
        }
        while (flag>0){
            tmp.append(1,flag%10+'0');
            flag/=10;
        }
        ans=Add(ans,tmp);
    }
    return ans;
}
int main(){
    int n;
    cin>>n;
    string l="1",ans="1";
    for(int i=2;i<=n;i++){
        l=Mul(l,toString(i));
        ans=Add(ans,l);
    }
    reverse(ans.begin(),ans.end());
    cout<<ans<<endl;
}