acdream 1427后缀数组

来源：互联网发布：网络利大于弊辩论赛编辑：程序博客网时间：2024/05/22 15:50

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

解法：后缀数组的分组解法实在是经典但是今天我明显非主流然后使用了后缀数组以后没有进行分组而是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 11111char 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;}

0 0