HDU 5558 后缀数组

思路：

这是一个错误的思路， 因为数据水才过= =

首先求出来后缀数组

把rank插到set里

每回差i两边离i近的rank值，更新

如果LCP相同，暴力左（右）继续更新sa的最小值

//By SiriusRen
#include <bits/stdc++.h>
using namespace std;
const int N=;
char s[N];
int n,T,base[N],f[N][],I=;
int cntA[N],cntB[N],A[N],B[N],sa[N],tsa[N],rk[N],ht[N];
void SA(){
    for(int i=;i<=;i++)cntA[i]=;
    for(int i=;i<=n;i++)cntA[s[i]]++;
    for(int i=;i<=;i++)cntA[i]+=cntA[i-];
    for(int i=n;i;i--)sa[cntA[s[i]]--]=i;
    rk[sa[]]=;
    for(int i=;i<=n;i++)rk[sa[i]]=rk[sa[i-]]+(s[sa[i]]!=s[sa[i-]]);
    for(int l=;rk[sa[n]]<n;l<<=){
        memset(cntA,,sizeof(cntA));
        memset(cntB,,sizeof(cntB));
        for(int i=;i<=n;i++)cntA[A[i]=rk[i]]++,cntB[B[i]=(i+l<=n?rk[i+l]:)]++;
        for(int i=;i<=n;i++)cntA[i]+=cntA[i-],cntB[i]+=cntB[i-];
        for(int i=n;i;i--)tsa[cntB[B[i]]--]=i;
        for(int i=n;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i];
        rk[sa[]]=;
        for(int i=;i<=n;i++)rk[sa[i]]=rk[sa[i-]]+(A[sa[i]]!=A[sa[i-]]||B[sa[i]]!=B[sa[i-]]);
    }
    for(int i=,j=;i<=n;i++){
        j=j?j-:;
        while(s[i+j]==s[sa[rk[i]-]+j])j++;
        ht[rk[i]]=j;
    }base[]=-;
    for(int i=;i<=n;i++)f[i][]=ht[i],base[i]=base[i>>]+;
    for(int j=;j<=;j++)
        for(int i=;i+(<<(j-))<=n;i++)
            f[i][j]=min(f[i][j-],f[i+(<<(j-))][j-]);
}
set<int>S;set<int>::iterator it;
int LCP(int x,int y){
    x++;int t=base[y-x+];
    return min(f[x][t],f[y-(<<t)+][t]);
}
int main(){
    for(scanf("%d",&T);I<=T;I++){
        printf("Case #%d:\n",I);
        scanf("%s",s+),n=strlen(s+);SA();
        S.clear();int t=;
        while(){
            it=S.upper_bound(rk[t]);int tmp=,rec=N;
            while(it!=S.end()){
                if(LCP(rk[t],*it)>=tmp)tmp=LCP(rk[t],*it),rec=min(sa[*it],rec);
                else break;
                it++;
            }it=S.upper_bound(rk[t]);
            while(it!=S.begin()){
                it--;int tt=LCP(*it,rk[t]);
                if(tt>tmp)tmp=tt,rec=sa[*it];
                else if(tt==tmp)rec=min(sa[*it],rec);
                else break;
            }
            if(!tmp){printf("-1 %d\n",s[t]);tmp=;}
            else printf("%d %d\n",tmp,rec-);
            for(int i=;i<=tmp;i++)S.insert(rk[t]),t++;
            if(t>n)break;
        }
    }
}