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 500003char 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;}