P1009 [NOIP1998 普及组] 阶乘之和

题目描述

用高精度计算出 S = 1! + 2! + 3! + \cdots + n!S=1!+2!+3!+⋯+n!（n \le 50n≤50）。

其中 ! 表示阶乘，定义为 n!=n\times (n-1)\times (n-2)\times \cdots \times 1n!=n×(n−1)×(n−2)×⋯×1。例如，5! = 5 \times 4 \times 3 \times 2 \times 1=1205!=5×4×3×2×1=120。

输入格式

一个正整数 nn。

输出格式

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

输入输出样例

输入 

3

输出 

9

说明/提示

【数据范围】

对于 100 \%100% 的数据，1 \le n \le 501≤n≤50。

大家都能想到这种代码

#include <iostream>

using namespace std;

int main()
{
    int n=0;
    cin>>n;
    int ret=1,sum=0;
    for(int j=0;j<=n;j++)
    {
     ret=1;
     for(int i=1;i<=j;i++)
     {
      ret*=i;
     }
     sum+=ret;
     }
    cout<<sum-1;
    return 0;
}

#include <iostream>

using namespace std;

int main()
{
    int n=0;
    cin>>n;
    int ret=1,sum=0;
    for(int j=0;j<=n;j++)
    {
     ret=1;
     for(int i=1;i<=j;i++)
     {
      ret*=i;
     }
     sum+=ret;
     }
    cout<<sum-1;
    return 0;
}

然后，就悲剧了

 50分？？？

50! > 2^64-1

超了...........

高精度代码如下

#include<bits/stdc++.h>

using namespace std;

const int t=101;
int a[t+1],b[t+1],n,s;
int main()
{
    a[0]=b[0]=1;
    cin>>n;
    for(int i=2;i<=n;i++)
    {
        for(int j=0;j<=t;j++)
            b[j]*=i;
        for(int j=0;j<=t;j++)
        {
            if(b[j]>=10)
            {
                b[j+1]+=b[j]/10;
                b[j]%=10;
            }
        }
        for(int z=0;z<=t;z++)
        {
            a[z]+=b[z];
           if(a[z]>=10)
           {
              a[z+1]+=a[z]/10;
              a[z]%=10;
            }
        }
    }
    for(int i=t-1;i>=0;i--)
        if(a[i]!=0)
        {
            s=i;
            break;
        }
        for(int j=s;j>=0;j--)
            cout<<a[j];
        return 0;
}