洛谷P1591 阶乘数码（高精度）//高精度加法、高精度乘法

原题链接：
https://www.luogu.com.cn/problem/P1591

//洛谷P1591 阶乘数码

#include
#include
const int N =10005 ; // 同一个数据出现多次，可抽象成一个变量
using namespace std;
void jia(int p[],int q[],int m[])
{
int a[N],b[N],c[N],x=0; //数组可能不够大
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
memset(c,0,sizeof©);
for(int i=0;i<N;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[N],b[N],c[N],x=0;
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
memset(c,0,sizeof©);
for(int i=0;i<N;i++)
{
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];
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];
}
}

int a[N],b[N],c[N];
int main()
{
int t;
cin>>t;
while(t–)
{
int n,s;
long long sum=0; //数据类型决定了数据的范围
cin>>n>>s;
memset(a,0,sizeof(a));//必须初始化
memset(b,0,sizeof(b));
memset(c,0,sizeof©);
a[0]=1;
a[1]=1;
b[0]=1;
b[1]=1;
c[0]=1;
c[1]=1;
for(int i=1;i<n;i++)
{
jia(b,c,b); //b++
multiple(a,b,a);
}
for(int i=1;i<=a[0];i++)
{
if(a[i]==s) sum++;
}
cout<<sum<<endl;

}
return 0;

}
//根据阶乘和的题解修改而来