uva12206 后缀数组

这题说的是给了一串字符 我们要将这个字符 中找出至少出现m次的最长字符串 一个字符课多次使用

利用后缀数组计算最长的lcp

这里有一个点 记得将后缀数组中加入一个空串 如果遇到全部相同的字符时 没办法 判断 倒数第二个和 第三个的大小 因此他们就会被遗漏

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <vector>
#include <string.h>
using namespace std;
const int maxn =;
struct SuffixArray
{
     int s[maxn];
     int sa[maxn];
     int t[maxn],t2[maxn],c[maxn];
     int rank[maxn],lcp[maxn];
     int n;
     void build(int m)
     {

         int *x=t,*y=t2,i;
         for(i=; i<m; i++ )c[i]=;
         for(i=; i<n; i++)c[ x[i]=s[i]  ]++;
         for(i=; i<m; i++)c[i]+=c[i-];
         for(i=n-; i>=; i--)sa[ --c[ x[i]  ]   ]=i;
         for(int k=; k<=n; k<<=)
            {
                int p=;
                 for(i=n-k; i<n; i++) y[ p++ ]=i;
                 for(i=; i<n; i++)if(sa[i]>=k) y[p++ ]= sa[i]-k;
                 for(i=; i<m; i++) c[ i ]=;
                 for(i=; i<n; i++) c[ x[ y[i] ]  ]++;
                 for(i=; i<m; i++)c[i]+=c[i-];
                 for(i=n-; i>=; i--)
                    sa[ -- c[ x[ y[ i ] ] ]   ]=y[i];
                 for(i=; i<n; i++)rank[ sa[i] ]=i;
                 p=;
                 swap(x,y);
                  x[ sa[] ] =;
                 for(int i=; i<n; i++)
                    x[ sa[i]  ]= y[ sa[i] ]==y[ sa[i-] ]&&y[ sa[i]+k ]==y[ sa[i-]+k  ] ?p-:p++;
                 if(p>=n)break;
                 m=p;
            }
     }
     void getLcp()
     {
          int i,k=;
          for(i=; i<n; i++)rank[ sa[ i ] ] =i;
          for(i=; i<n; i++)
            {
                if(k)k--;
                int j=sa[ rank[i]-  ];
                while(s[j+k]==s[i+k])k++;
                lcp[ rank[i]  ]=k;
            }
     }
}sa;
bool solve(int len,int m)
{
     int L=;
     for(int R=; R<=sa.n; R++)
     {
         if(sa.lcp[R]<len||R==sa.n)
            {
                if( R-L >= m )return true;
                L=R;
            }

     }
     return false;
}
int look(int L, int R)
{
      int ans=-;
      for(int i=L; i<R; i++)
        {
             ans=max(sa.sa[i],ans);
        }
      return ans;
}

int select(int len, int m)
{
     int L=,ad=-;

     for(int R=; R<=sa.n; R++)
     {
         if(sa.lcp[R]<len||R==sa.n)
            {
                if(R-L>=m)
                    {
                        ad=max(look(L,R),ad);
                    }

                   L=R;
            }

     }
     return ad;
}
char ss[maxn];
int main()
{
    int m;
    while(scanf("%d",&m)==)
    {

         if(m==)break;
        scanf("%s",ss);

        sa.n=strlen(ss);
        for(int i=; i<sa.n; i++)
            {
                 sa.s[i]= ss[i]-'a'+;
            }
        if(m==)
            {
               printf("%d %d\n",sa.n,); continue;
            }
        sa.s[sa.n]=;
        sa.n++;
        sa.build();
        sa.getLcp();

        if(solve(,m)==false)
            {
                puts("none");continue;
            }
        int L=,R=sa.n,ans;
        while(L<=R)
            {
               int mid=(R+L)/;
                if(solve(mid,m)){
                    ans=mid; L=mid+;
                }else R=mid-;
            }
        printf("%d %d\n",ans,select(ans,m));
    }
    return ;
}