bzoj2565: 最长双回文串 回文自动机

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新加数据一组

分析

首先我们有一个结论:(在WC2017被证明)
最长的双倍回文串中一定有一个回文串是不可拓展(最长的)的.
所以我们可以枚举取到最长的那个回文串,然后计算在剩下的字符中最长的回文串
用回文自动机可以在线性时间内求出以i字符结尾的最长回文串（其实manacher也可以）
那么对于这个以字符i为结尾的最长双回文串，答案就是$mx[i]+mx[i - mx[i]]$
左右都可能拓展，所以倒着再做一遍即可。

代码

/**************************************************************
    Problem: 2565
    User: 2014lvzelong
    Language: C++
    Result: Accepted
    Time:96 ms
    Memory:14608 kb
****************************************************************/

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#include<cmath>
using namespace std;
const int N = 110000;
int a[N], f[N], ch[N][26], s[N], len[N], mx[N], n, l, last, sz;
void readchar() {
    char ch = getchar();
    while(ch < 'a' || ch > 'z') ch = getchar();
    for(;ch >= 'a' && ch <= 'z'; ch = getchar()) a[++n] = ch - 'a';
}

void pre() {
    f[last = sz = 0] = 1;
    len[++sz] = s[l = 0] = -1;
    memset(ch[0], 0, sizeof(ch[0]));
    memset(ch[1], 0, sizeof(ch[1]));
}
void Extend(int c) {
    s[++l] = c; int p, x, np;
    for(p = last;s[l - len[p] - 1] != s[l]; p = f[p]) ;
    if(!ch[p][c]) {
        len[np = ++sz] = len[p] + 2; memset(ch[np], 0, sizeof(ch[np]));
        for(x = f[p]; s[l - len[x] - 1] != s[l]; x = f[x]) ;
        f[np] = ch[x][c]; ch[p][c] = np;
    }
    mx[l] = len[last = ch[p][c]]; 
}

int main() {
    readchar(); int ans = 0;
    pre(); for(int i = 1;i <= n; ++i) Extend(a[i]);
    for(int i = 1;i <= n; ++i) ans = max(ans, mx[i] + mx[i - mx[i]]); 
    pre(); for(int i = n; i; --i) Extend(a[i]);
    for(int i = 1;i <= n; ++i) ans = max(ans, mx[i] + mx[i - mx[i]]);
    printf("%d\n", ans);
    return 0;
}