【洛谷训练】阶乘之和

题目链接:阶乘之和

555555我太菜了,高精度已经要让我崩溃了。对于这道题,我想到要结合前两次的A+B和A*B的高精度算法,也做了融合,但是不知道为什么还是在大数上有问题。参考大佬的解法,发现方法其实更简单。(但是我还没有发现为什么自己过不了后两个点)

#include<iostream>
#include<cstring>
using namespace std;
int n,a[90],b[90],c[90],f[90],d=0,len_a,len_b=1,len_c=1,len_ans,m=1;
string s;
int main(){
    cin>>n;
    b[0]=1; //初始化
    for(int i=1;i<=n;i++){ //计算i的阶乘,已经算好了i-1的阶乘
        len_a=0; //i的长度
        int p=i;
        while(p>0){ //把i存进a数组
            a[len_a++]=p%10;
            p/=10;
        }
        for(int j=0;j<len_a;j++) //计算a*b(i*(i-1)的阶乘),即i的阶乘,看不懂的网上查,我也不知道为什么
            for(int k=0;k<=len_b;k++)
                c[j+k]+=a[j]*b[k];
        for(int j=0;j<len_c;j++) //需要进位的就进位
            if(c[j]>9) c[j+1]+=c[j]/10,c[j]%=10;
        if(c[len_c]) len_c++; //看最高位要不要进位
        len_ans=len_b,len_b=len_c,m=max(m,len_c); //把len_b赋值给len_ans,修改len_b的值,m为i阶乘的长度,看有没有进位
        for(int k=len_c-1;k>=0;k--) b[k]=c[k]; //把c存进b数组,即存进i的阶乘,下次循环b为i-1的阶乘
        len_c=len_a+len_ans;
        memset(c,0,sizeof(c)); //清零c数组,准备计算下个阶乘
        for(int j=0;j<m;j++){ //高精加,直接套模板
            f[j]+=b[j];
            if(f[j]>9) f[j+1]+=f[j]/10,f[j]%=10; //进位,注意不要写成f[j+1]++,f[j]-=10;就因为这里wa了一个点
        }
    }
    while(!f[m]&&m>0) m--; //去掉首导零
    for(int i=m;i>=0;i--) cout<<f[i]; //倒序输出
    return 0; //圆满结束
}

问了朋友之后,发现了自己的一些问题。朋友的解法其实和大佬的解法差不多,我要自己再打一遍。比较不一样的是加法那里,借鉴!
弄了好久终于打出来了…天哪我实在太笨了…一直漏东漏西。

#include<iostream>
#define MAXLINE 105
using namespace std;
int num1[MAXLINE]={1},num2[MAXLINE],num3[MAXLINE],num4[MAXLINE],sum[MAXLINE];
int len1=1,len2=0,len3=1,len=0;
void mul(int k)
{
	//change int(k) to int[] k-num2[]
	len2=0;
	while(1){
		num2[len2++]=k%10;
		k/=10;
		if(k==0) break;
	} 
	for(int i=0;i<len1;i++){
		for(int j=0;j<len2;j++){
			num3[i+j]+=num1[i]*num2[j];
		}
	}
	for(int i=0;i<len1+len2-1;i++){
		if(num3[i]>9){
			num3[i+1]+=num3[i]/10;
			num3[i]%=10;
		} 
	}
	len3=len1+len2;
	while(len3>1&&num3[len3-1]==0) len3--;
	len1=len3;
	for(int i=0;i<len3;i++)	num1[i]=num3[i];
	for(int i=0;i<MAXLINE;i++) num3[i]=0;
}
void add()
{
	for(int i=0;i<MAXLINE-1;i++){
		if(i<len1)	num4[i]+=num1[i];
		if(i<len)	num4[i]+=sum[i];
		num4[i+1]+=num4[i]/10;
		num4[i]%=10;
	}
	len=MAXLINE;
	while(len>0&&num4[len-1]==0) len--;
	for(int i=0;i<len;i++)	sum[i]=num4[i];
	for(int i=0;i<MAXLINE;i++) num4[i]=0;
}
int main(void)
{
	int num;
	cin>>num;
	for(int i=1;i<=num;i++){
		mul(i);
		add();
	}	
	for(int i=0;i<len;i++)	cout<<sum[len-i-1];
	cout<<endl;
	return 0;
} 
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值