JZOJ2682. 【WC2012选拔12.17】最长双回文串

Description

顺序和逆序读起来完全一样的串叫做回文串。比如acbca是回文串，而abc不是（abc的顺序为“abc”，逆序为“cba”，不相同）。输入长度为n的串S，求S的最长双回文子串T,即可将T分为两部分X，Y，（|X|,|Y|≥1）且X和Y都是回文串。

Input

一行由小写英文字母组成的字符串S。

Output

一行一个整数，表示最长双回文子串的长度。

Sample Input

baacaabbacabb

Sample Output

12

Data Constraint

对于10%的数据，2≤|S|≤103。
对于30%的数据，2≤|S|≤104。
对于100%的数据，2≤|S|≤105。

题解

题目要求的是从一个位置切开，
然后两边都是回文串。

所以就枚举哪一个位置被切开。
先用manacher求出以每个位置为回文中心的最长回文半径。
然后再对每一个点分别求出它左边和右边能覆盖到它的最远回文中心。
然后枚举切开的位置，
为了避免一个字母同时被划分到两个回文串里面，切开的位置只能是添加进去的字符，
利用它左边和右边的最远回文中心的位置，求出答案。

code

#include<queue>
#include<cstdio>
#include<iostream>
#include<algorithm>
#include <cstring>
#include <string.h>
#include <cmath>
#include <math.h>
#define ll long long
#define ull unsigned long long
#define N 100008
#define M 102
#define db double
#define P putchar
#define G getchar
#define mo 1000000007
#define zm 19260817
using namespace std;
char ch;
void read(int &n)
{
    n=0;
    ch=G();
    while((ch<'0' || ch>'9') && ch!='-')ch=G();
    ll w=1;
    if(ch=='-')w=-1,ch=G();
    while('0'<=ch && ch<='9')n=(n<<3)+(n<<1)+ch-'0',ch=G();
    n*=w;
}

db max(db a,db b){return a>b?a:b;}
int min(int a,int b){return a<b?a:b;}
ll abs(ll x){return x<0?-x:x;}
ll sqr(ll x){return x*x;}
void write(ll x){if(x>9) write(x/10);P(x%10+'0');}

char s[N*2],t[N*2];
int len,ans,p[N*2],l[N*2],r[N*2],mx,id;

void manacher()
{
    mx=0;
    for(int i=1;i<=len;i++)
    {
        if(mx>=i)p[i]=min(mx-i,p[id*2-i]);else p[i]=0;
        while(s[i-p[i]-1]==s[i+p[i]+1])p[i]++;
        if(mx<i+p[i])mx=p[i]+i,id=i;
    }
}

int main()
{
    scanf("%s",t+1);
    len=strlen(t+1);
    s[0]='#';
    for(int i=1;i<=len;i++)
        s[(i<<1)-1]='*',s[i<<1]=t[i];
    s[(len<<1)+1]='*';
    len=(len<<1)+1;
    s[len+1]='!';
    manacher();

    mx=0;
    for(int i=1;i<=len && mx<len;i++)
        if(i+p[i]>mx)
        {
            for(int j=mx+1;j<=i+p[i];j++)
                l[j]=i;
            mx=i+p[i];
        }

    mx=len+1;
    for(int i=len;i && mx;i--)
        if(i-p[i]<mx)
        {
            for(int j=i-p[i];j<mx;j++)
                r[j]=i;
            mx=i-p[i];
        }

    for(int i=1;i<=len;i++)
        if(s[i]=='*')ans=max(ans,(((i-l[i])<<1)+((r[i]-i)<<1)+1)>>1);
    write(ans);
}