SPOJ

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