spoj 1812 LCS2(SAM+DP)

【题目链接】

  http://www.spoj.com/problems/LCS2/en/

【题意】

  求若干个串的最长公共子串。

【思路】

  SAM+DP

  先拿个串建个SAM,然后用后面的串匹配,每次将所有的匹配长度记录在状态上取min,然后对所有状态取max即答案。

  需要更新fa,因为fa[p]一定比p更优,但匹配的时候可能只更新了p而没有更新fa[p],所以还需要递推一边。

  注意mn[p]初始化为l[p]

【代码】

 #include<cstdio>
#include<cstring>
using namespace std; const int N = *1e5+; char s[N];
int sz,last,root,l[N],ch[N][],fa[N],mn[N],mx[N];
int b[N],cnt[N];
void add(int x) {
int c=s[x]-'a';
int p=last,np=++sz; last=np;
l[np]=mn[np]=x+;
for(;p&&!ch[p][c];p=fa[p]) ch[p][c]=np;
if(!p) fa[np]=root;
else {
int q=ch[p][c];
if(l[p]+==l[q]) fa[np]=q;
else {
int nq=++sz; l[nq]=mn[nq]=l[p]+;
memcpy(ch[nq],ch[q],sizeof(ch[q]));
fa[nq]=fa[q];
fa[np]=fa[q]=nq;
for(;ch[p][c]==q;p=fa[p]) ch[p][c]=nq;
}
}
} int main() {
root=last=++sz;
scanf("%s",s);
int len=strlen(s);
for(int i=;i<len;i++) add(i);
for(int i=;i<=sz;i++) cnt[l[i]]++;
for(int i=;i<=len;i++) cnt[i]+=cnt[i-];
for(int i=;i<=sz;i++) b[cnt[l[i]]--]=i;
while(scanf("%s",s)==) {
int p=root; len=;
for(int i=;s[i];i++) {
int c=s[i]-'a';
if(ch[p][c]) { len++; p=ch[p][c]; }
else {
while(p&&!ch[p][c]) p=fa[p];
if(!p) { len=; p=root; }
else { len=l[p]+; p=ch[p][c]; }
}
if(len>mx[p]) mx[p]=len;
}
for(int i=sz;i;i--) {
p=b[i];
if(mx[p]<mn[p]) mn[p]=mx[p];
if(fa[p] && mx[fa[p]]<mx[p]) mx[fa[p]]=mx[p];
mx[p]=;
}
}
int ans=;
for(int i=;i<=sz;i++)
if(mn[i]>ans) ans=mn[i];
printf("%d",ans);
return ;
}
上一篇:C++设计模式——观察者模式


下一篇:Jquery-处理iframe的高度自适应