Power String poj2406(后缀数组+DC3算法)
来源:互联网 发布:子域名怎么看 编辑:程序博客网 时间:2024/06/06 05:48
这题其实应该不是用后缀数组做的。。。但是就当做练习吧。
倍增算法在这里会被卡掉,所以要用DC3算法。。这里我是套模板的。。
重点就在于,怎么找这个循环节。
做法是把height数组出来之后,枚举一下长度。
设长度为i,如果满足条件的话。
后缀0(整个串)和后缀i的LCP是 length-i
这个做法很巧妙,学习了。
//42704K 2829MS C++ 2709B 2013-12-18 13:31:43/* 题意: 给出一个由某个子串重复R次组成的字符串,求R的最大值 后缀数组: KMP应该会简单些,因为此处要练习后缀数组,故用后缀数组。首先考虑用DA的话会TLE, 因为其时间复杂度为O(nlgn),数据太大(n=2000000)。 此处用dc3,dc3的算法没研究,只是套用其模板,dc3算法会比DA快些,在此处勉勉强强 的过了。先用dc3求出后缀数组,然后再求height数组,再穷举字符串长度len。求k能整除len 且suffix(1) 与 suffix(k+1) 最长前缀等于len-k。 在求最长公共前缀时,由于suffix(1)是固定的,利用height数组的特性,求出height数组 中每一个数到height[rank[0]]间的最小值即可。*/#include<stdio.h>#include<string.h>#include<math.h>#include<algorithm>using namespace std;#define N 2000005#define F(x) ((x)/3+((x)%3==1?0:tb))#define G(x) ((x)<tb?(x)*3+1:((x)-tb)*3+2)int wa[N],wb[N],wv[N],ws[N];int rank[N],height[N];int sa[N],r[N];char c[N];int s[N];int cmp(int *y,int a,int b,int l){ return y[a]==y[b]&&y[a+l]==y[b+l];}int c0(int *y,int a,int b){ return y[a]==y[b]&&y[a+1]==y[b+1]&&y[a+2]==y[b+2];}int c12(int k,int *y,int a,int b){ if(k==2) return y[a]<y[b]||y[a]==y[b]&&c12(1,y,a+1,b+1); else return y[a]<y[b]||y[a]==y[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;}void dc3(int *r,int *sa,int n,int m){ 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; 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; for(i=0;i<tbc;i++) if(san[i]<tb) wb[ta++]=san[i]*3; if(n%3==1) wb[ta++]=n-1; sort(r,wb,wa,ta,m); for(i=0;i<tbc;i++) wv[wb[i]=G(san[i])]=i; 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;}void get_height(int *r,int n){ int i,j,k=0; for(i=0;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 rm[N];int sss[N];char arr[N];int main(){ while(~scanf("%s",arr)) {// memset(rm,0,sizeof(rm));// memset(height,0,sizeof(height)); int len=strlen(arr); if(len==1 && arr[0]=='.') break; for(int i=0;i<len;i++) r[i]=r[i]+1; r[len]=0; dc3(r,sa,len+1,256); get_height(r,len); int temp=height[rank[0]]; for(int i=rank[0]-1;i>=1;i--) { rm[i]=temp; temp=min(temp,height[i]); //cout<<rm[i]<<endl; } temp=0x3f3f3f3f; for(int i=rank[0]+1;i<=len;i++) { temp=min(temp,height[i]); rm[i]=temp; //cout<<rm[i]<<endl; } int ans=1; for(int i=1;i<=len;i++) { if(len%i!=0) continue; else { if(rm[rank[i]]==len-i) { ans=len/i; break; } } } printf("%d\n",ans); } return 0;}
0 0
- Power String poj2406(后缀数组+DC3算法)
- poj2406 Power Strings 后缀数组 DC3算法
- POJ2406:Power Strings(后缀数组DC3)
- 【*】后缀数组(dc3算法构造)
- POJ 2406 Power Strings <后缀数组(DC3) / KMP>
- POJ 2406Power Strings (dc3 后缀数组)
- poj 2406 Power Strings //DC3 后缀数组
- POJ2406——Power Strings(后缀数组)
- [KMP][后缀数组] poj2406 Power Strings
- POJ2406 Power Strings 【KMP 或 后缀数组】
- 构造后缀数组的DC3算法实现
- POJ 2406 Power Strings(KMP or 后缀数组 dc3板子)
- 后缀数组罗穗蹇模板(dc3)
- 后缀数组da+dc3
- 后缀数组 dc3 模板
- zoj 1905 Power String(后缀数组)
- 实现后缀数组的倍增算法和DC3算法
- 【kmp】【后缀数组】【连续重复子串】Power Strings POJ2406
- HDU 5616 Jam's balance
- UVA - 146 ID Codes
- Qt在Windows下的三种编程环境搭建
- UVA - 10098 Generating Fast
- UVa 1400
- Power String poj2406(后缀数组+DC3算法)
- 弟子规详解
- 蓝鸥 8.16
- Dijkstra、Bellman-Ford及Spfa算法思想对比
- 09 redis集群的搭建 以及遇到的问题
- UVA - 729 The Hamming Distance Problem
- 给RecyclerView最纯粹的下拉刷新和上拉加载更多
- [poj1511] Invitation Cards
- JSP 使用<%@include%>报Duplicate local variable path 错误 解决方法