这个也可以利用manacher算法来做
奇数长度的回文串是 从i开始的后缀 与 反转后字符串中从原来的i位子开始的后缀 的公共前缀
偶数长度类似
所以和最长公共子串类似的思路, 将原串S与反转后的S中间加入一个不会出现的字符拼接在一起, 对得到的字符串S’求高度数组
以原来S[i]为中心的奇数长度的回文串长度 = 后缀S’[i…]和后缀 S’[2n-i…]的公共前缀
以原来S[i-1]和S[i]为中心的回文串长度 = 后缀S’[i…]和后缀S’[2n-i+1…]的公共前缀
由于我这个程序高度数组可能某一个位置上可能为0,也就是排名最后的后缀,但是这样数据是aabaa会报错,进行更改后,现在如果数据是abcd会报成0。
#include<bits/stdc++.h>
using namespace std;
const int maxn=2e5+100;
char s[maxn];
int rk[maxn],tmp[maxn],sa[maxn],n,k,lcp[maxn];
bool cmp_sa(int i,int j)
{
if(rk[i]!=rk[j]){
return rk[i]<rk[j];
}else{
int ri=i+k<=n?rk[i+k]:-1;
int rj=j+k<=n?rk[j+k]:-1;
return ri<rj;
}
}
void construct_sa()
{
for(int i=0;i<=n;i++){
sa[i]=i;
rk[i]=i<n?s[i]:-1;
}
for(k=1;k<=n;k<<=1){
sort(sa,sa+n+1,cmp_sa);
tmp[sa[0]]=0;
for(int i=1;i<=n;i++){
tmp[sa[i]]=tmp[sa[i-1]]+(cmp_sa(sa[i-1],sa[i])?1:0);
}
copy(tmp,tmp+n+1,rk);
}
}
void construct_lcp()
{
for(int i=0;i<=n;i++){
rk[sa[i]]=i;
}
int h=0;
lcp[0]=0;
for(int i=0;i<n;i++){
int j=sa[rk[i]-1];
if(h>0){
h--;
}
for(;j+h<n&&i+h<n;h++){
if(s[j+h]!=s[i+h]){
break;
}
}
lcp[rk[i]-1]=h;
}
lcp[n]=lcp[n-1];
}
//ST表进行最值查询
int mi[maxn][30];
void init_rmq(int a[],int n)
{
for(int i=0;i<n;i++){
mi[i][0]=a[i];
}
for(int j=1;(1<<j)<=n;j++){
for(int i=0;i+(1<<j)<=n;i++){
mi[i][j]=min(mi[i][j-1],mi[i+(1<<(j-1))][j-1]);
}
}
}
int rmq(int l,int r)
{
int k=0;
while((1<<(k+1))<=r-l+1){
k++;
}
return min(mi[l][k],mi[r-(1<<k)+1][k]);
}
int main()
{
scanf("%s",s);
int len=strlen(s);
s[len]='$';
reverse_copy(s,s+len,s+len+1);
n=len*2+1;
s[n]=0;
for(int i=0;i<n;i++){
cout<<s[i];
}
cout<<endl;
construct_sa();
construct_lcp();
for(int i=0;i<=n;i++){
cout<<sa[i]<<" ";
}
cout<<endl;
for(int i=0;i<=n;i++){
cout<<rk[i]<<" ";
}
cout<<endl;
init_rmq(lcp,n+1);
for(int i=0;i<=n;i++){
cout<<i<<" ";
}
cout<<endl;
for(int i=0;i<=n;i++){
cout<<lcp[i]<<" ";
}
cout<<endl;
int ans=0;
if(len%2){
for(int i=0;i<len;i++){
int j=len*2-i;
int l=rmq(min(rk[i],rk[j]),max(rk[i],rk[j]));
cout<<i<<" "<<rk[i]<<" "<<j<<" "<<rk[j]<<" "<<l<<endl;
ans=max(ans,2*l-1);
}
}else{
for(int i=0;i<len;i++){
int j=len*2-i+1;
int l=rmq(min(rk[i],rk[j]),max(rk[i],rk[j]));
cout<<i<<" "<<rk[i]<<" "<<j<<" "<<rk[j]<<" "<<l<<endl;
ans=max(ans,2*l);
}
}
cout<<ans<<endl;
return 0;
}