题目大意:求两个字符串的最长公共子串长度

把两个串接在一起,中间放一个#,然后求出height

接下来还是老套路,二分出一个答案ans,然后去验证,如果有连续几个位置的h[i]>=ans,且存在sa[i]的最大值在第二个串里,最小值在第一个串里,说明答案成立

别再把后缀数组敲错了

  1. #include <cmath>
  2. #include <cstdio>
  3. #include <cstring>
  4. #include <algorithm>
  5. #define ll long long
  6. #define ull unsigned long long
  7. #define inf 0x3f3f3f3f
  8. #define N 205000
  9. using namespace std;
  10. //re
  11. int n,l1,l2,len;
  12. char s1[N],s2[N],str[N];
  13. int tr[N],rk[N],sa[N],hs[N],h[N];
  14. int gint()
  15. {
  16. int rett=,fh=;char c=getchar();
  17. while(c<''||c>''){if(c=='-')fh=-;c=getchar();}
  18. while(c>=''&&c<=''){rett=(rett<<)+(rett<<)+c-'';c=getchar();}
  19. return rett*fh;
  20. }
  21. bool check(int k,int x,int y){
  22. if(x+k>len||y+k>len) return ;
  23. else return (rk[x]==rk[y]&&rk[x+k]==rk[y+k])?:;
  24. }
  25. void get_sa()
  26. {
  27. int cnt=,i;
  28. for(i=;i<=len;i++) hs[str[i]]++;
  29. for(i=;i<=;i++) if(hs[i]) tr[i]=++cnt;
  30. for(i=;i<=;i++) hs[i]+=hs[i-];
  31. for(i=;i<=len;i++) rk[i]=tr[str[i]],sa[hs[str[i]]--]=i;
  32. for(int k=;cnt<len;k<<=)
  33. {
  34. for(i=;i<=cnt;i++) hs[i]=;
  35. for(i=;i<=len;i++) hs[rk[i]]++;
  36. for(i=;i<=cnt;i++) hs[i]+=hs[i-];
  37. for(i=len;i>=;i--) if(sa[i]>k) tr[sa[i]-k]=hs[rk[sa[i]-k]]--;
  38. for(i=;i<=k;i++) tr[len-i+]=hs[rk[len-i+]]--;
  39. for(i=;i<=len;i++) sa[tr[i]]=i;
  40. for(i=,cnt=;i<=len;i++) tr[sa[i]]=check(k,sa[i],sa[i-])?cnt:++cnt;
  41. for(i=;i<=len;i++) rk[i]=tr[i];
  42. }
  43. }
  44. void get_height()
  45. {
  46. for(int i=;i<=len;i++){
  47. if(rk[i]==) continue;
  48. for(int j=max(,h[rk[i-]]-);;j++)
  49. if(str[i+j-]==str[sa[rk[i]-]+j-]) h[rk[i]]=j;
  50. else break;
  51. }
  52. }
  53. int check(int ans)
  54. {
  55. for(int i=;i<=len;){
  56. if(h[i]<ans){i++;continue;}
  57. int mi=sa[i-],ma=sa[i-];
  58. for(;i<=len&&h[i]>=ans;i++)
  59. mi=min(mi,sa[i]),ma=max(ma,sa[i]);
  60. if(mi<=l1&&ma>l1+) return ;
  61. }return ;
  62. }
  63.  
  64. int main()
  65. {
  66. scanf("%s",s1+);
  67. scanf("%s",s2+);
  68. l1=strlen(s1+),l2=strlen(s2+);
  69. for(int i=;i<=l1;i++)
  70. str[i]=s1[i];
  71. str[l1+]='#';
  72. for(int i=;i<=l2;i++)
  73. str[i+l1+]=s2[i];
  74. len=l1+l2+;
  75. get_sa();
  76. get_height();
  77. int l=,r=min(l1,l2),ans;
  78. while(l<=r){
  79. int mid=(l+r)>>;
  80. if(check(mid)) ans=mid,l=mid+;
  81. else r=mid-;
  82. }
  83. printf("%d\n",ans);
  84. return ;
  85. }

