转载
题意(事实上未曾找到相关题目):
给定一字符串,求它的回文子序列个数。内容相同位置不同的子序列算不同的子序列。
例如字符串aba中,回文子序列为”a”, “a”, “aa”, “b”, “aba”,共5个。
字符串长度 s.length <= 50
思路:
- 区间型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]。这个自己应该也能分析出来,和上面类似。
- 求区间i到j时, 会用到d[i+1][j-1], d[i][j-1], d[i+1][j]。 这也正是为什么区间型DP一般都是从相距较小的区间开始, 然后不扩大区间。当下所求区间(i, j)所用的子区间,之前都已求完,so还是那句话,直接用就好啦。
那个地方意见不一致欢迎指出,互相学习!!
用dp[i][j]表示第i到第j个字符间的最长回文子序列的长度(i<=j),则状态转移方程为:
dp[i][j]=dp[i+1][j] + dp[i][j-1] - dp[i+1][j-1] , if(s[i]!=s[j])
dp[i][j]=dp[i+1][j] + dp[i][j-1] +1 , if (s[i]==s[j])
#include <iostream>
#include<string>
#define LL long long
using namespace std;
LL dp[50][50];
LL solve(string str) {
int l = str.length();
for (int j = 0; j<l; ++j) {
dp[j][j] = 1;
for (int i = j - 1; i >= 0; --i) {
if (str[i] == str[j]) dp[i][j] = dp[i + 1][j] + dp[i][j - 1] + 1;
else dp[i][j] = dp[i + 1][j] + dp[i][j - 1] - dp[i + 1][j - 1];
}
}
return dp[0][l - 1];
}
int main() {
string str;
cin >> str;
cout << solve(str) << endl;
return 0;
}