hdu4632(区间dp+容斥)

链接：点击打开链接

题意：给出一个串，求其中回文子序列的个数

代码：

#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
char s[1005];
long long dp[1005][1005];
const long long MOD=10007;
int main(){
    long long t,i,j,l,k,cas,len;
    scanf("%I64d",&t);
    for(cas=1;cas<=t;cas++){
        scanf("%s",s+1);
        len=strlen(s+1);
        memset(dp,0,sizeof(dp));
        for(i=1;i<=len;i++)
        dp[i][i]=1;
        for(i=1;i<=len-1;i++){
            if(s[i]==s[i+1])
            dp[i][i+1]=3;
            else
            dp[i][i+1]=2;
        }                                       //dp[i][j]表示i到j区间内的种数
        for(l=3;l<=len;l++){
            for(i=1;i<=len-l+1;i++){
                j=i+l-1;
                if(s[i]==s[j])                  //端点值相等则可以加上中间的种数
                dp[i][j]=(dp[i][j]+dp[i+1][j-1]+1)%MOD;
                dp[i][j]=(dp[i][j]+dp[i+1][j]+dp[i][j-1]-dp[i+1][j-1]+MOD)%MOD;
            }                                   //根据容斥原理减去中间重合的部分
        }
        printf("Case %I64d: %I64d\n",cas,dp[1][len]);
    }
    return 0;
}