hdu 5785 Interesting(manacher+前缀和)

题目链接：hdu 5785 Interesting

题意：

有一个长度为n的串（n<=10^6），对 1 <= i <= j < k <= length(s) . 如果[i，j]和[j+1,k]都是回文串。则对答案的贡献为 i*k ,求贡献和。

题解：

详细题解传送门

#include<bits/stdc++.h>
#define F(i,a,b) for(int i=a;i<=b;++i)
using namespace std;
typedef long long ll;

const int N=1e6+7,P=1e9+7;

struct Manacher{
    char str[N*2];
    int p[N*2],len,mx,id,tl,ans,i;
    void maxlen(char *s){
        len=strlen(s),mx=0,id=0,tl=0,str[tl++]='$',str[tl++]='#';
        for(i=0;i<len;i++)str[tl++]=s[i],str[tl++]='#';
        for(i=2,str[tl]=0,ans=0;i<tl;i++){
            p[i]=mx>i?min(p[(id<<1)-i],mx-i):1;
            while(str[i-p[i]]==str[i+p[i]])p[i]++;
            if(i+p[i]>mx)mx=i+p[i],id=i;
        }
    }
}M;
char s[N];
int dp[N*2][4],ans;

void add(int op,int l,int r,int v)
{
    if(l>r)return;
    dp[l][op]=(dp[l][op]+v)%P;
    dp[r+1][op]=(dp[r+1][op]-v)%P;
}

int main(){
    while(~scanf("%s",s))
    {
        memset(dp,0,sizeof(dp));
        M.maxlen(s);
        F(i,1,M.tl-2)
        {
            add(0,i,M.p[i]+i-1,i);
            add(1,i,M.p[i]+i-1,1);
            add(2,i-M.p[i]+1,i,i);
            add(3,i-M.p[i]+1,i,1);
        }
        F(j,1,M.tl)F(i,0,3)dp[j][i]=(dp[j][i]+dp[j-1][i]+P)%P;
        ans=0;
        for(int i=4;i<M.tl;i+=2)
        {
            ll a=(dp[i][2]-1ll*dp[i][3]*(i/2)%P+P)%P;
            ll b=(dp[i-2][0]-1ll*dp[i-2][1]*((i-2)/2)%P+P)%P;
            ans=(a*b%P+ans)%P;
        }
        printf("%d
",ans);
    }
    return 0;
}