BZOJ 2565 最长双回文串（回文自动机）

题意

给一个长度为N的字符串S。对于一个字符串AB，如果A和B都是回文串，那么称AB是一个双回文串。求问S最长双回文子串的长度？
N <= 100000

题解

正反双向构造回文自动机，得到某一个点为结尾和开始的最长回文串长度，记为f[i],g[i];

答案就是f[i]+g[i+1]的最大值。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const int N=100500;
int len[N],cnt[N],to[N][30],s[N],fail[N],last,tot,n;
int ans,L,p1[N],p2[N];
char s1[N],s2[N];
void init(){
    tot=1;n=0;
    len[0]=0;len[1]=-1;
    fail[0]=1;s[0]=-1;last=0;
    memset(to,0,sizeof(to));
}
void add(int c,int p[]){
    s[++n]=c;
    int cur,now,tmp;
    for(cur=last;s[n-len[cur]-1]!=c;cur=fail[cur]);
    if(!to[cur][c]){
        tot++;
        len[tot]=len[cur]+2;now=tot;
        for(tmp=fail[cur];s[n-len[tmp]-1]!=c;tmp=fail[tmp]);
        fail[now]=to[tmp][c];
        to[cur][c]=now;
    }
    last=to[cur][c];
    cnt[last]++;
    p[n]=len[last];
}
int main(){
    scanf("%s",s1+1);
    int L=strlen(s1+1);
    for(int i=1;i<=L;i++)s2[i]=s1[L-i+1];
    init();
    for(int i=1;i<=L;i++)add(s1[i]-'a',p1);
    init();
    for(int i=1;i<=L;i++)add(s2[i]-'a',p2);
    for(int i=1;i<=L-1;i++){
        ans=max(ans,p1[i]+p2[L-i]);
    } 
    printf("%d",ans);
    return 0; 
}