后缀自动机在复习

后缀自动机，处理所有子串。

注意：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;
}