BZOJ 1566: [NOI2009]管道取珠

啊 太蠢了这道题没看出来
$\sum ai^2$ 其实就是两个人分别取 取出相同的方案数啊！
为什么我这题的代码这么慢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;
}