1. 程式人生 > >●SPOJ 1811 Longest Common Substring

●SPOJ 1811 Longest Common Substring

scan spa ostream class first longest open main rule

題鏈:

http://poj.org/problem?id=2774

題解:

求兩個字符串(S,T)的最長公共子串。
對 S串建後綴自動機。
接下來就用這個自動機去求出能和 S串匹配的 T的每一個前綴的最長的後綴。
最終答案就是對每個 T的前綴得到的答案取最大值就好了。

代碼:

#include<cstdio>
#include<cstring>
#include<iostream>
#define MAXN 250050
#define filein(x) freopen(#x".in","r",stdin);
#define fileout(x) freopen(#x".out","w",stdout);
using namespace std;
struct SAM{
	int size,last,q,p,nq,np;
	int pre[MAXN*2],step[MAXN*2],ch[MAXN*2][26];
	int Newnode(int a,int b){
		step[size]=a; memcpy(ch[size],ch[b],sizeof(ch[b]));
		return size++;
	}
	void Extend(int x){
		p=last; last=np=Newnode(step[p]+1,0);
		while(p&&!ch[p][x]) ch[p][x]=np,p=pre[p];
		if(!p) pre[np]=1;
		else{
			q=ch[p][x];
			if(step[q]!=step[p]+1){
				nq=Newnode(step[p]+1,q);
				pre[nq]=pre[q]; pre[q]=pre[np]=nq;
				while(p&&ch[p][x]==q) ch[p][x]=nq,p=pre[p];
			}
			else pre[np]=q;
		}
	}
	void Build(char *S){
		memset(ch[0],0,sizeof(ch[0]));
		size=1; last=Newnode(0,0); pre[last]=0;
		for(int i=0;S[i];i++) Extend(S[i]-‘a‘);
	}
	int Match(char *T){
		p=1; int ans=0;
		for(int j=0,x,match=0;T[j];j++){ //與 T的每個前綴匹配 
			x=T[j]-‘a‘;
			if(ch[p][x]) match++,p=ch[p][x];
			else{
				while(p&&!ch[p][x]) p=pre[p];
				if(!p) match=0,p=1;
				else match=step[p]+1,p=ch[p][x];
			}
			ans=max(ans,match);
		}
		return ans;
	}
}suf;
char S[MAXN],T[MAXN];
int main()
{
	scanf(" %s %s",S,T);
	suf.Build(S);
	printf("%d",suf.Match(T));
	return 0;
}

●SPOJ 1811 Longest Common Substring