根据height数组的定义,和当前后缀串i最长的相同串的长度就是max(height[i],height[i+1]),这个后缀贡献的最短不同串长度就是len=max(height[i],height[i+1])

考虑这个贡献的范围,对于ii+len,贡献是len,对于w=i+len+1n可以更新w-len+1,也就是这长为len的串再向右延伸,可以证明这也是不重子串

这两个东西用两个线段树维护,w-len+1这个直接把-len+1打在线段树上,查询的时候再加回w即可

  1. #include<iostream>
  2. #include<cstdio>
  3. #include<cstring>
  4. using namespace std;
  5. const int N=500005;
  6. int n,wa[N],wb[N],wv[N],wsu[N],sa[N],rk[N],he[N];
  7. char s[N];
  8. struct qwe
  9. {
  10. int l,r,mn[2],lz[2];
  11. }t[N<<2];
  12. bool cmp(int r[],int a,int b,int l)
  13. {
  14. return r[a]==r[b]&&r[a+l]==r[b+l];
  15. }
  16. void saa(char r[],int n,int m)
  17. {
  18. int *x=wa,*y=wb;
  19. for(int i=0;i<=m;i++)
  20. wsu[i]=0;
  21. for(int i=1;i<=n;i++)
  22. wsu[x[i]=r[i]]++;
  23. for(int i=1;i<=m;i++)
  24. wsu[i]+=wsu[i-1];
  25. for(int i=n;i>=1;i--)
  26. sa[wsu[x[i]]--]=i;
  27. for(int j=1,p=1;j<=n&&p<n;j<<=1,m=p)
  28. {
  29. p=0;
  30. for(int i=n-j+1;i<=n;i++)
  31. y[++p]=i;
  32. for(int i=1;i<=n;i++)
  33. if(sa[i]>j)
  34. y[++p]=sa[i]-j;
  35. for(int i=1;i<=n;i++)
  36. wv[i]=x[y[i]];
  37. for(int i=0;i<=m;i++)
  38. wsu[i]=0;
  39. for(int i=1;i<=n;i++)
  40. wsu[wv[i]]++;
  41. for(int i=1;i<=m;i++)
  42. wsu[i]+=wsu[i-1];
  43. for(int i=n;i>=1;i--)
  44. sa[wsu[wv[i]]--]=y[i];
  45. swap(x,y);
  46. x[sa[1]]=1;
  47. p=1;
  48. for(int i=2;i<=n;i++)
  49. x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p:++p;
  50. }
  51. for(int i=1;i<=n;i++)
  52. rk[sa[i]]=i;
  53. for(int i=1,j,k=0;i<=n;he[rk[i++]]=k)
  54. for(k?k--:0,j=sa[rk[i]-1];r[i+k]==r[j+k];k++);
  55. }
  56. void pd(int w,int ro)
  57. {
  58. if(t[ro].lz[w]!=1e9)
  59. {
  60. t[ro<<1].lz[w]=min(t[ro<<1].lz[w],t[ro].lz[w]);
  61. t[ro<<1].mn[w]=min(t[ro<<1].mn[w],t[ro].lz[w]);
  62. t[ro<<1|1].lz[w]=min(t[ro<<1|1].lz[w],t[ro].lz[w]);
  63. t[ro<<1|1].mn[w]=min(t[ro<<1|1].mn[w],t[ro].lz[w]);
  64. t[ro].lz[w]=1e9;
  65. }
  66. }
  67. void build(int ro,int l,int r)
  68. {
  69. t[ro].l=l,t[ro].r=r;
  70. t[ro].mn[0]=t[ro].mn[1]=1e9;
  71. t[ro].lz[0]=t[ro].lz[1]=1e9;
  72. if(l==r)
  73. return;
  74. int mid=(l+r)>>1;
  75. build(ro<<1,l,mid);
  76. build(ro<<1|1,mid+1,r);
  77. }
  78. void update(int w,int ro,int l,int r,int v)
  79. {
  80. if(t[ro].l==l&&t[ro].r==r)
  81. {
  82. t[ro].lz[w]=min(t[ro].lz[w],v);
  83. t[ro].mn[w]=min(t[ro].mn[w],v);
  84. return;
  85. }
  86. pd(w,ro);
  87. int mid=(t[ro].l+t[ro].r)>>1;
  88. if(r<=mid)
  89. update(w,ro<<1,l,r,v);
  90. else if(l>mid)
  91. update(w,ro<<1|1,l,r,v);
  92. else
  93. update(w,ro<<1,l,mid,v),update(w,ro<<1|1,mid+1,r,v);
  94. t[ro].mn[w]=min(t[ro<<1].mn[w],t[ro<<1|1].mn[w]);
  95. }
  96. int ques(int ro,int p)
  97. {
  98. if(t[ro].l==t[ro].r)
  99. return min(t[ro].mn[0],t[ro].mn[1]+p);
  100. pd(0,ro);
  101. pd(1,ro);
  102. int mid=(t[ro].l+t[ro].r)>>1;
  103. if(p<=mid)
  104. return ques(ro<<1,p);
  105. else
  106. return ques(ro<<1|1,p);
  107. }
  108. int main()
  109. {
  110. scanf("%s",s+1);
  111. n=strlen(s+1);
  112. saa(s,n,200);
  113. build(1,1,n);
  114. for(int i=1;i<=n;i++)
  115. {
  116. int len=max(he[rk[i]],he[rk[i]+1]);
  117. if(i+len<=n)
  118. update(0,1,i,i+len,len+1);
  119. if(i+len<n)
  120. update(1,1,i+len+1,n,1-i);
  121. }
  122. for(int i=1;i<=n;i++)
  123. printf("%d\n",ques(1,i));
  124. return 0;
  125. }

