LG3804 【模板】后缀自动机

题意

给定一个只包含小写字母的字符串S,

请你求出 SS 的所有出现次数不为 1 的子串的出现次数乘上该子串长度的最大值。

对于100%的数据,|S|<=10^6

分析

学习资料:OI wikiMenci's Blog,张天杨《后缀自动机及其应用》。

我们令叶子节点的size=1.暴力建出parent树然后dfs,求出每个节点的right集合size,然后求len×size的最大值就行了。

时间复杂度:\(O(|S|)\)

代码

#include<bits/stdc++.h>
#define rg register
#define il inline
#define co const
template<typename T>il T read()
{
    rg T data=0,w=1;
    rg char ch=getchar();
    while(!isdigit(ch))
    {
        if(ch=='-') w=-1;
        ch=getchar();
    }
    while(isdigit(ch))
    {
        data=data*10+ch-'0';
        ch=getchar();
    }
    return data*w;
}
template<typename T>il T read(rg T&x)
{
    return x=read<T>();
}
typedef long long ll;

co int N=2e6;
namespace SAM
{
    int tot,last;
    int ch[N][26],fail[N],len[N],siz[N];
    void init()
    {
        tot=last=0;
        fail[0]=-1,len[0]=0;
    }
    void extend(int k)
    {
        int cur=++tot;
        len[cur]=len[last]+1,siz[cur]=1;
        int p=last;
        while(~p&&!ch[p][k])
        {
            ch[p][k]=cur;
            p=fail[p];
        }
        if(p==-1)
            fail[cur]=0;
        else
        {
            int q=ch[p][k];
            if(len[p]+1==len[q])
                fail[cur]=q;
            else
            {
                int clone=++tot;
                std::copy(ch[q],ch[q]+26,ch[clone]);
                fail[clone]=fail[q],len[clone]=len[p]+1;
                while(~p&&ch[p][k]==q)
                {
                    ch[p][k]=clone;
                    p=fail[p];
                }
                fail[q]=fail[cur]=clone;
            }
        }
        last=cur;
    }
    void ins(char*s,int n)
    {
        for(int i=0;i<n;++i)
            extend(s[i]-'a');
    }
    
    int nx[N],to[N];
    ll ans;
    void build()
    {
        for(int i=1;i<=tot;++i)
            nx[i]=to[fail[i]],to[fail[i]]=i;
    }
    void dfs(int x)
    {
        for(int i=to[x];i;i=nx[i])
        {
            dfs(i);
            siz[x]+=siz[i];
        }
        if(siz[x]>1)
            ans=std::max(ans,(ll)siz[x]*len[x]);
    }
    void solve()
    {
        build();
        dfs(0);
        printf("%lld\n",ans);
    }
}
char buf[N];

int main()
{
//  freopen(".in","r",stdin);
//  freopen(".out","w",stdout);
    SAM::init();
    scanf("%s",buf);
    SAM::ins(buf,strlen(buf));
    SAM::solve();
    return 0;
}

猜你喜欢

转载自www.cnblogs.com/autoint/p/10345193.html