思路:设dp[i][j] 为i到j内回文子串的个数。先枚举所有字符串区间。再依据容斥原理。
那么状态转移方程为
dp[i][j] = dp[i][j-1] + dp[i+1][j] - dp[i+1][j-1]
如果 a[i] = a[j] , dp[i][j] += (dp[i+1][j-1] + 1)。
或者
dp[i][j] =dp[i][j-1] + dp[i+1][j]
如果 a[i] != a[j] , dp[i][j] -= dp[i+1][j-1]。
需要注意的是 假如区间两端元素相等,那么区间内任一回文串加上这两个元素组成的新字符串也是回文串。
1 #include<cstdio> 2 #include<iostream> 3 #include<cstring> 4 5 using namespace std; 6 const int maxn = 1005; 7 const int mod = 10007; 8 9 char a[maxn]; 10 int n,dp[maxn][maxn]; 11 12 int solve(){ 13 memset(dp,0,sizeof(dp)); 14 for(int i = n-1 ; i >= 0 ;i--){ 15 for(int j = i ;j < n ; j++){ 16 dp[i][j] = (dp[i+1][j] + dp[i][j-1] - dp[i+1][j-1] + mod) % mod; 17 if(a[i] == a[j]) dp[i][j] = (dp[i][j] + dp[i+1][j-1] + 1) % mod; 18 } 19 } 20 return dp[0][n-1]; 21 } 22 23 int main(){ 24 int cas; 25 scanf("%d",&cas); 26 for(int T = 1 ; T <= cas; T++){ 27 scanf("%s",a); 28 n = strlen(a); 29 printf("Case %d: %d ",T,solve()); 30 } 31 32 return 0; 33 }