【BZOJ】2342: [Shoi2011]双倍回文（Manacher）

题目

传送门：QWQ

分析

（sb如我写了发不知道什么东西在洛谷上竟然水了84分

嗯咳

设$ i $为双重回文的中心

如果$ j~i $ 可以被算作答案，只有满足如下两式：

• $ p[j]+j \geq i $
• $ 2*(i-j) \leq p[j] $

计算时我们先做一次马拉车，然后按照 $ p[j]+j \geq i $排序，保证它的单调，接着把满足$ 2*(i-j) \leq p[j] $扔进set里询问。

代码

#include <bits/stdc++.h>
using namespace std;
const int maxn=;
char s2[maxn],s[ maxn];
int p[ maxn], len;
set<int> t;
void Manacher(){
    len=strlen(s2+);
    for(int i=;i<=len;i++){
        s[i*-]='#'; s[i*]=s2[i];
    }
     s[len=len*+]='#';
     int right=, pos=;
     for(int i=;i<=len;i++){
         if(i<right){ p[i]=min(p[*pos-i],right-i); } else p[i]=;
        while(i+p[i]<=len &&i-p[i]> && s[i+p[i]]==s[i-p[i]]) p[i]++;
        if(i+p[i]>right){
            right=i+p[i]; pos=i;
        }
     }
}
int q[maxn], f[maxn];
bool cmp(int a,int b){ return (a-f[a])<(b-f[b]); }
int main(){
    int n;
    scanf("%d%s",&n,s2+);
    Manacher();
    for(int i=;i<=n;i++) q[i]=i, f[i]=(p[i*+]-)/;
    sort(q+,q++n,cmp);
    int now=,ans=;
    for(int i=;i<=n;i++){
        while(now<=n&&q[now]-f[q[now]]<=i) {
            t.insert(q[now]);
            now++;
        }
        set<int>::iterator tmp=t.upper_bound(i+f[i]/);
        if(tmp!=t.begin ()){
            ans=max(ans,(*--tmp - i));
        }
    }
    printf("%d\n",ans*);
    return ;
}
/*
17
qwertyuaabbaabbaa
*/