P1009 [NOIP1998 普及组] 阶乘之和

#include<bits/stdc++.h>
using namespace std;
int n,a[100],b[100],c[100],result[100],len_a,len_b=1,len_c=1,len_ans,m=1;
// a数组存储i   b数组存储（i-1)的阶乘  c数组存储i的阶层  result数组存储累加的和  m存储result数组的长度
int main()
{
    cin>>n;
    b[0]=1;
    for(int i=1;i<=n;i++){ //计算i的阶乘，已经算好了i-1的阶乘
        len_a=0;
        int temp=i;
        while (temp>0){ //把i存进a数组
            a[len_a++]=temp%10;
            temp/=10;
        }
        for(int j=0;j<len_a;j++)  //计算a*b（i*（i-1）的阶乘），即i的阶乘
            for(int k=0;k<len_b;k++)
                c[j+k]+=a[j]*b[k]; // 累加之前的结果 这里wa了一下
        for(int j=0;j<len_c;j++){
            if(c[j]>9){
                c[j+1]+=c[j]/10;
                c[j]=c[j]%10;
            }
        }
        len_ans=len_a+len_b; //用len_ans暂存len_a+len_b的值
        while (c[len_c]==0&&len_c>0)len_c--; //去除数字首多余的0
        if(c[len_c]) len_c++; //看最高位要不要进位
        len_b=len_c;m=max(m,len_c);len_c=len_ans;
        for(int j=len_b-1;j>=0;j--) b[j]=c[j];
        memset(c,0,sizeof(c));
        for(int j=0;j<m;j++){ //高精度加法
            result[j]+=b[j];
            if(result[j]>9){
                result[j+1]+=result[j]/10;
                result[j]%=10;
            }
        }
        if(result[m]!=0) m++; //有进位则修改长度
    }
    while(!result[m]&&m>0) m--;
    for(int i=m;i>=0;i--) cout<<result[i]; //倒序输出
    return 0;
}