2018.12.15 spoj Substrings（后缀自动机）

传送门

后缀自动机基础题。

求长度为iii的子串出现次数的最大值。

---

对原串建出samsamsam，然后用sizsizsiz更新每个maxlenmaxlenmaxlen的答案。

然后由于后缀链接将其转化成了一种树形结构，因此直接在上面树形递推即可。

代码

#include<bits/stdc++.h>
#define ri register int
using namespace std;
const int N=5e5+5;
int n;
char s[N];
struct SAM{
	int f[N],siz[N],cnt[N],rk[N],len[N],link[N],son[N][26],tot,rt,last;
	SAM(){last=rt=tot=1,len[0]=-1,fill(son[0],son[0]+26,1);}
	inline void expend(int x){
		int p=last,np=++tot;
		siz[last=np]=1,len[np]=len[p]+1;
		while(p&&!son[p][x])son[p][x]=np,p=link[p];
		if(!p){link[np]=rt;return;}
		int q=son[p][x],nq;
		if(len[q]==len[p]+1){link[np]=q;return;}
		len[nq=++tot]=len[p]+1,memcpy(son[nq],son[q],sizeof(son[q])),link[nq]=link[q];
		while(p&&son[p][x]==q)son[p][x]=nq,p=link[p];
		link[np]=link[q]=nq;
	}
	inline void topsort(){
		for(ri i=1;i<=tot;++i)++cnt[len[i]];
		for(ri i=1;i<=last;++i)cnt[i]+=cnt[i-1];
		for(ri i=1;i<=tot;++i)rk[cnt[len[i]]--]=i;
		for(ri i=tot;i;--i)siz[link[rk[i]]]+=siz[rk[i]];
		for(ri i=1;i<=tot;++i)f[len[i]]=max(f[len[i]],siz[i]);
		for(ri i=n-1;i;--i)f[i]=max(f[i],f[i+1]);
		for(ri i=1;i<=n;++i)cout<<f[i]<<'\n';
	}
}sam;
int main(){
	freopen("lx.in","r",stdin);
	scanf("%s",s+1),n=strlen(s+1);
	for(ri i=1;i<=n;++i)sam.expend(s[i]-'a');
	sam.topsort();
	return 0;
}