码迷,mamicode.com
首页 > 其他好文 > 详细

codevs 3160 最长公共子串(SAM)

时间:2018-07-21 16:52:41      阅读:135      评论:0      收藏:0      [点我收藏+]

标签:i++   char   xtend   复杂度   name   def   code   自动   oid   

题解:

因为父亲节点是第一个right集合不同的后缀

所以我们用ac自动机匹配的方法向下找,不符合了就往父亲方向跳

时间复杂度O(n)

代码:

#include <bits/stdc++.h>
#define ll long long
#define rint register int
#define rep(i,h,t) for (rint i=h;i<=t;i++)
#define dep(i,t,h) for (rint i=t;i>=h;i--)
using namespace std;
const int N=3e6;
char s[N],s1[N];
int size[N],len[N],ch[N][26];
int lst=1,node=1,t[N],a[N],fa[N],ans;
void extend(int c)
{
  int f=lst,p=++node; lst=p;
  len[p]=len[f]+1; size[p]=1;
  while (f&&!ch[f][c]) ch[f][c]=p,f=fa[f];
  if (!f) { fa[p]=1; return;};
  int x=ch[f][c],y=++node;
  if (len[f]+1==len[x]) {fa[p]=x; node--;return;};
  len[y]=len[f]+1; fa[y]=fa[x]; fa[x]=fa[p]=y;
  memcpy(ch[y],ch[x],sizeof(ch[x]));
  while (f&&ch[f][c]==x) ch[f][c]=y,f=fa[f];
}
void js()
{
  int l=strlen(s1),now=1,t=0;
  rep(i,1,l)
  {
    int x=s1[i-1]-a;
    if (ch[now][x]) now=ch[now][x],++t;
    else 
    {
      while (now&&!ch[now][x]) now=fa[now];
      if (!now) now=1,t=0;
      else t=len[now]+1,now=ch[now][x];
    }
    ans=max(ans,t);
  }
}
int main()
{
  freopen("1.in","r",stdin);
  freopen("1.out","w",stdout);
  ios::sync_with_stdio(false);
  cin>>s;
  int l=strlen(s);
  rep(i,1,l) extend(s[i-1]-a);
  cin>>s1;
  js();
  cout<<ans<<endl;
  return 0;
}

 

codevs 3160 最长公共子串(SAM)

标签:i++   char   xtend   复杂度   name   def   code   自动   oid   

原文地址:https://www.cnblogs.com/yinwuxiao/p/9346667.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!