CF 245H Queries for Number of Palindromes

题意：给出长度不超过50000的字符串，询问区间$[l,r]$内有多少回文子串。
思路：因为字符串长度比较长，首先我们要解决的就是判断$s_l...s_r$是不是回文字符串。
这个可以递推得到。
设p[l][r]为起始位置为l，结束位置为r的子串是否是回文串。1表示是，0表示不是
递推方程为：

$$p[l][r] = 1 \iff p[l+1][r-1] = 1 \land s_l =s_r$$
边界条件是
$$p[l][l] = p[l][l+1] = 1$$
统计回文子串同样可以递推得到
设dp[l][r]为其实位置为l,结束位置为r的子串中含有的回文子串的个数。
递推方程为：
$$dp[l][r] = dp[l+1][r] + dp[l][r-1] - dp[l+1][r-1] + p[l][r]$$
边界条件是
$$dp[l][l] =dp[l][l+1] = 1$$
整体的算法的复杂度就是 $\varTheta(n^2)$
代码如下：

#include <cstdio>
#include <algorithm>
#include <cstring>

using namespace std;

char s[5010];
int p[5010][5010],dp[5010][5010];

int pail(int l, int r)
{
    if(p[l][r] >= 0)
        return p[l][r];
    if(l >= r)
        return p[l][r] = 1;
    if(s[l] == s[r] && pail(l+1,r-1))
        return p[l][r] = 1;
    else
        return p[l][r] = 0;
}

int sum(int l, int r)
{
    if(dp[l][r] >= 0)
        return dp[l][r];
    if(l == r)
        return dp[l][r] = 1;
    if(l > r)
        return dp[l][r] = 0;
    else
        return dp[l][r] = sum(l,r-1) + sum(l+1,r) - sum(l+1,r-1) + pail(l,r);
}



int main(void)
{
    //freopen("input.txt","r",stdin);
    int T;
    scanf("%s",s);
    scanf("%d",&T);
    memset(p,-1,sizeof(p));
    memset(dp,-1,sizeof(dp));
    for(int i = 0; i < T; ++i){
        int l,r;
        scanf("%d%d",&l,&r);
        printf("%d\n",sum(--l,--r));
    }
    return 0;
}