HN2015集训 月宫的符卡序列

可能看到这个题有的人第一反应是回文自动机之类的，但是注意到这个题是跟回文串中心相关，而回文自动机则是与回文串末尾相关，可能无法处理。
那么与回文串中心相关的算法还有什么，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;
    }
}