POJ 1743 (后缀数组)

思路:
首先构建后缀数组,这里要注意theme可以相差k,如abcdefgh中,d-a == e-b因此可以将原数组相邻元素做差,s[i] = s[i] - s[i+1],这样就转化为大小为n-1的数组建立后缀树组。
然后利用二分法搜索可能的长度,判断该长度是不是无重叠重复子串的长度。对于每一个长度d,将height中大于等于d的保留,小于的去除,这样就将height分为几组,因为求最大重复子串就是求某个区间中最小的height,因此最长的无重叠子串一定在同一组内,判断是否重叠就是,利用sa数组看它们在实际数组中的位置间隔是否比d大,如果大则不重叠。
后缀数组是套用的吉林大学的模板(网上大多数的解题报告都是套的国家集训队的模板),真心给这些大牛跪了,感觉比国家集训队论文给的代码效率要高,并且简洁到令人发指的地步,只有望其项背的份了,不过简洁的代码太尼玛难懂了,看了整整一天才明白啊有木有!!

#include <iostream>
using namespace std;

const int N = 20050;
int s[N];                         //store numbers
int n, sa[N], height[N], rank[N], sum[N], tsa[N];

inline int min(int a, int b){return a < b ? a : b;}
inline int max(int a, int b){return a > b ? a : b;}

void makesa(int na){                       // O(N * log N) 
	int i, j, len; 
	//n = strlen(s) + 1;
	//na = (n < 256 ? 256 : n); 
	memset(sum, 0, na * sizeof(int)); 
	for (i = 0; i < n ; i++) sum[ rank[i] = s[i]]++; 
	for (i = 1; i < na; i++) sum[i] += sum[i - 1]; 
	for (i = 0; i < n; i++) sa[ --sum[ rank[i] ] ] = i; //第一关键字排序 
	for (len = 1; len < n; len <<= 1) { 
		//第二关键字排序:首先将sa看做第二关键字,j就是第一关键字,将第一关键字按照第二关键字
		//出现的顺序放入到tsa中,也就完成了双关键字的排序
		for (i = 0; i < n; i++) { 
			j = sa[i] - len; if (j < 0) j += n; 
			tsa[ sum[ rank[j] ]++ ] = j; 
		} 
		sa[ tsa[ sum[0] = 0 ] ] = j = 0; //j代表排名
		for (i = 1; i < n; i++) {
			//如果当前项的两个关键字中,第一个和第二个都和排在前面项的相同,则代表
			//该项和前一项排名相同,j不变
			if (rank[ tsa[i] ] != rank[ tsa[i-1] ] || 
				rank[ tsa[i]+len ]!=rank[ tsa[i-1]+len ]) 
				sum[++j] = i; //这里sum记录了相同项中,第一次出现的项的排名,也就是相同项中排名最高者 
			sa[ tsa[i] ] = j; //这里将rank值放入sa,sa这里只是作为一个缓冲
		} 
		memcpy(rank, sa , n * sizeof(int)); 
		memcpy(sa  , tsa, n * sizeof(int)); 
		if (j >= n - 1) break; 
	} 
}

void lcp(){                          // O(4 * N) 
	int i, j, k; 
	j = rank[0];
	height[0] = 0;
	for (i = k = 0; i < n - 1; i++, k++) 
		while (k >= 0 && s[i] != s[ sa[j-1] + k ]) 
			height[j] = (k--), j = rank[ sa[j] + 1 ]; 
}
/*
利用二分法搜索可能的长度,然后判断是不是无重叠重复子串的长度。
对于每一个长度d,将height中大于等于d的保留,小于的去除,这样就将height非为几组,
因为求最大重复子串就是求某个区间中最小的height,因此最长的无重叠子串一定在同一组内,
判断是否重叠就是,利用sa数组看它们在实际数组中的位置间隔是否比d大,如果da则不重叠。
*/
bool check(int mid)
{
	int low = sa[1], high = sa[1];
	for (int i = 2; i < n; ++i){
		if (height[i] >= mid){
			low = min(low, sa[i]);
			high = max(high, sa[i]);
			if (high - low >= mid) return true;
		}
		else{
			low = sa[i];
			high = sa[i];
		}
	}
	return false;
}

void solve()
{
	int l = 1, r = n, mid, i, ans = 0;
	while (l < r){
		mid = (l + r) / 2;
		if (check(mid)) 
			ans = mid, l = mid + 1;
		else r = mid;
	}
	if (ans >= 4) printf("%d\n", ans+1);
	else printf("0\n");
}

int main()
{
	int i;
	while (scanf("%d", &n) && n){
		for (i = 0; i < n; ++i)
			scanf("%d", &s[i]);
		for (i = 0; i < n-1; ++i)
			s[i] = s[i] - s[i+1] + 100;
		s[i] = 0;
		makesa(200);
		lcp();
		solve();
	}
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值