字符串算法模板整理

1.字符串哈希：

ll H(ll* h,int l,int r) {return h[r+1]-h[l]*p[r-l+1];}
void getH(char* s,ll* h,int n) {h[0]=0; for(int i=1; i<=n; ++i)h[i]=h[i-1]*M+s[i-1];}

2.KMP：

void build(char* s,int n) {
    fa[0]=-1;
    for(int i=2,j=0; i<=n; ++i) {
        for(; ~j&&s[j]!=s[i-1]; j=fa[j]);
        fa[i]=++j;
    }
}

3.字典树/AC自动机：

int newnode() {int u=++tot; pre[u]=ed[u]=0; memset(go[u],0,sizeof go[u]); return u;}
void ins(char* s,int x) {
    int n=strlen(s),u=1;
    for(int i=0; i<n; u=go[u][id[s[i]]],++i)
        if(!go[u][id[s[i]]])go[u][id[s[i]]]=newnode();
    ed[u]|=1<<x;
}
void build() {
    queue<int> q;
    for(int i=0; i<4; ++i) {
        if(go[1][i])pre[go[1][i]]=1,q.push(go[1][i]);
        else go[1][i]=1;
    }
    while(q.size()) {
        int u=q.front();
        q.pop();
        ed[u]|=ed[pre[u]];
        for(int i=0; i<4; ++i) {
            if(go[u][i])pre[go[u][i]]=go[pre[u]][i],q.push(go[u][i]);
            else go[u][i]=go[pre[u]][i];
        }
    }
}
void init() {tot=0,newnode();}

4.后缀自动机：

int newnode(int l) {int u=++tot; mxl[u]=l,memset(go[u],0,sizeof go[u]); return u;}
void add(int ch) {
    int p=last,np=last=newnode(mxl[p]+1);
    for(; p&&!go[p][ch]; p=fa[p])go[p][ch]=np;
    if(!p)fa[np]=1;
    else {
        int q=go[p][ch];
        if(mxl[q]==mxl[p]+1)fa[np]=q;
        else {
            int nq=newnode(mxl[p]+1);
            memcpy(go[nq],go[q],sizeof go[q]);
            fa[nq]=fa[q],fa[q]=fa[np]=nq;
            for(; p&&go[p][ch]==q; p=fa[p])go[p][ch]=nq;
        }
    }
}

5.回文树：

struct PAM {
    int tot,last,go[N][M],fa[N],mxl[N],siz[N],rt[N],cc[N],nc;
    int newnode(int l) {int u=++tot; fa[u]=0,siz[u]=0,mxl[u]=l,memset(go[u],0,sizeof go[u]); return u;}
    int getfa(int u) {for(; u&&cc[nc]!=cc[nc-mxl[u]-1]; u=fa[u]); return u;}
    void add(int ch) {
        cc[++nc]=ch;
        int u=getfa(last);
        if(!go[u][ch]) {
            int v=newnode(mxl[u]+2),p=getfa(fa[u]);
            fa[v]=go[p][ch]?go[p][ch]:2;
            go[u][ch]=v;
        }
        rt[nc-1]=last=go[u][ch];
        siz[go[u][ch]]++;
    }
    void build(char* s,int n) {
        tot=nc=0,cc[0]=-1,last=newnode(-1),newnode(0),fa[2]=1;
        for(int i=0; i<n; ++i)add(s[i]-'a');
        for(int i=tot; i>=1; --i)siz[fa[i]]+=siz[i];
    }
    int getmxl(int r) {return mxl[rt[r]];}
};

6.Manacher：

