博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
BZOJ1396 : 识别子串
阅读量:5339 次
发布时间:2019-06-15

本文共 1328 字,大约阅读时间需要 4 分钟。

枚举左端点$i$,那么可行的右端点$j$的最小值单调不下降,可以通过双指针求出,检验可以通过在后缀数组里检查相邻height值做到$O(1)$。

那么左端点为$i$,右端点在$[j,n]$,它对前面一段的贡献为定值,对后面一段的贡献为等差数列,线段树维护即可。

时间复杂度$O(n\log n)$。

 

#include
#include
const int N=100010,M=262150;int n,i,j,rk[N],sa[N],height[N],tmp[N],cnt[N],ta[M],tb[M];char s[N];void suffixarray(int n,int m){ int i,j,k;n++; for(i=0;i
=n-1)break; } for(j=rk[height[i=k=0]=0];i
=y-x+1||height[k+1]>=y-x+1;}inline void up(int&a,int b){if(a>b)a=b;}void build(int x,int a,int b){ ta[x]=tb[x]=N; if(a==b)return; int mid=(a+b)>>1; build(x<<1,a,mid),build(x<<1|1,mid+1,b);}void tag(int*v,int x,int a,int b,int c,int d,int p){ if(c<=a&&b<=d){up(v[x],p);return;} int mid=(a+b)>>1; if(c<=mid)tag(v,x<<1,a,mid,c,d,p); if(d>mid)tag(v,x<<1|1,mid+1,b,c,d,p);}void dfs(int x,int a,int b){ if(a==b){ up(ta[x],tb[x]+a); printf("%d\n",ta[x]+1); return; } int mid=(a+b)>>1; up(ta[x<<1],ta[x]),up(ta[x<<1|1],ta[x]); up(tb[x<<1],tb[x]),up(tb[x<<1|1],tb[x]); dfs(x<<1,a,mid),dfs(x<<1|1,mid+1,b);}int main(){ gets(s);n=strlen(s); suffixarray(n,128); build(1,0,n-1); for(i=0;i
=n)break; tag(ta,1,0,n-1,i,j,j-i); tag(tb,1,0,n-1,j,n-1,-i); } dfs(1,0,n-1); return 0;}

  

转载于:https://www.cnblogs.com/clrs97/p/5277707.html

你可能感兴趣的文章
Django之Models
查看>>
CSS 透明度级别 及 背景透明
查看>>
Linux 的 date 日期的使用
查看>>
PHP zip压缩文件及解压
查看>>
SOAP web service用AFNetWorking实现请求
查看>>
Java变量类型,实例变量 与局部变量 静态变量
查看>>
mysql操作命令梳理(4)-中文乱码问题
查看>>
Python环境搭建(安装、验证与卸载)
查看>>
一个.NET通用JSON解析/构建类的实现(c#)
查看>>
Windows Phone开发(5):室内装修 转:http://blog.csdn.net/tcjiaan/article/details/7269014
查看>>
详谈js面向对象 javascript oop,持续更新
查看>>
关于这次软件以及pda终端的培训
查看>>
jQuery上传插件Uploadify 3.2在.NET下的详细例子
查看>>
如何辨别一个程序员的水平高低?是靠发量吗?
查看>>
新手村之循环!循环!循环!
查看>>
正则表达式的用法
查看>>
线程安全问题
查看>>
SSM集成activiti6.0错误集锦(一)
查看>>
下拉刷新
查看>>
linux的子进程调用exec( )系列函数
查看>>