bzoj 1396: 识别子串 && bzoj 2865: 字符串识别【后缀数组+线段树】的更多相关文章

  1. BZOJ 5496: [2019省队联测]字符串问题 (后缀数组+主席树优化建图+拓扑排序)

    题意 略 分析 考场上写了暴力建图40分溜了-(结果只得了30分) 然后只要优化建边就行了 首先给出的支配关系无法优化,就直接A向它支配的B连边. 考虑B向以B作为前缀的所有A连边,做一遍后缀数组,两 ...

  2. BZOJ1396 识别子串 和 BZOJ2865 字符串识别

    字符串识别 2865: 字符串识别 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 839  Solved: 261[Submit][Status][D ...

  3. 【BZOJ4556】[TJOI2016&HEOI2016] 字符串(后缀自动机+线段树合并+二分)

    点此看题面 大致题意: 给你一个字符串\(s\),每次问你一个子串\(s[a..b]\)的所有子串和\(s[c..d]\)的最长公共前缀. 二分 首先我们可以发现一个简单性质,即要求最长公共前缀,则我 ...

  4. BZOJ 1396: 识别子串( 后缀数组 + 线段树 )

    这道题各位大神好像都是用后缀自动机做的?.....蒟蒻就秀秀智商写一写后缀数组解法..... 求出Height数组后, 我们枚举每一位当做子串的开头. 如上图(x, y是height值), Heigh ...

  5. BZOJ 2865 字符串识别 | 后缀数组 线段树

    集训讲字符串的时候我唯一想出正解的题-- 链接 BZOJ 2865 题面 给出一个长度为n (n <= 5e5) 的字符串,对于每一位,求包含该位的.最短的.在原串中只出现过一次的子串. 题解 ...

  6. BZOJ.1396.识别子串(后缀自动机/后缀数组 线段树)

    题目链接 SAM:能成为识别子串的只有那些|right|=1的节点代表的串. 设这个节点对应原串的右端点为r[i],则如果|right[i]|=1,即\(s[\ [r_i-len_i+1,r_i-le ...

  7. BZOJ 2865 字符串识别(后缀数组+线段树)

    很容易想到只考虑后缀长度必须为\(max(height[rk[i]],height[rk[i]+1])+1\)(即\([i,i+x-1]\)代表的串只出现过一次)然后我正着做一遍反着做一遍,再取一个\ ...

  8. BZOJ4556 [Tjoi2016&Heoi2016]字符串 【后缀数组 + 主席树 + 二分 + ST表】

    题目 佳媛姐姐过生日的时候,她的小伙伴从某东上买了一个生日礼物.生日礼物放在一个神奇的箱子中.箱子外边写了 一个长为n的字符串s,和m个问题.佳媛姐姐必须正确回答这m个问题,才能打开箱子拿到礼物,升职 ...

  9. 【BZOJ4556】字符串(后缀数组,主席树)

    [BZOJ4556]字符串(后缀数组,主席树) 题面 BZOJ 题解 注意看题: 要求的是\([a,b]\)的子串和[c,d]的\(lcp\)的最大值 先来一下暴力吧 求出\(SA\)之后 暴力枚举\ ...

随机推荐

  1. winform窗体继承泛型类时,设计器无法使用解决办法

    当我们使用winform程序时,winform窗体程序和控件程序都是可以通过设计器进行控件拖拽的,但如果我们继承了一个带有泛型的form基类.那么设计器是无法使用的. 目前我没有找到根本解决的办法,但 ...

  2. CloudStack系统虚拟机启动但不在Host表中

    有网友问到这个问题,CloudStack中,系统虚拟机正常启动,Running状态. 可是在host表中没有对应项,上传下载模板等功能也不正常. 原因:系统虚拟机启动之后,会通过管理网段主动连接man ...

  3. 当电视沦为“情怀”,5G能不能拯救它?(zz)

    文|佘凯文 来源|智能相对论(aixdlun) 现阶段,智能家居行业极度期待5G的到来,甚至超过手机.行业对于颠覆性的升级的欲望极其强烈,纵观整个智能家居行业,除了像智能音箱外的偶尔单品能够“引爆”市 ...

  4. C++ Primer 学习笔记与思考_7 void和void*指针的使用方法

    (一)void的含义 void的字面意思是"无类型",void差点儿仅仅有"凝视"和限制程序的作用,由于从来没有人会定义一个void变量,让我们试着来定义: v ...

  5. GIF Movie Gear逆向实战+注册代码+补丁

    GIF Movie Gear逆向实战+注册代码+补丁 准备 我是在windows 8.1 x64上进行的操作.有不足之处,还望大虾指出. 获取资源 网站下载:http://www.gamani.com ...

  6. (转) Universal-Image-Loader使用大全(史上最屌)

    转载自http://blog.csdn.net/zenjj11/article/details/38728481 项目介绍: Android上最让人头疼的莫过于从网络获取图片.显示.回收,不论什么一个 ...

  7. 关于用HOOK拦截键盘的一些问题

    因为MSDN上说要这样做,所以我就这样做的,读懂MSDN是关键,下面来仔细阅读一下MSDN,看它到底是怎样描述的.阅读的时候我先给出原文,再进行自己的一些翻译或描述. 先看回调函数KeyboardPr ...

  8. spark通信原理

    https://github.com/apache/spark/tree/master/core/src/main/scala/org/apache/spark/network https://git ...

  9. sqlldr trailing nullcols

    由于要导入到tmp_content表的一些列(列:要导入的源文件txt or csv文件)为空,也,按理讲我当时另存为(在windows处理)csv,以,分隔.就是这个列没有内容,也该显示, ,之类的 ...

  10. Vue.js新手入门指南

    最近在逛各大网站,论坛,以及像SegmentFault等编程问答社区,发现Vue.js异常火爆,重复性的提问和内容也很多,楼主自己也趁着这个大前端的热潮,着手学习了一段时间的Vue.js,目前用它正在 ...