阶乘和

描述

用高精度计算出S=1!+2!+3!+…+n!（n≤50）,其中“!”表示阶乘，例如：5!=5×4×3×2×1。

输入正整数n，输出计算结果S。

格式

输入格式

一个正整数n。

输出格式

计算结果S。

样例

输入样例

5

输出样例

153

#include <iostream>
#include <stdio.h>
#include <cstring>  
#include <stdio.h>
#include <algorithm>
#include <cmath>
using namespace std;
//#define N 10001
char a[101];
char b[101];
char c[501];
char d[501];
char abc[501];
int cnt;
int acnt;

void add(char a[],char b[], int len){
		
	int flag = 0;
	cnt = 0;
	
	for(int i=0; i<len; ++i){
		if(a[i] == '\0'){
			a[i] = 0;
		}
		if(b[i] == '\0'){
			b[i] = 0;
		}
		
		int sum = (a[i]  + b[i] ) + flag;
		flag = sum/10;
		sum = sum%10;
		d[cnt++] = sum;
	}
	if(flag>=1){
		d[cnt++] += flag;
	}
//	printf("add = ");
//	for(int i=cnt-1; i>=0; --i)	{
//		printf("%d", d[i]);
//	}
//	printf("\n");
}

void Reverse(char arr[], int len){
	
	for(int i=0; i<len/2; ++i){
		char ch = arr[i];
		arr[i] = arr[len - i - 1];
		arr[len - i - 1] = ch;		
	}	
}



int chen(char a[],char b[], int lena, int lenb){
	
	for(int i=0; i<500; ++i){
		d[i] = 0;
	}
	
	int flag = 0, i, j;
//	printf("%d %d\n",lena, lenb);
	for(j=0; j<lenb; ++j){
		
		flag = 0;
		int cnt = 0;
		for(i=0; i<lena; ++i){
			
//			printf("%c %c %d\n", a[i], b[j], flag);
			int t = ( a[i] - '0' ) * (b[j] - '0') + flag;
			
			flag = t/10;
			t = t%10;
			c[cnt++] = t;
//			printf("t = %d %d\n",t,flag);
			
		}
		if(flag >= 1){
			c[cnt++] = flag;
		}
		
		
		for(int k=cnt-1 + j; k>=j; --k){
			c[k] = c[k-j];
		}
		for(int k=0; k<j; ++k){
			c[k] = 0;
		}
		cnt += j;
		
		add(c, d, cnt);

//		printf("res = ");
//		for(int k=cnt-1; k>=0; --k){
//			printf("%d ", c[k]);
//		}
//		printf("\n");
		
	}
	
}

void fun(int x){
		
//	printf("x = %d\n",x);
	int lena;
	int lenb;
	if(x<10){
		lenb = 1;
	}else{
		lenb = 2;
	}
	int y = 0;
	while(x){
		int t = x%10;
		x = x/10;
		b[y++] = t + '0';
	}
	
	lena = cnt;
	if(lena == 0){
		lena = 1;
	}
//	for(int i=0; i<lena; ++i){
//		printf("%c",a[i]);
//	}
//	printf("\n");
//	for(int i=0; i<lenb; ++i){
//		printf("%c",b[i]);
//	}
//	printf("\n");
	
	chen(a, b, lena, lenb);
	for(int i=0; i<cnt; ++i){
		a[i] = d[i] + '0';
	}
	
	
}

int main(){	
	int n;
	char res[500];
	for(int i=0; i<500; ++i){
		abc[i] = 0;
	}
	
	cin>>n;
	a[0] = '1';
	for(int i=1; i<=n; ++i){
		fun(i);
		int j;
		for(j=0; j<cnt; ++j){
			res[j] = d[j];
		}
		res[j] ='\0';
		add(abc, res, cnt);
		for(j=0; j<cnt; ++j){
			abc[j] = d[j];
		}
		
//		printf("res = ");
//		for(int i=cnt-1; i>=0; --i)	{
//			printf("%d", d[i]);
//		}
//		printf("\n");		
//		printf("res_1 = ");		
	}	
	for(int i=cnt-1; i>=0; --i)	{
		printf("%d", abc[i]);
	}
	printf("\n");
	
	
	
	return 0;

}