后缀自动机模板
入門推薦博客:https://www.luogu.com.cn/blog/Kesdiael3/hou-zhui-zi-dong-ji-yang-xie?
另一個講的很清楚的博客:https://www.cnblogs.com/zjp-shadow/p/9218214.html
簡單介紹一下各個參數:構造SAM時強制在線,一個一個字母按照順序依次插入
fa(i)結點與 i 結點的關系:
len(fa(i))+1=minlen(i)//當前節點最小長度為minlen,最大長度為len,可以通過后綴鏈接求出minlen
一般如果需要從子節點向上更新信息的話,需要先對長度進行基數排序,達到滿足拓撲序
const int N=3e5+100;int tot,last,id[N<<1],tong[N<<1];struct Node {int ch[26];int fa,len; }st[N<<1];inline int newnode() {tot++;for(int i=0;i<26;i++)st[tot].ch[i]=0;st[tot].fa=st[tot].len=0;return tot; }void add(int x) {int p=last,np=last=newnode();st[np].len=st[p].len+1;while(p&&!st[p].ch[x])st[p].ch[x]=np,p=st[p].fa;if(!p)st[np].fa=1;else{int q=st[p].ch[x];if(st[p].len+1==st[q].len)st[np].fa=q;else{int nq=newnode();st[nq]=st[q]; st[nq].len=st[p].len+1;st[q].fa=st[np].fa=nq;while(p&&st[p].ch[x]==q)st[p].ch[x]=nq,p=st[p].fa;//向上把所有q都替換成nq}} }void radix_sort() {memset(tong,0,sizeof(tong));for(int i=1;i<=tot;i++)tong[st[i].len]++;for(int i=1;i<=tot;i++)tong[i]+=tong[i-1];for(int i=1;i<=tot;i++)id[tong[st[i].len]--]=i; }void init() {last=1;tot=0;newnode(); }總結
- 上一篇: HDU - 5394 Trie in T
- 下一篇: SPOJ - LCS Longest C