poj 3693 重复次数最多的子串
来源:互联网 发布:淘宝 海外买手现场 编辑:程序博客网 时间:2024/06/08 07:42
#include<stdio.h>#include<iostream>#include<string.h>#include<algorithm>using namespace std;const int nMax=100005;int r[nMax],a[nMax];int sa[nMax], rank[nMax], height[nMax];int wa[nMax], wb[nMax], wv[nMax], wd[nMax];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 n, int m){ // 倍增算法 r为待匹配数组 n为总长度(包含结尾的0) m为字符范围 int i, j, p, *x = wa, *y = wb, *t; for(i = 0; i < m; i ++) wd[i] = 0; for(i = 0; i < n; i ++) wd[x[i]=r[i]] ++; for(i = 1; i < m; i ++) wd[i] += wd[i-1]; for(i = n-1; i >= 0; i --) sa[-- wd[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 ++) wd[i] = 0; for(i = 0; i < n; i ++) wd[wv[i]] ++; for(i = 1; i < m; i ++) wd[i] += wd[i-1]; for(i = n-1; i >= 0; i --) sa[-- wd[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 ++; } }}void calHeight(int *r, int n){ // 求height数组。 int i, j, k = 0; for(i = 1; i <= n; i ++) rank[sa[i]] = i; // 1->n for(i = 0; i < n; i++){ for(k ? k -- : 0, j = sa[rank[i]-1]; r[i+k] == r[j+k]; k ++); height[rank[i]] = k; }}int Log[nMax],dp[nMax][20];void initRMQ(int n){ // 求log2(x) Log[0]=-1; for(int i=1;i<=nMax;i++) Log[i]=(i&(i-1))?Log[i-1]:Log[i-1]+1; for(int i=1;i<=n;i++) dp[i][0]=height[i]; for(int j=1;j<=Log[n];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 lcp(int x,int y){ int a=rank[x],b=rank[y]; if(a>b) swap(a,b); a++; int k=Log[b-a+1]; return min(dp[a][k],dp[b-(1<<k)+1][k]);}int main(){// freopen("test.txt","r",stdin); char s[nMax]; int num=0; while(scanf("%s",s) && s[0]!='#') { int n=strlen(s);num++; for(int i=0;i<n;i++) r[i]=s[i]-'a'+1; r[n]=0; da(r,n+1,27); calHeight(r,n); initRMQ(n); int ret=0,cnt=0,flag=0; for(int len=1;len<=n;len++) { for(int j=len;j<=n;j+=len) { int lc=lcp(j-len,j),m=len-lc%len; if(j>len && (lc%len)) lc=max(lc,lcp(j-len-m,j-m)); lc=lc/len+1; if(lc>ret) { ret=lc; cnt=0; a[cnt++]=len; } else if(lc==ret && a[cnt-1]!=len) a[cnt++]=len; } } for(int i=1;i<=n;i++) { for(int j=0;j<cnt;j++) if(lcp(sa[i],sa[i]+a[j])>=a[j]*(ret-1)) { flag=1; printf("Case %d: ",num); for(int k=0;k<a[j]*ret;k++) printf("%c",s[k+sa[i]]); printf("\n"); break; } if(flag) break; } } return 0;}
- poj 3693 重复次数最多的子串
- POJ 3693 求重复次数最多的子串
- 连续重复次数最多的重复子串 hdu2459 poj 3693
- POJ 3693 重复次数最多的连续重复子串 后缀数组
- POJ - 3693 Maximum repetition substring(后缀数组求重复次数最多的连续重复子串)
- 后缀数组(重复次数最多的最长重复子串)—— POJ 3693
- poj 3693 后缀数组 重复次数最多的连续重复子串
- poj 3693 重复次数最多的连续重复子串(后缀数组+RMQ)
- POJ 3693 Maximum repetition substring(后缀数组[重复次数最多的连续重复子串])
- poj 3693 求出现次数最多的连续重复子串(具体的串)黑盒
- POJ 3693 Maximum repetition substring (求重复次数最多的连续子串,4级)
- POJ 3693 重复次数最多的子串是什么?【后缀数组】
- POJ 3693 Maximum repetition substring(重复次数最多的连续子串 字典序最小)
- POJ 3693 后缀数组 重复次数最多的连续重复子串 倍增法以及D3法
- poj 3693/hdu 2459 Maximum repetition substring spoj 687. Repeats ( 后缀数组 重复次数最多的连续重复子串)
- Maximum repetition substring+POj+后缀数组之求重复次数最多的连续重复子串
- 后缀数组(重复次数最多的连续重复子串)
- POJ-3693--后缀数组求字典序最小重复次数最多子串
- 埃及分数,迭代加深搜索
- Web Application的invoker servlet
- 使用PL/SQL 实现Oracle Streams AQ
- 挂载大文件作为新磁盘
- urlrewrite(伪静态技术运用)
- poj 3693 重复次数最多的子串
- 在ArcEngine中使用INewLineFeedback实现 距离测量(面积测量)的正确方法
- C#使用System.Timers.Timer类实现一个Button锁定功能(windows form)
- mysql中文乱码总结
- 深入C++的new
- 论坛里搜集的正则 干货留存
- 使用List集合做数据源,并且使用AspNetPager来分页
- 巧用记事本记录时间
- JavaScript正则表达式的使用篇-更新中~