BZOJ 2565 最长双回文串 （Manacher）

2565: 最长双回文串

Time Limit: 10 Sec Memory Limit: 128 MB
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

题目大意：
给出一个长度为n(2≤n≤105)的字符串，求最长双回文子串T（双回文指可将T分为非空的连续两部分X,Y，且X和Y都是回文串）的长度。

思路：
考虑枚举X，Y的交界点，那么我们只需要求出以每个位置为边界最多向左和向右延伸的最远回文子串长度。
以向左延伸为例。实际上就是要对每个位置i 找到最靠左的能包含到i的回文中心。于是我们从左到右扫描回文中心，同时记录一下最靠右的右边界r，扫描到一个位置时，将r+1 到i+p[i]-1 这一段的信息设为i，并且将右边界r 设为i+p[i]-1。
这样就找到了每个位置最远能向左延伸到的回文子串长度，向右同理。
这样O(n)的预处理之后找出了每个分隔符之下S和T的中心的最佳位置, 最后枚举一下分隔符的位置即可找到最优解。

#include <iostream> 
#include <algorithm>  
#include <cstdio>  
#include <cstring>  
#include <cmath>  
#define LL long long 
using namespace std;  

const int N = 200020;  

char ss[N], s[N];   
int tot=0, p[N];  
int l[N], r[N];  

//l[i]表示第i位的位置作为两个回文串的分割字符时候前面的回文串的中心在l[i]位置时最优 
//r[i]表示第i位的位置作为两个回文串的分割字符的时候后面的回文串的中心在r[i]位置时最优

void manacher(){
    memset(p, 0, sizeof(p)); 
    int id = 0, mx = 0;//正向
    for(int i=1; i<=tot; i++){
        if(mx > i) p[i] = min(p[id*2-i], mx-i+1);
        else p[i] = 1;
        while(s[i+p[i]] == s[i-p[i]]) p[i]++;
        if(i+p[i]-1 > mx){
            for(int j=mx+1; j<=i+p[i]-1; j++)  
                l[j] = i; 
                //因为i单调递增，所以第一次mx更新的时候一定是最优的（此时中心点i最靠左） 
                //右端点固定了，中心点越左，回文串越长  
              mx = i+p[i]-1;
              id = i;
        }
    }
    memset(p, 0, sizeof(p)); 
    id = mx = tot + 1;//反向 
    for(int i=tot; i>=1; i--){
        if(mx < i) p[i] = min(p[id*2-i], i-mx-1);
        else p[i] = 1;
        while(s[i-p[i]] == s[i+p[i]]) p[i]++;
        if(i-p[i]+1 < mx){
            for(int j=mx-1; j>=i-p[i]+1; j--)  
                r[j] = i;//同上 
              mx = i-p[i]+1;
              id = i;
        }
    }
}

int main(){
    scanf("%s", ss);  
    int len = strlen(ss);  
    s[0] = '+';  
    for(int i=0; i<len; i++){
        s[++tot] = '#'; 
        s[++tot] = ss[i];
    }
    s[++tot] = '#'; s[tot+1] = '-';  
    manacher();  
    int ans = 0;  
    for(int i=1; i<=tot; i+=2){//分隔符只能是插入的'#'字符
        int len = (2 * i - 2 * l[i] + 1) >> 1;
        len += (2 * r[i] - 2 * i + 1) >> 1;
        ans = max(ans, len);
    }  
    printf("%d\n", ans);  
    return 0;  
}