acdream 1427后缀数组

题意：求一个串中不重叠出现两次以上的子串的数量

解法：后缀数组的分组解法实在是经典 但是今天我明显非主流 然后使用了 后缀数组以后没有进行分组 而是rmq与处理后判断所有后缀的height最多可以有多少成为该串的做法 这当然还需要字符串的hash 但是WA了 然后就没有然后了 

以下是老老实实的分组做法：

#include<map>
#include<math.h>
#include<cstdio>
#include<string.h>
#include<string>
#include<algorithm>
using namespace std;
#define ll unsigned long long
#define mod 131
#define maxn 11111
char txt[maxn];
int sa[maxn],no_1[maxn],no_2[maxn],ex[maxn],he[maxn],rk[maxn],n,limi;
int mi[25][maxn],lg[maxn],len;
void _sa(char *st){
    int i,k,p,*x=no_1,*y=no_2;
    for(i=0; i<limi; i++) ex[i]=0;
    for(i=0; i<n; i++) ex[x[i]=st[i]]++;
    for(i=1; i<limi; i++) ex[i]+=ex[i-1];
    for(i=n-1; i>=0; i--)
        sa[--ex[x[i]]]=i;
    for(k=1,p=1; p<n; k<<=1,limi=p){
        for(p=0,i=n-k; i<n; i++) y[p++]=i;
        for(i=0; i<n; i++) if(sa[i]>=k) y[p++]=sa[i]-k;
        for(i=0; i<limi; i++) ex[i]=0;
        for(i=0; i<n; i++) ex[x[y[i]]]++;
        for(i=1; i<limi; i++) ex[i]+=ex[i-1];
        for(i=n-1; i>=0; i--) sa[--ex[x[y[i]]]]=y[i];
        for(swap(x,y),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++;
    }
}
void _he(char *st){
    int i,j,k=0;
    for(i=0;i<n;i++) rk[sa[i]]=i;
    for(i=0;i<n-1;i++)
    {
        if(k) k--;
        j=sa[rk[i]-1];
        while(st[i+k]==st[j+k]) k++;
        he[rk[i]]=k;
    }
}int q,x,y;
//char now[1111];

void rmq_init(){
    for(int i=0;i<n;i++){mi[0][i]=he[i];}
    for(int i=1;i<=lg[n];i++)
        for(int j=0;j+(1<<i)-1<n;j++){
            mi[i][j]=min(mi[i-1][j],mi[i-1][j+(1<<(i-1))]);
        }
}
int rmq_min(int l,int r){
    if(l>r){
        l--;r++;swap(l,r);
    }
    int tmp=lg[r-l+1];
    return min(mi[tmp][l],mi[tmp][r-(1<<tmp)+1]);
}
int abs(int x){return x>0?x:-x;}
map<ll,int>mp;
int main(){
    lg[0]=-1;for(int i=1;i<maxn;i++)lg[i]=lg[i>>1]+1;
    while(~scanf("%s",txt)){
        mp.clear();
        limi=255,n=(int)strlen(txt);n++;
        _sa(txt);_he(txt);
        rmq_init();
        n--;int ans=0;
        
        //        for(int i=1;i<=n;++i){
        //
        //            printf("%d\n",rmq_min(i,i));
        //
        //        }
        
        //        printf("%d\n",rmq_min(2,3));
        int mi=0,ma=0,fl=0;
        for(int z=1;z<=n;++z){
            fl=0;
            for(int i=1;i<=n;++i){
                //                for(int a=sa[i];a<=n;++a)printf("%c",txt[a]);
                //                printf("\n");
                //                printf("%d %d\n",sa[i],he[i]);
                if(he[i]>=z){
                    if(!fl){
                        fl=1;ma=max(sa[i],sa[i-1]);mi=min(sa[i],sa[i-1]);
                    }else{
                        ma=max(sa[i],ma);mi=min(sa[i],mi);
                    }
                }else{
                    if(!fl){
                        //                        continue;
                    }else{
                        if(ma-mi>=z)ans++;
                        fl=0;
                    }
                }
                //               printf("len:%d i:%d %d %d %d %d\n\n",z,i,ma,mi,ans,fl);
            }
            if(ma-mi>=z&&fl)ans++;
        }
        printf("%d\n",ans);
    }
    return 0;
}