2018.12.15 spoj Longest Common Substring II(后缀自动机)

传送门

后缀自动机基础题。

给出10个串求最长公共子串。


我们对其中一个建一个samsamsam,然后用剩下九个去更新范围即可。

代码:

#include<bits/stdc++.h>
#define ri register int
using namespace std;
const int N=2e5+5;
int T=0,n;
char s[N];
struct SAM{
	int rt,tot,last,len[N],lz[N],cnt[N],rk[N],val[N],link[N],son[N][26],mn[N],mx[N];
	SAM(){rt=tot=last=1,len[0]=-1,fill(son[0],son[0]+26,1);}
	inline void expend(int x){
		int p=last,np=++tot;
		last=np,len[np]=len[p]+1;
		while(p&&!son[p][x])son[p][x]=np,p=link[p];
		if(!p){link[np]=rt;return;}
		int q=son[p][x],nq;
		if(len[q]==len[p]+1){link[np]=q;return;}
		len[nq=++tot]=len[p]+1,memcpy(son[nq],son[q],sizeof(son[q])),link[nq]=link[q];
		while(p&&son[p][x]==q)son[p][x]=nq,p=link[p];
		link[np]=link[q]=nq;
	}
	inline void topsort(){
		for(ri i=1;i<=tot;++i)mn[i]=len[i];
		for(ri i=1;i<=tot;++i)++cnt[len[i]];
		for(ri i=1;i<=last;++i)cnt[i]+=cnt[i-1];
		for(ri i=1;i<=tot;++i)rk[cnt[len[i]]--]=i;
	}
	inline void update(){
		int p=1,nowlen=0;
		for(ri i=1;i<=tot;++i)mx[i]=0;
		for(ri i=1,x;i<=n;++i){
			x=s[i]-'a';
			if(son[p][x])p=son[p][x],++nowlen;
			else{
				while(!son[p][x])p=link[p];
				nowlen=len[p]+1,p=son[p][x];
			}
			mx[p]=max(mx[p],nowlen);
		}
		for(ri p,i=tot;i;--i)p=rk[i],mn[p]=min(mn[p],mx[p]),mx[link[p]]=min(len[link[p]],max(mx[link[p]],mx[p]));
	}
	inline void query(){
		int ans=0;
		for(ri i=1;i<=tot;++i)ans=max(ans,mn[i]);
		if(ans<2)ans=0;
		cout<<ans;
	}
}sam;
int main(){
	while(~scanf("%s",s+1)){
		++T,n=strlen(s+1);
		if(T==1){for(ri i=1;i<=n;++i)sam.expend(s[i]-'a');sam.topsort();}
		else sam.update();
	}
	sam.query();
	return 0;
}
上一篇:Tencent Cloud Developers Conference(2018.12.15)


下一篇:IntelliJ Idea 常用快捷键列表(转)