Queries for Number of Palindromes（O(n^2)求所有子区间回文串个数）

题目链接：https://vjudge.net/problem/Gym-226149C
题意：给定字符串s(长度≤5e3)，询问n(≤1e6)个连续子串中回文串的个数

思路1：回文树

因为n只有5e3，所以可以提前预处理出每一段l,r的答案。时间复杂度O(n²+Q)

//回文树复杂度：设字符集大小为m,母串长度为n,则空间复杂度为O(nm),时间复杂度为O(nlogm)[也可以说是O(n),因为logm极小]
//原谅我都没有想到可以先预处理…heihei

#include <stdio.h>
#include <math.h>
#include <string.h>
#include <algorithm>
#include <iostream>
#include <vector>
#include <queue>
using namespace std;
//#define int long long
#define LL long long
const int manx=5e3+10;
const int max_len=2e6+10;
char s[manx];
int ans[manx][manx];
struct pal_tree
{
    int net[manx][26],fail[manx],num[manx],len[manx],S[manx];//cnt[manx];
    int p,last,n;
    int newnode(int l)
    {
        for(int i=0;i<26;i++)
            net[p][i]=0;
        num[p]=0;
        len[p]=l;
        return p++;
    }
    void init()
    {
        p=0;
        newnode(0);
        newnode(-1);
        fail[0]=1;
        S[0]=-1;
        last=n=0;
    }
    int getfail(int x)
    {
        while(S[n-len[x]-1]!=S[n])x=fail[x];
        return x;
    }
    int add(int c)
    {
        c-='a';
        S[++n]=c;
        int cur=getfail(last);
        if(!net[cur][c])
        {
            int now=newnode(len[cur]+2);
            fail[now]=net[getfail(fail[cur])][c];
            net[cur][c]=now;
            num[now]=num[fail[now]]+1;
        }
        last=net[cur][c];
        return num[last];
    }
}tree;
int main()
{
    int n,ss,e;
    scanf("%s",s);
    for(int i=0;s[i]!='\0';i++)//枚举起点
    {
        tree.init();
        for(int j=i;s[j]!='\0';j++)//枚举终点
            ans[i][j]=(i==j?0:ans[i][j-1])+tree.add(s[j]);
    }
    scanf("%d",&n);
    while(n--)
    {
        scanf("%d%d",&ss,&e);
        printf("%d\n",ans[ss-1][e-1]);
    }
    return 0;
}

思路2：区间dp 容斥定理

我们令dp[i][j]表示区间[i,j]中回文串的个数，那么长度为len的子串就可以由它的长度为len-1的子串来更新。
转移方程：dp[i][j]=dp[i+1][j]+dp[i][j-1]-dp[i+1][j-1]+【i到j为一个回文串?1:0 】

和另外一个区间dp的题很相似：HDU - 4632 - Palindrome subsequence，但这个题中所求的回文串不要求是连续的子序列。题解

#include <stdio.h>
#include <math.h>
#include <string.h>
#include <algorithm>
#include <iostream>
#include <vector>
#include <queue>
using namespace std;
//#define int long long
#define LL long long
const int manx=5e3+10;
const int max_len=2e6+10;
char s[manx];
int dp[manx][manx],is[manx][manx];
int main()
{
    int n,ss,e;
    scanf("%s",s);
    int len=strlen(s);
    for(int l=1;l<=len;l++)
    {
        for(int i=0;i+l-1<len;i++)
        {
            int j=i+l-1;
            if(l==1)dp[i][j]=is[i][j]=1;
            else if(l==2)is[i][j]=(s[i]==s[j]?1:0);
            else if(s[i]==s[j]&&is[i+1][j-1])is[i][j]=1;
            else is[i][j]=0;
        }
    }
    for(int l=2;l<=len;l++)
        for(int i=0;i+l-1<len;i++)
        {
            int j=i+l-1;
            dp[i][j]=dp[i+1][j]+dp[i][j-1]-dp[i+1][j-1]+is[i][j];
        }
    scanf("%d",&n);
    while(n--)
    {
        scanf("%d%d",&ss,&e);
        printf("%d\n",dp[ss-1][e-1]);
    }
    return 0;
}

还有一种dp的思路：https://blog.csdn.net/hhaile/article/details/8223547