标签:
时间限制:2000ms
单点时限:1000ms
内存限制:256MB
描述
给定字符串,求它的回文子序列个数。回文子序列反转字符顺序后仍然与原序列相同。例如字符串aba中,回文子序列为”a”, “a”, “aa”, “b”, “aba”,共5个。内容相同位置不同的子序列算不同的子序列。
输入
第一行一个整数T,表示数据组数。之后是T组数据,每组数据为一行字符串。
输出
对于每组数据输出一行,格式为”Case #X: Y”,X代表数据编号(从1开始),Y为答案。答案对100007取模。
数据范围
1 ≤ T ≤ 30
小数据
字符串长度 ≤ 25
大数据
字符串长度 ≤ 1000
样例输入
5
aba
abcbaddabcba
12111112351121
ccccccc
fdadfa
样例输出
Case #1: 5
Case #2: 277
Case #3: 1333
Case #4: 127
Case #5: 17
分析:
1. 区间型DP, 先来说一下状态方程: 如求i到j这个区间共有多少回文子序列 ,两种情况①当s[i] == s[j]时,例如“a….a”,d[i][j] = d[i][j-1] + d[i+1][j] + 1(由d[i][j] = 2*d[i+1][j-1] + d[i][j-1] + d[i+1][j] + 1 - 2*d[i+1][j-1] 化
简来的。 解释:已知 i+1 到 j-1 有d[i+1][j-1]个回文子序列, 又有s[i] == s[j],
那么可与中间(串i+1到j-1)已知的回文子序列再构成d[i+1][j-1]个回文子序列,再加上原来中间串所包含的回文子序列共2*d[i+1]
[j-1]个, 两个a组合”aa”也是回文, 所以再加1, 再分别计算左、右两边的a和中间串所构成的回文子序列,
但是这个时候还没完,注意:在分别计算左右两边a和中间串时,又算了两遍中间串包含的回文子序列,
所以在减去2*d[i+1][j-1]个);②当s[i] != s[j]时,d[i][j] = d[i][j-1]+d[i+1][j] - d[i+1][j-1]。这个自己应该也能分析出来,和上面类似。
2. 求区间i到j时, 会用到d[i+1][j-1], d[i][j-1], d[i+1][j]。 这也正是为什么区间型DP一般都是从相距较小的区间开始, 然后不扩大区间。当下所求区间(i, j)所用的子区间,之前都已求完,so还是那句话,直接用就好啦。
那个地方意见不一致欢迎指出,互相学习!!
#include<iostream>
#include<cstdio>
#include<string.h>
#include<cstring>
using namespace std;
const int MOD = 100007;
int t, len, d[1010][1010];
char s[1010];
void dp()
{
for(int i = 1; i < len; i++)
{
for(int j = 0; j < (len - i); j++)
{
if(s[j] == s[j+i])
d[j][i+j] = (d[j][i+j-1] + d[j+1][i+j] + 1) % MOD;
else
d[j][i+j] = (d[j][i+j-1] + d[j+1][i+j] - d[j+1][i+j-1] + MOD) % MOD;
}
}
}
int main()
{
cin >> t;
for(int q = 1; q <= t; q++)
{
memset(s, 0, sizeof(s));
memset(d, 0, sizeof(d));
cin >> s;
len = strlen(s);
for(int i = 0; i < len; i++) d[i][i] = 1;
dp();
printf("Case #%d: %d\n", q, d[0][len-1] % MOD);
}
return 0;
}
标签:
原文地址:http://www.cnblogs.com/xinshui/p/4455815.html