【编程之美挑战赛】回文字符序列（区间dp）

给定字符串，求它的回文子序列个数。回文子序列反转字符顺序后仍然与原序列相同。例如字符串aba中，回文子序列为"a", "a", "aa", "b", "aba"，共5个。内容相同位置不同的子序列算不同的子序列。

输入

第一行一个整数T，表示数据组数。之后是T组数据，每组数据为一行字符串。

输出

对于每组数据输出一行，格式为"Case #X: Y"，X代表数据编号（从1开始），Y为答案。答案对100007取模。

数据范围

1 ≤ T ≤ 30

小数据

字符串长度 ≤ 25

大数据

字符串长度 ≤ 1000

样例输入

5
aba
abcbaddabcba
12111112351121
ccccccc
fdadfa

样例输出

Case #1: 5
Case #2: 277
Case #3: 1333
Case #4: 127
Case #5: 17  

hint:5:

1.f 
2.d 
3.a 
4.f 
5.d 
6.a 
7.ff 
8.dd 
9.aa 
10.fdf 
11.faf 
12.fdf 
13.dad 
14.ada 
15.afa 
16.fddf 
17.fdadf

http://hihocoder.com/contest/msbop2015qual/problem/2

大神是用dp做的。。

用dp[i][j]表示这一段里有多少个回文串，那首先dp[i][j]=dp[i+1][j]+dp[i][j-1]，但是dp[i+1][j]和dp[i][j-1]可能有公共部分，所以要减去dp[i+1][j-1]。
如果str[i]==str[j]的话，还要加上dp[i+1][j-1]+1。

以后保险起见还是所有减法的操作都+mod吧。。

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <string>
#include<cmath>
#define mod 100007  
using namespace std;
int dp[1005][1005];
int main(){
	int t,cnt=0;
	cin>>t;
	while(t--){
		string x;
		cin>>x; 
		x='0'+x;
		int n=x.length()-1;
		memset(dp,0,sizeof(dp));
		ll s=0;
		for(int i=1;i<=n;++i){
			for(int j=i;j>=1;--j){
				dp[j][i]+=dp[j+1][i];
				dp[j][i]+=dp[j][i-1];
				dp[j][i]-=dp[j+1][i-1]; //减去公共部分（加了两次） 
				if(x[i]==x[j])
					dp[j][i]+=dp[j+1][i-1]+1; //(+1)是新生成回文串【x[i]x[j]】 
				
				dp[j][i]+=100007; //并不理解为什么上面的减法操作会导致负数
				dp[j][i]%=100007; 
			}
		}
		printf("Case #%d: %d\n",++cnt,dp[1][n]);
	}
	return 0;
}