题意:多个字符串的最长公共子串
解法:SAM
- 将第一个串建立自动机
- 后面的每个字符串就行匹配,维护每个位置的值,然后所有串中此位置取最小值
- 拓扑序,对父节点进行更新
- /*************************************************************************
- > File Name: lcs2.cpp
- > Author: cy
- > Mail: 1002@qq.com
- > Created Time: 14/11/25 20:48:22
- ************************************************************************/
- #include<iostream>
- #include<cstring>
- #include <algorithm>
- #include<cstdlib>
- #include<vector>
- #include<cmath>
- #include<stdlib.h>
- #include<iomanip>
- #include<list>
- #include<deque>
- #include<map>
- #include <stdio.h>
- #include <queue>
- #define inf 0x3f3f3f3f
- #define INF 0x3FFFFFFFFFFFFFFFLL
- #define rep(i,n) for(i=0;i<n;i++)
- #define reP(i,n) for(i=1;i<=n;i++)
- #define ull unsigned long long
- #define ll long long
- #define cle(a) memset(a,0,sizeof(a))
- #define MAXN 500010+5
- #define MAXC 26
- using namespace std;
- char ch[MAXN];
- int rk[MAXN],sz[MAXN],ans[MAXN];
- struct Suffix_Automaton{
- int child[MAXN][MAXC];
- int f[MAXN];//指向最近的接受态
- int l[MAXN];//此节点代表的最长串的长度
- int head,tail,tot;//头节点,最后的接受态,总节点
- void clear(){
- memset(child,0,sizeof(child));
- memset(f,0,sizeof(f));
- memset(l,0,sizeof(l));
- head=tail=tot=1;
- }
- void add(int c){
- int p=tail,np=++tot;l[np]=l[p]+1;tail=np;
- for(;p&&!child[p][c];p=f[p]) child[p][c]=np;
- if(!p) f[np]=head;
- else if(l[child[p][c]]==l[p]+1) f[np]=child[p][c];
- else{
- int q=child[p][c],r=++tot;
- memcpy(child[r],child[q],MAXC*sizeof(int));
- f[r]=f[q];l[r]=l[p]+1;f[q]=f[np]=r;
- for(;p&&child[p][c]==q;p=f[p]) child[p][c]=r;
- }
- }
- void insert(char*ch){
- int l=strlen(ch);
- for(int i=0;i<l;i++) add(ch[i]-'a');
- }
- int ws[MAXN];
- void tpu(){
- int len=strlen(ch);
- for(int i=0;i<=len;i++)ws[i]=0;
- for(int i=1;i<=tot;i++)ws[l[i]]++;
- for(int i=1;i<=len;i++)ws[i]+=ws[i-1];
- for(int i=tot;i>0;i--)rk[ws[l[i]]--]=i;
- }
- int ans[MAXN];
- //为什么要沿着父节点更新?
- int now[MAXN];
- void solve(){
- for(int i=1;i<=tot;i++)//最长可能性
- {
- now[i]=l[i];
- }
- while(scanf("%s",&ch)!=EOF){
- cle(ans);
- int len=strlen(ch);
- int p=head;
- int temp=0;
- for(int i=0;i<len;i++)
- {
- int tt=ch[i]-'a';
- if(child[p][tt])temp++,p=child[p][tt];
- else{
- while(true){
- if(!p)break;
- if(child[p][tt])break;
- p=f[p];
- }
- if(!p){
- temp=0;
- p=head;
- }
- else{
- temp=l[p]+1;
- p=child[p][tt];
- }
- }
- ans[p]=max(temp,ans[p]);
- }
- for(int i=tot;i>0;i--)
- {
- ans[f[rk[i]]]=max(ans[f[rk[i]]],ans[rk[i]]);
- }
- for(int i=1;i<=tot;i++)
- {
- now[i]=min(now[i],ans[i]);
- }
- }
- int ret=0;
- for(int i=1;i<=tot;i++){
- ret=max(now[i],ret);
- }
- printf("%d\\n",ret);
- }
- }SAM;
- int main()
- {
- #ifndef ONLINE_JUDGE
- freopen("in.txt","r",stdin);
- //freopen("out.txt","w",stdout);
- #endif
- scanf("%s",&ch);
- SAM.clear();
- SAM.insert(ch);
- SAM.tpu();
- SAM.solve();
- return 0;
- }