BZOJ 1566: [NOI2009]管道取珠

14 篇文章 0 订阅

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

没错已经很久没写了

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

5

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;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值