[回文自动机 Manacher] BZOJ4166: 月宫的符卡序列

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;
}