poj 3261 Milk Patterns 后缀数组最长重复子串

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

给一串数组，数组最少含有k个相同子串，可重叠，求这样子串的最长长度。

后缀数组求出 height[]，若连续k个height[]都大于mid，就可以了。当然要找最大的mid，这里可用二分查找.。

关于后缀数组可以参考：http://hi.baidu.com/fhnstephen/blog/item/4b20757c37245d0429388a76.html

#include <iostream>

#include<stdio.h>

#include<stdlib.h>

using namespace std;

const int maxn = 1000005;

const int maxm = 1000001;

int wa[maxn], wb[maxn];

int r[maxn], b[1000001];

int rank[maxn];

int num[maxn], sa[maxn];

int height[200001];

int cmp (int *r, int a, int b, int l) {

return r[a] == r[b] && r[a+l] == r[b+l];

}

void make_sa(int n, int m){

int i, j, p, *x = wa, *y = wb;

for (i = 0; i < m; i++) b[i] = 0;

for (i = 0; i < n; i++) b[x[i] = num[i]]++;

for (i = 1; i < m; i++) b[i] += b[i-1];

for (i = n-1; i >= 0; i--) sa[--b[x[i]]] = i;

for (j = 1, p=1; p < n; j*=2, 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;

for (i = 0; i < n; i++) r[i] = x[y[i]];

for (i = 0; i < m; i++) b[i] = 0;

for (i = 0; i < n; i++) b[r[i]]++;

for (i = 1; i < m; i++) b[i] += b[i-1];

for (i = n-1; i>= 0; i--) sa[--b[r[i]]] = y[i];

int *t = x; x = y; y = t;

for (i = p = 1, x[sa[0]] = 0; i < n; i++)

x[sa[i]] = cmp(y, sa[i-1], sa[i], j) ? p-1:p++;

}

return;

}

void calheight(int len)

{

int i,j,k=0;

for (i =1; i<len+1; i++) rank[sa[i]] = i;

for(i=0;i<len;height[rank[i++]]=k)

for(k?k--:0,j=sa[rank[i]-1];num[i+k]==num[j+k];k++);

return ;

}

int k;

bool judge(int key,int len)

{

int ssum=0;

for(int i=1;i<len;i++)

{

if(height[i]>=key)

{

ssum++;

if(ssum+1>=k)

{

return true;

}

else{

ssum=0;

}

return false;

}

int main()

{

int n,i;

while(scanf("%d%d",&n,&k)!=EOF)

{

for(i=0;i<n;i++)

{

scanf("%d",&num[i]);

}

num[n]=0;

make_sa(n+1,maxm);

calheight(n);

int low=1,hight=n;

while(low<=hight)

{

int mid=(low+hight)>>1;

if(judge(mid,n+1))

{

low=mid+1;

}

else{

hight=mid-1;

}

if(hight<n)

{

printf("%d\n",hight);

}

return 0;

}

转载于:https://www.cnblogs.com/zxj015/archive/2011/03/19/2740273.html

poj 3261 Milk Patterns 后缀数组最长重复子串相关推荐

试题算法训练后缀数组——最长重复子串
资源限制时间限制:100ms 内存限制:256.0MB 问题描述给定一个长度为n的数串,求至少出现k 次的最长重复子串的长度,这k 个子串可以重叠.保证有子串出现至少k次. 输入格式第一行:两个 ...
POJ - 3261 Milk Patterns(二分+后缀数组)
题目链接:点击查看题目大意:给出一个字符串,以及一个k,现在求出现次数大于等于k次的最大可重叠子串的长度题目分析:可以说是后缀数组的模板题目了吧..直接跑出height数组,因为height数组代 ...
HDU1403(后缀数组--最长公共子串)
题目:Longest Common Substring 看代码注释请戳这里题意:判断给定的两个串中,最长的公共串. 思路:将它们合并为一个串,然后利用后缀数组求解. 首先是二倍增算法:时间复杂度为O ...
后缀数组求最长重复子串
问题描述给定一个字符串,求出其最长重复子串例如:abcdabcd 最长重复子串是 abcd,最长重复子串可以重叠例如:abcdabcda,这时最长重复子串是 abcda,中间的 a 是被重叠的. ...
后缀数组--(可重叠最长重复子串问题)
问题描述:给定一个字符串,求最长重复子串,这两个子串可以重叠. 其实问题可以转化为height数组的最大值.至于为什么是这样,我可以这样解释: 求可重叠最长重复子串等价于求两个后缀的最长公共前缀的最大 ...
#1407 : 后缀数组二·重复旋律2 (不可重叠最长重复子串问题)
题目链接思路求不可重叠最长重复子串,也可以利用height[i]height[i]height[i]. 二分枚举答案KKK 当height[i]≥Kheight[i] \ge Kheight[i] ...
最长公共子序列|最长公共子串|最长重复子串|最长不重复子串|最长回文子串|最长递增子序列|最大子数组和...
最长公共子序列|最长公共子串|最长重复子串|最长不重复子串|最长回文子串|最长递增子序列|最大子数组和文章作者:Yx.Ac 文章来源:勇幸|Thinking (http://www.ahathi ...
最长重复子串和最长不重复子串求解
最长重复子串和最长不重复子串求解本文内容框架: §1 最长重复子串基本方法.KMP算法求解.后缀数组求解 §2 最长不重复子串基本方法.动态规划.动态规划+Hash §3 小结 §1最长重复子串 ...
POJ1743——不可重迭的最长重复子串
题意:给定一系列的整数作为音阶,旋律为相邻音阶之差.问最长的主旋律是多长,主旋律需满足3个条件: 1. 长度至少为5: 2. 至少重复出现2次: 3. 主旋律各不重迭. 对源进行变换以后,就是问不重迭 ...

poj 3261 Milk Patterns 后缀数组最长重复子串

poj 3261 Milk Patterns 后缀数组最长重复子串相关推荐

最新文章

热门文章

poj 3261 Milk Patterns 后缀数组 最长重复子串

poj 3261 Milk Patterns 后缀数组 最长重复子串相关推荐

最新文章

热门文章

poj 3261 Milk Patterns 后缀数组最长重复子串

poj 3261 Milk Patterns 后缀数组最长重复子串相关推荐