BZOJ 1566 管道取珠（DP）

题目链接：http://61.187.179.132/JudgeOnline/problem.php?id=1566

题意：

思路：假设得到同一个序列S 的方案有ai=3，x,y,z，那么ai^2就是(x,x)(x,y)(x,z)(y,x)(y,y)(y,z)(z,x)(z,y)(z,z)，因此我 们可以将原来的一套装置看做两套。f[i][j][p][q]表示从第一套装置的上面取i个下面取j个、从第二套装置的上面取p个下面取q个，且两套装置 得到的两个序列相同的方案数，则答案就是f[n][m][n][m]。由于任何时候我们都是同时枚举两个装置中的，因此i+j=p+q在任意时刻都相等， 因此只要三维就够了。

 

int f[N][N][N],n,m;
char s1[N],s2[N];


void up(int &x,int y)
{
    x+=y;
    if(x>=mod) x-=mod;
}


int main()
{
    RD(n,m);
    RD(s1+1); RD(s2+1);
    reverse(s1+1,s1+n+1);
    reverse(s2+1,s2+m+1);
    f[0][0][0]=1;
    int i,j,k,t;
    for(i=0;i<=n;i++) for(j=0;j<=m;j++) for(k=0;k<=n;k++)
    {
        t=i+j-k;
        if(t<0||t>m) continue;
        if(s1[i+1]==s1[k+1]) up(f[i+1][j][k+1],f[i][j][k]);
        if(s1[i+1]==s2[t+1]) up(f[i+1][j][k],f[i][j][k]);
        if(s2[j+1]==s1[k+1]) up(f[i][j+1][k+1],f[i][j][k]);
        if(s2[j+1]==s2[t+1]) up(f[i][j+1][k],f[i][j][k]);
    }
    PR(f[n][m][n]);
}