HDU2459 后缀数组+RMQ

题目大意：

在原串中找到一个拥有连续相同子串最多的那个子串

比如dababababc中的abababab有4个连续的ab，是最多的

如果有同样多的输出字典序最小的那个

这里用后缀数组解决问题：

枚举连续子串的长度l ， 那么从当前位置0出发每次递增l，拿 i 和 i+l 开头的后缀求一个前缀和val ， 求解依靠RMQ 得到区间 rank(i),rank(i+l)

那么连续的子串个数应该是val/l+1

但是由于你不一定是从最正确的位置出发，那么我们就需要不断将这个i往前推l位，直到某一位字符不匹配，推移的过程中，可能与形成连续串多出的

部分形成一个新的子串，那么个数应该加1，且不断更新推移过程中的rank值，尽量取到rank值小的开头的字符串

 #include <cstdio>
 #include <cstring>
 #include <vector>
 #include <cmath>
 #include <algorithm>
 #include <iostream>
 using namespace std;
 typedef long long ll;
 const int N = ;
 int r[N] , sa[N] , _rank[N] , height[N];
 int wa[N] , wb[N] , wv[N] , wsf[N];
 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 *sa , int n , int m)
 {
     int i,j,p,*x=wa,*y=wb,*t;
     for(i=;i<m;i++)wsf[i]=;
     for(i=;i<n;i++)wsf[x[i]=r[i]]++;
     for(i= ; i<m ; i++) wsf[i]+=wsf[i-];
     for(i=n-;i>=;i--) sa[--wsf[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++) wsf[i]=;
         for(i=;i<n;i++) wsf[wv[i]]++;
         for(i=;i<m;i++) wsf[i]+=wsf[i-];
         for(i=n-;i>=;i--) sa[--wsf[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 callHeight(int *r , int *sa , int n)
 {
     int i,j,k=;
     for(i=;i<=n;i++) _rank[sa[i]]=i;
     for(i=;i<n;height[_rank[i++]]=k)
         for(k?k--:,j=sa[_rank[i]-];r[i+k]==r[j+k];k++);
     return;
 }

 int dp[N<<][] , n;
 char s[N];
 void ST()
 {
     memset(dp , 0x3f , sizeof(dp));
     for(int i= ; i<=n ; i++)dp[i][]=height[i];
     for(int k= ; (<<k)<=n ; k++){
         for(int i= ; i<=n ; i++){
             dp[i][k] = min(dp[i][k-] , dp[i+(<<(k-))][k-]);
         }
     }
 }
 int RMQ(int s , int t)
 {
     int d = t-s+;
     int k = (int)log2(d*1.0);
    // cout<<s<<" "<<t<<" "<<k<<" "<<t-(1<<k)+1<<endl;
     return min(dp[s][k] , dp[t-(<<k)+][k]);
 }
 int RMQPOS(int x , int y)
 {
   //  cout<<"cal pos: "<<x<<" "<<y<<" ";
     x = _rank[x] , y = _rank[y];
     if(x>y){int t=x;x=y,y=t;}
    // cout<<x<<" "<<y<<" "<<RMQ(x+1,y)<<" ";
     return RMQ(x+,y);
 }
 void solve(int &mxTime , int &ansl , int &ansLf)
 {
     ansl = ansLf = mxTime = ;
     for(int l= ; l<=n/ ; l++) {//最外层循环节长度
         for(int i= ; i+l<n ; i+=l){
             int mxl = RMQPOS(i , i+l);
             int time = mxl/l+;

             int del = time*l-mxl , curpos=i , mxRank=_rank[i];
             int t;
             for(t= ; t<l ; t++){
                 if(i<t || s[i-t]!=s[i+l-t]) break;
                 if(_rank[i-t]<mxRank){
                     mxRank = _rank[i-t];
                     curpos = i-t;
                 }
                 if(t==del){
                     mxRank = _rank[i-t];
                     curpos = i-t , time++;
                 }
             }

             if(mxTime<time||(mxTime==time&&mxRank<_rank[ansLf])) mxTime=time,ansl=time*l,ansLf=curpos;

         }
     }
 }
 int main()
 {
     //freopen("a.in"  , "r" , stdin);
     int cas = ;
     while(scanf("%s" , s))
     {
         n = strlen(s);
         if(n== && s[]=='#') break;
         printf("Case %d: ",++cas);
         for(int i= ; i<n ; i++) r[i] = s[i]-'a'+;
         r[n]=;
         da(r,sa,n+,);
         callHeight(r,sa,n);
         ST();

         int mxTime , ansl , ansLf;
         solve(mxTime , ansl , ansLf);
         if(mxTime==){
             char minc='a';
             for(int i= ; i<n ; i++){
                 minc=min(minc , s[i]);
             }
             printf("%c\n" , minc);
             continue;
         }
         for(int i= , j=ansLf ; i<ansl ; i++,j++) printf("%c" , s[j]);
         puts("");
     }
     return ;
 }