POJ

传送门：POJ - 3693  

题意：给你一个字符串，求重复次数最多的连续重复子串，如果有一样的，取字典序小的字符串。

题解：

比较容易理解的部分就是枚举长度为L，然后看长度为L的字符串最多连续出现几次。既然长度为L的串重复出现，那么str[0],str[l],str[2*l]……中肯定有两个连续的出现在字符串中。

那么就枚举连续的两个，然后从这两个字符前后匹配，看最多能匹配多远。即以str[i*l],str[i*l+l]前后匹配，这里是通过查询suffix(i*l),suffix(i*l+l)的最长公共前缀。通过rank值能找到i*l,与i*l+l的排名，我们要查询的是这段区间的height的最小值，通过RMQ预处理达到查询为0(1)的复杂度

设LCP长度为M, 则答案显然为M / L + 1, 但这不一定是最好的, 因为答案的首尾不一定再我们枚举的位置上. 我的解决方法是, 我们考虑M % L的值的意义, 我们可以认为是后面多了M % L个字符, 但是我们更可以想成前面少了(L - M % L)个字符! 所以我们求后缀j * L - (L - M % L)与后缀(j+ 1) * L - (L - M % L)的最长公共前缀。即把之前的区间前缀L-M%L即可。

然后把可能取到最大值的长度L保存，由于 题目要求字典序最小，通过sa数组进行枚举，取到的第一组，肯定是字典序最小的。

题解的出处：https://blog.csdn.net/acm_cxlove/article/details/7941205 

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

const int maxn = 1e5+10;
int wa[maxn],wb[maxn],wsf[maxn],wv[maxn],sa[maxn];
int rnk[maxn],height[maxn];
char s[maxn];
int r[maxn];

//sa:字典序中排第i位的起始位置在str中第sa[i]  sa[1~n]为有效值

//rnk:就是str第i个位置的后缀是在字典序排第几 rnk[0~n-1]为有效值

//height：字典序排i和i-1的后缀的最长公共前缀  height[2~n]为有效值，第二个到最后一个

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

void getsa(int *r,int *sa,int n,int m)//n为添加0后的总长
{
    int i,j,p,*x=wa,*y=wb,*t;
    for(i=0; i<m; i++)  wsf[i]=0;
    for(i=0; i<=n; i++)  wsf[x[i]=r[i]]++;
    for(i=1; i<m; i++)  wsf[i]+=wsf[i-1];
    for(i=n; i>=0; i--)  sa[--wsf[x[i]]]=i;
    p=1;
    j=1;
    for(; p<=n; j*=2,m=p){
        for(p=0,i=n+1-j; i<=n; i++)  y[p++]=i;
        for(i=0; i<=n; i++)  if(sa[i]>=j)  y[p++]=sa[i]-j;
        for(i=0; i<=n; i++)  wv[i]=x[y[i]];
        for(i=0; i<m; i++)  wsf[i]=0;
        for(i=0; i<=n; i++)  wsf[wv[i]]++;
        for(i=1; i<m; i++)  wsf[i]+=wsf[i-1];
        for(i=n; i>=0; i--)  sa[--wsf[wv[i]]]=y[i];
        swap(x,y);
        x[sa[0]]=0;
        for(p=1,i=1; i<=n; i++)
            x[sa[i]]=cmp(y,sa[i-1],sa[i],j)? p-1:p++;
    }
}

void getheight(int *r,int n)//n为添加0后的总长
{
    int i,j,k=0;
    for(i=1; i<=n; i++)  rnk[sa[i]]=i;
    for(i=0; i<n; i++){
        if(k)
            k--;
        else
            k=0;
        j=sa[rnk[i]-1];
        while(r[i+k]==r[j+k])
            k++;
        height[rnk[i]]=k;
    }
}

int dp[maxn][40];

void rmq_init(int n){

    int m=floor(log(n+0.0)/log(2.0));
    for(int i=1;i<=n;i++)dp[i][0]=height[i];
    for(int j=1;j<=m;j++){
        for(int i=n;i;i--){
            dp[i][j]=dp[i][j-1];
            if(i+(1<<(j-1))<=n){
                dp[i][j]=min(dp[i][j],dp[i+(1<<(j-1))][j-1]);
            }
        }
    }
}

int rmq(int l,int r){

    int a=rnk[l],b=rnk[r];
    if(a>b)swap(a,b);
    a++;
    int k=floor(log(b-a+1.0)/log(2.0));
    return min(dp[a][k],dp[b-(1<<k)+1][k]);
}

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int t=1;
    while(cin>>s){
        if(s[0]=='#') break;
        int len=strlen(s);
        for(int i=0;i<len;i++) r[i]=s[i]-'a'+1;
        r[len]=0;
        getsa(r,sa,len,150);
        getheight(r,len);
        rmq_init(len);
        int ans=0;
        int pos=0,p=0;
        for(int k=1;k<len;k++){     //枚举长度
            for(int i=0;i+k<len;i+=k){  //第i段
                int n=rmq(i,i+k);   //每一段的公共前缀最小值
                n--;
                for(int j=0;j<=k-1;j++){    //枚举每一段的起点
                    int now=i-j;
                    if((now<0||s[now]!=s[now+k])&&j) break;
                    n++;
                    int sum=n/k+1;
                    if(sum>ans||(sum==ans&&rnk[now]<rnk[pos])){
                        ans=sum;
                        pos=now;
                        p=k;
                    }
                }
            }
        }
        cout<<"Case "<<t++<<": ";
        if(ans<=1){
            char tmp='z';
            for(int i=0;i<len;i++) tmp=min(tmp,s[i]);
            cout<<tmp<<endl;
        }
        else{
            for(int i=0;i<ans*p;i++){
                cout<<s[i+pos];
            }
            cout<<endl;
        }
    }
    return 0;
}