hdu4632:区间dp+容斥原理+记忆化搜索

hdu4632:区间dp+容斥原理+记忆化搜索

题意：

求给出字符串str中子序列（注：不是子字符串）是回文串的个数。

思路：

1. 状态设计：dp[s][e]表示str[s]~str[e]区间的回文子序列个数
2. 初始化/递归结束条件：
   1）s>e:表示空串，dp[s][e]为零；
   2）s==e：单个字符构成一个回文串，返回1；
3. 状态转移：dp[s][e]的上一个状态是dp[s][e-1]或dp[s+1][e]，将它们加起来，区间dp[s+1][e-1]中的回文子序列是被重复计数的，因此需要减去这个区间的dp。若str[s]、str[e]相等，则他们构成一个子序列（+1），同时，他们与区间dp[s+1][e-1]中每一个回文子序列都构成一个回文子序列，（+dp[s+1][e-1])。
   1）str[s]!=str[e]：dp[s][e]=(dp[s][e-1]+dp[s+1][e]-dp[s+1][e-1]+mod)%mod
   2)str[s]==str[e]：dp[s][e]=(dp[s][e-1]+dp[s+1][e]+1)%mod

代码：

using namespace std;
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
const int mod = 10007;
const int maxn = 1e3+5;
int dp[maxn][maxn];
char str[maxn];

int dpDfs(int s, int e)//记忆化搜索
{
    if(s > e)return 0;//空字符串无回文
    if(s == e)return dp[s][e] = 1;//单个字符的字符串
    if(dp[s][e]) return dp[s][e];//已经搜索过，直接返回
    //容斥原理
    if(str[s] != str[e]){
        return dp[s][e] = (dpDfs(s, e-1) + dpDfs(s+1, e) - dpDfs(s+1, e-1) + mod) % mod;
    }
    else{
        return dp[s][e] = (dpDfs(s,e-1) + dpDfs(s+1,e) + 1) % mod;
    }
}

int main()
{
    //freopen("in.txt", "r", stdin);
    int t;
    cin >> t;
    for(int Case = 1; Case <= t; Case++){
        cin >> str;
        int n = strlen(str);
        //init
        memset(dp, 0, sizeof dp);
        //dpDfs
        printf("Case %d: %d\n", Case, dpDfs(0,n-1));
    }
    return 0;
}