【POJ 2774】Long Long Message 最长公共子串

最新推荐文章于 2020-01-17 21:38:30 发布

转载最新推荐文章于 2020-01-17 21:38:30 发布 · 140 阅读

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/abclzr/p/5418374.html

本文深入探讨了一种高效的字符串匹配算法，通过使用特定的数据结构和排序技巧来加速字符串的比较过程。该算法首先对输入字符串进行预处理，创建辅助数组并进行排序，然后计算每个位置的最长公共前缀长度，最后通过比较特定条件下的高度值来确定最大匹配长度。文章提供了详细的代码实现，包括初始化数组、排序、计算排名和高度等关键步骤。

还是模板啊，手残&&打成||查错查了1h+TAT

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
const int N = 2000003;

int t1[N], t2[N], c[N];
void st(int *x, int *y, int *sa, int n, int m) {
	int i;
	for(i = 0; i < m; ++i) c[i] = 0;
	for(i = 0; i < n; ++i) ++c[x[y[i]]];
	for(i = 1; i < m; ++i) c[i] += c[i - 1];
	for(i = n - 1; i >= 0; --i) sa[--c[x[y[i]]]] = y[i];
}
void mkhz(int *a, int *sa, int n, int m) {
	int i, j, p, *t, *x = t1, *y = t2;
	for(i = 0; i < n; ++i) x[i] = a[i], y[i] = i;
	st(x, y, sa, n, m);
	for(p = 1, j = 1; p < n; j <<= 1, m = p) {
		for(p = 0, i = n - j; i < n; ++i) y[p++] = i;
		for(i = 0; i < n; ++i) if (sa[i] >= j) y[p++] = sa[i] - j;
		st(x, y, sa, n, m);
		for(t = x, x = y, y = t, p = 1, x[sa[0]] = 0, i = 1; i < n; ++i)
			x[sa[i]] = y[sa[i]] == y[sa[i - 1]] && y[sa[i] + j] == y[sa[i - 1] + j] ? p - 1 : p++;
	}
}
void mkh(int *r, int *sa, int *rank, int *h, int n) {
	int i, j, k = 0;
	for(i = 1; i <= n; ++i) rank[sa[i]] = i;
	for(i = 1; i <= n; h[rank[i++]] = k)
		for(k ? --k : 0, j = sa[rank[i] - 1]; r[i + k] == r[j + k]; ++k);
}

char s[N];
int a[N], sa[N], h[N], rank[N], n;
int main() {
	while(~scanf("%s", s + 1)) {
		int tmp = strlen(s + 1);
		s[tmp + 1] = '$';
		scanf("%s", s + tmp + 2);
		int n = strlen(s + 1);
		for(int i = 1; i <= n; ++i) a[i] = s[i];
		mkhz(a, sa, n + 1, 130);
		mkh(a, sa, rank, h, n);
		int ans = 0;
		++tmp;
		for(int i = 2; i <= n; ++i)
			if (((tmp < sa[i - 1] && tmp > sa[i]) || (tmp < sa[i] && tmp > sa[i - 1])) && h[i] > ans)
				ans = h[i];
		printf("%d\n", ans);
	}
	return 0;
}

Orz

转载于:https://www.cnblogs.com/abclzr/p/5418374.html