SPOJ 8222 Substrings（后缀自动机）

【题目链接】 http://www.spoj.com/problems/NSUBSTR/

【题目大意】

　　给一个字符串S，令F(x)表示S的所有长度为x的子串中，出现次数的最大值。

　　求出所有的F。

【题解】

　　在SAM中，一个串出现的次数就是|Right(s)|，我们按长度从小到大分配内存单位，
　　从后往前计算可以获得Right值大小，用所有的Right去更新相应长度的答案即可。

【代码】

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=500005;
char s
;
struct sam{
int p,q,np,nq,cnt,last,a
[26],l
,f
;
sam(){cnt=0;last=++cnt;}
void extend(int c){
p=last;np=last=++cnt;l[np]=l[p]+1;
while(!a[p][c]&&p)a[p][c]=np,p=f[p];
if(!p)f[np]=1;
else{
q=a[p][c];
if(l[p]+1==l[q])f[np]=q;
else{
nq=++cnt;l[nq]=l[p]+1;
memcpy(a[nq],a[q],sizeof(a[q]));
f[nq]=f[q]; f[np]=f[q]=nq;
while(a[p][c]==q)a[p][c]=nq,p=f[p];
}
}
}int b
,x
,r
;
void build(){
scanf("%s",s+1);
int len=strlen(s+1);
for(int i=1;i<=len;i++)extend(s[i]-'a');
for(int i=1;i<=cnt;i++)b[l[i]]++;
for(int i=1;i<=len;i++)b[i]+=b[i-1];
for(int i=1;i<=cnt;i++)x[b[l[i]]--]=i;
for(int i=p=1;i<=len;i++){p=a[p][s[i]-'a'];r[p]++;}
for(int i=cnt;i;i--)r[f[x[i]]]+=r[x[i]];
}
int F
;
void CalF(){
int len=strlen(s+1);
for(int i=1;i<=cnt;i++)F[l[i]]=max(F[l[i]],r[i]);
for(int i=1;i<=len;i++)printf("%d\n",F[i]);
}
}sam;
int main(){
sam.build();
sam.CalF();
return 0;
}