hash被卡…
本来以为是回文自动机裸题
发现fail树上一条链的节点表示的回文子串的中点是不一样的…
不过回文树上的链是一样的
那么用建出回文树(我用回文自动机建的,manacher建不知道为什么WA了),然后找到以每个点为中点的最大回文子串,这个用manacher找
在对应节点加上贡献就行了
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>using namespace std;typedef unsigned long long ull;
typedef long long ll;const int N=1000010,P1=19980403,P2=998244353,base1=131,base2=129;int t,n,p,cnt=1,nxt[N][30],fa[N],len[N],fail[N],val[N];
char a[N];ull hs1[N],pw1[N],hs2[N],pw2[N];struct Hst{int G[P1+10],nxt[N<<1],pos[N<<1],cnt;ull val[N<<1];void clear(){ memset(G,0,sizeof(G)); cnt=0; }void insert(ull x,int p){ull u=x>>32,v=x-(u<<32);//for(int i=G[u];i;i=nxt[i])// if(val[i]==x) return ;++cnt; val[cnt]=v; nxt[cnt]=G[u]; G[u]=cnt; pos[cnt]=p;}int find(ull x){ull u=x>>32,v=x-(u<<32);for(int i=G[u];i;i=nxt[i])if(val[i]==v) return pos[i];}
}H;inline ull Getsub(int l,int r){ull a=(hs1[r]+P1-hs1[l-1]*pw1[r-l+1]%P1)%P1;ull b=(hs2[r]+P2-hs2[l-1]*pw2[r-l+1]%P2)%P2;return a<<32|b;
}inline void extend(int c,int ps){while(a[ps-len[p]-1]!=a[ps]) p=fail[p];if(!nxt[p][c]){int np=++cnt,k=fail[p]; len[np]=len[p]+2;while(a[ps-len[k]-1]!=a[ps]) k=fail[k];fail[np]=nxt[k][c]; fa[np]=p;nxt[p][c]=np;H.insert(Getsub(ps-len[np]+1,ps),np);}p=nxt[p][c];
}char b[N<<1];int r[N<<1];inline void manacher(){int t=0;for(int i=1;i<=n;i++)b[++t]=a[i],b[++t]='%';int mxr=0,p=0;for(int i=1;i<=t;i++){int j=0;if(i<=mxr)j=min(i+r[2*p-i],mxr+1);j=max(j,i+1);while(2*i-j>0 && j<=t && b[2*i-j]==b[j]) j++;r[i]=j-i;int L=i-r[i]+1,R=i+r[i]-1;if(b[L]=='%') L++,R--;if(L<=R) val[H.find(Getsub(L+1>>1,R+1>>1))]^=(i+1>>1)-1;if(i+r[i]-1>mxr) mxr=i+r[i]-1,p=i;}
}inline char nc(){static char buf[100000],*p1=buf,*p2=buf;return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}inline void read(int &x){char c=nc(); x=0;for(;c>'9'||c<'0';c=nc());for(;c>='0'&&c<='9';x=x*10+c-'0',c=nc());
}inline int read(char *a){char c=nc(); int x=0;for(;c>'z'||c<'a';c=nc());for(;c>='a'&&c<='z';a[++x]=c,c=nc()); return x;
}int main(){read(t);while(t--){n=read(a); cnt=1; p=0;memset(nxt,0,sizeof(nxt));memset(val,0,sizeof(val));H.clear();len[1]=-1; fail[0]=fail[1]=1; pw1[0]=pw2[0]=1;for(int i=1;i<=n;i++){hs1[i]=(1LL*hs1[i-1]*base1+a[i])%P1,pw1[i]=1LL*pw1[i-1]*base1%P1;hs2[i]=(1LL*hs2[i-1]*base2+a[i])%P2,pw2[i]=1LL*pw2[i-1]*base2%P2;}for(int i=1;i<=n;i++) extend(a[i]-'a',i);int ans=0;manacher();for(int i=cnt;i;i--) val[fa[i]]^=val[i];for(int i=2;i<=cnt;i++) ans=max(ans,val[i]);printf("%d\n",ans);}return 0;
}