题解

#include<iostream>
#include<algorithm>
using namespace std;
long long f[60][150],sum[150];
void stoh(long long h[],string s){
	h[0]=s.size();
	reverse(s.begin(),s.end());
	for(int i=1;i<=s.size();i++){
		h[i]=s[i-1]-'0';
	}
}
void mul(long long a[],long long b,long long c[]){
	c[0]=a[0];
	for(int i=1;i<=c[0];i++){
		c[i]+=a[i]*b;
		c[i+1]+=c[i]/10;
		c[i]%=10;
	}
	while(c[c[0]+1]>0){
		c[0]++;
		c[c[0]+1]=c[c[0]]/10;
		c[c[0]]%=10;
	}
}
void add(long long a[],long long c[]){
	c[0]=max(c[0],a[0]);
	for(int i=1;i<=c[0];i++){
		c[i]+=a[i];
		c[i+1]+=c[i]/10;
		c[i]%=10;
	}
	while(c[c[0]]==0&&c[0]>1){
		c[0]--;
	}
	while(c[c[0]+1]>0){
		c[0]++;
		c[c[0]+1]=c[c[0]]/10;
		c[c[0]]%=10;
	}
}
void printh(long long sum[]){
	for(int i=sum[0];i>0;i--){
		cout<<sum[i];
	}
}
int main(){
	int n;
	cin>>n;
	stoh(f[0],"1");
	stoh(sum,"0");
	for(int i=1;i<=n;i++){
		mul(f[i-1],i,f[i]);
		add(f[i],sum);
	}
	printh(sum);
	return 0;
}