POJ 3261 Milk Patterns

原创 2012年03月24日 08:56:43

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

题意:N头牛,每头牛有一个编号,求可重叠的至少出现K次的最大子串。

思路: 后缀数组+二分 。 二分答案,然后将后缀分成若干组。不同的是,这里要判断的是有没有一个组的后缀个数不小于k。如果有,那么存在k 个相同的子串满足条件,否则不存在。这个做法的时间复杂度为O(nlogn)。

代码: 

#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;
#define MAXN 20010
#define MAXN1 1000010
int num[MAXN] ;
int sa[MAXN] , rank[MAXN] ,height[MAXN] ;  
int wa[MAXN1] , wb[MAXN1] ,wv[MAXN1],wd[MAXN1] ;   
int N,M;  
int cmp(int *r , int a , int b , int l){  
    return r[a] == r[b] && r[a+l] == r[b+l] ;     
}  

void DA(int *r,int n,int m){        //O(NlogN)  
    int i, j , p , *x=wa, *y=wb,*t ;  
    for( i = 0 ; i < m ; i++ )   wd[i] = 0 ;  
    for( i = 0 ; i < n ; i++ )       wd[x[i]=r[i]] ++ ;  
    for( i = 1 ; i < m ; i++ )       wd[i] += wd[i-1] ;  
    for( i = n-1 ;i >= 0 ; i-- )     sa[--wd[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++)    wv[i] = x[y[i]] ;  
        for( i = 0 ; i < m ; i++)    wd[i] = 0 ;  
        for( i = 0 ; i < n ; i++)    wd[wv[i]] ++ ;  
        for( i = 1 ; i < m ; i++)    wd[i] += wd[i-1] ;  
        for( i = n-1 ; i >= 0 ; i--) sa[ --wd[wv[i]]] = y[i] ;  
        for( t = x , x = y , y = t , p = 1 , x[ sa[0] ] = 0,i = 1;i < n ; i++){  
            x[sa[i]] = cmp( y ,sa[i-1] ,sa[i] , j ) ? p-1: p++;   
        }     
    }  
}   
  
void calHeight(int *r , int n){  
    int i , j , k = 0 ;  
    for( i = 1 ; i <= n ; i++)   rank[sa[i]] = i ;  
    for( i = 0 ; i < n ; height[ rank[i++]]=k){  
        for( k ? k-- : 0 , j=sa[rank[i]-1]; r[i+k]==r[j+k] ; k++) ;   
    }     
}
bool is_ok(int len){
	int i = 1 ;
	while(i <= N){
		while(i <= N && height[i]<len)	i++ ;
		int c = 1 ;
		while(i <= N && height[i]>=len){		//以height[]数组分组
			i++ ;
			c ++ ;	
		}
		if(c >= M)	return true ;
	}	
	return false ;
}

int main(){
	while(scanf("%d %d",&N,&M) == 2){
		for(int i=0;i<N;i++){
			scanf("%d",&num[i]);	
			num[i]++ ;
		}		
		num[N] = 0 ;
		DA(num,N+1,MAXN1);
		calHeight(num, N) ;
		int low, high , mid ;
		low = 1 ; high = N ;
		while(low < high){
			mid = (low + high + 1) >> 1;
			if(is_ok(mid))	
				low = mid ;
			else 
				high = mid - 1;
		}
		printf("%d\n",low);
	}
	return 0 ;	
}


相关文章推荐

【POJ3261】Milk Patterns 后缀数组

水题不好意思说题解。 说说题意吧: 给一个字符串(数字串),然后求最长k次重复子串。 即某串在字符串中重复了至少k次,求这种串的最长长度。 代码: #include #in...
  • Vmurder
  • Vmurder
  • 2015年01月04日 22:19
  • 742

【POJ】3261 Milk Patterns 【后缀数组】

传送门:【POJ】3261 Milk Patterns 题目分析:构好后缀数组,然后二分答案就好了~_~ 代码如下: #include #include #include using...

poj-3261-Milk Patterns

后缀数组。。 二分答案。。 #include #include #include #include #include using namespace std; #define...

POJ 3261 Milk Patterns (离散化+后缀数组 可重叠k次最长重复子串)

题意:给出一串长度为n的字符,再给出一个k值,要你求重复次数大于等于k次的最长子串长度 思路:首先离散化,然后采用后缀数组的常用处理方法,二份答案,按照二份值k将height数组分组,对于k是否可行...

poj 3261 Milk Patterns 【后缀数组】

给定一个字符串,求至少出现 k 次的最长重复子串,这 k 个子串可以重叠。 算法分析: 这题的做法和上一题差不多,也是先二分答案,然后将后缀分成若干组。不同的是,这里要判断的是有没有一个组的后缀个...

POJ 3261 Milk Patterns

刚开始看到最大值要取到1000000,当时就吓尿了,心想难道要改排序算法? 赶紧看了下讨论,发现数据很水,200000就能过了 思路:就是答案二分,然后长度分组,满足条件返回true #include...
  • lj94093
  • lj94093
  • 2015年04月14日 11:03
  • 380

poj 3261 Milk Patterns 后缀数组+二分

题意是求可以重复的出现次数超过k次的子串,最长是多少。 样例输入 8 2 1 2 3 2 3 2 3 1 输出 4 二分验证是否有连续的k-1个height大于所验证的答案...
  • Loi_a
  • Loi_a
  • 2016年12月27日 08:03
  • 154

POJ 3261 Milk Patterns 最长出现k次的子串长度(后缀数组)

Milk PatternsTime Limit: 5000MS Memory Limit: 65536K DescriptionFarmer John has noticed that t...
  • w4149
  • w4149
  • 2017年07月27日 19:36
  • 73

poj 3261 Milk Patterns (后缀数组)

poj 3261 Milk Patterns 题意: 给定一个字符串S,求至少出现k次的最长重复子串,这k个子串可以重叠。 限制: 1 思路: 二分答案长度用lcp判断 ...
  • whai362
  • whai362
  • 2015年08月14日 14:20
  • 291

POJ-3261-Milk Patterns

利用后缀数组的height数组做,具体可以参考罗大神的论文~ 代码: #include #include #include using namespace std; const int max...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:POJ 3261 Milk Patterns
举报原因:
原因补充:

(最多只允许输入30个字)