void Manacher(char* s,int n) {
    Ma[0]='$',Ma[1]='#',rd[0]=0;
    for(int i=0; i<n; ++i)Ma[(i+1)<<1]=s[i],Ma[((i+1)<<1)+1]='#';
    for(int i=1,p=0; i<(n<<1)+2; ++i) {
        rd[i]=i<p+rd[p]?min(p+rd[p]-i,rd[(p<<1)-i]):1;
        for(; Ma[i+rd[i]]==Ma[i-rd[i]]; ++rd[i]);
        if(i+rd[i]>p+rd[p])p=i;
    }
    for(int i=1; i<(n<<1)+2; ++i)L[i]=R[i]=0;
    for(int i=1; i<(n<<1)+2; ++i) {
        L[i+rd[i]-1]=max(L[i+rd[i]-1],rd[i]*2-1);
        R[i-rd[i]+1]=max(R[i-rd[i]+1],rd[i]*2-1);
    }
    for(int i=n<<1; i>=1; --i)L[i]=max(L[i],L[i+1]-2);
    for(int i=2; i<(n<<1)+2; ++i)R[i]=max(R[i],R[i-1]-2);
}

7.EXKMP:

void getlen(char* a,char* b,int na,int nb,int* ex,int* len) {
    for(ex[0]=0; ex[0]<na&&ex[0]<nb&&a[ex[0]]==b[ex[0]]; ++ex[0]);
    for(int i=1,p=0; i<na; ++i) {
        ex[i]=i<=p+ex[p]?min(p+ex[p]-i,len[i-p]):0;
        for(; i+ex[i]<na&&ex[i]<nb&&a[i+ex[i]]==b[ex[i]]; ++ex[i]);
        if(i+ex[i]>p+ex[p])p=i;
    }
}
void exkmp(char* a,char* b,int na,int nb) {
    len[0]=nb;
    getlen(b+1,b,nb-1,nb,len+1,len);
    getlen(a,b,na,nb,ex,len);
}

8.后缀数组：

void Sort(int* x,int* y,int m) {
    for(int i=0; i<m; ++i)c[i]=0;
    for(int i=0; i<n; ++i)++c[x[i]];
    for(int i=1; i<m; ++i)c[i]+=c[i-1];
    for(int i=n-1; i>=0; --i)sa[--c[x[y[i]]]]=y[i];
}
void da(int* s,int n,int m=1000) {
    int *x=buf1,*y=buf2;
    x[n]=y[n]=-1;
    for(int i=0; i<n; ++i)x[i]=s[i],y[i]=i;
    Sort(x,y,m);
    for(int k=1; k<n; k<<=1) {
        int p=0;
        for(int i=n-k; i<n; ++i)y[p++]=i;
        for(int i=0; i<n; ++i)if(sa[i]>=k)y[p++]=sa[i]-k;
        Sort(x,y,m),p=1,y[sa[0]]=0;
        for(int i=1; i<n; ++i)y[sa[i]]=x[sa[i-1]]==x[sa[i]]&&x[sa[i-1]+k]==x[sa[i]+k]?p-1:p++;
        if(p==n)break;
        swap(x,y),m=p;
    }
}
void getht() {
    for(int i=0; i<n; ++i)rnk[sa[i]]=i;
    ht[0]=0;
    for(int i=0,k=0; i<n; ++i) {
        if(k)--k;
        if(!rnk[i])continue;
        for(; s[i+k]==s[sa[rnk[i]-1]+k]; ++k);
        ht[rnk[i]]=k;
    }
}
void initST() {
    for(int i=1; i<n; ++i)ST[i][0]=ht[i];
    for(int j=1; (1<<j)<=n; ++j)
        for(int i=1; i+(1<<j)-1<n; ++i)
            ST[i][j]=min(ST[i][j-1],ST[i+(1<<(j-1))][j-1]);
}
int lcp(int l,int r) {
    if(l==r)return n-sa[l];
    if(l>r)swap(l,r);
    l++;
    int k=Log[r-l+1];
    return min(ST[l][k],ST[r-(1<<k)+1][k]);
}

9.序列自动机：

void build() {
    for(int i=0; i<M; ++i)go[n][i]=go[n+1][i]=n+1;
    for(int i=n-1; i>=0; --i)memcpy(go[i],go[i+1],sizeof go[i]),go[i][s[i]-'a']=i+1;
}