为啥没人写hash或者sa?不明觉厉一脸懵逼…………
于是我也跟风写一发SAM吧…>_<…
把第一个串建SAM,后面的每个放到SAM里跑,匹配到对应的节点,就与当前长度取max(len[p]=max(len[p],tmp)),表示这个串在这个节点能匹配的最大长度,然后每个串取min,意味着取出来的是公共字串(ans[i]=min(ans[i],len[i])),最后所有ans取max,表示所有公共字串里最长的那个。注意每个串搞定之后要toposort一遍…>_<…
这个取max取min取max我也是无语了…………
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
//by:MirrorGray
using namespace std;
const int N=4111;
char s[N];
int n,cnt,root,last;
int gs[N],q[N],len[N],ans[N];
struct SAM{
int nxt[26],fa,mx;
}tr[N];
int nwn(int l){
tr[++cnt].mx=l;
return cnt;
}
int extend(int x){
int np=nwn(tr[last].mx+1);
int p=last;last=np;
for(;p && !tr[p].nxt[x];p=tr[p].fa)tr[p].nxt[x]=np;
if(!p)return tr[np].fa=root,233;
int q=tr[p].nxt[x];
if(tr[q].mx==tr[p].mx+1)return tr[np].fa=q,233;
int nq=nwn(tr[p].mx+1);
tr[nq].fa=tr[q].fa;tr[q].fa=tr[np].fa=nq;
memcpy(tr[nq].nxt,tr[q].nxt,sizeof(tr[q].nxt));
for(;p && tr[p].nxt[x]==q;p=tr[p].fa)tr[p].nxt[x]=nq;
}
void toposort(){
for(int i=1;i<=n;i++)gs[i]=0;
for(int i=1;i<=cnt;i++)gs[ tr[i].mx ]++;
for(int i=1;i<=n;i++)gs[i]+=gs[i-1];
for(int i=cnt;i;i--)q[ gs[tr[i].mx]-- ]=i;
}
int main(){
int m;scanf("%d",&m);m--;
scanf("%s",s+1);n=strlen(s+1);
root=last=nwn(0);
for(int i=1;s[i];i++)extend(s[i]-'a');
toposort();
for(int i=1;i<=cnt;i++)ans[i]=tr[i].mx;
while(m--){
memset(len,0,sizeof(len));
scanf("%s",s+1);
int p=root,tmp=0;
for(int i=1;s[i];i++){
int c=s[i]-'a';
while(p!=root && !tr[p].nxt[c])p=tr[p].fa,tmp=tr[p].mx;
if(tr[p].nxt[c])tmp++,p=tr[p].nxt[c];
len[p]=max(len[p],tmp);
}
for(int i=cnt;i;i--)len[ tr[q[i]].fa ]=max(len[ tr[q[i]].fa ],len[ q[i] ]);
for(int i=1;i<=cnt;i++)ans[i]=min(ans[i],len[i]);
}
int pr=0;
for(int i=1;i<=cnt;i++)pr=max(pr,ans[i]);
printf("%d\n",pr);
return 0;
}