信息学奥赛一本通:1173:阶乘和

1173:阶乘和


时间限制: 1000 ms         内存限制: 65536 KB
提交数: 17291     通过数: 8783 

【题目描述】

用高精度计算出S=1!+2!+3!+…+n!(n≤100),其中“!!”表示阶乘,例如:5!=5×4×3×2×1。

输入正整数n,输出计算结果S。

【输入】

一个正整数n。

【输出】

计算结果S。

【输入样例】

5

【输出样例】

153

第一种:

在a数组中做阶乘,然后加在b数组中。

#include<iostream>
#include<cmath> 
#include<cstring>
int a[100000005],b[100000005],n;
using namespace std;
int main()
{
	int i,j,k,m,t;
     a[0]=1;
     k=1;
	scanf("%d",&m);
	for(t=1;t<=m;t++)//从1的阶乘开始到输入的数。 
	{
		
	for(i=1;i<=t;i++)//每数的阶乘放在a中 
	{
		for(j=0;j<k;j++)//从0开始计算第个数位乘i 
		{
		//printf("%d\n ",m);
			a[j]=a[j]*i;
		}
		for(j=0;j<k;j++)//进位算,大于等于10时进位,取剩下的模。 
		{
			if(a[j]>=10)
			a[j+1]+=a[j]/10;
			a[j]=a[j]%10;
		}
		while(a[k]>0)//看最后进位了吗 
		{
		a[k+1]=a[k]/10;
		a[k]=a[k]%10;
			k++;
	     }
    }
	
		for(int x=0;x<k;x++)//从一的阶乘放加在b数组中,把a[0]放1,
//其它位清零,为下次阶乘准备。 
		{
		b[x]+=a[x];
		if(x==0)
		a[x]=1;
		else
		a[x]=0;
		}
	}
	    for(j=0;j<k;j++)//b中的数有进位的吗? 
		{
			if(b[j]>=10)
			b[j+1]+=b[j]/10;
			b[j]=b[j]%10;
		}
	   while(b[k]>0)//处理最后的数位
		{
	
		b[k+1]=b[k]/10;
		b[k]=b[k]%10;
			k++;
	    }
	for(i=k-1;i>=0;i--)//从高位输出每位上的数
	printf("%d",b[i]);
	return 0;
}

第二种:从网上找的,没注释啊,我不是很明白那

#include<iostream>
#include<cmath> 
#include<cstring>
using namespace std;
const int maxn = 1e5 + 5;
int n;
int a[maxn] = {0,1},s[maxn];
int num = 1,alen = 1,jin = 1,slen = 1;
int i,j;
int main(){
    ios::sync_with_stdio(0);
    cin >> n;
    for(i = 1; i <= n; i++){
        for(j = 1; j <= alen; j++)
            a[j] *= num;
        for(j = 1; j <= alen; j++) {
            if (a[j] >= 10) {
                a[j + 1] += a[j] / 10;
                a[j] %= 10;
                jin++;
            }
        }
        alen++;
        while(a[alen] >= 10){
            a[alen + 1] = a[alen] / 10;
            a[alen] %= 10;
            alen++;
            jin++;
        }
        num++;
        slen = max(slen,alen);
        for(j = 1; j <= slen; j++){
            s[j] += a[j];
            if(s[j] >= 10){
                if(j == slen) slen++;
                s[j + 1] += s[j] / 10;
                s[j] %= 10;
            }
        }
        alen = jin;
    }
    for( i = slen; i > 1; i--){
        if(s[i])
            break;
    }
    for(j = i; j>= 1; j--)
        cout << s[j];
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值