编程之美 资格赛B 字符串回文子序列



时间限制: 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

典型的字符串dp

显然的状态 dp[i][j] 表示闭合字串区间i-j的最大回文子序列数

那么怎么转移呢?

考虑到 在求dp[i][j]的时候 它的子区间 (不包含本身)全部都已经求出来了

那么 显然可以递推到更小的区间上求解(分治)

递推的时候 有两种决策, 如果sl == sr 那么dp[i][j] = dp[i+1][j] + dp[i][j-1] + 1

为什么+1呢: 比如字符串abba 在dp[2][4] + dp[1][3]的时候 多算了dp[2][3]区间 少算了 以sl开头 和以sr结尾的回文串个数.但是 以sl开头且sr结尾的回文串个数恰好等于dp[2][3] + 1 (显而易见, aa 原来是3 但是 1在前面dp[1][3]算了 还有一个1在dp[2][4]算了,所以只剩下1, 然后用两端构造回文,其个数等价于不用这两端构造)

当sl != sr的时候 还是 dp[i][j] = dp[i+1][j] + dp[i][j-1] - dp[i+1][j-1];

因为 多算了(sl sr不可构造) 所以减去

代码如下//
//  Created by Running Photon on 2015-04-17
//  Copyright (c) 2015 Running Photon. All rights reserved.
//
#include <algorithm>
#include <cctype>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iomanip>
#include <iostream>
#include <map>
#include <queue>
#include <string>
#include <set>
#include <vector>
#define ll long long
#define CLR(x) memset(x, 0, sizeof x)
using namespace std;
const int inf = 0x3f3f3f3f;
const int MOD = 100007;
const int maxn = 1e3 + 10;

int t;
char s[maxn];
int dp[maxn][maxn];
int main() {
#ifdef LOCAL
 freopen("in.txt", "r", stdin);
// freopen("out.txt","w",stdout);
#endif
 ios_base::sync_with_stdio(0);
 scanf("%d", &t); int cas = 0;
 while(t--){
  scanf("%s", s + 1);
  CLR(dp); cas++;
  int len = strlen(s + 1);
  for(int i = 1; i <= len; i++) dp[i][i] = 1;
  for(int i = len - 1; i > 0; i--) for(int j = i + 1; j <= len; j++){
   if(s[i] == s[j]){
    dp[i][j] = (dp[i+1][j] + dp[i][j-1] + 1) % MOD;
   }
   else{
    dp[i][j] = (dp[i+1][j] + dp[i][j-1] - dp[i+1][j-1]) % MOD;
   }
  }
  printf("Case #%d: %d\n", cas, dp[1][len]);
 }
 return 0;
}