2565: 最长双回文串 - BZOJ

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

HINT

样例说明

　　从第二个字符开始的字符串aacaabbacabb可分为aacaa与bbacabb两部分，且两者都是回文串。

数据规模及限制

　　对于10%的数据，2≤|S|≤10^3。

　　对于30%的数据，2≤|S|≤10^4。

　　对于100%的数据，2≤|S|≤10^5。

其实很简单

我们先用manacher算法求出以每个点为中心的最长回文串

然后求出l[i]和r[i]，分别表示向前和向后的最远的 回文串能覆盖到i 的中心（当然i是manacher算法里面加入的‘#’）因为只有这个才可以做分割点

 const
     maxn=;
 var
     c:array[..maxn]of char;
     p,ll,rr:array[..maxn]of longint;
     n:longint;

 function min(x,y:longint):longint;
 begin
     if x<y then exit(x);
     exit(y);
 end;

 function max(x,y:longint):longint;
 begin
     if x>y then exit(x);
     exit(y);
 end;

 procedure init;
 var
     i,id,r:longint;
 begin
     n:=;
     while not eoln do
       begin
         inc(n);
         read(c[n]);
         inc(n);
       end;
     c[]:='$';
     c[n+]:='#';
     id:=;
     r:=;
     for i:= to n do
       begin
         if r>=i then p[i]:=min(p[id<<-i],r-i+)
         else p[i]:=;
         while c[i+p[i]]=c[i-p[i]] do
           inc(p[i]);
         if i+p[i]->r then
         begin
           id:=i;
           r:=i+p[i]-;
         end;
       end;
 end;

 procedure work;
 var
     i,k,ans:longint;
 begin
     k:=;
     for i:= to n do
       if i and = then
       begin
         while k+p[k]-<i do
           inc(k);
         ll[i]:=i-k;
       end;
     k:=n;
     for i:=n downto  do
       if i and = then
       begin
         while k-p[k]+>i do
           dec(k);
         rr[i]:=k-i;
       end;
     ans:=;
     for i:= to n do
       if i and = then ans:=max(ans,ll[i]+rr[i]);
     writeln(ans);
 end;

 begin
     init;
     work;
 end.