参考http://www.cnblogs.com/lidaxin/p/5198507.html
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 = 5*1e5+10;
char s[N];
int sz,last,root,l[N],ch[N][26],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+1;
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]+1==l[q]) fa[np]=q;
else {
int nq=++sz; l[nq]=mn[nq]=l[p]+1;
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=0;i<len;i++) add(i);
for(int i=1;i<=sz;i++) cnt[l[i]]++;
for(int i=1;i<=len;i++) cnt[i]+=cnt[i-1];
for(int i=1;i<=sz;i++) b[cnt[l[i]]--]=i;
while(scanf("%s",s)==1) {
int p=root; len=0;
for(int i=0;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=0; p=root; }
else { len=l[p]+1; 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]=0;
}
}
int ans=0;
for(int i=1;i<=sz;i++)
if(mn[i]>ans) ans=mn[i];
printf("%d",ans);
return 0;
}