hdu 2459 Maximum repetition substring(后缀数组)

题目链接：hdu 2459 Maximum repetition substring

题意：

让你找一个重复最多的子串，并且输出。

题解：

这个是论文题，看的cxlove的题解，不是很理解为什么这样就能完全找完，当作结论使吧。

#include<bits/stdc++.h>
#define F(i,a,b) for(int i=a;i<=b;++i)
using namespace std;
namespace suffixarray{    
    #define FN(n) for(int i=0;i<n;i++)
    const int N =1E5+7;//字符串长度
    int rnk[N],sa[N],height[N],c[N];char s[N];
    void getsa(int n,int m,int *x=rnk,int *y=height){
        FN(m)c[i]=0;FN(n)c[x[i]=s[i]]++;FN(m)c[i+1]+=c[i];
        for(int i=n-1;i>=0;i--)sa[--c[x[i]]]=i;
        for(int k=1,p;p=0,k<=n;k=p>=n?N:k<<1,m=p){
            for(int i=n-k;i<n;i++)y[p++]=i;
            FN(n)if(sa[i]>=k)y[p++]=sa[i]-k;
            FN(m)c[i]=0;FN(n)c[x[y[i]]]++;FN(m)c[i+1]+=c[i];
            for(int i=n-1;i>=0;i--)sa[--c[x[y[i]]]]=y[i];
            swap(x,y),p=1,x[sa[0]]=0;
            for(int 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++;
        }
        FN(n)rnk[sa[i]]=i;
        for(int i=0,j,k=0;i<n-1;height[rnk[i++]]=k)
        for(k=k?k-1:k,j=sa[rnk[i]-1];s[i+k]==s[j+k];k++);
    }
}
using namespace suffixarray;

int f[N][20],ans[N],cas;
void rmq(int *a,int n)
{
    for(int i=1;i<=n;i++)f[i][0]=a[i];
    for(int j=1;1<<j<n;j++)for(int i=1;i<=n;i++)
    if(i+(1<<j)-1<=n)f[i][j]=min(f[i][j-1],f[i+(1<<j-1)][j-1]);
    else break;
}
inline int find(int l,int r)
{
    l=rnk[l],r=rnk[r];
    if(l>r)swap(l,r);
    l++;
    int k=31-__builtin_clz(r-l+1);
    return min(f[l][k],f[r-(1<<k)+1][k]);
}

int main()
{
    while(scanf("%s",s),s[0]!='#')
    {
        int len=strlen(s);
        getsa(len+1,130);
        rmq(height,len);
        int cnt=0,mx=0;
        F(l,1,len-1)for(int i=0;i+l<len;i+=l)
        {
            int r=find(i,i+l);
            int step=r/l+1,k=i-(l-r%l);
            if(k>=0&&r%l)step+=(find(k,k+l)>=r);
            if(step>mx)
            {
                mx=step,cnt=0;
                ans[cnt++]=l;
            }else if(step==mx)ans[cnt++]=l;
        }
        int tmp=-1,st;
        for(int i=1;i<=len&&tmp==-1;i++)F(j,0,cnt-1)
        {
            int l=ans[j];
            if(find(sa[i],sa[i]+l)>=(mx-1)*l)
            {
                tmp=l,st=sa[i];
                break;
            }
        }
        printf("Case %d: ",++cas);  
        for(int i=st,j=0;j<tmp*mx;j++,i++) printf("%c",s[i]);  
        printf("
");
    }
    return 0;
}