poj 3693 Maximum repetition substring
来源:互联网 发布:windows桌面图标没了 编辑:程序博客网 时间:2024/06/05 01:12
题目链接:http://poj.org/problem?id=3693
题目思路:这个题目和前面那个题目的本质是一样的,只是要求最小字典序,由于一个解中对应几个可行串,所以要暴力查找,最坏复杂度很高,由于数据水,还是勉强过了。
还有一种最小字典序的处理方法,那就是在第一步中记录能够达到最大重复数的重复长度的集合,然后按照sa[1],sa[2],sa[3],……的顺序去暴力枚举,每到一个位置就从小到大地判重复长度的集合中有没有可行的,如果有一个可行的就立即输出结果并退出。
#include<stdio.h>#include<stdlib.h>#include<string.h>#include<string>#include<queue>#include<algorithm>#include<vector>#include<stack>#include<list>//#include<iostream>#include<map>using namespace std;#define inf 0x3f3f3f3f#define M 110000int max(int a,int b){return a>b?a:b;}int min(int a,int b){return a<b?a:b;}int sa[M],rank[M],height[M],r[M];int wa[M],wb[M],wv[M],ws[M],dp[20][M],mm[M];bool cmp(int *y,int a,int b,int l){ return y[a]==y[b]&&y[a+l]==y[b+l];}void da(int n,int m){ int i,j,p; int *x=wa,*y=wb; 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) { p=0; for(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<m;i++) ws[i]=0; for(i=0;i<n;i++) wv[i]=x[y[i]]; 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]; swap(x,y); x[sa[0]]=0; p=1; for(i=1;i<n;i++) { if(cmp(y,sa[i-1],sa[i],j)) x[sa[i]]=p-1; else x[sa[i]]=p++; } }}void calh(int n){ int i,k,tmp; for(i=1;i<=n;i++) rank[sa[i]]=i; k=0; for(i=0;i<n;i++) { tmp=sa[rank[i]-1]; for(;r[i+k]==r[tmp+k];k++) ; height[rank[i]]=k; k?--k:0; }}void init(int n){ int i,j; mm[0]=-1; for(i=1;i<=n;i++) { if((i&(i-1))==0) mm[i]=mm[i-1]+1; else mm[i]=mm[i-1]; } for(i=1;i<=n;i++) dp[0][i]=height[i]; for(i=1;i<=mm[n];i++) { int tmp=1<<(i-1); for(j=1;j+tmp<=n;j++) dp[i][j]=min(dp[i-1][j],dp[i-1][j+tmp]); }}int lcp(int a,int b){ a=rank[a],b=rank[b]; if(a>b) swap(a,b); a++; int l=mm[b-a+1]; b-=(1<<l)-1; return min(dp[l][a],dp[l][b]);}char s[M];int main(){ int i,j,g,t,n,num,k,ans,rec1,rec2,st,ed,count=1; while(scanf("%s",s)!=EOF) { if(s[0]=='#') break; n=strlen(s); for(i=0;i<n;i++) { r[i]=s[i]; } r[n]=0; da(n+1,128); calh(n); init(n); // for(i=1;i<=n;i++) // { // printf("i %d sa %d height %d\n",i,sa[i],height[i]); // } // puts("akkk"); ans=1; st=ed=0; for(i=1;i<=n/2;i++) { for(j=0;j+i<n;j+=i) { k=lcp(j,j+i); num=k/i+1; rec1=j; rec2=j+num*i; int tmp=j-(i-k%i); g=j-1; if(tmp>=0&&lcp(tmp,tmp+i)>=i) { num++; rec1=tmp; rec2=tmp+num*i; g=tmp-1; } for(;g>=0&&g>j-i;g--) { if(lcp(g,g+i)>=j-g) { if(rank[g]<rank[rec1]) { rec1=g; rec2=g+num*i; } } else break; } if(ans<num) { ans=num; st=rec1; ed=rec2; } else if(ans==num) { //printf("i %d ans %d st %d rec1 %d\n",i,ans,st,rec1); tmp=lcp(st,rec1); if(tmp<rec2-rec1) { if(rank[rec1]<rank[st]) { st=rec1; ed=rec2; } } else { if(rec2-rec1<ed-st) { st=rec1; ed=rec2; } } } } } //printf("%d %d\n",st,ed); printf("Case %d: ",count++); for(i=st;i<ed;i++) { printf("%c",r[i]); } puts(""); }}
- poj 3693 Maximum repetition substring
- POJ 3693 Maximum repetition substring
- POJ 3693 Maximum repetition substring
- poj 3693 Maximum repetition substring
- POJ 3693 Maximum repetition substring
- POJ 3693 Maximum repetition substring
- POJ 3693Maximum repetition substring
- poj 3693 Maximum repetition substring
- Maximum repetition substring POJ
- Maximum repetition substring POJ
- poj 3693 Maximum repetition substring //后缀数组
- poj 3693 Maximum repetition substring (后缀数组)
- poj 3693 Maximum repetition substring(后缀数组)
- poj 3693 Maximum repetition substring(后缀数组)
- 【后缀数组】【poj 3693】Maximum repetition substring
- POJ 3693 Maximum repetition substring(好题)
- POJ 3693 Maximum Repetition Substring 后缀数组
- 「POJ 3693」Maximum repetition substring
- 深入探索ScrollWindow
- T_SQL 2008 语言基础学习笔记
- Android 通过软引用实现图片缓存,防止内存溢出
- Linux下java环境变量设置方法
- hdu2531 Catch him
- poj 3693 Maximum repetition substring
- java.lang.NoClassDefFoundError: org.apache.commons.lang.StringEscapeUtils
- Android ViewPager多页面滑动切换以及动画效果
- linux 命令大全
- Oracle数据库11g新特性:数据库重放
- 胜者树与败者树
- VLC学习文档(载)
- ADO对象连接数据库
- 异常