BZOJ2565:最长双回文串(Manacher)

Description

顺序和逆序读起来完全一样的串叫做回文串。比如acbca是回文串，而abc不是（abc的顺序为“abc”，逆序为“cba”，不相同）。
输入长度为n的串S，求S的最长双回文子串T,即可将T分为两部分X，Y，（|X|,|Y|≥1）且X和Y都是回文串。

Input

一行由小写英文字母组成的字符串S。

Output

一行一个整数，表示最长双回文子串的长度。

Sample Input

baacaabbacabb

Sample Output

12

HINT

样例说明

从第二个字符开始的字符串aacaabbacabb可分为aacaa与bbacabb两部分，且两者都是回文串。

对于100%的数据，2≤|S|≤10^5

Solution

在manacher的时候顺便统计一下以每个位置为结束/开始的最长回文串最后枚举统计答案就好了

Code

#include<iostream>
#include<cstring>
#include<cstdio>
#define N (200000+1000)
using namespace std;

int n,tot,len[N],ans,E[N],S[N];
char s[N],a[N];

void Manacher()
{
    int x,mid=0,maxn=0;
    for (int i=1; i<=tot; ++i)
    {
        if (i>maxn) x=1;
        else x=min(maxn-i+1,len[mid*2-i]);
        S[i-x+1]=max(S[i-x+1],x*2-1);
        E[i+x-1]=max(E[i+x-1],x*2-1);//这里一开始忘加了 
        while (s[i+x]==s[i-x])
        {
            ++x;
            S[i-x+1]=max(S[i-x+1],x*2-1);
            E[i+x-1]=max(E[i+x-1],x*2-1);
        }
        len[i]=x;
        if (i+x-1>maxn) maxn=i+x-1,mid=i; 
    }
}

int main()
{
    scanf("%s",a);
    n=strlen(a);
    s[++tot]='@'; s[++tot]='#';
    for (int i=0; i<n; ++i)
        s[++tot]=a[i], s[++tot]='#';
    s[++tot]='$';
    Manacher();
    
    for (int i=4; i<tot-2; i+=2)
        ans=max(ans,(S[i]+E[i]-1)/2);
    printf("%d",ans);
}