管道取珠是小X很喜欢的一款游戏。在本题中,我们将考虑该游戏的一个简单改版。游戏画面如图1所示:
(图1)
游戏初始时,左侧上下两个管道分别有一定数量的小球(有深色球和浅色球两种类型),而右侧输出管道为空。每一次操作,可以从左侧选择一个管道,并将该管道中最右侧的球推入右边输出管道。
例如:我们首先从下管道中移一个球到输出管道中,将得到图2所示的情况。
(图2)
假设上管道中有n个球, 下管道中有m个球,则整个游戏过程需要进行n+m次操作,即将所有左侧管道中的球移入输出管道。最终n+m个球在输出管道中从右到左形成输出序列。
爱好数学的小X知道,他共有C(n+m,n)种不同的操作方式,而不同的操作方式可能导致相同的输出序列。举个例子,对于图3所示的游戏情形:
(图3)
我们用A表示浅色球,B表示深色球。并设移动上管道右侧球的操作为U,移动下管道右侧球的操作为D,则共有C(2+1,1)=3种不同的操作方式,分别为UUD,UDU,DUU;最终在输出管道中形成的输出序列(从右到左)分别为BAB,BBA,BBA。可以发现后两种操作方式将得到同样的输出序列。
假设最终可能产生的不同种类的输出序列共有K种,其中:第i种输出序列的产生方式(即不同的操作方式数目)有ai个。聪明的小X早已知道,
Σai=C(n+m,n)
因此,小X希望计算得到:
Σ(ai)^2
你能帮助他计算这个值么?由于这个值可能很大,因此只需要输出该值对1024523的取模即可(即除以1024523的余数)。
说明:文中C(n+m,n)表示组合数。组合数C(a,b)等价于在a个不同的物品中选取b个的选取方案数。
输入输出格式
输入格式:
输入文件中的第一行为两个整数n,m,分别表示上下两个管道中球的数目。
第二行中为一个AB字符串,长度为n,表示上管道中从左到右球的类型。其中:A表示浅色球,B表示深色球。
第三行中为一个AB字符串,长度为m,表示下管道中的情形。
输出格式:
输出文件中仅一行为一个整数,即为 除以1024523的余数。
输入输出样例
输入样例#1:
2 1
AB
B
输出样例#1:
5
说明
【样例说明】
样例即为文中(图3)。共有两种不同的输出序列形式,序列BAB有1种产生方式,而序列BBA有2种产生方式,因此答案为5。
【数据规模和约定】
对于30%的数据,满足:m,n<=12;
对于100%的数据,满足:m,n<=500。
Σ(ai)^2:选两次,两次一样的序列有多少种配对方式。(取两次, 形成序列相同的方案数).
dp[i][j][k][t]:表示第一次上管道选取了前i个,第一次下管道选取了前j个,第二次上管道选取了前k个,第二次下管道选取了前jt个.
//第一次选的i+j个和第二次选的k+t个形成相同序列
因为t = i + j - k.所以只需要三维。
初值dp[0][0][0] = 1;
滚动数组降低空间复杂的。//否则会超时。
//洛谷开o2,大牛分站可过
#include<iostream>
#include<cstdio>
using namespace std;
int n,m,dp[2][505][505];
char a[505],b[505];
const int mod=1024523;
int main()
{
freopen("ha.in","r",stdin);
scanf("%d%d",&n,&m);
scanf("%s%s", a + 1, b + 1);
dp[0][0][0]=1;
int pas = 1, now = 0;
for(register int i=0;i<=n;i++)
{
for(register int j=0;j<=m;j++)
for(register int k=0;k<=min(i+j,n);k++)
{
if(i != 0)dp[now][j][k] = 0;
int t=i+j-k;
if(t<0||t>m) continue;
if(i!=0&&k!=0&&a[i]==a[k])
{
dp[now][j][k]+=dp[pas][j][k-1]; if(dp[now][j][k] >= mod)dp[now][j][k] -= mod;
}
if(i!=0&&t!=0&&a[i]==b[t])
{
dp[now][j][k]+=dp[pas][j][k]; if(dp[now][j][k] >= mod) dp[now][j][k] -= mod;
}
if(j!=0&&k!=0&&b[j]==a[k])
{
dp[now][j][k]+=dp[now][j-1][k-1]; if(dp[now][j][k] >= mod) dp[now][j][k]-= mod;
}
if(j!=0&&t!=0&&b[j]==b[t])
{
dp[now][j][k]+=dp[now][j-1][k]; if(dp[now][j][k] >= mod) dp[now][j][k]-= mod;
}
}
now ^= 1; pas ^= 1;
}
printf("%d",dp[pas][m][n]);
return 0;
}