后綴自動機，處理所有子串。

注意：nq節點的建立，是因為新加了一個字母使原本為一個狀態的東西必須分為兩個狀態，例如ba到baa，a的出現次數變多，a狀態與ba狀態分離。baa的新pnt不能是ba，因為空節點可以通過兩種方式到達ba狀態，它的pnt的狀態就不唯一了。換句話說，ｂａ不能由ｂａａ截斷前綴得到，不符合ｐｎｔ的節點的后綴與兒子相同，則跳ｐｎｔ指針時會出現問題。所以建立ｎｑ節點，把是ｎｐ節點ｐｎｔ的狀態分離出來，讓ｎｑ＝ｎｐ的ｐｎｔ。因為ｎｑ節點比ｂａ節點短，所以ｂａ的ｐｎｔ＝ｎｐ的ｐｎｔ＝ｎｑ。

每個節點的轉移他們的ｐｎｔ都應該具備

題目大意

給定一個字符串，求每個長度的子串出現次數的最大值

tyvj 4123用pnt樹size從大往小更新，用并查集維護是否更新過

模板代碼

#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> using namespace std; #define maxn 200020struct node{int next[26],val;int pnt; }sam[maxn * 2]; struct node2{int sz,num;bool operator < (const node2 a) const{return sz > a.sz; } }data[maxn * 2]; char a[maxn]; int n,tot,last,root; int fa[maxn]; int num[maxn],ans[maxn],b[maxn];inline void add(int x){int np = ++tot,p = last;sam[np].val = sam[p].val + 1;data[np].sz = 1;while ( !sam[p].next[x] && p ) sam[p].next[x] = np, p = sam[p].pnt;int q = sam[p].next[x];if ( !q ) sam[np].pnt = root,sam[root].next[x] = np;else if ( q && sam[p].val + 1 == sam[q].val ) sam[np].pnt = q;else{int nq = ++tot;sam[nq].val = sam[p].val + 1;sam[nq].pnt = sam[q].pnt;sam[np].pnt = sam[q].pnt = nq;memcpy(sam[nq].next,sam[q].next,sizeof(sam[q].next));while ( p && sam[p].next[x] == q ) sam[p].next[x] = nq , p = sam[p].pnt;if ( sam[p].next[x] == q ) sam[p].next[x] = nq; //特判p跳到根后的情況}last = np; } inline void getsize(){for (int i = 0 ; i <= tot ; i++) num[sam[i].val]++;for (int i = 1 ; i <= n ; i++) num[i] += num[i - 1]; //從小到大for (int i = tot ; i >= 0 ; i--) b[--num[sam[i].val]] = i;for (int i = tot ; i >= 0 ; i--) if ( sam[b[i]].pnt != b[i] ) data[sam[b[i]].pnt].sz += data[b[i]].sz; } inline int getfa(int x){if ( x == fa[x] ) return x;return fa[x] = getfa(fa[x]); } int main(){scanf("%s",a + 1);n = (int) strlen(a + 1);for (int i = 1 ; i <= n ; i++) add(a[i] - 'a');getsize();for (int i = 1 ; i <= n + 2 ; i++) fa[i] = i;for (int i = 1 ; i <= tot ; i++) data[i].num = i;sort(data + 1,data + tot + 1);for (int i = 1 ; i <= tot ; i++){int x = data[i].num;for (int j = getfa(sam[sam[x].pnt].val + 1) ; j <= sam[x].val ; j = getfa(j + 1)){ans[j] = max(data[i].sz,ans[j]);fa[j] = fa[j + 1];}}for (int i = 1 ; i <= n ; i++) printf("%d\n",ans[i]);return 0; }

　　

?

轉載于:https://www.cnblogs.com/zqq123/p/5201312.html

總結

以上是生活随笔為你收集整理的后缀自动机在复习的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。