后缀自动机(SAM)模板

 struct SAM{
     int ch[maxn][],fa[maxn],len[maxn],cnt,last;
     void Init()
     {
         memset(ch,,sizeof(ch));
         memset(fa,,sizeof(fa));
         last=cnt=;
     }
     void Add(int c)
     {
         int p=last,np=last=++cnt;
         len[np]=len[p]+;
         while(!ch[p][c]&&p){
             ch[p][c]=np;p=fa[p];
         }
         if(p==)
             fa[np]=;
         else{
             int q=ch[p][c];
             if(len[p]==len[q]-){
                 fa[np]=q;
             }
             else{
                 int nq=++cnt;len[nq]=len[p]+;
                 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=;
 char s[N];
 int fa[N],pos[N],vis[N],sa[N],rank[N];
 int son[N][],end[N],rht[N],lcp[N];
 int ch[N][],len[N],id[N],tot;
 int od[N],wv[N],lst,cnt;

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

 void Insert(int c){
     int p=lst,np=lst=++cnt;end[lst]=;
     id[len[np]=len[p]+]=np;rht[np]=;
     while(p&&!ch[p][c])ch[p][c]=np,p=fa[p];
     if(!p)fa[np]=;
     else{
         int q=ch[p][c],nq;
         if(len[q]==len[p]+)fa[np]=q;
         else{
             len[nq=++cnt]=len[p]+;
             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=;i<=cnt;i++)wv[len[i]]++;
     for(int i=;i<=cnt;i++)wv[i]+=wv[i-];
     for(int i=;i<=cnt;i++)od[wv[len[i]]--]=i;
     for(int i=cnt;i>=;i--)rht[fa[od[i]]]+=rht[od[i]];
 }

 void Build_Tree(){
     int l=strlen(s+);
     for(int i=l;i>=;i--)Insert(s[i]-'a');
     for(int i=l;i>=;i--)
         for(int x=id[i],p=l+;x&&!pos[x];x=fa[x])
             p-=len[x]-len[fa[x]],pos[x]=p;
     for(int x=;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]+]=++tot]=l-len[x]+;
     for(int i=;i<;i++)if(son[x][i])DFS(son[x][i],l);
 }

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

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

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