SPOJ-LCS 后缀自动机

题目
求字符串s,p的lcs
将s加入sam中
遍历p,当前字符u
初始结点v=1,长度len=0;
如果该结点有u的出边,匹配,++len
否则通过fa回跳到能匹配的位置。
如果一个也不能匹配,v=1,len=0
板子

#include
#include
#include
#include
using namespace std;
typedef long long ll;
const int N=2e5+5e4+10;
struct SAM{
	int t[N<<1][26],len[N<<1]={-1},fa[N<<1],sz=2,last=1;
	void init(){memset(t,0,sizeof(t));sz=2;last=1;}
	void add(int ch){
		int p=last,np=last=sz++;
		len[np]=len[p]+1;
		for(;p&&!t[p][ch];p=fa[p])t[p][ch]=np;
		if(!p){fa[np]=1;return;}
		int q=t[p][ch];
		if(len[p]+1==len[q])fa[np]=q;
		else{
			int nq=sz++;len[nq]=len[p]+1;
			memcpy(t[nq],t[q],sizeof t[0]);
			fa[nq]=fa[q],fa[np]=fa[q]=nq;
			for(;t[p][ch]==q;p=fa[p])t[p][ch]=nq;
		}
	}
}sf;
char s[N],p[N];
int main(){
	sf.init();
	scanf("%s%s",s,p);
	int n=strlen(s),m=strlen(p);
	for(int i=0;i<n;++i)sf.add(s[i]-'a');
	int ans=0,v=1,len=0;
	for(int i=0;i<m;++i,ans=max(ans,len)){
		int u=p[i]-'a';
		if(sf.t[v][u]){
			v=sf.t[v][u],++len;
			continue;
		}
		while(v&&!sf.t[v][u])v=sf.fa[v];
		if(v)len=sf.len[v]+1,v=sf.t[v][u];
		else v=1,len=0;
	}
	cout<<ans;
	return 0;
}

你可能感兴趣的:(spoj,ACM_字符串)