洛谷P1009 [NOIP1998 普及组] 阶乘之和（高精度）//高精度乘法、高精度加法

//学习记录

//洛谷P1009 [NOIP1998 普及组] 阶乘之和

#include
#include
using namespace std;
void jia(int p[],int q[],int m[])
{
int a[205],b[205],c[205],x=0;
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
memset(c,0,sizeof©);
for(int i=0;i<205;i++)
{
a[i]=p[i];
b[i]=q[i];
c[i]=m[i];
}
c[0]=max(a[0]+1,b[0]+1);
for(int i=1;i<=c[0];i++)
{
c[i]=a[i]+b[i]+x;
x=c[i]/10;
c[i]%=10;
}
while(c[c[0]]==0&&c[0]>1) c[0]–;
for(int i=0;i<=c[0];i++)
{
m[i]=c[i];
}
}
void multiple(int p[],int q[],int m[])
{
int a[205],b[205],c[205],x=0;
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
memset(c,0,sizeof©);//乘法未初始化c数组，几个小时卡在这个点没看出来
for(int i=0;i<205;i++)//i=0
{
a[i]=p[i];
b[i]=q[i];
}
c[0]=a[0]+b[0]+1;
for(int i=1;i<=a[0];i++)
{
for(int j=1;j<=b[0];j++)
{
c[i+j-1]=a[i]*b[j]+x+c[i+j-1]; //若未初始化c，则此处出错
x=c[i+j-1]/10;
c[i+j-1]%=10;
}
if(x!=0)
{
c[i+b[0]]+=x;
x=0;
}
}
while(c[c[0]]==0&&c[0]>1) c[0]–;
for(int i=0;i<=c[0];i++)
{
m[i]=c[i];
}
}

void put(int r[])
{
for(int i=r[0];i>=1;i–) cout<<r[i];
cout<<endl;
}

int s[205],a[205],b[205];
int main()
{
int n;
cin>>n;
a[0]=1;
a[1]=1;
s[0]=1;
for(int i=1;i<n;i++)
{
b[1]=(i+1)%10;
b[2]=(i+1)/10;
if(b[2]==0) b[0]=1;
else b[0]=2;
jia(s,a,s);
multiple(a,b,a);
}
//put(s);
jia(s,a,s);
put(s);
return 0;
}
//理解函数的封闭性与封装、接口价值和作用
//如无必要，函数调用最好不改变原来传入函数的数据