poj 2406 后缀数组dc3 板子 计算出现次数最多的循环节 dc3 板子
来源:互联网 发布:eclipse java界面设计 编辑:程序博客网 时间:2024/05/21 19:37
题意:给你一个串,该串可以由其一个子串重复k’次形成,让你求的k的最大值
正解应该是 kmp
题意:给你一个串,该串可以由其一个子串重复k’次形成,让你求的k的最大值
思路:首先得思想是枚举k然后变为判断性问题,接下来就是如何确定成立条件了
1.k必须要能被串长整除,即len%k==0成立
2.sa[0],sa[k]代表的两个串的最长公共前缀的长度应该==len-k,
3.*还有人说需要满足rank[k]-rank[0]=1;但我感觉此条件是k成立的必要条件而已,不要也可以
还有需要注意的一点是,该题的数据量太大,如果用倍增法的话就要tle了,所以要用时间复杂度更低的DC3算法
#include <cstdio>#include <algorithm>#include <cstring>using namespace std;const int N = 1e6+10;#define F(x)((x)/3+((x)%3==1?0:tb)) //F(x)求出原字符串的suffix(x)在新的字符串中的起始位置#define G(x)((x)<tb?(x)*3+1:((x)-tb)*3+2) //G(x)是计算新字符串的suffix(x)在原字符串中的位置,和F(x)为互逆运算int wa[N],wb[N],wv[N],WS[N];int sa[N*3] ; //第i小的后缀,起始位置在源字符串的位置int ra[N],height[N]; //rank 以i为起始位置的后缀在后缀排列中的名次int s[N*3]; //承接字符串,用来计算char str[N]; int c0(int *r,int a,int b) { return r[a]==r[b] && r[a+1]==r[b+1] && r[a+2]==r[b+2];}int c12(int k,int *r,int a,int b) { if(k==2) return r[a]<r[b] || ( r[a]==r[b] && c12(1,r,a+1,b+1) ); else return r[a]<r[b] || ( r[a]==r[b] && wv[a+1]<wv[b+1] );}void sort(int *r,int *a,int *b,int n,int m) { int i; for(i=0; i<n; i++) wv[i]=r[a[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--) b[--WS[wv[i]]]=a[i]; return;}//注意点:为了方便下面的递归处理,r数组和sa数组的大小都要是3*nvoid dc3(int *r,int *sa,int n,int m) { //rn数组保存的是递归处理的新字符串,san数组是新字符串的sa int i , j , *rn = r+n , *san = sa+n , ta = 0 ,tb = (n+1)/3 , tbc = 0 , p; r[n] = r[n+1] = 0; for(i=0; i<n; i++) { if(i%3!=0) wa[tbc++]=i; //tbc表示起始位置模3为1或2的后缀个数 } sort(r+2,wa,wb,tbc,m); sort(r+1,wb,wa,tbc,m); sort(r,wa,wb,tbc,m); for(p=1,rn[F(wb[0])]=0,i=1; i<tbc; i++) rn[F(wb[i])]=c0(r,wb[i-1],wb[i])?p-1:p++; if(p<tbc) dc3(rn,san,tbc,p); else { for(i=0; i<tbc; i++) san[rn[i]]=i; }//对所有起始位置模3等于0的后缀排序 for(i=0; i<tbc; i++) { if(san[i]<tb) wb[ta++]=san[i]*3; } if(n%3==1) //n%3==1,要特殊处理suffix(n-1) wb[ta++]=n-1; sort(r,wb,wa,ta,m); for(i=0; i<tbc; i++) wv[wb[i] = G(san[i])]=i;//合并所有后缀的排序结果,保存在sa数组中 for(i=0,j=0,p=0; i<ta&&j<tbc; p++) sa[p]=c12(wb[j]%3,r,wa[i],wb[j])?wa[i++]:wb[j++]; for(; i<ta; p++) sa[p]=wa[i++]; for(; j<tbc; p++) sa[p]=wb[j++]; return;}//height[i]=suffix(sa[i-1])和suffix(sa[i])的最长公共前缀,也就是排名相邻的两个后缀的最长公共前缀void calheight(int *r,int *sa,int n) { int i,j,k=0; for(i=1; i<=n; i++) ra[sa[i]]=i; for(i=0; i<n; height[ra[i++]]=k) for(k?k--:0,j=sa[ra[i]-1]; r[i+k]==r[j+k]; k++);}int main(){ while(scanf("%s",str)!=EOF) { if(str[0]=='.') break; int len=strlen(str); for(int i=0;i<len;i++) s[i]=str[i]-'a'+1; s[len]=0; dc3(s,sa,len+1,105); calheight(s,sa,len); int flag=0; for(int i=1;i<=len;i++) { if(len % i == 0 && ra[0] == ra[i] + 1 && height[ra[0]] == len - i) { printf("%d\n",len/i ); flag=1; break; } } if(!flag) printf("1\n"); }}
阅读全文
0 0
- poj 2406 后缀数组dc3 板子 计算出现次数最多的循环节 dc3 板子
- POJ 2406 Power Strings(KMP or 后缀数组 dc3板子)
- poj 2406 Power Strings //DC3 后缀数组
- 统计每个长度 出现的最多次数 后缀自动机 板子
- 重复次数最多的连续字串 后缀数组板子 黑盒子
- 后缀数组的DC3模版【后缀数组】
- POJ 2406 Power Strings <后缀数组(DC3) / KMP>
- POJ 2406Power Strings (dc3 后缀数组)
- 后缀数组da+dc3
- 后缀数组 dc3 模板
- 构造后缀数组的DC3算法实现
- 后缀数组的DC3构造法
- dc3 后缀数组构建 code
- POJ 2774 Long Long Message <后缀数组(DC3)>
- poj 2406 连续重复子串(后缀数组DC3/kmp)
- 实现后缀数组的倍增算法和DC3算法
- 后缀数组的倍增模板与DC3模板
- 【*】后缀数组(dc3算法构造)
- 聊聊你不知道的建造者设计模式
- C# 委托 事件 观察者模式 总结
- 缓存无关算法
- git的一些基础命令
- BZOJ 2440: [中山市选2011]完全平方数 莫比乌斯反演 +容斥原理
- poj 2406 后缀数组dc3 板子 计算出现次数最多的循环节 dc3 板子
- Linux基本命令入门
- 数组小练习
- GCD HDU
- jvm相关参数及其含义
- hadoop启动不正常问题解决
- 解决oracle11g的最大内存占用
- 使用ffmpeg循环推流(循环读取视频文件)推送EasyDSS RTMP流媒体服务器的方法
- ubuntu 安装软件后将软件添加进启动器