luogu3804 【模板】后缀自动机

http://www.elijahqi.win/archives/2992
题目描述
给定一个只包含小写字母的字符串
S ,
请你求出
S 的所有出现次数不为
1 的子串的出现次数乘上该子串长度的最大值。
输入输出格式
输入格式：
一行一个仅包含小写字母的字符串
S
输出格式：
一个整数，为 所求答案
输入输出样例
输入样例#1： 复制
abab
输出样例#1： 复制
4
说明
对于
10% 的数据，
∣S∣<=1000
对于
100% 的数据，
∣S∣<=106
拓扑排序统计出每个子串出现的个数 然后直接跑一遍统计即可

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 2200000
#define ll long long
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline void read_s(char *s){
    int op=0;char ch=gc();
    while(ch<'a'||ch>'z') ch=gc();
    while(ch<='z'&&ch>='a') s[op++]=ch,ch=gc();
}
char s[N>>1];ll ans;
int len[N],fa[N],ch[N][26];
int cnt=1,root=1,last=1,c[N>>1],size[N],rk[N];
inline void insert1(int x){
    int np=++cnt,p=last;len[np]=len[p]+1;
    for (;p&&!ch[p][x];p=fa[p]) ch[p][x]=np;
    if (!p) fa[np]=root;else{
        int q=ch[p][x];if(len[p]+1==len[q]) fa[np]=q;else{
            int nq=++cnt;fa[nq]=fa[q];memcpy(ch[nq],ch[q],sizeof(ch[q]));
            fa[q]=fa[np]=nq;len[nq]=len[p]+1;
            for (;ch[p][x]==q;p=fa[p]) ch[p][x]=nq;
        }
    }
    last=np;size[np]=1;
}
int main(){
//  freopen("luogu3804.in","r",stdin);
    read_s(s);int ln=strlen(s);
    for (int i=0;i<ln;++i) insert1(s[i]-'a');
    for (int i=1;i<=cnt;++i) ++c[len[i]];
    for (int i=1;i<=ln;++i) c[i]+=c[i-1];
    for (int i=1;i<=cnt;++i) rk[c[len[i]]--]=i;
    for (int i=cnt;i;--i){
        int p=rk[i];
        size[fa[p]]+=size[p];
        if (size[p]>1) ans=max(ans,(ll)size[p]*len[p]);
    }printf("%lld\n",ans);
    return 0;
}