18.7.29 luogu P1009 阶乘之和[高精加法&乘法]

题目描述

用高精度计算出 S=1!+2!+3!+…+n! (n≤50)S=1!+2!+3!+…+n!(n≤50)

其中“!”表示阶乘，例如： 5!=5 imes 4 imes 3 imes 2 imes 15!=5×4×3×2×1 。

输入输出格式

输入格式：

 

一个正整数 NN 。

 

输出格式：

 

一个正整数 SS ，表示计算结果。

 

输入输出样例

输入样例#1： 

3

输出样例#1： 

9

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <math.h>
#include <iostream>

using namespace std;

struct num {
    int a[100];
    num(){ memset(a, 0, sizeof(a)); }
    num(int x) {
        memset(a, 0, sizeof(a));
        int i = 0;
        while (x) {
            i++;
            a[i] = x % 10;
            x /= 10;
        }
        a[0] = i;
    }
    num operator +(num p) {
        int l = max(a[0], p.a[0]);
        num ans;
        for (int i = 1; i <= l; i++) {
            ans.a[i] += a[i] + p.a[i];
            if (ans.a[i] > 9) {
                ans.a[i] -= 10;
                ans.a[i + 1]++;
            }
        }
        if (ans.a[l + 1] != 0)l++;
        ans.a[0] = l;
        return ans;
    }
    num operator *(num p) {
        num ans;
        for(int i=1;i<=a[0];i++)
            for (int j = 1; j <= p.a[0]; j++) {
                ans.a[i + j - 1] += a[i] * p.a[j];
                ans.a[i + j ] += ans.a[i + j - 1] / 10;
                ans.a[i + j - 1] %= 10;
                ans.a[0] = max(ans.a[0], i + j - 1);
                if (ans.a[i + j ] != 0)
                    ans.a[0] = max(ans.a[0], i + j );
            }
        return ans;
    }
    void print() {
        printf("%d", a[a[0]]);
        for (int i = a[0] - 1; i >= 1; i--)
            printf("%d", a[i]);
        printf("
");
    }
};

int main()
{
    int n;
    scanf("%d", &n);
    num tmp(1),ans(0);
    for (int i = 1; i <= n; i++) {
        tmp=tmp*i;
        ans = ans + tmp;
    }
    ans.print();
    return 0;
}

可当模板