SPOJ - NSUBSTR

#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中间是有很长一段不管的。所以有些长度出现最多的可能更新不到。。哎呀感觉没有图说不清楚。网上很多解释是说什么长的子串一定比短的多出现。其实是有一定道理的。