2015编程之美资格赛：回文子序列数 DP

2015编程之美资格赛第二题：回文子序列数

题目链接：http://hihocoder.com/contest/msbop2015qual/problem/2

分析：动态规划。用一个数组dp（j，i）表示在给定字符串str中第j个元素到第i个元素之间回文子序列的个数，在j==i时dp（j，i）=1；j和i相邻时判断str[ j ]和str[ i ]是否相等，分别赋值给dp（i，j）3和2。然后是一般情况：对于dp（j，i），我们先让它的值加上一个dp（j，i-1）+1 （即，先加上在j到i-1区间的回文子序列个数，然后再加上第i个元素这个单独的元素的回文子序列个数1），然后用一个循环k，从i-1循环到j，当str[ k ]=str[ i ]的时候，更新dp（j,i）的值：dp（j，i）+=dp（k+1，e-1）+1，即如果i-1到j中某个字符和i相等，那么在区间（j+1，i-1）中的每对回文序列都可以再分别在两头加上str[ k ]和str[ i ]，再构成一个新的回文序列，所以我们要拿区间（j+1，i-1）内的回文子序列的个数再加上由str[ k ]和str[ i ]单独构成回文序列数1，然后把这个值加到dp（j，i）上。最后别忘了在每求出一步dp值时对mod=100007取模。

实现代码如下：

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
int dp[1010][1010];
char str[1010];
int len;
void solve(int s,int e)
{
    if(s==e)
    {
        dp[s][e]=1;
        return ;
    }
    if(e-s==1)
    {
        if(str[e]==str[s]) dp[s][e]=3;
        else dp[s][e]=2;
        return ;
    }
    /*dp[s][e]+=dp[s][e-1]+1;
    for(int i=e-1;i>=s;i--)
      if(str[i]==str[e])
        dp[s][e]+=dp[i+1][e-1]+1;*/
    dp[s][e]=dp[s][e-1]+dp[s+1][e]-dp[s+1][e-1];
    if(str[s]==str[e]) dp[s][e]+=dp[s+1][e-1]+1;
    dp[s][e]+=100007;
    dp[s][e]%=100007;
}
int main()
{
    int t,T=1;
    cin>>t;
    while(t--)
    {
        memset(dp,0,sizeof(dp));
        scanf("%s",str);
        len=strlen(str);
        for(int i=0;i<len;i++)
          for(int j=i;j>=0;j--)
            solve(j,i);
        printf("Case #%d: %d\n",T++,dp[0][len-1]);
    }
    return 0;
}

代码2：

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
int dp[1010][1010];
char str[1010];
int len;
void solve(int s,int e)
{
    if(s==e)
    {
        dp[s][e]=1;
        return ;
    }
    if(str[e]==str[s])
      dp[s][e]=dp[s][e-1]+dp[s+1][e]+1;
    else dp[s][e]=dp[s][e-1]+dp[s+1][e]-dp[s+1][e-1];
    dp[s][e]+=100007;
    dp[s][e]%=100007;
}
int main()
{
    int t,T=1;
    cin>>t;
    while(t--)
    {
        memset(dp,0,sizeof(dp));
        scanf("%s",str);
        len=strlen(str);
        for(int i=0;i<len;i++)
          for(int j=i;j>=0;j--)
            solve(j,i);
        printf("Case #%d: %d\n",T++,dp[0][len-1]);
    }
    return 0;
}