bzoj 2342: [Shoi2011]双倍回文

Description

Input

输入分为两行，第一行为一个整数，表示字符串的长度，第二行有个连续的小写的英文字符，表示字符串的内容。

 

Output

输出文件只有一行，即：输入数据中字符串的最长双倍回文子串的长度，如果双倍回文子串不存在，则输出0。

Sample Input

16
ggabaabaabaaball

Sample Output

12

HINT

N<=500000

 

题解：

先manacher求出f数组，然后显然我们枚举中点x如果把x-y作为第三段，那么 y-f[y]<=x && y<=x+f[x]/2 即可满足

于是我们以y-f[y]排序一个数组，维护一个set，然后每次放入y-f[y]<=i

显然 y越靠近x+f[x]/2答案越大，我们就找x+f[x]/2的前驱即可

 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<set>
using namespace std;
const int N=1000005;
char S[N],s[N];int f[N],d[N];
struct node{
    int id,x;
    bool operator<(const node &p)const{return x<p.x;}
}a[N];
set<int>t;
int main()
{
    //freopen("pp.in","r",stdin);
    int n;
    scanf("%d",&n);
    scanf("%s",S+1);
    int l=0;
    for(int i=1;i<=n;i++){
        s[++l]='#';
        s[++l]=S[i];
    }
    s[++l]='#';
    int id=1;
    for(int i=1;i<=l;i++){
        if(i<id+f[id])f[i]=min(id+f[id]-i,f[(id<<1)-i]);
        else f[i]=1;
        while(s[i+f[i]]==s[i-f[i]])f[i]++;
        if(f[i]>f[id])id=i;
    }
    n=0;
    for(int i=3;i<=l;i+=2)d[++n]=(f[i]-1)>>1;
    for(int i=1;i<=n;i++)a[i].id=i,a[i].x=i-d[i];
    sort(a+1,a+n+1);
    int p=0,ans=0,x;
    for(int i=1;i<=n;i++){
        while(p<n && a[p+1].x<=i)p++,t.insert(a[p].id);
        x=*--t.upper_bound(i+(d[i]>>1));
        if(x==*t.begin())continue;
        if(x-i>ans)ans=x-i; 
    }
    printf("%d",ans<<2);
    return 0;
}