SPOJ 8222 Substrings（后缀自动机）

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

题意：给定一个串S，定义函数f(x)表示长度为x的所有子串中出现次数最多的子串出现的次数。输出f(i)(1<=i<=len(S))。

思路：建立串S的后缀自动机。然后就是对于每个节点u，从u开始到达接收状态的路径的个数就是u节点对应串的出现次数。

const int KIND=26;
struct SAM
{
    SAM *son[KIND],*pre;
    int len,cnt;

    void init()
    {
        clr(son,NULL);
        pre=NULL;
        len=0;
    }
};

SAM sam[N],*head,*last,*b[N];
int d[N],cnt,f[N],len;
char s[N];

void initSAM()
{
    sam[0].init();
    head=last=&sam[0];
    cnt=1;
}


void insert(int x)
{
    SAM *p=&sam[cnt++],*u=last;
    p->len=last->len+1;
    last=p;
    for(;u&&!u->son[x];u=u->pre) u->son[x]=p;
    if(!u) p->pre=head;
    else if(u->son[x]->len==u->len+1) p->pre=u->son[x];
    else
    {
        SAM *r=&sam[cnt++],*q=u->son[x];
        *r=*q; r->len=u->len+1;
        p->pre=q->pre=r;
        for(;u&&u->son[x]==q;u=u->pre) u->son[x]=r;
    }
}

int main()
{
    RD(s);
    len=strlen(s);
    initSAM();
    int i;
    FOR0(i,len) insert(s[i]-'a');
    FOR0(i,cnt) d[sam[i].len]++;
    FOR1(i,len) d[i]+=d[i-1];
    FOR0(i,cnt) b[--d[sam[i].len]]=&sam[i];
    last=head;
    FOR0(i,len)
    {
        last=last->son[s[i]-'a'];
        last->cnt=1;
    }
    FORL1(i,cnt-1)
    {
        f[b[i]->len]=max(f[b[i]->len],b[i]->cnt);
        b[i]->pre->cnt+=b[i]->cnt;
    }
    FORL1(i,len-1) f[i]=max(f[i],f[i+1]);
    FOR1(i,len) PR(f[i]);
    return 0;

}