大整数类实现阶乘

#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
const int maxn = 10000;

struct BigInt
{
	BigInt(int x = 0)
	{
		memset(arr, 0, sizeof(arr));//初始化全为0
		for (len = 1; x; len++)
		{
			arr[len] = x % 10;//对x逆序保存
			x /= 10;
		}
		len--;
	}

	int& operator[](int i)
	{
		return arr[i];//用x[i]代替x.arr[i]
	}

	void flatten(int L)//把不是一位数的位处理成一位数,故称展平,L要不小于有效长度
	{
		len = L;
		for (int i = 1; i <= len; ++i)
		{
			arr[i + 1] += arr[i] / 10;
			arr[i] %= 10;
		}
		while (!arr[len])//消除前置0
		{
			len--;
		}
	}

	void print()//反向输出
	{
		for (int i = max(1, len); i >= 1; --i)
		{
			cout << arr[i];
		}
		cout << endl;
	}

	int len;
	int arr[maxn];
};

BigInt operator+(BigInt &a,BigInt &b)
{
	BigInt c;
	int len = max(a.len, b.len);
	for (int i = 1; i <= len; ++i)
	{
		c[i] += a[i] + b[i];
	}
	c.flatten(len + 1);//最大长度不超过len+1
	return c;
}

BigInt operator*(BigInt& a, BigInt& b)
{
	BigInt c;
	int len = a.len + b.len;
	for (int i = 1; i <= a.len; ++i)
	{
		for (int j = 1; j <= b.len; ++j)
		{
			c[i + j - 1] += a[i] * b[j];//a[i] * b[j]的结果产生在i+j-1位上
		}
	}
	c.flatten(len );//最大长度不超过len+1
	return c;
}

BigInt operator*(BigInt& a, int& b)
{
	BigInt c;
	int len = a.len;
	for (int i = 1; i <= len; ++i)
	{
		c[i] += a[i] *b;
	}
	c.flatten(len + 11);//int类型最长10位,所以可以这样展平
	return c;
}

int main()
{
	int n;
	cin >> n;
	BigInt fac(1);
	BigInt ans(0);
	for (int i = 1; i <= n; ++i)
	{
		fac = fac * i;
		ans = ans + fac;
	}
	ans.print();

	return 0;
}