【hdu3518】Boring counting

题意：找出一个字符串中至少重复出现两次的字串的个数（重复出现时不能重叠）。

后缀数组

枚举字串长度h,对于每一次的h，利用height数组，找出连续的height大于等于h的里面最左端和最右端得为之l和r。
如果l+h-1<r的话，说明没有重叠，答案加1.

 #include<algorithm>
 #include<cstdlib>
 #include<cstring>
 #include<cstdio>
 #include<cmath>
 using namespace std;

 #define N 1010

 int wa[N],wb[N],ws[N],wv[N];
 int sa[N],rank1[N],height[N];
 int r[N];

 char s[N];

 int n,ans;

 int cmp(int *r,int a,int b,int l)
 {
     return r[a]==r[b] && r[a+l]==r[b+l];
 }

 void da(int *r,int n,int m)
 {
     int i,j,p,*x=wa,*y=wb,*t;
     for (i=;i<m;i++) ws[i]=;
     for (i=;i<n;i++) ws[x[i]=r[i]]++;
     for (i=;i<m;i++) ws[i]+=ws[i-];
     for (i=n-;i>=;i--) sa[--ws[x[i]]]=i;
     for (j=,p=;p<n;j<<=,m=p)
     {
         for (p=,i=n-j;i<n;i++) y[p++]=i;
         for (i=;i<n;i++) if (sa[i]>=j) y[p++]=sa[i]-j;
         for (i=;i<n;i++) wv[i]=x[y[i]];
         for (i=;i<m;i++) ws[i]=;
         for (i=;i<n;i++) ws[wv[i]]++;
         for (i=;i<m;i++) ws[i]+=ws[i-];
         for (i=n-;i>=;i--) sa[--ws[wv[i]]]=y[i];
         for (t=x,x=y,y=t,p=,x[sa[]]=,i=;i<n;i++)
             x[sa[i]]=cmp(y,sa[i-],sa[i],j) ? p- : p++;
     }
     return ;
 }

 void calheight(int *r,int n)
 {
     int i,j,k=;
     for (i=;i<=n;i++) rank1[sa[i]]=i;
     for (i=;i<n;height[rank1[i++]]=k)
         for (k ? k-- : ,j=sa[rank1[i]-];r[i+k]==r[j+k];k++);
     return ;
 }

 int main()
 {
     while (scanf("%s",s)!=EOF && s[]!='#')
     {
         ans=;
         n=strlen(s);
         for (int i=;i<n;i++)
             r[i]=s[i]-'a'+;
         r[n]=;
         da(r,n+,);
         calheight(r,n);
         for (int i=;i<=n/;i++)
         {
             int l=n+,r=;
             for (int j=;j<=n;j++)
             {
                 if (height[j]>=i)
                 {
                     r=max(r,sa[j-]);
                     r=max(r,sa[j]);
                     l=min(l,sa[j-]);
                     l=min(l,sa[j]);
                 }
                 else
                 {
                     if (r-l>=i)
                         ans++;
                     r=;
                     l=n+;
                 }
             }
             if (r-l>=i)
                 ans++;
         }
         printf("%d\n",ans);
     }
     return ;
 }