POJ 2774 Long Long Message (后缀数组+二分)的更多相关文章

  1. poj 2774 Long Long Message 后缀数组基础题

    Time Limit: 4000MS   Memory Limit: 131072K Total Submissions: 24756   Accepted: 10130 Case Time Limi ...

  2. POJ 2774 Long Long Message 后缀数组

    Long Long Message   Description The little cat is majoring in physics in the capital of Byterland. A ...

  3. poj 2774 Long Long Message 后缀数组LCP理解

    题目链接 题意:给两个长度不超过1e5的字符串,问两个字符串的连续公共子串最大长度为多少? 思路:两个字符串连接之后直接后缀数组+LCP,在height中找出max同时满足一左一右即可: #inclu ...

  4. POJ 2774 Long Long Message 后缀数组模板题

    题意 给定字符串A.B,求其最长公共子串 后缀数组模板题,求出height数组,判断sa[i]与sa[i-1]是否分属字符串A.B,统计答案即可. #include <cstdio> #i ...

  5. POJ - 2774 Long Long Message (后缀数组/后缀自动机模板题)

    后缀数组: #include<cstdio> #include<algorithm> #include<cstring> #include<vector> ...

  6. POJ 2774 Long Long Message ——后缀数组

    [题目分析] 用height数组RMQ的性质去求最长的公共子串. 要求sa[i]和sa[i-1]必须在两个串中,然后取height的MAX. 利用中间的字符来连接两个字符串的思想很巧妙,记得最后还需要 ...

  7. POJ 3261 Milk Patterns(后缀数组+二分答案)

    [题目链接] http://poj.org/problem?id=3261 [题目大意] 求最长可允许重叠的出现次数不小于k的子串. [题解] 对原串做一遍后缀数组,二分子串长度x,将前缀相同长度超过 ...

  8. POJ 3294 Life Forms(后缀数组+二分答案)

    [题目链接] http://poj.org/problem?id=3294 [题目大意] 求出在至少在一半字符串中出现的最长子串. 如果有多个符合的答案,请按照字典序输出. [题解] 将所有的字符串通 ...

  9. POJ 1743 Musical Theme(后缀数组+二分答案)

    [题目链接] http://poj.org/problem?id=1743 [题目大意] 给出一首曲子的曲谱,上面的音符用不大于88的数字表示, 现在请你确定它主旋律的长度,主旋律指的是出现超过一次, ...

随机推荐

  1. 如何构建通用 api 中间层

    零.问题的由来 开门见山地说,这篇文章是一篇安利软文~,安利的对象就是最近搞的 tua-api. 顾名思义,这就是一款辅助获取接口数据的工具. 发请求相关的工具辣么多,那我为啥要用你呢? 理想状态下, ...

  2. vue 2.0 购物车小球抛物线

    备注:此项目模仿 饿了吗.我用的是最新的Vue, 视频上的一些写法已经被废弃了. 布局代码 <div class="ball-container"> <trans ...

  3. Python学习笔记(1)Pycharm基本认识

    Pycharm的使用 集成开发环境(IDE,集成开发环境) VIM  #经典的linux下的文本编辑器 Pycharm,主要用于python开发的ide PyCharm常用快捷键 Ctrl + / # ...

  4. 【Codeforces Round #502 (in memory of Leopoldo Taravilse, Div. 1 + Div. 2) D】The Wu

    [链接] 我是链接,点我呀:) [题意] 给你n个字符串放在multiset中. 这些字符串都是长度为m的01串. 然后给你q个询问 s,k 问你set中存在多少个字符串t 使得∑(t[i]==s[i ...

  5. BA-深化设计流程

    本文分三个层次描述了常规BA系统深化设计应包含了步骤和文件,第一个.第二个为转载别人的步骤,第三个为本人写的步骤. 深化系统大体步骤: 1.认真阅读招标文件,明确招标方需求. 2.仔细查阅图纸,确定被 ...

  6. six.moves的用法

    six是用来兼容python 2 和 3的,我猜名字就是用的2和3的最小公倍数. six.moves 是用来处理那些在2 和 3里面函数的位置有变化的,直接用six.moves就可以屏蔽掉这些变化 P ...

  7. tomcat内存大小设置

    tomcat内存大小设置 如果安装为windows服务,需要进行内存设置的时候,选择configure...界面, 在Java Tab页面内可以进行内存参数的设置. 学习了:http://elf884 ...

  8. Routh-Hurwitz Criterion 劳斯稳定判据

    Routh-Hurwitz Criterion 为什么仅仅要有一个极点在右半平面,那么系统就不会稳定? 比如H(s) =( 1/(s+1) ) *  ( 1/(s+3) ) * ( 1/(s-2) ) ...

  9. HTML5实现歌词同步

    开篇 HTML5的最强大之处莫过于对媒体文件的处理,如利用一个简单的vedio标签就能够实现视频播放.相似地,在HTML5中也有相应的处理音频文件的标签,那就是audio标签 在线Demo audio ...

  10. DNS递归查询和迭代查询的差别

    转载地址:http://blog.csdn.net/wuchuanpingstone/article/details/6720723 递归查询和迭代查询的差别 (1)递归查询 递归查询是一种DNS s ...