Bzoj2342 [Shoi2011]双倍回文

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 1915  Solved: 713
[Submit][Status][Discuss]

Description

 

Input

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

 

Output

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

Sample Input

16
ggabaabaabaaball

Sample Output

12

HINT

N<=500000

http://hzwer.com/5375.html

算法太神以至于看不懂，先存下代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<set>
using namespace std;
const int mxn=800010;
int p[mxn],q[mxn];
int len,ans;
char s[mxn];
set<int>t;
int cmp(int a,int b){
    return (a-p[a])<(b-p[b]);
}
void manacher(){
    int id=0,mx=0;
    for(int i=1;i<=len;i++){
        if(mx>=i)p[i]=min(p[id]+id-i,p[2*id-i]);
        else p[i]=0;
        while(s[i+1+p[i]]==s[i-p[i]])p[i]++;
        //字符串长度是偶数，对称轴在i和i+1之间，
        //那么左边是i-p[i]，右边对称位置是i+1+p[i] 
        if(p[i]+i>mx)mx=p[i]+i,id=i;
    }
    return;
}
int main(){
    scanf("%d %s",&len,s+1);
    s[0]='$';
    int i,j;
    manacher();
    for(i=1;i<=len;i++)q[i]=i;
    sort(q+1,q+len+1,cmp);
    int now=1;
    for(i=1;i<=len;i++){
        while(now<=len && q[now]-p[q[now]]<=i){
            t.insert(q[now]);
            now++;
        }
        set<int>::iterator tmp=t.upper_bound(i+p[i]/2);
        if(tmp!=t.begin())ans=max(ans,(*--tmp-i)*4);
    }
    printf("%d
",ans);
    return 0;
}