UOJ35 后缀排序 后缀数组

大家都很强, 可与之共勉 。

清嗓子,之前的后缀数组的板子都有问题……

  特此补一个

# include <bits/stdc++.h>

# define N 100010

class Suffix_Array  {
private :
    int n, a [N], b [N], c [N], sa [N], height [N], rank [N] ;
public :
    void build ( char* s, int n, int m )  {
        int *x = a, *y = b, *t ;  this -> n = n ;

        for ( register int i = 0 ; i < m ; c [i ++] = 0 ) ;
        for ( register int i = 0 ; i < n ; ++ c [x [i] = s [i]], ++ i ) ;     
        for ( register int i = 1 ; i < m ; ++ i ) c [i] += c [i - 1] ;
        for ( register int i = n - 1 ; ~ i ; sa [-- c [x [i]]] = i, -- i );

        for ( register int i, k = 1, p = 1 ; p < n && k <= n ; k <<= 1, m = p )  {
            for ( p = 0, i = n - 1 ; i >= n - k ; y [p ++] = i -- ) ;
            for ( i = 0 ; i < n ; sa [i] >= k ? y [p ++] = sa [i ++] - k : ++ i ) ;

            for ( i = 0 ; i < m ; c [i ++] = 0 ) ;
            for ( i = 0 ; i < n ; ++ c [x [y [i ++]]] ) ; 
            for ( i = 1 ; i < m ; ++ i )  c [i] += c [i - 1] ;
            for ( i = n - 1 ; ~ i ; sa [-- c [x [y [i]]]] = y [i], -- i ) ;

            for ( t = x, x = y, y = t, p = 1, x [sa [0]] = 0, i = 1 ; i < n ; ++ i )
                x [sa [i]] = ( y [sa [i - 1]] == y [sa [i]] && y [sa [i - 1] + k] == y [sa [i] + k] ) ? p - 1 : p ++ ;

        }

        for ( register int i = 0 ; i < n ; rank [sa [i]] = i, ++ i ) ;
        for ( register int i = 0, k = 1, j; i < n ; height [rank [i ++]] = k )
            for ( k ? -- k : 0, j = sa [rank [i] - 1] ; s [i + k] == s [j + k] ; ++ k ) ; // s [len - 1] must be '$' ;
    }

    void print ( )  {
        for ( int i = 1 ; i < n ; ++ i )  printf ( "%d%c", sa [i] + 1, ( i == n - 1 ) ? '\n' : ' ' ) ;
        for ( int i = 2 ; i < n ; ++ i )  printf ( "%d%c", height [i], ( i == n - 1 ) ? '\n' : ' ' ) ;
    }
} sa ;

int main ( )  {
    static char s [N] ;
    scanf ( "%s", s ) ;
    int len = strlen ( s ) ;
    s [len] = '$' ;
    sa.build ( s, len + 1, 500 ) ;
    sa.print ( ) ;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值