码迷,mamicode.com
首页 > 其他好文 > 详细

●SPOJ LCS2Longest Common Substring II

时间:2018-03-10 22:00:50      阅读:204      评论:0      收藏:0      [点我收藏+]

标签:span   --   else   body   new   16px   spoj   substr   匹配   

题链:

http://www.spoj.com/problems/LCS2/
题解:

后缀自动机。
对第一个串建立后缀自动机,
然后把后面的每个串分别与该串的自动机去匹配,求出相应的数组val*[s]:
表示第*个串与第一个串的自动机的状态s的最大匹配长度,

(求法就是两个串用后缀自动机求LCS这一过程):
对于当前已经匹配的子串T,长度为now,此刻在状态s,现在要匹配第i个字符x,
若trans(s,x)!=0,则s=trans(s,x),now++,i++,并更新val*[s]=max(val*[s],now);
否则,s=parent[s],直到trans(s,x)!=0或者s=0。
由于我们在每一步更新s时,没有更新到其祖先,所以最后再桶排之后去依次更新父亲。

而最后的状态s和每个串的最大匹配长度就是min(min(val*[s]),maxs[s](这个是状态s所允许的最大长度));


代码:

 

#include<bits/stdc++.h>
#define MAXN 100005
#define INF 0x3f3f3f3f
using namespace std;
struct SAM{
	int size,last,lens;
	int maxs[MAXN*3],trans[MAXN*3][26],parent[MAXN*3],minmatch[MAXN*3];
	int Newnode(int a,int b){
		++size; maxs[size]=a; minmatch[size]=INF;
		memcpy(trans[size],trans[b],sizeof(trans[b]));
		return size;
	}
	void Extend(int x){
		static int p,np,q,nq;
		p=last; last=np=Newnode(maxs[p]+1,0);
		for(;p&&!trans[p][x];p=parent[p]) trans[p][x]=np;
		if(!p) parent[np]=1;
		else{
			q=trans[p][x];
			if(maxs[p]+1!=maxs[q]){
				nq=Newnode(maxs[p]+1,q);
				parent[nq]=parent[q];
				parent[q]=parent[np]=nq;
				for(;p&&trans[p][x]==q;p=parent[p]) trans[p][x]=nq;
			}
			else parent[np]=q;
		}
	}
	void Build(char *S){
		memset(trans[0],0,sizeof(trans[0]));
		size=0; last=Newnode(0,0); lens=strlen(S);
		for(int i=0;i<lens;i++) Extend(S[i]-‘a‘);
	}
	void Update(int *val){
		static int tmp[MAXN],order[MAXN*3];
		memset(tmp,0,sizeof(tmp));
		for(int p=1;p<=size;p++) tmp[maxs[p]]++;
		for(int i=1;i<=lens;i++) tmp[i]+=tmp[i-1];
		for(int p=1;p<=size;p++) order[tmp[maxs[p]]--]=p;
		for(int i=size,p;i;i--)
			p=order[i],val[parent[p]]=max(val[parent[p]],val[p]);
	}
	void Match(char *T){
		static int val[MAXN*3],p,i,len,now;
		for(p=1;p<=size;p++) val[p]=0;
		p=1; i=0; now=0; len=strlen(T);
		while(i<len){
			if(!p) p=1,now=0,i++;
			else if(!trans[p][T[i]-‘a‘]) p=parent[p],now=maxs[p];
			else now++,p=trans[p][T[i]-‘a‘],i++;
			val[p]=max(val[p],now);
		}
		Update(val);
		for(p=1;p<=size;p++) minmatch[p]=min(minmatch[p],val[p]);
	}
}SUF;
int main(){
	static char S[MAXN];
	scanf("%s",S); SUF.Build(S);
	while(~scanf("%s",S)) SUF.Match(S);
	int ans=0;
	for(int p=1;p<=SUF.size;p++)
		ans=max(ans,min(SUF.minmatch[p],SUF.maxs[p]));
	printf("%d\n",ans);
	return 0;
}

 

  

 

●SPOJ LCS2Longest Common Substring II

标签:span   --   else   body   new   16px   spoj   substr   匹配   

原文地址:https://www.cnblogs.com/zj75211/p/8541876.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!