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

【wikioi】3160 最长公共子串(后缀自动机)

时间:2014-12-03 23:14:17      阅读:276      评论:0      收藏:0      [点我收藏+]

标签:des   blog   http   io   ar   os   sp   for   on   

http://codevs.cn/problem/3160/

sam的裸题。。。(之前写了spoj上另一题sam的题目,但是spoj被卡评测现在还没评测完QAQ打算写那题题解时再来详细介绍sam的。。。。那就再等等吧。

求两个串的lcs话,就是先建立a串的sam,然后用b的字串去匹配a中。

因为sam中每个状态的len对应最长子串,因此自动机不断trans匹配时,如果没找到下一个点,那么在parent树的祖先中找是否还有子串可以更新(因为祖先的max比这个节点小,且都包含当前状态的right,所以祖先能匹配的串当前状态一定能,所以当前状态不能了,就去找祖先的子串)

所以不断找然后更新即可。

如果还不懂详看clj论文

#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>
#include <set>
#include <map>
using namespace std;
typedef long long ll;
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << (#x) << " = " << (x) << endl
#define error(x) (!(x)?puts("error"):0)
#define rdm(x, i) for(int i=ihead[x]; i; i=e[i].next)
inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<‘0‘||c>‘9‘; c=getchar()) if(c==‘-‘) k=-1; for(; c>=‘0‘&&c<=‘9‘; c=getchar()) r=r*10+c-‘0‘; return k*r; }

struct sam {
	static const int N=200005;
	int c[N][26], l[N], f[N], root, last, cnt;
	sam() { cnt=0; root=last=++cnt; }
	void add(int x) {
		int now=last, a=++cnt; last=a;
		l[a]=l[now]+1;
		for(; now && !c[now][x]; now=f[now]) c[now][x]=a;
		if(!now) { f[a]=root; return; }
		int q=c[now][x];
		if(l[q]==l[now]+1) { f[a]=q; return; }
		int b=++cnt;
		memcpy(c[b], c[q], sizeof c[q]);
		l[b]=l[now]+1;
		f[b]=f[q];
		f[q]=f[a]=b;
		for(; now && c[now][x]==q; now=f[now]) c[now][x]=b;
	}
	void build(char *s) {
		int len=strlen(s);
		rep(i, len) add(s[i]-‘a‘);
	}
	int find(char *s) {
		int ret=0, len=strlen(s), now=1, t=0;
		rep(i, len) {
			int x=s[i]-‘a‘;
			if(c[now][x]) now=c[now][x], ++t;
			else {
				while(now && !c[now][x]) now=f[now];
				if(!now) now=root, t=0;
				else t=l[now]+1, now=c[now][x];
			}
			ret=max(ret, t);
		}
		return ret;
	}
}a;
const int N=100005;
char s[N];
int main() {
	scanf("%s", s);
	a.build(s);
	scanf("%s", s);
	printf("%d\n", a.find(s));
	return 0;
}

  

 


 

 

题目描述 Description

给出两个由小写字母组成的字符串,求它们的最长公共子串的长度。

输入描述 Input Description

读入两个字符串

输出描述 Output Description

输出最长公共子串的长度

样例输入 Sample Input
yeshowmuchiloveyoumydearmotherreallyicannotbelieveit
yeaphowmuchiloveyoumydearmother
样例输出 Sample Output

27

数据范围及提示 Data Size & Hint

单个字符串的长度不超过100000

【wikioi】3160 最长公共子串(后缀自动机)

标签:des   blog   http   io   ar   os   sp   for   on   

原文地址:http://www.cnblogs.com/iwtwiioi/p/4141438.html

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