BZOJ——2342双倍回文

传送门：点我

Description

Input

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

Output

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

Sample Input

16
ggabaabaabaaball

Sample Output

12

HINT

N<=500000

思路：

以每个#为双倍回文中心枚举，找出左端或右端的回文中心（由于对称找一端即可），以左端为例，判断小回文中心+小回文半径是否>=双倍回文中心，并且总长度一定是4的倍数。

贴代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 500003
char s[N],t[N<<1];
int p[N<<1];
int main()
{
    int n,i,j,mid;
    while(~scanf("%d",&n))
    {
        scanf("%s",s);
        for(i=1; i<=n; i++)
            t[2*i]=s[i-1],t[2*i+1]='#';
        t[0]='!',t[1]='#';
        int len=2*n+1;
        int ans=0,id=0,mx=0;
        for(i=1; i<=len; i++)
        {
            if(mx>i) p[i]=min(p[2*id-i],mx-i);
            else p[i]=1;
            while(t[i-p[i]]==t[i+p[i]]) p[i]++;
            if(i+p[i]>mx) mx=i+p[i],id=i;
        }
        for(i=3; i<=len; i+=2)//以#(i)为分割中心
        {
            int R=p[i]-1;
            while(R>3)
            {
                if(R%4==0)//双倍回文串长度是4的倍数
                {
                    if(R<ans) break;//小于答案就不用算了
                    j=i-R;//扩充后双倍回文的最左端下标
                    mid=(i+j)>>1;//左边回文的中心下标
                    if(mid+p[mid]-1>=i)
                        ans=max(ans,R);
                }
                R-=2;//如果没有找到一直缩短枚举
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}