这一题的状态方程不是特别难想,还可以用滚动数组优化,就是高精度略嫌麻烦,我把答案转换成字符串来输出了
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
using namespace std;
const int xmax=10000+10;
const int zmax=100+10;
const int maxn=30;
const int base=10000;
char x[xmax],z[zmax];
int f[xmax][maxn],last[xmax][maxn];//last记录上一次的d[x][z]
int n;
void add_bignum(int a[],int b[],int c[])
{
for(int i=maxn-1;i>=1;i--)
{
int sum=b[i]+c[i];
a[i]+=sum%base;//本位和
a[i-1]+=sum/base;//进位
}
a[0]+=b[0]+c[0];
}
void print_ans(int d[])
{
int i,len=0,rlen=0;
char s[120],str[5];//把有效位转换成字符串;
for(i=0;i<maxn;i++)
if(d[i])//找到第一个不为零的数组下标i,并把它转换成字符串插入s
{
sprintf(str,"%d",d[i]);
rlen=strlen(str);
for(int j=0;j<rlen;j++)
s[len++]=str[j];
s[len]='\0';
break;
}
for(int j=i+1;j<maxn;j++)
{
sprintf(str,"%d",d[j]);
rlen=strlen(str);
for(int k=0;k<4-rlen;k++)//处理前导零
s[len++]='0';
for(int m=0;m<rlen;m++)
s[len++]=str[m];
}
s[len]='\0';
printf("%s\n",s);
}
int main()
{
//freopen("in.txt","r",stdin);
//freopen("myout.txt","w",stdout);
scanf("%d",&n);
while(n--)
{
scanf("%s %s",x,z);
int xlen=strlen(x),zlen=strlen(z);
int count=0;
memset(f,0,sizeof(f));
memset(last,0,sizeof(last));
for(int j=xlen-1;j>=0;j--)
{
if(z[zlen-1]==x[j])
{
count++;
}
f[j][maxn-1]=count;
}
// for(int j=0;j<xlen;j++)
// printf("%d ",f[j]);
// printf("\n");
for(int i=zlen-2;i>=0;i--)
for(int j=xlen-1;j>=0;j--)
{
if(z[i]==x[j])//j=xlen-1,i<zlen-1?
{
//last[j]=f[j];滚动数组,如果没有高精度的话,这样就可以了,实际上f[j]=d[z+1][x]
memcpy(last[j],f[j],maxn*sizeof(int));
//f[j]=last[j+1]+f[j+1];//f[j]=d[z+1][x+1]+d[z][x+1]
//f[j]=[0];
memset(f[j],0,maxn*sizeof(int));
// for(int k=0;k<maxn;k++)
// printf("%d ",f[j][k]);
// printf("\n");
add_bignum(f[j],last[j+1],f[j+1]);
}
else
{
//last[j]=f[j];//这个last真不知道为什么会有、、、这几天看一下
memcpy(last[j],f[j],maxn*sizeof(int));
//f[j]=f[j+1];//d[z][x+1]
memcpy(f[j],f[j+1],maxn*sizeof(int));
}
}
// for(int i=0;i<maxn;i++)
// printf("%d ",f[0][i]);
// printf("\n");
print_ans(f[0]);
}
return 0;
}