P1009 阶乘之和

P1009 阶乘之和

• 题目提供者洛谷OnlineJudge
• 标签数论(数学相关)高精1998NOIp提高组NOIp普及组
• 难度普及-
• 通过/提交1139/3791

提交该题 讨论 题解 记录

题目描述

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

其中“！”表示阶乘，例如：5！=5*4*3*2*1。

输入输出格式

输入格式：

一个正整数N。

输出格式：

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

输入输出样例

输入样例#1：

3

输出样例#1：

9

#include<iostream>
#include<cstring>
using namespace std;
int a[],n,i,j,x[],y[];
void add()//处理进位
{
    memset(x,,sizeof(x));
    x[]=max(y[],a[]);
    for(i=;i<=x[];i++){x[i]+=y[i]+a[i],x[i+]=x[i]/,x[i]%=;}
    while(x[x[]+]>){x[x[]+]=x[x[]+]/,x[x[]+]%=,x[]++;}
    y[]=x[];
    for(i=;i<=x[];i++){y[i]=x[i];}
}
int main()
{
    cin>>n;
    a[]=,a[]=,y[]=,y[]=;
    for(j=;j<=n;j++)
    {
        memset(x,,sizeof(x));//清空阶乘数组
        x[]=a[];
        for(i=;i<=a[];i++){x[i]+=a[i]*j,x[i+]=x[i]/,x[i]%=;}//模拟阶乘
        while(x[x[]+]>){x[x[]+]=x[x[]+]/,x[x[]+]%=,x[]++;}//加上阶乘
        for(i=;i<=x[];i++){a[i]=x[i];}
        a[]=x[];
        add();
    }
    for(i=y[];i>=;i--){cout<<y[i];}//倒序存储就要倒序输出
    return ;
}