[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() ;
}