关键词:区间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;
}