KMP
略
Manacher
inline void manacher(){
int len=strlen(s);
now[0]='!',now[1]='$',n=1;
for(re int i=0;i<len;i++){
now[++n]=s[i],now[++n]='$';
}
now[n+1]='@';
mx=mid=1;
for(re int i=1;i<=n;i++){
if(i<mx)p[i]=min(mx-i,p[2*mid-i]);
else p[i]=1;
while(now[i-p[i]]==now[i+p[i]])p[i]++;
if(mx<i+p[i])mx=i+p[i],mid=i;
}
}
PAM
namespace PAM{
int nxt[N][27],fail[N],tot,len[N],s[N],cnt,last;
inline void init(){
tot=1,fail[0]=1,len[1]=-1,s[0]=-1;
}
inline int find(int p,int k){
while(s[k-len[p]-1]!=s[k])p=fail[p];
return p;
}
inline void insert(int c){
s[++cnt]=c;int p=find(last,cnt);
if(!nxt[p][c]){
fail[++tot]=nxt[find(fail[p],cnt)][c];
nxt[p][c]=tot,len[tot]=len[p]+2;
}
last=nxt[p][c];
}
}
AC自动机
namespace AC{
int fail[N],end[N],nxt[N][26],tot;
inline void insert(char c[],int p){
int len=strlen(c),u=0;
for(int i=0;i<len;i++){
if(!nxt[u][c[i]-'a'])nxt[u][c[i]-'a']=++tot;
u=nxt[u][c[i]-'a'];
}
end[u]=p;
}
inline void buildfail(){
queue<int>q;
for(int i=0;i<26;i++){
int v=nxt[0][i];
if(v)fail[v]=0,q.push(v);
}
while(!q.empty()){
int u=q.front();q.pop();
for(int i=0;i<26;i++){
int v=nxt[u][i];
if(v)fail[v]=nxt[fail[u]][i],q.push(v);
else nxt[u][i]=nxt[fail[u]][i];
}
}
}
}
统计时记得沿 f a i l fail fail跳所有子串
后缀数组+ST表
二分大法吼啊
namespace SA{
int n,m,sa[N],rk[N],cnt[N],sa2[N],st[22][N],ht[N],lg[N];
inline void bucsort(){
for(int i=1;i<=m;i++)cnt[i]=0;
for(int i=1;i<=n;i++)cnt[rk[sa2[i]]]++;
for(int i=1;i<=m;i++)cnt[i]+=cnt[i-1];
for(int i=n;i>=1;i--)sa[cnt[rk[sa2[i]]]--]=sa2[i];//n->1
}
inline void buildsa(){
m=128;//
for(int i=1;i<=n;i++)rk[i]=a[i],sa2[i]=i;
bucsort();
for(int i=1,pos;i<=n&&pos<n;i<<=1,pos=0){
pos=0;
for(int j=n-i+1;j<=n;j++)sa2[++pos]=j;//
for(int j=1;j<=n;j++)if(sa[j]>i)sa2[++pos]=sa[j]-i;
bucsort();
swap(rk,sa2);
rk[sa[1]]=1,pos=1;
for(int j=2;j<=n;j++)
rk[sa[j]]=((sa2[sa[j]]==sa2[sa[j-1]])&&(sa2[sa[j]+i]==sa2[sa[j-1]+i])?pos:++pos);//[j-1]
m=pos;
}
for(int i=1;i<=n;i++)rk[sa[i]]=i;
for(int i=1,k=0,j;i<=n;ht[rk[i++]]=k)
for(k?k--:0,j=sa[rk[i]-1];a[j+k]==a[i+k];k++);//a,rk[i]-1
}
void buildst(){lg[0]=-1;
for(int i=1;i<=n;i++)lg[i]=lg[i>>1]+1;
for(int i=1;i<=n;i++)st[0][i]=ht[i];
for(int i=1;(1<<i)<=n;i++)
for(int j=1;j+(1<<i)-1<=n;j++)
st[i][j]=min(st[i-1][j],st[i-1][j+(1<<(i-1))]);
}
inline int lcp(int x,int y){
x=rk[x],y=rk[y];
if(x>y)swap(x,y);x++;//x+1
int t=lg[y-x+1];
return min(st[t][x],st[t][y-(1<<t)+1]);
}
}
SAM
namespace SAM{
int nxt[N][26],len[N],link[N],siz[N],tot,last;
inline void init(){
last=tot=1;
}
inline void insert(int c){
int cur=++tot,p=last;last=cur;
siz[cur]=1,len[cur]=len[p]+1;//p!=last
for(;p&&!nxt[p][c];p=link[p])nxt[p][c]=cur;
if(!p)link[cur]=1;
else{
int q=nxt[p][c];
if(len[q]==len[p]+1)link[cur]=q;
else {
int clo=++tot;
len[clo]=len[p]+1;link[clo]=link[q];
memcpy(nxt[clo],nxt[q],sizeof(nxt[q]));
for(;p&&nxt[p][c]==q;p=link[p])nxt[p][c]=clo;
link[q]=link[cur]=clo;
}
}
}
int A[N],B[N];
inline void getsiz(){//维护endpos大小
for(int i=1;i<=tot;i++)A[len[i]]++;
for(int i=1;i<=tot;i++)A[i]+=A[i-1];
for(int i=1;i<=tot;i++)B[A[len[i]]--]=i;
for(int i=tot;i>=1;i--)siz[link[B[i]]]+=siz[B[i]];
}
}
理解 p a r e n t t r e e parent\ tree parent tree和 D A G DAG DAG意义