啊 太蠢了这道题没看出来
∑ai2 其实就是两个人分别取 取出相同的方案数啊!
为什么我这题的代码这么慢QAQ
于是乎我发现一个可以快很多的地方 见代码
#include<bits/stdc++.h>
#define me(a,x) memset(a,x,sizeof a)
using namespace std;
const int N=502,mod=1024523;
char s1[N],s2[N]; int f[2][N][N];
int main()
{
int n,m,i,j,k,t=0;
scanf("%d%d%s%s",&n,&m,s1+1,s2+1);
for(i=1;i<=n/2;i++)swap(s1[n-i+1],s1[i]);
for(i=1;i<=m/2;i++)swap(s2[m-i+1],s2[i]);
f[0][0][0]=1;
for(i=0;i<n+m;i++){
int u1=i<n?i:n;
for(j=0;j<=u1;j++)for(k=0;k<=u1;k++)if(f[t][j][k]>0){//这个if神剪枝。。
int u=f[t][j][k];
if(k<n && j<n && s1[j+1]==s1[k+1])(f[1-t][j+1][k+1]+=u)%=mod;
if(k<n && i-j<m && s2[i-j+1]==s1[k+1])(f[1-t][j][k+1]+=u)%=mod;
if(j<n && i-k<m && s1[j+1]==s2[i-k+1])(f[1-t][j+1][k]+=u)%=mod;
if(i-j<m && i-k<m && s2[i-j+1]==s2[i-k+1])(f[1-t][j][k]+=u)%=mod;
}
me(f[t],0); t=1-t;
}
printf("%d\n",f[t][n][n]);
return 0;
}