题目链接:阶乘之和
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;
}