2015编程之美 资格赛 B题 区间dp

关键词：区间dp，统计回文子串数量
求长度为n的字符串的回文子串的个数（不一定连续！！！）
数据较大 mod100007输出

dp[i][j]：[i,j]区间内回文子串个数
dp[i][j]=dp[i][j-1]+sum{(dp[k+1][j-1]+1)|i<=k<=j-1&&a[k]==a[j]}
//部分细节需要单独处理
复杂度：O(n^3)

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#define mem(a,b) memset(a,sizeof(a),b)
#define mod 100007
using namespace std;

const int maxn = 1000+10;

int t;
char ch[maxn];
int dp[maxn][maxn];
int head[maxn];
struct Edge{
    int next;
}edge[maxn];

int main(){
    //freopen("a.txt","r",stdin);
    scanf("%d",&t);
    for(int cas=1;cas<=t;cas++){
        scanf("%s",ch);
        int n=strlen(ch);
        memset(head,-1,sizeof(head));
        for(int i=0;i<n;i++){
            int id=ch[i];
            edge[i].next=head[id],head[id]=i;
        }
        for(int i=0;i<n;i++) dp[i][i]=1;
        for(int k=1;k<=n;k++){
            for(int i=0;i+k<n;i++){
                int j=i+k;
                dp[i][j]=(dp[i][j-1]+1)%mod;
                int id=ch[j];
                for(int l=head[id];l!=-1&&l>=i;l=edge[l].next){
                    if(l<j){
                        if(l==j-1) dp[i][j]=(dp[i][j]+1)%mod;
                        else dp[i][j]=(dp[i][j]+dp[l+1][j-1]+1)%mod;
                    }
                }
            }
        }
        printf("Case #%d: %d\n",cas,dp[0][n-1]%mod);
    }
    return 0;
}

通过合理利用状态信息，状态转移方程可以进行简化！
dp[i][j]={dp[i][j-1]+dp[i+1]][j]-dp[i+1][j-1]}+(a[i]==a[j])*(dp[i+1][j-1]+1)
第一部分是不包含同时以a[i],a[j]结尾的回文子串个数
第二部分是同时以a[i]和a[j]结尾的回文子串个数
复杂度O(n^2)

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<limits.h>
#include<queue>
#define mem(a,b) memset(a,sizeof(a),b)
#define ll long long
#define INF 0x3f3f3f3f
using namespace std;

const int maxn = 1000+10 ;

int dp[maxn][maxn];
char a[maxn];
int t;

int main(){
    //freopen("a.txt","r",stdin);
    scanf("%d",&t);
    int cas=0;
    while(t--){
        scanf("%s",a);
        int n=strlen(a);
        memset(dp,0,sizeof(dp));
        for(int i=0;i<n;i++) dp[i][i]=1;
        for(int j=1;j<n;j++){
            for(int i=j-1;i>=0;i--){
                dp[i][j]=dp[i+1][j]+dp[i][j-1]-dp[i+1][j-1]+(a[i]==a[j])*(dp[i+1][j-1]+1);
            }
        }
        printf("Case #%d: %d\n",++cas,dp[0][n-1]);
    }
    return 0;
}