[分块 回文自动机] LOJ#6070. 「2017 山东一轮集训 Day4」基因

回文自动机那套理论可以在 翁文涛的集训队论文 里看

每个串的回文自动机形态是固定的，一个子串的回文自动机也就是一个子图

那么就可以用记忆化的方式来降低复杂度……我感觉好像就是论文里提到的“不基于势能分析的插入方法“

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <vector>
#include <map>
#include <cstring>

using namespace std;

typedef unsigned long long ll;

const int N=100010,M=350,P=31;

int n,ct,q,S,it;
int a[N];
int vis[N];

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 rea(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());
}

int p1,p2,cnt,len[N],fail[N];
int pos[M][N],pp[M][N],nxt[N][30],trans[N][30];

inline void PushBack(int l,int r){
    int c=a[r]; 
    if(a[r-len[p1]-1]!=c || r-len[p1]-1<l) p1=trans[p1][c];
    if(!nxt[p1][c]){
        int cur=++cnt,k=fail[p1]; len[cur]=len[p1]+2;
        if(a[r-len[k]-1]==c) fail[cur]=nxt[k][c];
        else fail[cur]=nxt[trans[k][c]][c];
        memcpy(trans[cur],trans[fail[cur]],sizeof(trans[cur]));
        trans[cur][a[r-len[fail[cur]]]]=fail[cur];
        nxt[p1][c]=cur;
    }
    p1=nxt[p1][c];
    if(len[p1]==r-l+1) p2=p1;
}

inline void PushFront(int l,int r){
    int c=a[l];
    if(a[l+len[p2]+1]!=c || l+len[p2]+1>r) p2=trans[p2][c];
    if(!nxt[p2][c]){
        int cur=++cnt,k=fail[p2]; len[cur]=len[p2]+2;
        if(a[l+len[k]+1]==c) fail[cur]=nxt[k][c];
        else fail[cur]=nxt[trans[k][c]][c];
        memcpy(trans[cur],trans[fail[cur]],sizeof(trans[p2]));
        trans[cur][a[l+len[fail[cur]]]]=fail[cur];
        nxt[p2][c]=cur;
    }
    p2=nxt[p2][c];
    if(len[p2]==r-l+1) p1=p2;
}

int ans[M][N];

int main(){
    rea(ct);
    rea(n); rea(q); S=sqrt(n);
    fail[0]=fail[1]=1; len[1]=-1; cnt=1;
    for(int i=0;i<26;i++)
        trans[0][i]=trans[0][i]=1;
    a[0]=-1;
    for(int i=1;i<=n;i++)
        while((a[i]=nc()-'a')>25 || a[i]<0);
    int m=0;
    memset(pos,0x7f,sizeof(pos));
    for(int i=0;i*S<n;i++,m++){
        ++it; p1=p2=0;
        for(int j=i*S+1;j<=n;j++){
            PushBack(i*S+1,j),pp[i][j]=p2;
            ans[i][j]=ans[i][j-1];
            if(vis[p1]^it){
                vis[p1]=it,ans[i][j]++;
                pos[i][p1]=j;
            }
        }
    }
    int lst=0;
    while(q--){
        int l,r; rea(l); rea(r); ++it;
        l^=ct*lst; r^=ct*lst;
        if(l>r) swap(l,r);
        if((l-1)/S==(r-1)/S){
            lst=0; p1=p2=0;
            for(int i=l;i<=r;i++){
                PushBack(l,i);
                if(vis[p1]^it)
                    vis[p1]=it,lst++;
            }
            printf("%d\n",lst);
            continue;
        }
        int L=(l-1)/S+1;
        lst=ans[L][r];
        p2=pp[L][r];
        for(int i=L*S;i>=l;i--){
            PushFront(i,r);
            if(vis[p2]^it){
                vis[p2]=it;
                if(pos[L][p2]>r) lst++;
            }
        }
        printf("%d\n",lst);
    }
    return 0;
}