【SPOJ - LCS2】Longest Common Substring II【SAM】

时间:2023-03-09 14:53:20
【SPOJ - LCS2】Longest Common Substring II【SAM】

题意

求出多个串的最长公共子串。

分析

 刚学SAM想做这个题的话最好先去做一下那道codevs3160。求两个串的LCS应该怎么求?把一个串s1建自动机,然后跑另一个串s2,然后找出s2每个前缀的最长公共后缀。那么多个的时候,我们也用这种类似的方法,但是我们求最长公共后缀的时候要求第一个串的。我们把其中一个串建SAM,然后把其他的串都在上面跑,维护两个值,Max[u]和Min[u]。自动机中每个状态u的Right存的是结尾集合。那么对于一个字符串,我们可以求出他和自动机中每个状态的最长公共后缀。然后,我们通过Max[fa[u]]=max(Max[fa[u]],Max[u])来确定左右状态的最长公共后缀,然后更新Min[o]。

下面的代码没有AC,但是我找了好久BUG没找到··如果有人看完并且找到了bug麻烦跟我说一下万分感谢!

  

 #include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm> using namespace std;
const int maxn=+;
char s[maxn];
struct state{
int len,link;
int next[];
}st[*maxn];
int n,last,cur,sz;
int Min[*maxn],Max[*maxn],c[*maxn]; void init(){
sz=;
cur=last=;
st[].link=-;
st[].len=;
} void build_sam(int c){
cur=sz++;
st[cur].len=st[last].len+;
Min[cur]=st[cur].len;
int p;
for(p=last;p!=-&&st[p].next[c]==;p=st[p].link)
st[p].next[c]=cur;
if(p==-)
st[cur].link=;
else{
int q=st[p].next[c];
if(st[q].len==st[p].len+)
st[cur].link=q;
else{
int clone=sz++;
st[clone].len=st[p].len+;
Min[clone]=st[clone].len;
st[clone].link=st[q].link;
for(int i=;i<;i++)
st[clone].next[i]=st[q].next[i];
for(;p!=-&&st[p].next[c]==q;p=st[p].link)
st[p].next[c]=clone;
st[cur].link=st[q].link=clone;
}
}
last=cur;
} int cmp(int a,int b){
return st[a].len<st[b].len;
} int main(){
scanf("%s",s);
n=strlen(s);
init(); for(int i=;i<n;i++)
build_sam(s[i]-'a');
for(int i=;i<sz;i++)
c[i]=i;
sort(c,c+sz,cmp); while(scanf("%s",s)!=EOF){
n=strlen(s);
int u=,len=;
for(int i=;i<n;i++){
int c=s[i]-'a';
if(u!=-&&st[u].next[c]==)
u=st[u].link,len=st[u].len;
if(u==-)
u=,len=;
else{
u=st[u].next[c];
len++;
Max[u]=max(Max[u],len);
}
} for(int i=sz-;i>=;i--){
int o=c[i];
Min[o]=min(Min[o],Max[o]);
if(st[o].link!=-){
Max[st[o].link]=max(Max[st[o].link],Max[o]);
}
Max[o]=;
}
}
int ans=;
for(int i=;i<sz;i++){
ans=max(ans,Min[i]);
}
printf("%d\n",ans); return ;
}