【bzoj 1566】管道取珠（DP）

最新推荐文章于 2019-01-04 20:43:59 发布

zP1nG

最新推荐文章于 2019-01-04 20:43:59 发布

阅读量421

点赞数

分类专栏： zP1nG的bzoj

本文链接：https://blog.csdn.net/zP1nG/article/details/78084169

版权

zP1nG的bzoj 专栏收录该内容

183 篇文章 4 订阅

订阅专栏

传送门biu~
ai是一种情况的方案数，那么把ai^2看成两个人分别取珠，情况相同的方案数。
f[i][j][k]表示取了i个珠子，第一个人在第一个管道取j个，第二个人在第一个管道取k个，且结果相同的方案数。由f[0][0][0]=1开始dp最终解f[n+m][n][n]。

#include <bits/stdc++.h>
using namespace std;
const int modd=1024523;
int n,m;
int f[2][505][505];
int a[505],b[505];
char ac[505],bc[505];
int main(){
    scanf("%d%d",&n,&m);
    scanf("%s%s",ac,bc);
    for(int i=1;i<=n;++i)   a[i]=ac[n-i]-'A';
    for(int i=1;i<=m;++i)   b[i]=bc[m-i]-'A';
    f[0][0][0]=1;
    for(int i=0;i<n+m;++i){
        int now=i&1;
        for(int j=0;j<=i && j<=n;++j)
            for(int k=0;k<=i && k<=n;++k)
                if(f[now][j][k]){
                    if(a[j+1]==a[k+1])      (f[now^1][j+1][k+1]+=f[now][j][k])%=modd;
                    if(a[j+1]==b[i-k+1])    (f[now^1][j+1][k]+=f[now][j][k])%=modd;
                    if(b[i-j+1]==a[k+1])    (f[now^1][j][k+1]+=f[now][j][k])%=modd;
                    if(b[i-j+1]==b[i-k+1])  (f[now^1][j][k]+=f[now][j][k])%=modd;
                    f[now][j][k]=0;
                }
    }
    printf("%d",f[(n+m)&1][n][n]);
    return 0;
}