BZOJ 1566: [NOI2009]管道取珠

来源：互联网发布：百度开放云绑定域名编辑：程序博客网时间：2024/05/06 05:21

/****************失踪人口的回归****************/

没错已经很久没写了

Description

Input

第一行包含两个整数n, m，分别表示上下两个管道中球的数目。第二行为一个AB字符串，长度为n，表示上管道中从左到右球的类型。其中A表示浅色球，B表示深色球。第三行为一个AB字符串，长度为m，表示下管道中的情形。

Output

仅包含一行，即为 Sigma(Ai^2) i从1到k 除以1024523的余数。

Sample Input

2 1
AB
B

Sample Output

HINT

样例即为文中(图3)。共有两种不同的输出序列形式，序列BAB有1种产生方式，而序列BBA有2种产生方式，因此答案为5。
【大致数据规模】
约30%的数据满足 n, m ≤ 12；
约100%的数据满足n, m ≤ 500。

呣唔我一开始以为那个a^2很难搞

然而一看大神题解让我茅厕茅塞顿开

所谓a^2次相当于取俩次珠，俩次结果相同的方案数

就是输出序列相同的方案数

于是就成了一道大水题。求相同序列有多少个

#include<cstdio>#include<algorithm>#define mod 1024523#define g getchar()#define ll long long#define inf 0x3f3f3f3fusing namespace std;inline ll read(){ll x=0,f=1;char ch=g;for(;ch<'0'||ch>'9';ch=g)if(ch=='-')f=-1;for(;ch>='0'&&ch<='9';ch=g)x=x*10+ch-'0';return x*f;}inline void out(ll x){int a[25],wei=0;if(x<0)putchar('-'),x=-x;for(;x;x/=10)a[++wei]=x%10;if(wei==0){puts("0");return;}for(int j=wei;j>=1;--j)putchar('0'+a[j]);putchar('\n');}int n,m;char s1[510],s2[510];int f[2][510][510];    //空间怕不够用滚动数组int main(){//freopen("","r",stdin);//freopen("","w",stdout);scanf("%d%d",&n,&m);scanf("%s%s",s1+1,s2+1);f[0][0][0]=1;for(int k=0;k<=n+m-1;++k)for(int i=max(0,k-m);i<=min(k,n);++i)for(int j=max(0,k-m);j<=min(k,n);++j){int x=f[k%2][i][j];f[k%2][i][j]=0;if(i<n&&j<n&&s1[n-i]==s1[n-j])f[(k+1)%2][i+1][j+1]=(x+f[(k+1)%2][i+1][j+1])%mod;if(i<n&&k-j<m&&s1[n-i]==s2[m-k+j])f[(k+1)%2][i+1][j]=(x+f[(k+1)%2][i+1][j])%mod;if(k-i<m&&k-j<m&&s2[m-k+i]==s2[m-k+j])f[(k+1)%2][i][j]=(x+f[(k+1)%2][i][j])%mod;if(k-i<m&&j<n&&s2[m-k+i]==s1[n-j])f[(k+1)%2][i][j+1]=(x+f[(k+1)%2][i][j+1])%mod;}out(f[(n+m)%2][n][n]);return 0;}

0 0