bzoj4829 [TJOI2017]dna(后缀数组)
来源:互联网 发布:数据预测模型有哪些 编辑:程序博客网 时间:2024/06/06 07:22
luogu3763有题面可以去看的说。。大牛分站开o2优化了(没错我就是这样才过的。)后缀数组处理一下lcp,然后枚举三次就好了。bzoj上我的总是T。嗯,一定是姿势不对。。不过我看神犇们都是用SAM写的,贼快。我也要去学习一下。。
#include <cstdio>#include <cstring>#define N 200005char s[N];int tst,n,len1,len2,m,st[N][19],log[N];int rank[N<<1],rank1[N],sa[N],h[N],tmp[N],count[N];inline int swap(int &x,int &y){x^=y;y^=x;x^=y;}inline int min(int x,int y){return x<y?x:y;}inline int lcp(int x,int y){int l=rank[x],r=rank[y];if(l>r) swap(l,r);++l;int k=log[r-l+1];return min(st[l][k],st[r-(1<<k)+1][k]);}int main(){//freopen("dna.in","r",stdin);scanf("%d",&tst);log[0]=-1;for(int i=1;i<=N-5;++i) log[i]=log[i>>1]+1;while(tst--){int ans=0;scanf("%s",s+1);len1=strlen(s+1);m='Z'+1;s[len1+1]='Z';scanf("%s",s+len1+2);n=strlen(s+1);len2=strlen(s+len1+2);for(int i=1;i<=m;++i) count[i]=0;for(int i=1;i<=n;++i) count[s[i]]=1;for(int i=1;i<=m;++i) count[i]+=count[i-1];for(int i=1;i<=n;++i) rank[i]=count[s[i]];int k=0;for(int p=1;k!=n;p<<=1,m=k){for(int i=1;i<=m;++i) count[i]=0;for(int i=1;i<=n;++i) count[rank[i+p]]++;for(int i=1;i<=m;++i) count[i]+=count[i-1];for(int i=n;i>=1;--i) tmp[count[rank[i+p]]--]=i;for(int i=1;i<=m;++i) count[i]=0;for(int i=1;i<=n;++i) count[rank[tmp[i]]]++;for(int i=1;i<=m;++i) count[i]+=count[i-1];for(int i=n;i>=1;--i) sa[count[rank[tmp[i]]]--]=tmp[i];memcpy(rank1,rank,sizeof(rank1));rank[sa[1]]=k=1;for(int i=2;i<=n;++i){if(rank1[sa[i]]!=rank1[sa[i-1]]||rank1[sa[i]+p]!=rank1[sa[i-1]+p]) ++k;rank[sa[i]]=k;}}k=0;for(int i=1;i<=n;++i){if(rank[i]==1){h[1]=0;continue;}if(i==1||h[rank[i-1]]<=1) k=0;if(k) --k;while(s[i+k]==s[sa[rank[i]-1]+k]) ++k;h[rank[i]]=k;}for(int i=1;i<=n;++i) st[i][0]=h[i];for(int i=1;i<=log[n];++i)for(int j=1;j<=n;++j)if(j+(1<<i-1)<=n) st[j][i]=min(st[j][i-1],st[j+(1<<i-1)][i-1]);for(int i=1;i<=len1-len2+1;++i){int x=lcp(i,len1+2);if(x>=len2-3){ans++;continue;}++x;x+=lcp(i+x,len1+2+x);if(x>=len2-2){ans++;continue;}++x;x+=lcp(i+x,len1+2+x);if(x>=len2-1){ans++;continue;}++x;x+=lcp(i+x,len1+2+x);if(x>=len2){ans++;continue;}}printf("%d\n",ans);}return 0;}
阅读全文
0 0
- bzoj4829 [TJOI2017]dna(后缀数组)
- UVa760 - DNA Sequencing(后缀数组)
- UVA 题目760 DNA Sequencing (后缀数组求两个串最长公共子串,字典序输出)
- 后缀数组(倍增)
- 后缀数组(详细)
- 后缀数组(详细)
- 后缀数组(转载)
- 后缀数组(不完善)
- 后缀数组(poj3729)
- suffix_array(后缀数组)
- poj3581_Sequence(后缀数组)
- 后缀数组(一)
- bzoj1031(后缀数组)
- POJ1743(后缀数组)
- 后缀数组(一)
- 后缀数组(SA)
- 后缀数组(倍增)
- 后缀数组(UVa11107)
- nyoj-114 某种序列
- failure during conversion to COFF: file invalid or corrupt
- python核心高级学习总结8------动态性、__slots__、生成器、迭代器、装饰、闭包
- URI和URL
- NYOJ 455 黑色帽子
- bzoj4829 [TJOI2017]dna(后缀数组)
- How to excute jquery by nodejs?
- Python 初始化 dict 方法
- 大白话讲解Promise(三)——搞懂jquery中的Promise
- MySQL恢复之update忘加where条件误操作后数据恢复(提前条件binlog为row行格式)
- HDU1021
- css_eg02_css画一个倒三角
- 基于YCbCr的肤色检测OpenCV实现
- spring学习笔记2--bean装配