2018.12.15 spoj1812 Longest Common Substring（后缀自动机）

传送门

后缀自动机模板题。

题意简述：求两个字串的最长公共子串长度。

---

对其中一个构建后缀自动机，用另外一个在上面跑即可。

代码：

#include<bits/stdc++.h>
#define ri register int
using namespace std;
const int N=5e5+5;
int n;
char s[N];
struct SAM{
	int tot,last,rt,len[N],son[N][26],link[N];
	SAM(){tot=last=rt=1,len[0]=-1;fill(son[0],son[0]+26,1);}
	inline void expand(int x){
		int p=last,np=++tot;
		last=np,len[np]=len[p]+1;
		for(;p&&!son[p][x];p=link[p])son[p][x]=np;
		if(!p){link[np]=rt;return;}
		int q=son[p][x],nq;
		if(len[p]+1==len[q]){link[np]=q;return;}
		len[nq=++tot]=len[p]+1,memcpy(son[nq],son[q],sizeof(son[q])),link[nq]=link[q];
		for(;p&&son[p][x]==q;p=link[p])son[p][x]=nq;
		link[q]=link[np]=nq;
	}
	inline void query(){
		int p=1,nowlen=0,mxlen=0;
		for(ri i=1;i<=n;++i){
			if(son[p][s[i]-'a']){p=son[p][s[i]-'a'],mxlen=max(mxlen,++nowlen);continue;}
			while(!son[p][s[i]-'a'])p=link[p];
			mxlen=max(mxlen,nowlen=len[p]+1),p=son[p][s[i]-'a'];
		}
		cout<<mxlen;
	}
}sam;
int main(){
	freopen("lx.in","r",stdin);
	scanf("%s",s+1),n=strlen(s+1);
	for(ri i=1;i<=n;++i)sam.expand(s[i]-'a');
	scanf("%s",s+1),n=strlen(s+1),sam.query();
	return 0;
}