【bzoj 1566】管道取珠(DP)

来源:互联网 发布:软件行业目标市场分析 编辑:程序博客网 时间:2024/04/30 05:15

传送门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;}

【看数据大小貌似不需要滚动数组QAQ

原创粉丝点击