[hihocoder1449]后缀自动机

hihocoder1449

题目描述：

​ 小Hi平时的一大兴趣爱好就是演奏钢琴。我们知道一个音乐旋律被表示为一段数构成的数列。

现在小Hi想知道一部作品中所有长度为K的旋律中出现次数最多的旋律的出现次数。但是K不是固定的，

小Hi想知道对于所有的K的答案。

分析：

• 求所有长度为K的子串中出现次数最多的串出现的次数。注意，我们构造出后缀自动机以后，一个点代表一个等价类，同时我们还知道这个等价类的长度在minlen~maxlen之间。假如我们知道每一个等价类出现的次数，就可以来更新答案。
• 其实就是在parent树上dfs一下。这样所有的等价类次数都可以求出来。
• 现在，如何求ans数组呢？我们如果直接把每个等价类长度的区间答案尝试更新，复杂度最坏到$n^2$。其实，我们直接把每个等价类的最大长度尝试更新，然后长度小的最值可以由长度大的最值更新过来。这样一边循环就好了。
• 什么？你问我为什么？因为答案具有单调性啊233……

Coding

#include<bits/stdc++.h>
#define ll long long
#define rint register int
using namespace std;
const int N=1e6+10;
int n,tot=1,la=1,ne[N*2][30],fa[N*2],len[N*2],siz[N*2],c[N*2],A[N*2];
char s[N];int ans[N];
void add(int c){
	int p=la,np=la=++tot;
	len[np]=len[p]+1;siz[np]=1;
	for(;p&&!ne[p][c];p=fa[p]) ne[p][c]=np;
	if(!p) fa[np]=1;
	else{
		int q=ne[p][c];
		if(len[p]+1==len[q]) fa[np]=q;
		else{
			int nq=++tot;
			memcpy(ne[nq],ne[q],sizeof(ne[q]));
			fa[nq]=fa[q];fa[q]=fa[np]=nq;
			len[nq]=len[p]+1;
			for(;p&&ne[p][c]==q;p=fa[p]) ne[p][c]=nq;
		}
	}
}
int main(){
	scanf("%s",s+1);
	n=strlen(s+1);
	for(int i=1;i<=n;++i) add(s[i]-'a');
	for(int i=1;i<=tot;++i) c[len[i]]++;
	for(int i=1;i<=tot;++i) c[i]+=c[i-1];
	for(int i=1;i<=tot;++i) A[c[len[i]]--]=i;
	for(int i=tot;i>0;--i){
		int p=A[i];
		siz[fa[p]]+=siz[p];
	}
	for(int i=1;i<=tot;++i) ans[len[i]]=max(ans[len[i]],siz[i]);
	for(int i=n-1;i>=1;--i) ans[i]=max(ans[i],ans[i+1]);
	for(int i=1;i<=n;++i) cout<<ans[i]<<endl;
	return 0;
}