后缀自动机——Longest Common Substring

题解:

让我们匹配两个串的最长公共子串。求这个长度。我们把a串建立后缀自动机,然后按照b串依次去匹配,每次记录匹配值,并且取最大。在匹配的时候,如果当前状态可以转移 则答案++,否则的话我们就去找后缀链,看能不能匹配上,如果有代表有这个边的出边,则答案就是lenmax(now)+1,now是当前状态,如果最后都没有匹配上,答案为0,并状态从1重新开始匹配。

#include <bits/stdc++.h>
using namespace std;
const int maxn = 2e6+10;
char a[maxn],b[maxn];
struct SuffixAutoMation
{
    int last=1,cnt=1;
    int trans[maxn<<1][26],slink[maxn<<1],l[maxn<<1];
    int t[maxn<<1],a[maxn<<1];
    int num[maxn<<1];
    inline void add(int x)
    {
        int p = last,np = ++cnt;last = np;l[np] = l[p]+1;
        for(;p && !trans[p][x];p=slink[p]) trans[p][x] = np;
        if(!p) slink[np] = 1;
        else{
            int q = trans[p][x];
            if(l[p]+1 == l[q]) slink[np] = q;
            else
            {
                int nq = ++cnt;l[nq] = l[p]+1;
                memcpy(trans[nq],trans[q],sizeof(trans[q]));
                slink[nq] = slink[q];
                slink[q] = slink[np] = nq;
                for(;trans[p][x] == q; p = slink[p]) trans[p][x] = nq;
            }
        }
        num[np] = 1;
    }
    void solve()
    {
        int res=0,ans=0,len=strlen(b+1),u=1;
        for(int i=1;i<=len;i++){
            if(trans[u][b[i]-'a']){
                ans++;
                u=trans[u][b[i]-'a'];
            }else{
                while(u&&!trans[u][b[i]-'a']) u=slink[u];
                if(!u) u=1,ans=0;
                else{
                    ans=l[u]+1;
                    u=trans[u][b[i]-'a'];
                }
            }
            res=max(res,ans);
        }
        cout<<res<<endl;
    }
}sam;
signed main()
{
    scanf("%s",a+1);
    scanf("%s",b+1);
    int len=strlen(a+1);
    for(int i=1;i<=len;i++) sam.add(a[i]-'a');
    sam.solve();
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值