求字符串回文子序列的个数

转载
题意（事实上未曾找到相关题目）：

给定一字符串，求它的回文子序列个数。内容相同位置不同的子序列算不同的子序列。
例如字符串aba中，回文子序列为”a”, “a”, “aa”, “b”, “aba”，共5个。
字符串长度 s.length <= 50
思路：

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还是那句话，直接用就好啦。
   那个地方意见不一致欢迎指出，互相学习！！

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