SPOJ NSUBSTR 后缀自动机+DP

题意

给一个长度为S的字符串，问长度为1-S的子串最多出现了多少次。

题解

DP其实是非常简单的。。。主要就是一个后缀自动机的模板，敲对了就行。至于算法，认真观察后缀自动机，对后缀自动机理解透彻了就发现这是一道模板题。对于这道题，我发现的规律就是，对于一个点我们很容易算出来出现次数，然后的话我们可以看一下这个点对应的字符串最长长度是多少，那么小于等于这个长度的所有字符串出现次数最少为这个点的出现次数（RT数组）。

代码

#include <bits/stdc++.h>
#define UP(i,l,h) for(int i=l;i<h;i++)
#define DOWN(i,h,l) for(int i=h-1;i>=l;i--)
#define W(t) while(t)
#define MEM(a,b) memset(a,b,sizeof(a))
#define MAXN 350010
using namespace std;
static const int NODE=MAXN<<1,C=26;
char s[MAXN];
int dp[MAXN];
struct SuffixAuto {
    int allc,n,last,par[NODE],len[NODE],tran[NODE][C],c[NODE],id[NODE],rt[NODE];
    int newNode() {
        int now=++allc;
        MEM(tran[now],0);
        return now;
    }
    void init() {
        MEM(c,0);
        MEM(id,0);
        MEM(rt,0);
        allc=0;
        last=newNode();
        par[last]=len[last]=0;
    }
    void extend(int c) {
        int p=last,np=newNode();
        len[np]=len[last]+1;
        rt[np]=1;
        for(; p&&!tran[p][c]; p=par[p]) tran[p][c]=np;
        if(!p) par[np]=1;
        else {
            int q=tran[p][c];
            if(len[q]==len[p]+1) par[np]=q;
            else {
                int nq=++allc;
                par[nq]=par[q];
                len[nq]=len[p]+1;
                memcpy(tran[nq],tran[q],sizeof(tran[q]));
                par[np]=par[q]=nq;
                for(tran[p][c]=nq,p=par[p]; p&&tran[p][c]==q; p=par[p]) tran[p][c]=nq;
            }
        }
        last=np;
    }
    void topsort() {
        UP(i,1,allc+1) c[len[i]]++;
        UP(i,1,n+1) c[i]+=c[i-1];
        UP(i,1,allc+1) id[c[len[i]]--]=i;
        DOWN(i,allc+1,1) rt[par[id[i]]]+=rt[id[i]];
    }
};
SuffixAuto sa;

int main() {
    W(~scanf("%s",s)) {
        MEM(dp,0);
        sa.init();
        int len=strlen(s);
        sa.n=len;
        UP(i,0,len) {
            sa.extend(s[i]-'a');
        }
        sa.topsort();
        UP(i,1,sa.allc+1) {
            int x=sa.id[i];
            dp[sa.len[x]]=max(dp[sa.len[x]],sa.rt[x]);
//            cout<<sa.rt[x]<<endl;
        }
        DOWN(i,len,1) dp[i]=max(dp[i+1],dp[i]);
        UP(i,1,len+1) printf("%d\n",dp[i]);
    }
}