luogu_1009 阶乘之和

#include<bits/stdc++.h>

using namespace std;

int a[1010],n,sum[1010],m[3],z[1010];

void mul(){

memset(z,0,sizeof(z));

z[0]=a[0]+m[0];

for(int i=1;i<=a[0];i++)

for(int j=1;j<=m[0];j++)z[i+j-1]+=a[i]*m[j];

for(int i=1;i<=z[0];i++){z[i+1]+=z[i]/10; z[i]%=10;}

while(!z[z[0]] && z[0]>1)z[0]--;

for(int i=0;i<=z[0];i++)a[i]=z[i];

}

void add(){

sum[0]=max(sum[0],z[0]);

for(int i=1;i<=sum[0];i++){

sum[i]+=z[i]; sum[i+1]+=sum[i]/10; sum[i]%=10;

}

while(sum[sum[0]+1])sum[0]++;

}

int main(){

cin>>n;

a[0]=sum[0]=1; a[1]=sum[1]=1;

for(int i=2;i<=n;i++){

int x=i; m[0]=0;

for(int j=1;x;j++){m[j]=x%10; x/=10; m[0]++;}

mul(); add();

}

for(int i=sum[0];i;i--)printf("%d",sum[i]);

puts("");

return 0;

}