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 ;

}