JZOJ 2682 【WC2012选拔12.17】最长双回文串

Description:

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

题解：

manacher的裸题。
先用马拉车求出r数组。
显然可以先枚举X部分的回文中心i，对于每一个i，目标就是找到最大的j,满足(j>i+r[i]-1)and(j - r[j]<=i+r[i]-1)。
倒着枚举i，用个树状数组维护一下即可。

Code：

#include<cstdio>
#include<cstring>
#define fo(i, x, y) for(int i = x; i <= y; i ++) 
#define fd(i, x, y) for(int i = x; i >= y; i --)
#define min(a, b) ((a) < (b) ? (a) : (b))
#define max(a, b) ((a) > (b) ? (a) : (b))
#define low(x) ((x) & -(x))
using namespace std;

const int N = 100005 * 2;

int n, r[N], f[N];
char s[N];

void Build() {
    int p = 0, mx = 0;
    fo(i, 1, n) {
        r[i] = mx > i ? min(r[p * 2 - i], mx - i + 1) : 1;
        while (i > r[i] && s[i - r[i]] == s[i + r[i]]) r[i] ++;
        if(i + r[i] - 1 > mx) p = i, mx = i + r[i] - 1;
    }
}

void add(int x, int y) {
    while(x <= n) {
        f[x] = max(f[x], y);
        x += low(x);
    }
}

int find(int x) {
    int s = 0;
    while(x) {
        s = max(s, f[x]);
        x -= low(x);
    }
    return s;
}

int main() {
    freopen("a.in", "r", stdin);
    scanf("%s", s + 1);
    n = strlen(s + 1);
    fd(i, n, 1) s[i * 2] = s[i];
    fo(i, 1, n) s[i * 2 - 1] = '@';
    n = n * 2 + 1; s[n] = '@';
    Build();
    int ans = 0;
    fd(i, n, 1) {
        int s = i + r[i] > n ? find(i + r[i] - 1) : find(i + r[i]);
        if(s > 0) ans = max(ans, s - i);
        add(i - r[i] + 1, i);
    }
    printf("%d", ans);
}