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两部分,且两者都是回文串。
对于100%的数据,2≤|S|≤10^5
思路
首先用manacher求出每个点所能扩展的最长回文串,然后对于每个点记录l[i]r[i]代表以当前点结尾或者开头的最长回文串长度…
对于lr的求法有两种…一种是我这种比较low的求法…一种是LTY大佬边做manacher边求lr的做法…
代码
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int N=200000+5;
int ans,len,p[N],l[N],r[N];
char s[N],str[N];
void prework()
{
int i;
for (i=0;s[i];i++)
str[2*i+1]='#',str[(i+1)*2]=s[i];
len=2*i+1;str[0]='$';str[len]=str[len+1]='#';
}
void manacher()
{
int id,mx=0;
for (int i=1;i<len;i++)
{
p[i]=i<mx?min(p[id*2-i],mx-i):1;
while(str[i-p[i]]==str[i+p[i]]) p[i]++;
if (p[i]+i>mx) mx=p[i]+i,id=i;
}
}
int main()
{
scanf("%s",s);ans=0;
prework();manacher();
for (int i=1;i<=len;i++)
l[i-p[i]+1]=max(l[i-p[i]+1],p[i]-1),r[i+p[i]-1]=max(r[i+p[i]-1],p[i]-1);
for (int i=1;i<=len;i++)
{
if (i&1) l[i]=max(l[i],l[i-1]-2);
else l[i]=max(l[i],l[i-1]);
}
for (int i=1;i<len;i++)
if (l[i]&&r[i]) ans=max(ans,l[i]+r[i]);
printf("%d\n",ans);
return 0;
}