BZOJ 1566: [NOI2009]管道取珠

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

没错已经很久没写了

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 0x3f3f3f3f
using 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;
}