可能看到这个题有的人第一反应是回文自动机之类的,但是注意到这个题是跟回文串中心相关,而回文自动机则是与回文串末尾相关,可能无法处理。
那么与回文串中心相关的算法还有什么,manacher(马拉车)。
学过的同学都知道吗,本质不同的回文串最多n个,也就是拉车到i的时候,本来已有f[i]所表示的回文串是之前已经出现过的,暴力新扩展出来的才是新出现的回文串。
那么就构造出了一套拓扑序,最后反拓扑序异或回去就行了。
还有一些处理技巧,双hash,用map会T成hmp,于是改用hash_table,快的抠脚。
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<map>
#include<ext/pb_ds/assoc_container.hpp>
#define N 2000005
#define ull unsigned long long
using namespace std;
const int seed1=131;
const int seed2=129;
const int mod=20000808;
template <typename T>
inline void _read(T& x){
char ch=getchar();bool sign=true;
while(!isdigit(ch)){if(ch=='-')sign=false;ch=getchar();}
for(x=0;isdigit(ch);ch=getchar())x=x*10+ch-'0';
if(!sign)x=-x;
}
int T,n,m,tot;
char s[N],S[N<<1];
ull pow1[N],pow2[N],_hash1[N],_hash2[N];
int fa[N],f[N],sum[N];
void compile1(){
pow1[0]=1;_hash1[0]=0;
for(int i=1;i<=n;i++){
pow1[i]=pow1[i-1]*seed1%mod;
_hash1[i]=(_hash1[i-1]*seed1+s[i])%mod;
}
}
void compile2(){
pow2[0]=1;_hash2[0]=0;
for(int i=1;i<=n;i++){
pow2[i]=pow2[i-1]*seed2%mod;
_hash2[i]=(_hash2[i-1]*seed2+s[i])%mod;
}
}
ull substr1(int l,int r){
return (_hash1[r]+mod-((_hash1[l-1]*pow1[r-l+1])%mod))%mod;
}
ull substr2(int l,int r){
return (_hash2[r]+mod-((_hash2[l-1]*pow2[r-l+1])%mod))%mod;
}
ull substr(int l,int r){
return (substr1(l,r)<<32)|(substr2(l,r));
}
__gnu_pbds::gp_hash_table<ull,int> MAP;
void manacher(){
int i,j,k,len;
len=1;
S[0]='$';S[1]='#';
for(i=1;i<=n;i++){
S[++len]=s[i];
S[++len]='#';
}
S[len+1]='\0';
int maxn=0,id=0,last=0;
for(i=1;i<=len;i++){
if(maxn>i)f[i]=min(maxn-i,f[(id<<1)-i]);
else f[i]=1;
if(f[i]==1)last=0;
else last=MAP[substr(((i-f[i])>>1)+1,(i+f[i]-1)>>1)];
while(S[i+f[i]]==S[i-f[i]]){
f[i]++;
int& temp=MAP[substr(((i-f[i])>>1)+1,(i+f[i]-1)>>1)];
if(!temp){
temp=++tot;fa[tot]=last;sum[tot]=0;
}
last=temp;
}
sum[last]^=(i>>1)-1;
if(i+f[i]>maxn){maxn=i+f[i];id=i;}
}
}
int main(){
cin>>T;
while(T--){
tot=0;
MAP.clear();
memset(fa,0,sizeof(fa));
memset(f,0,sizeof(f));
memset(sum,0,sizeof(sum));
int i,j,k;
scanf("%s",s+1);
n=strlen(s+1);
compile1();
compile2();
manacher();
int ans=0;
for(i=tot;i;i--){
ans=max(ans,sum[i]);
sum[fa[i]]^=sum[i];
}
cout<<ans<<endl;
}
}