后缀自动机模板

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
const int N=2000010;

int tot=1, last=1;
struct Node{
    int len, fa;
    int ch[26];
}node[N];
char str[N];
ll f[N], ans;
int head[N], ver[N], nxt[N], idx=1;

void extend(int c){
    int p=last, np=last=++tot;
    f[tot]=1;
    node[np].len=node[p].len+1;
    for(; p && !node[p].ch[c]; p=node[p].fa)
        node[p].ch[c]=np;
    if(!p) node[np].fa=1;
    else {
        int q=node[p].ch[c];
        if(node[q].len==node[p].len+1) 
            node[np].fa=q;
        else {
            int nq=++tot;
            node[nq]=node[q], node[nq].len=node[p].len+1;
            node[q].fa=node[np].fa=nq;
            for( ; p && node[p].ch[c]==q; p=node[p].fa)
                node[p].ch[c]=nq;
        }
    }
}

void add(int a, int b){
    ver[++idx]=b, nxt[idx]=head[a], head[a]=idx;
}

void dfs(int x){
    for(int i=head[x]; i; i=nxt[i]){
        dfs(ver[i]);
        f[x]+=f[ver[i]];
    }
    if(f[x]>1) ans=max(ans, f[x]*node[x].len);
}

int main(){
    scanf("%s", str);
    for(int i=0; str[i]; i++) extend(str[i]-'a');
    for(int i=2; i<=tot; i++) add(node[i].fa, i);
    dfs(1);
    cout<<ans<<endl;
    return 0;
}