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
分析:
建两棵回文树
一棵从前往后 记录前缀回文长度
一棵从后往前 记录后缀回文长度
然后遍历一遍计算答案
code:
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#define ll long long
using namespace std;
const int maxm=1e5+5;
struct Pam{
int nt[maxm][26];//下一个节点
int fail[maxm];//失配指针
int cnt[maxm];//节点i表示的本质不同的串的个数(最后cal完才是准确的)
int num[maxm];//节点i表示的回文串的回文后缀个数
int len[maxm];//节点i表示的回文串长度
int s[maxm];//存放添加的字符
int last;//新添加一个字母后所形成的最长回文串对应节点
int n;//添加的字符数,也是s数组的长度
int p;//节点指针,表示节点数量
int newnode(int x){//创建新节点
for(int i=0;i<26;i++){
nt[p][i]=0;
}
cnt[p]=num[p]=0;
len[p]=x;//长度为x
return p++;
}
void init(){//初始化
p=n=last=0;
newnode(0);//偶数长度树
newnode(-1);//奇数长度树
s[n]=-1;//设置第一个字符为-1(也可以是其他不可能被匹配的字符)
fail[0]=1;
}
int getfail(int x){//利用fail找到匹配的节点
while(s[n-len[x]-1]!=s[n])x=fail[x];
return x;
}
void add(int c){
c-='a';
s[++n]=c;
int cur=getfail(last);
if(!nt[cur][c]){//如果没有就新建一个节点
int now=newnode(len[cur]+2);//从cur拓展而来
fail[now]=nt[getfail(fail[cur])][c];
nt[cur][c]=now;
num[now]=num[fail[now]]+1;
}
last=nt[cur][c];
cnt[last]++;
}
void cal(){
for(int i=p-1;i>=0;i--){//累加cnt,如果fail[v]=u,则u一定是v的子回文串
cnt[fail[i]]+=cnt[i];
}
}
}p,pp;
char s[maxm];
int len1[maxm],len2[maxm];
int main(){
p.init();
pp.init();
scanf("%s",s+1);
int n=strlen(s+1);
for(int i=1;i<=n;i++){
p.add(s[i]);
pp.add(s[n-i+1]);
len1[i]=p.len[p.last];
len2[n-i+1]=pp.len[pp.last];
}
int ans=0;
for(int i=1;i<=n-1;i++){
ans=max(ans,len1[i]+len2[i+1]);
}
printf("%d\n",ans);
return 0;
}