回文自动机那套理论可以在 翁文涛的集训队论文 里看
每个串的回文自动机形态是固定的,一个子串的回文自动机也就是一个子图
那么就可以用记忆化的方式来降低复杂度……我感觉好像就是论文里提到的“不基于势能分析的插入方法“
#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;
}