ABC141 E - Who Says a Pun?（SAM）

题意：

给定长度为n的串S，
现在你要在S串中找到两个相同的子串，满足这两个子串互不相交，
问能找到的串的最大长度是多少。

数据范围：n<=5e3

解法：

对串S建立后缀自动机，
拓扑排序求出记录每个节点的最左边出现的位置L和最右边出现的位置R。

对于节点x：
1.如果R-L>=maxlen[x]，那么该节点长度为[minlen[x],maxlen[x]]的串出现了至少两次，
那么可以用maxlen[x]更新答案。
2.如果R-L<maxlen[x]，那么该节点只有长度为[minlen[x]，R-L]的串出现了至少两次，
那么可以用R-L更新答案。

对上述情况取max就是最后的答案。

code：

#include<bits/stdc++.h>
using namespace std;
const int maxm=2e6+5;
char s[maxm];
int n;
struct SAM{
    int ch[maxm][26];
    int fa[maxm],l[maxm];//l[]是等价类的最长字符串长度len
    int last=1,tot=1;//tot是节点数量
    //
    int lc[maxm],rc[maxm];
    int cnt[maxm];
    int idx[maxm];
    int solve(){
        //基数排序
        for(int i=0;i<=tot;i++)cnt[i]=0;
        for(int i=1;i<=tot;i++)cnt[l[i]]++;
        for(int i=1;i<=tot;i++)cnt[i]+=cnt[i-1];
        for(int i=1;i<=tot;i++)idx[cnt[l[i]]--]=i;
        //逆向拓扑
        for(int i=tot;i>=1;i--){
            int x=idx[i];
            if(!lc[fa[x]]||lc[x]<lc[fa[x]]){
                lc[fa[x]]=lc[x];
            }
            if(!rc[fa[x]]||rc[x]>rc[fa[x]]){
                rc[fa[x]]=rc[x];
            }
        }
        //
        int ans=0;
        for(int i=1;i<=tot;i++){
            if(rc[i]-lc[i]>=l[i]){
                ans=max(ans,l[i]);
            }else{
                ans=max(ans,rc[i]-lc[i]);
            }
        }
        return ans;
    }
    //
    void add(int c,int pos){//传入s[i]-'a'而不是s[i]
        int p=last,np=++tot;last=np;l[np]=l[p]+1;
        for(;p&&!ch[p][c];p=fa[p])ch[p][c]=np;
        if(!p)fa[np]=1;
        else{
            int q=ch[p][c];
            if(l[p]+1==l[q])fa[np]=q;
            else{
                int nq=++tot;l[nq]=l[p]+1;
                memcpy(ch[nq],ch[q],sizeof ch[q]);
                fa[nq]=fa[q];fa[q]=fa[np]=nq;
                for(;ch[p][c]==q;p=fa[p])ch[p][c]=nq;
            }
        }
        //记录最左和最右位置
        lc[np]=rc[np]=pos;
    }
    void init(){
        for(int i=1;i<=tot;i++){
            memset(ch[i],0,sizeof ch[i]);
            fa[i]=l[i]=0;
        }
        last=tot=1;
    }
}S;
signed main(){
    scanf("%d",&n);
    scanf("%s",s+1);
    S.init();
    for(int i=1;i<=n;i++){
        S.add(s[i]-'a',i);
    }
    int ans=S.solve();
    cout<<ans<<endl;
    return 0;
}

---