0分求助

#include<bits/stdc++.h>
using namespace std;
const int N=5001;
int a[N]={0,1},b[N]={0,1};
int main()
{
	int n;
	scanf("%d",&n);
	for(int i=2;i<=n;i++)
	{
		for(int j=1;j<=N;j++)
			b[j]*=i;
		for(int j=1;j<N;j++)
			b[j+1]+=b[j]/10,b[j]%=10;
		for(int j=1;j<=N;j++)
		{
			a[j]+=b[j];
			a[j+1]+=a[j]/10,a[j]%=10;
		}
	}
	for(int i=N;i>=1;i--)
	{
		if(a[i]!=0)
		{
			for(int j=i;j>=1;j--)
				printf("%d",a[i]);
			break;
		}
	}
	return 0;
}