#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cstdlib>
using namespace std;
struct suffix_automaton{
int fa,son[28],l,key;
} a[610000];
int i,j,k,l,n,m,len,root,d[300000],rnk[600000],f[300000];
char c[300000];
void build(){
int r,p,q,np,last=root,x;
for(i=1;i<=l;++i){
x=c[i]-96;
np=++len; a[np].l=i; a[np].key=1;
for(p=last;p&&!a[p].son[x];p=a[p].fa)a[p].son[x]=np;
last=np;
if(!p)a[np].fa=root;
else{
q=a[p].son[x];
if(a[p].l+1==a[q].l)a[np].fa=q;
else{
r=++len;
for(j=1;j<=26;++j)a[r].son[j]=a[q].son[j];
a[r].fa=a[q].fa;
a[r].l=a[p].l+1;
a[q].fa=a[np].fa=r;
for(;p&&a[p].son[x]==q;p=a[p].fa)a[p].son[x]=r;
}
}
}
}
int main(){
root=++len;
scanf("%s",c+1);
l=strlen(c+1);
build();
for(i=1;i<=len;++i)d[a[i].l]++;
for(i=1;i<=l;++i)d[i]+=d[i-1],f[i]=1;
for(i=1;i<=len;++i)rnk[d[a[i].l]--]=i;//长度rank[x]对应哪个节点
for(j=len;j>=1;--j){ //按长度更新
i=rnk[j];
f[a[i].l]=max(f[a[i].l],a[i].key);
a[a[i].fa].key+=a[i].key;
}
for(i=l-1;i>=1;--i)f[i]=max(f[i],f[i+1]);//看解释
for(i=1;i<=l;++i)printf("%d\n",f[i]);
}
这题的最后一步没有理解好,就是之前已经遍历每个点, 在从长到短序的所有子串更新完答案F数组后。然而最后要从大道小再用f[i+1]更新f[i],本来觉得sam中每个节点不是已经能代表所有子串了吗,还要更新干嘛。然而有个问题,np的父亲是p,就是从st[np].len到st[p].len中间是有很长一段不管的。所以有些长度出现最多的可能更新不到。。哎呀感觉没有图说不清楚。网上很多解释是说什么长的子串一定比短的多出现。其实是有一定道理的。