noip1998 高精度计算 （高精度）

A1115. 高精度计算

时间限制：1.0s   内存限制：256.0MB  
总提交次数：653  
AC次数：250   平均分：51.49

将本题分享到：

      
   

查看未格式化的试题   提交   试题讨论

试题来源

　　NOIP1998 普及组

问题描述

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

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

　　输入正整数N，输出计算结果S。

输入格式

　　输入包含一个整数N（N<=50）

输出格式

　　输出为一个整数，为计算结果S。

解析：高精度。

代码：

#include<cstdio>
#include<algorithm>
using namespace std;

const long long len=1e8;
const int maxn=200;
long long a[maxn],b[maxn];

void multip(int x)
{
long long i,last=0;
for(i=1;i<=b[0];i++)
{
b[i]=b[i]*x+last;
last=b[i]/len,b[i]%=len;
}
if(last>0)b[++b[0]]=last;
}

void add()
{
long long i,k,last=0;
k=min(a[0],b[0]);
for(i=1;i<=a[0];i++)
{
a[i]=a[i]+b[i]+last;
last=a[i]/len,a[i]%=len;
}
for(i=k+1;i<=b[0];i++)
a[i]=b[i]+last,last=a[i]/len,a[i]%=len;
for(i=k+1;i<=a[0] && last>0;i++)
a[i]+=last,last=a[i]/len,a[i]%=len;
a[0]=max(a[0],b[0]);
if(last>0)a[++a[0]]=last;
}

int main()
{
int i,n;
scanf("%d",&n);
a[0]=a[1]=b[0]=b[1]=1;
for(i=2;i<=n;i++)multip(i),add();

printf("%I64d",a[a[0]]);
for(i=a[0]-1;i>=1;i--)printf("%08I64d",a[i]);
return 0;
}