给定字符串,求它的回文子序列个数。回文子序列反转字符顺序后仍然与原序列相同。例如字符串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
2.解题思路:本题利用容斥原理解决。本题是一个计数问题,要求找一个串中有多少个回文子串。本题的难点就在于子串的选取可以是不连续的。因此应该想到运用递推法来解决。定义dp(i,j)表示串[i,j]之间的回文子串的个数,那么可以得到如下递推式:
dp(i,j)=dp(i+1,j)+dp(i,j-1)+res (j-i>0)
其中res要分情况讨论:(1)如果s[i]==s[j]那么中间重叠部分dp(i+1,j-1)不需要减去,只需要再加上一个空串的情况即可,即res=1反之,根据容斥原理,需要减去中间重叠计算的部分,即res=-dp(i+1,j-1)。
当i==j时,即只有一个字符,那就是1,这样便不难通过递推式求出最终的答案了,最终的答案是dp(1,len+1)。len表示输入的串的长度,起点从1开始。
3.代码:
#define _CRT_SECURE_NO_WARNINGS #include<iostream> #include<algorithm> #include<string> #include<sstream> #include<set> #include<vector> #include<stack> #include<map> #include<queue> #include<deque> #include<cstdlib> #include<cstdio> #include<cstring> #include<cmath> #include<ctime> #include<functional> using namespace std; #define mod 100007 int dp[1005][1005]; char s[1005]; int F(int l, int r) { if (dp[l][r] != -1) return dp[l][r]; if (l > r) return dp[l][r] = 0; if (l == r) return dp[l][r] = 1; int& ret = dp[l][r]; ret = (F(l + 1, r) + F(l, r - 1)) % mod; if (s[l] == s[r]) ++ret;//首尾相等时,多加上一个中间为空的情况 else ret -= F(l + 1, r - 1);//否则,减去重叠部分 ret = (ret + mod) % mod; return ret; } int main() { //freopen("t.txt", "r", stdin); int cas = 0, t; scanf("%d", &t); while (t--) { memset(dp, -1, sizeof(dp)); scanf("%s", s + 1); int ans = F(1, strlen(s + 1)); printf("Case #%d: %d\n", ++cas, ans); } return 0; }
原文:http://blog.csdn.net/u014800748/article/details/45148441