[BZOJ1717]-[Usaco2006 Dec]Milk Patterns-后缀数组Height模板题

说在前面

回血第二题=w=
之后就看着某学习笔记的题表写题了hhhhh
（话说光之子的桌面真的很好看呢!!!）

---

题目

BZOJ1717传送门

题面

农夫John发现他的奶牛产奶的质量一直在变动。经过细致的调查，他发现：虽然他不能预见明天产奶的质量，但连续的若干天的质量有很多重叠。我们称之为一个“模式”。 John的牛奶按质量可以被赋予一个0到1000000之间的数。并且John记录了N(1<=N<=20000)天的牛奶质量值。他想知道最长的出现了至少K(2<=K<=N)次的模式的长度。比如1 2 3 2 3 2 3 1 中 2 3 2 3出现了两次。当K=2时，这个长度为4。

输入输出格式

输入格式：
第一行包含两个整数 N,K，含义如题。
接下来N行，每行一个整数表示当天的质量值。

输出格式：
输出一行一个整数，表示N天中最长的出现了至少K次的模式的长度

---

解法

知道Height大概是怎么用的，然后去看一点点博客，这就是一道模板题了=w=

就是先求出Height，然后二分长度，如果有连续K个Height[i]大于等于二分的mid，就check成功，大概就是这样吧= =

---

下面是自带大常数的代码

#include <cstdio>#include <cstring>#include <algorithm>using namespace std ;int N , K , a[20005] , uninum ;struct Unique_Data{    int num , id ;    bool operator < ( const Unique_Data &A ) const {        return num < A.num ;    }}U[20005] ;void Unique_(){    sort( U + 1 , U + N + 1 ) ;    for( int i = 1 , last = -1 ; i <= N ; i ++ ){        if( last != U[i].num ){            uninum ++ ;            last = U[i].num ;        }        a[ U[i].id ] = uninum ;    }}int RA[20005] , SA[20005] , He[20005] ;int tax[20005] , tp[20005] ;bool equ( int *f , int x , int y , int len ){    int xl = ( x + len > N ? 0 : f[x+len] ) ,        yl = ( y + len > N ? 0 : f[y+len] ) ;    return f[x] == f[y] && xl == yl ;}void Rsort( int m ){    for( int i = 1 ; i <= m ; i ++ ) tax[i] = 0 ;    for( int i = 1 ; i <= N ; i ++ ) tax[ RA[tp[i]] ] ++ ;    for( int i = 1 ; i <= m ; i ++ ) tax[i] += tax[i-1] ;    for( int i = N ; i ; i -- ) SA[ tax[RA[tp[i]]]-- ] = tp[i] ;}void getSA(){    for( int i = 1 ; i <= N ; i ++ )        RA[i] = a[i] , tp[i] = i ;    int m = uninum ; Rsort( m ) ;    for( int w = 1 , p = 0 ; p != N ; w += w , m = p ){        p = 0 ;        for( int i = N - w + 1 ; i <= N ; i ++ ) tp[++p] = i ;        for( int i = 1 ; i <= N ; i ++ )            if( SA[i] > w ) tp[++p] = SA[i] - w ;        Rsort( m ) ;        memcpy( tp , RA , sizeof( RA ) ) ;        RA[ SA[1] ] = p = 1 ;        for( int i = 2 ; i <= N ; i ++ )            RA[ SA[i] ] = ( equ( tp , SA[i] , SA[i-1] , w ) ? p : ++p ) ;    }    for( int i = 1 , k = 0 , j ; i <= N ; He[ RA[i] ] = k , i ++ )        for( k = k ? k-1 : k , j = SA[RA[i]-1] ; a[i+k] == a[j+k] ; k ++ ) ;}bool check( int mid ){    for( int i = 1 , len = 0 ; i <= N ; i ++ ){        if( He[i] < mid ) len = 0 ;        else if( (++len) == K - 1 ) return true ;    }    return false ;}void solve(){    getSA() ;    /*for( int i = 1 ; i <= N ; i ++ )        printf( "%d " , a[i] ) ;    puts( "" ) ;    for( int i = 1 ; i <= N ; i ++ )        printf( "%d " , RA[i] ) ; */    int lf = 1 , rg = N , ans = 0 ;    while( lf <= rg ){        int mid = ( lf + rg ) >> 1 ;        if( check( mid ) ){            ans = mid ;            lf = mid + 1 ;        } else rg = mid - 1 ;    }    printf( "%d" , ans ) ;}int main(){    scanf( "%d%d" , &N , &K ) ;    for( int i = 1 ; i <= N ; i ++ ){        scanf( "%d" , &a[i] ) ;        U[i].num = a[i] , U[i].id = i ;    }    Unique_() ;    solve() ;}