bzoj2565 最长双回文串

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

2015.4.25新加数据一组

 

正解：回文自动机。

正反构造两个回文自动机，顺便记录以当前点为终点的最长回文后缀。然后直接枚举端点，取和的$max$即可。

 

//It is made by wfj_2048~
#include <algorithm>
#include <iostream>
#include <complex>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <vector>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <set>
#define inf (1<<30)
#define N (200010)
#define il inline
#define RG register
#define ll long long
#define File(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)

using namespace std;

int ans;

struct PAM{

    int ch[N][26],fa[N],l[N],len[N],n,la,sz;
    char s[N];

    il void pre(){ l[++sz]=-1,fa[0]=1; return; }

    il void add(RG int c,RG int n){
    RG int x=la; while (s[n-l[x]-1]!=s[n]) x=fa[x];
    if (!ch[x][c]){
        RG int v=++sz,k=fa[x]; l[v]=l[x]+2;
        while (s[n-l[k]-1]!=s[n]) k=fa[k];
        fa[v]=ch[k][c],ch[x][c]=v;
    }
    la=ch[x][c],len[n]=l[ch[x][c]]; return;
    }
    
}tr1,tr2;

il void work(){
    scanf("%s",tr1.s+1),tr1.n=strlen(tr1.s+1),tr2.n=tr1.n,tr1.pre();
    for (RG int i=1;i<=tr2.n;++i) tr2.s[i]=tr1.s[tr1.n-i+1]; tr2.pre();
    for (RG int i=1;i<=tr1.n;++i)
    tr1.add(tr1.s[i]-97,i);
    for (RG int i=1;i<=tr2.n;++i) tr2.add(tr2.s[i]-97,i);
    for (RG int i=1;i<tr1.n;++i) ans=max(ans,tr1.len[i]+tr2.len[tr2.n-i]);
    printf("%d\n",ans); return;
}

int main(){
    File("PAM");
    work();
    return 0;
}