后缀自动机(SAM)模板

struct SAM{
    int ch[maxn][26],fa[maxn],len[maxn],cnt,last;
    void Init()
    {
        memset(ch,0,sizeof(ch));
        memset(fa,0,sizeof(fa));
        last=cnt=1;
    }
    void Add(int c)
    {
        int p=last,np=last=++cnt;
        len[np]=len[p]+1;
        while(!ch[p][c]&&p){
            ch[p][c]=np;p=fa[p];
        }
        if(p==0)
            fa[np]=1;
        else{
            int q=ch[p][c];
            if(len[p]==len[q]-1){
                fa[np]=q;
            }
            else{
                int nq=++cnt;len[nq]=len[p]+1;
                memcpy(ch[nq],ch[q],sizeof(ch[q]));
                fa[nq]=fa[q];fa[q]=fa[np]=nq;
                while(ch[p][c]==q&&p){
                    ch[p][c]=nq;
                    p=fa[p];
                }
            }
        }
    }
};

 　　然后是自整理最全SAM模版，超级大杀器！

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int N=200010;
char s[N];
int fa[N],pos[N],vis[N],sa[N],rank[N];
int son[N][26],end[N],rht[N],lcp[N];
int ch[N][26],len[N],id[N],tot;
int od[N],wv[N],lst,cnt;

void Init(){
    memset(ch,0,sizeof(ch));
    memset(end,0,sizeof(end));
    memset(son,0,sizeof(son));
    memset(vis,0,sizeof(vis));
    lst=cnt=1;tot=0;
}

void Insert(int c){
    int p=lst,np=lst=++cnt;end[lst]=1;
    id[len[np]=len[p]+1]=np;rht[np]=1;
    while(p&&!ch[p][c])ch[p][c]=np,p=fa[p];
    if(!p)fa[np]=1;
    else{
        int q=ch[p][c],nq;
        if(len[q]==len[p]+1)fa[np]=q;
        else{
            len[nq=++cnt]=len[p]+1;
            fa[nq]=fa[q];fa[q]=fa[np]=nq;
            memcpy(ch[nq],ch[q],sizeof(ch[q]));
            while(ch[p][c]==q)ch[p][c]=nq,p=fa[p];
        }
    }
}

void Get_Right(){
    for(int i=1;i<=cnt;i++)wv[len[i]]++;
    for(int i=1;i<=cnt;i++)wv[i]+=wv[i-1];
    for(int i=1;i<=cnt;i++)od[wv[len[i]]--]=i;
    for(int i=cnt;i>=1;i--)rht[fa[od[i]]]+=rht[od[i]];
}

void Build_Tree(){
    int l=strlen(s+1);
    for(int i=l;i>=1;i--)Insert(s[i]-'a');
    for(int i=l;i>=1;i--)
        for(int x=id[i],p=l+1;x&&!pos[x];x=fa[x])
            p-=len[x]-len[fa[x]],pos[x]=p;
    for(int x=2;x<=cnt;x++)son[fa[x]][s[pos[x]]-'a']=x;    
}

void DFS(int x,int l){
    if(end[x])sa[rank[l-len[x]+1]=++tot]=l-len[x]+1;
    for(int i=0;i<26;i++)if(son[x][i])DFS(son[x][i],l);
}

void Build_SA(){
    int l=strlen(s+1),k=0;DFS(1,l);
    for(int i=1,j;i<=l;lcp[rank[i++]]=k)
        for(k?k--:k,j=sa[rank[i]-1];s[i+k]==s[j+k];k++);
}

int main(){
    //freopen();
    //freopen();
    Init();
    scanf("%s",s+1);
    int l=strlen(s+1);
    Build_Tree();Build_SA();
    for(int i=1;i<=l;i++)printf("%d ",sa[i]);printf("
");
    for(int i=2;i<=l;i++)printf("%d ",lcp[i]);printf("
");
    return 0;
}

 　　警告：这里的SAM都是naive-sam，无法跑trie和多串。