poj 3693 Maximum repetition substring

呵呵呵呵呵呵呵呵呵呵，sb（神犇）题看了一天，还是不懂

题目要求的是最多重复的，那么就来找重复的，可以先枚举一个重复的单元（比如ababab，就枚举ab）的长度，

然后再原串中，会有ch[0],ch[L],c[2*L],,,,,这些，如果重复的话肯定是会有ch[i*L]==ch[(i+1)*L]的，那么我们枚举这个东西。

在找到之后，需要做的就是找ch[(i+1)*L]之后可以有多长重复  L1，以及ch[i*L]之前有多少重复  L2，这个问题用后缀数组的height处理出的RMQ解决（呵呵呵），这一步貌似是叫求两个后缀的最长公共前缀。。。（对余ch[i*L]往前的，就把串反过来搞），然后答案就是(L1+L2)/L+1。

然而这个恶心题还要求最小字典序。。。。那么考虑，另外寻找重复部分的开头，我们知道，已经找出了重复串的长度，那么在这里面的任意位置都是满足重复这个性质的，所以考虑(L1+L2)%L（长度不及L的部分，可以画图看一下），然后在搞一个ST表来记录区间的rank，然后在 i*L-L2 到 i*L-L2+(L1+L2)%L查询最小的rank就好。。

（本蒟蒻就会这么多，大神勿喷2333）

 #include<cstdio>
 #include<cstring>
 #include<iostream>
 #define N 100005
 #define LL long long
 #define inf 0x3f3f3f3f
 using namespace std;
 int n,mx,ans,ansl,ansr;
 int bin[],Log[];
 int mn[][];
 char ch[];
 void rmq(int mn[N][], int *a)
 {
     for (int i=; i<=n; i++) mn[i][]=a[i];
     for (int i=; i<=Log[n]; i++)
         for (int j=; j<=n; j++)
             if (j+bin[i]-<=n)
                 mn[j][i]=min(mn[j][i-],mn[j+bin[i-]][i-]);
             else break;
 }
 struct Suffix{
     int k,p,q;
     int rk[][N],sa[][N],v[N],a[N];
     int h[N],mn[N][];
     void clear()
     {
         memset(a,,sizeof(a));
         memset(v,,sizeof(v));
         memset(rk,,sizeof(rk));
     }
     void cal_sa(int *sa, int *rk, int *SA, int *RK)
     {
         for (int i=; i<=n; i++) v[rk[sa[i]]]=i;
         for (int i=n; i>=; i--)
             if (sa[i]>k) SA[v[rk[sa[i]-k]]--]=sa[i]-k;
         for (int i=n-k+; i<=n; i++) SA[v[rk[i]]--]=i;
         for (int i=; i<=n; i++)
             RK[SA[i]]=RK[SA[i-]]+(rk[SA[i-]]!=rk[SA[i]] || rk[SA[i-]+k]!=rk[SA[i]+k]);
     }
     void get_sa()
     {
         p=,q=;
         for (int i=; i<=n; i++) v[a[i]]++;
         for (int i=; i<=; i++) v[i]+=v[i-];
         for (int i=; i<=n; i++) sa[p][v[a[i]]--]=i;
         for (int i=; i<=n; i++)
             rk[p][sa[p][i]]=rk[p][sa[p][i-]]+(a[sa[p][i]]!=a[sa[p][i-]]);
         for (k=; k<n; k<<=,swap(p,q)) cal_sa(sa[p],rk[p],sa[q],rk[q]);
     }
     void get_height()
     {
         k=;
         for (int i=; i<=n; i++)
         {
             if (rk[p][i]==) h[]=;
             else{
                 int j=sa[p][rk[p][i]-];
                 while (a[i+k]==a[j+k]) k++;
                 h[rk[p][i]]=k;
                 if (k>) k--;
             }
         }
     }
     void pre()
     {
         get_sa(); get_height(); rmq(mn,h);
     }
     int lcp(int x, int y)
     {
         x=rk[p][x],y=rk[p][y];
         if (x>y) swap(x,y);
         x++;
         int t=Log[y-x+];
         return min(mn[x][t],mn[y-bin[t]+][t]);
     }
 }c[];
 int query(int x, int y)
 {
     int t=Log[y-x+];
     return min(mn[x][t],mn[y-bin[t]+][t]);
 }
 void solve(int L)
 {
     int l=,r=,t;
     for (int i=; i+L<=n; i+=L)
         if (ch[i]==ch[i+L])
         {
             r=c[].lcp(i,i+L),l=c[].lcp(n-i+,n-i-L+);
             if ((l+r)/L+>mx)
                 mx=(l+r)/L+,ans=inf;
             if ((l+r)/L+==mx)
             {
                 t=query(i-l,i-l+(l+r)%L);
                 if (t<ans)
                 {
                     ans=t;
                     ansl=c[].sa[c[].p][t],ansr=ansl+mx*L-;
                 }
             }
         }
 }
 int main()
 {
     bin[]=; for (int i=; i<=; i++) bin[i]=bin[i-]<<;
     Log[]=-; for (int i=; i<=; i++) Log[i]=Log[i/]+;
     int txt=;
     while (scanf("%s",ch+))
     {
         if (ch[]=='#') break;
         printf("Case %d: ",++txt);
         c[].clear(); c[].clear();
         n=strlen(ch+);
         for (int i=; i<=n; i++) c[].a[i]=ch[i]-'a'+;
         for (int i=; i<=n; i++) c[].a[i]=ch[n-i+]-'a'+;
         c[].pre(); c[].pre();
         rmq(mn,c[].rk[c[].p]);
         mx=; ans=inf;
         for (int i=; i<=n; i++)
             if (c[].rk[c[].p][i]<ans)
                 ans=c[].rk[c[].p][i],ansl=ansr=i;
         for (int i=; i<=n; i++) solve(i);
         for (int i=ansl; i<=ansr; i++)
             putchar(ch[i]);
         puts("");
     }
     return ;
 }