poj 3693 Maximum repetition substring（08合肥 RMQ+后缀数组）

来源：互联网发布：淘宝直通车如何出价编辑：程序博客网时间：2024/05/17 03:56

传送门：http://poj.org/problem?id=3693

题目：给出一个串，求重复次数最多的连续重复子串；

分析：

枚举重复单元的长度，然后理所当然的枚举起点。利用rmq处理，后缀i，i+l的最长前缀。

lcp/l+1，为当前相邻l长度单元的串的重复次数，但是由于i+=l，提高了效率，但是i不一定刚好是重复串的起点，所以如果r%l！=0，把串往前移l-r%l个单位。找到最大的重复次数，然后利用sa数组的排名（1-n），，可以字典序输出。

详细链接：http://blog.csdn.net/acm_cxlove/article/details/7941205

代码：

#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>#include<vector>#include<string>#include<queue>#include<deque>#include<stack>#include<map>#include<set>#define INF 0x7fffffff#define SUP 0x80000000#define mem(a,b) memset(a,b,sizeof(a))using namespace std;typedef long long LL;const int N=100007;/******************auther:ACsorry******************///以下为倍增算法求后缀数组int wa[N],wb[N],wv[N],Ws[N];int cmp(int *r,int a,int b,int l){    return r[a]==r[b]&&r[a+l]==r[b+l];}void da(const char *r,int *sa,int n,int m){int i,j,p,*x=wa,*y=wb,*t;for(i=0;i<m;i++) Ws[i]=0;for(i=0;i<n;i++) Ws[x[i]=r[i]]++;for(i=1;i<m;i++) Ws[i]+=Ws[i-1];for(i=n-1;i>=0;i--) sa[--Ws[x[i]]]=i;for(j=1,p=1;p<n;j*=2,m=p){for(p=0,i=n-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++) Ws[i]=0;for(i=0;i<n;i++) Ws[wv[i]]++;for(i=1;i<m;i++) Ws[i]+=Ws[i-1];for(i=n-1;i>=0;i--) sa[--Ws[wv[i]]]=y[i];for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;}return;}int sa[N],Rank[N],height[N];//求height数组void calheight(const char *r,int *sa,int n){int i,j,k=0;for(i=1;i<=n;i++) Rank[sa[i]]=i;for(i=0;i<n;height[Rank[i++]]=k)for(k?k--:0,j=sa[Rank[i]-1];r[i+k]==r[j+k];k++);return;}int dp[N][20];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=1;i+(1<<j)-1<=n;i++){            dp[i][j]=min(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);        }    }}int Rmq_query(int L,int R){    int a=Rank[L],b=Rank[R];    if(a>b) swap(a,b);    a++;//这个地方++，应为height[a]处理的是 a，a-1的公共串    int m=floor((log(b-a+1.0)/log(2.0)));    return min(dp[a][m],dp[b-(1<<m)+1][m]);}char str[N];int main(){    int cas=0;    while(scanf("%s",str)==1&&str[0]!='#')    {        int n=strlen(str);        da(str,sa,n+1,130);        calheight(str,sa,n);        Rmq_init(n);        int cnt=0,maxx=0,a[N];        for(int l=1;l<n;l++){            for(int i=0;i+l<n;i+=l){    //j+=l.复杂度的需要                int r=Rmq_query(i,i+l);                int step=r/l+1;                int k=i-(l-r%l);                if(k>=0&&r%l){                    if(Rmq_query(k,k+l)>=r)                        step++;                }                if(step>maxx)                {                    maxx=step;                    cnt=0;                    a[cnt++]=l;                }                else if(step==maxx)                    a[cnt++]=l;            }        }        int len=-1,st;        for(int i=1;i<=n&&len==-1;i++){            for(int j=0;j<cnt;j++){                if(Rmq_query(sa[i],sa[i]+a[j])/a[j]+1==maxx)                {                    len=a[j];                    st=sa[i];                    break;                }            }        }        str[st+maxx*len]=0;        printf("Case %d: %s\n",++cas,str+st);    }    return 0;}

//宁愿精彩的活，也不愿平庸一辈子。

0 1