BZOJ 2342 SHOI 2011 双倍回文 Manacher

在Manacher处理出所有极长回文子串后。
枚举双倍回文子串的中心i，考虑快速查询出中心左侧的其他回文中心j，并满足： j+pj≥i,j−pj≥i−pi
即 j+pj≥i,j≥i−12pi
枚举的中心肯定都在#上。
发现由于i是单调增的（我们这么枚举），而 j+pj 总会在某个i时不满足，因此此时j可以被剔除。

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define FOR(i,j,k) for(i=j;i<=k;++i)
const int N = 1000002;
int fa[N], p[N]; char s[N], a[N];
int find(int x) { return x == fa[x] ? x : fa[x] = find(fa[x]); }
int main() {
    int ans, i, j, mx = 0, id, n, m;
    scanf("%d%s", &n, s + 1);
    m = 2 * n + 1;
    for (i = 1; i <= n; ++i) {
        a[i * 2] = s[i];
        a[i * 2 + 1] = '#';
    }
    a[0] = '&'; a[m + 1] = '^'; a[1] = '#';
    for (i = 1; i <= m; ++i) {
        if (mx > i) p[i] = min(mx - i, p[2 * id - i]);
        else p[i] = 1;
        while (a[i - p[i]] == a[i + p[i]]) ++p[i];
        if (i + p[i] > mx) mx = i + p[i], id = i;
    }
    FOR(i,1,m) fa[i] = a[i] == '#' ? i : i + 1;
    FOR(i,3,m) if (a[i] == '#') {
        for(j = find(max(i - p[i] / 2, 1));
            j < i && j + p[j] <= i;
            j = fa[j] = find(j + 1));
        if (j < i) ans = max(ans, i - j << 1);
    }
    printf("%d", ans);
    return 0;
}