SPOJ 1812 LCS2 [后缀自动机 DP]

题意:

求多个串<=10的最长连续子串


一个串建SAM,然后其他串在上面走

每个状态记录所有串在这个状态的公共子串的最小值

一个串在上面走的时候记录与每个状态公共子串的最大值,注意出现次数向父亲传递,一个状态能到达说明了Suffix Link指向的状态可以取到最大子串,这一步对val后基数排序然后倒着更新就行了

注意两点:

1.空间要开两倍,基数排序用的东西也要两倍哦

2.答案不能用mn[1]更新!!!!因为mn[1]没有意义啊root状态一个子串也没有还可能有某些bug

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=2e5+;
int n,m;
char s[N];
struct State{
int ch[],par,val;
}t[N<<];
int sz,root,last;
inline int nw(int _){t[++sz].val=_;return sz;}
void iniSAM(){sz=;root=last=nw();}
void extend(int c){
int p=last,np=nw(t[p].val+);
while(p&&t[p].ch[c]==) t[p].ch[c]=np,p=t[p].par;
if(p==) t[np].par=root;
else{
int q=t[p].ch[c];
if(t[q].val==t[p].val+) t[np].par=q;
else{
int nq=nw(t[p].val+);
memcpy(t[nq].ch,t[q].ch,sizeof(t[q].ch));
t[nq].par=t[q].par;
t[q].par=t[np].par=nq;
while(p&&t[p].ch[c]==q) t[p].ch[c]=nq,p=t[p].par;
}
}
last=np;
}
int mx[N<<],mn[N<<];
void walk(char s[]){
int sum=,u=root,n=strlen(s+);
for(int i=;i<=n;i++){
int c=s[i]-'a';
if(t[u].ch[c]) u=t[u].ch[c],mx[u]=max(mx[u],++sum);
else{
while(u&&!t[u].ch[c]) u=t[u].par;
if(!u) u=root,sum=;
else sum=t[u].val,u=t[u].ch[c],mx[u]=max(mx[u],++sum);
}
}
} int c[N<<],a[N<<],ans;
void RadixSort(int len){
for(int i=;i<=sz;i++) c[t[i].val]++;
for(int i=;i<=len;i++) c[i]+=c[i-];
for(int i=sz;i>=;i--) a[c[t[i].val]--]=i;
for(int i=;i<=sz;i++) mn[i]=t[i].val;
}
void solve(){
iniSAM();
scanf("%s",s+);
int len=strlen(s+);
for(int i=;i<=len;i++) extend(s[i]-'a');
RadixSort(len);
while(scanf("%s",s+)!=EOF){
walk(s);
for(int j=sz;j>=;j--){
int u=a[j];
mn[u]=min(mn[u],mx[u]);
if(mx[u]&&t[u].par) mx[t[u].par]=t[t[u].par].val;
mx[u]=;
}
}
for(int i=;i<=sz;i++) ans=max(ans,mn[i]);
printf("%d",ans);
}
int main(){
//freopen("in","r",stdin);
solve();
}
上一篇:SQL Server,MySql,Oracle数据库的默认端口号


下一篇:mysql之系统默认数据库