P6640-[BJOI2020]封印【SAM,二分】
正题
题目链接:https://www.luogu.com.cn/problem/P6640
题目大意
给出两个字符串s,ts,ts,t。qqq次给出l,rl,rl,r询问sl∼rs_{l\sim r}sl∼r与ttt的最长公共子串。
解题思路
对于末尾的条件很好做,直接上SAMSAMSAM就好了,设fif_ifi表示一个最大的数满足si−fi+1,is_{i-f_i+1,i}si−fi+1,i是ttt的子串。
然后对于一个询问l,rl,rl,r就是相当于求max{min{fl+i−1,i}}(i∈[1,r−l+1])max\{\ min\{f_{l+i-1},i\}\ \}(i\in[1,r-l+1])max{ min{fl+i−1,i} }(i∈[1,r−l+1])这个东西。
考虑二分一个答案xxx,那么在[l,l+x−2][l,l+x-2][l,l+x−2]这个范围内一定不会有答案,在剩下的[l+x−1,r][l+x-1,r][l+x−1,r]内取一个最大值和xxx做一下对比即可去掉minminmin的条件。
用一个STSTST维护一下就好了,时间复杂度O(n+qlogn)O(n+q\log n)O(n+qlogn)
code
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=4e5+10,T=19;
int n,m,q,last,cnt,lg[N],f[N][T];
int ch[N][26],len[N],fa[N];
char s[N],t[N];
void Insert(int c){int p=last,np=last=++cnt;len[np]=len[p]+1;for(;p&&!ch[p][c];p=fa[p])ch[p][c]=np;if(!p)fa[np]=1;else{int q=ch[p][c];if(len[p]+1==len[q])fa[np]=q;else{int nq=++cnt;len[nq]=len[p]+1;memcpy(ch[nq],ch[q],sizeof(ch[nq]));fa[nq]=fa[q];fa[q]=fa[np]=nq;for(;p&&ch[p][c]==q;p=fa[p])ch[p][c]=nq;}}return;
}
int Ask(int l,int r){int z=lg[r-l+1];return max(f[l][z],f[r-(1<<z)+1][z]);
}
int main()
{scanf("%s",t+1);m=strlen(t+1);scanf("%s",s+1);n=strlen(s+1);scanf("%d",&q);last=cnt=1;for(int i=1;i<=n;i++)Insert(s[i]-'a');for(int i=1,p=1,l=0;i<=m;i++){int c=t[i]-'a';while(!ch[p][c])p=fa[p],l=len[p];if(!p)p=1,l=0;else p=ch[p][c],l++;f[i][0]=l;}for(int i=2;i<=n;i++)lg[i]=lg[i>>1]+1;for(int j=1;(1<<j)<=n;j++)for(int i=1;i+(1<<j)-1<=n;i++)f[i][j]=max(f[i][j-1],f[i+(1<<j-1)][j-1]);while(q--){int L,R;scanf("%d%d",&L,&R);int l=1,r=R-L+1;while(l<=r){int mid=(l+r)>>1;if(Ask(L+mid-1,R)>=mid)l=mid+1;else r=mid-1;}printf("%d\n",r);}return 0;
}
P6640-[BJOI2020]封印【SAM,二分】相关推荐
- [学习笔记]后缀自动机
解决大部分字符串问题的大杀器 给一下clj课件:戳我 SAM是一个博大精深的算法.虽然没有像网络流家族,Tarjan家族一样拉帮结派,但是自身包含很多方法. 一.前言 字符串常见问题:各种匹配 1.L ...
- 基础省选+NOI-第6部分 字符串
1.Kmp KMP算法 KMP算法_哔哩哔哩_bilibili KMP算法计算next函数值(教材版,超简单!) KMP算法计算next函数值(教材版,超简单!)_哔哩哔哩_bilibili KMP算 ...
- Before NOIP 2018
目录 总结 刷题 2018 - 9 - 24 2018 - 9 - 25 2018 - 9 - 26 2018 - 9 - 27 2018 - 9 - 28 2018 - 9 - 29 2018 - ...
- Luogu4022 CTSC2012 熟悉的文章 广义SAM、二分答案、单调队列
传送门 先将所有模板串扔进广义SAM.发现作文的\(L0\)具有单调性,即\(L0\)更小不会影响答案,所以二分答案. 假设当前二分的值为\(mid\),将当前的作文放到广义SAM上匹配. 设对于第\ ...
- [BZOJ]4180: 字符串计数 SAM+矩阵乘法+二分
Description SD有一名神犇叫做Oxer,他觉得字符串的题目都太水了,于是便出了一道题来虐蒟蒻yts1999. 他给出了一个字符串T,字符串T中有且仅有4种字符 'A', 'B', 'C', ...
- UVA 11419 SAM I AM(最大二分匹配最小点覆盖:König定理)
题意:在方格图上打小怪,每次可以清除一整行或一整列的小怪,问最少的步数是多少,又应该在哪些位置操作(对输出顺序没有要求). 分析:最小覆盖问题 这是一种在方格图上建立的模型:令S集表示"行& ...
- BZOJ 2780: [Spoj]8093 Sevenk Love Oimaster( 后缀数组 + 二分 + RMQ + 树状数组 )
全部串起来做SA, 在按字典序排序的后缀中, 包含每个询问串必定是1段连续的区间, 对每个询问串s二分+RMQ求出包含s的区间. 然后就是求区间的不同的数的个数(经典问题), sort queries ...
- BZOJ 2806 Luogu P4022 [CTSC2012]Cheat (广义后缀自动机、DP、二分、单调队列)
题目链接: (bzoj) https://www.lydsy.com/JudgeOnline/problem.php?id=2806 (luogu) https://www.luogu.org/pro ...
- BZOJ.4180.字符串计数(后缀自动机 二分 矩阵快速幂/倍增Floyd)
题目链接 先考虑 假设S确定,使构造S操作次数最小的方案应是:对T建SAM,S在SAM上匹配,如果有S的转移就转移,否则操作数++,回到根节点继续匹配S.即每次操作一定是一次极大匹配. 简单证明:假设 ...
最新文章
- UVA10110 Light, more light
- 南邮java实验一报告_南邮JAVA程序设计实验1 综合图形界面程序设计
- mysql与php驱动,用PHP和MySQL构建一个数据库驱动的网站(四)-PHP教程,PHP基础
- Codeforces 1181
- 虚拟存储体系由()两级存储器构成。
- 【转】一、用于VS2019的QT配置
- 四管前级怎么去掉高低音音调_一些歌曲音调太高怎么才能唱上去??
- drupal_Drupal,物联网和开放式硬件的交集
- bzoj1189 [HNOI2007]紧急疏散EVACUATE spfa+网络流+二分
- GitHub#algorithm#:《剑指offer》 的50道面试题
- 从数据平台到报表开发,我靠这个零代码报表工具,转行后月薪3W
- 引入方式之外部样式表(CSS、HTML)
- 2017 ACM/ICPC Asia Regional Beijing Online 记录
- 没有Angular 3,下一个Angular主版本将是Angular 4
- css line-height多种用法与之间的区别
- GIS实验之房价数据可视化分析
- utc时间 单位换算_utc时间(utc时间转换北京时间)
- 怎样将UltraISO做的启动U盘还原成原来的样子
- 迭代器 iter()函数与next()函数 for...in...循环的本质
- 操作的基本原则,每日必读【不断更新中】
热门文章
- java跨平台的特性_【简答题】什么是跨平台特性?Java怎样实现跨平台特性?
- python程序设计实践教程陈东_Python
- latex 1图加标题_使用VsCode编译latex心得
- html 存储登录状态,Vue中保存用户登录状态实例代码
- 罗马数字转换成数字java_C趣味编程百例(31)将阿拉伯数字转换为罗马数字
- c语言抬起按键,c语言获得键盘的按键
- python网页动图_python,tensorflow线性回归Django网页显示Gif动态图
- 7-3 最小生成树-kruskal (10 分)(思路+详解+并查集详解+段错误超时解决)宝 Come
- [蓝桥杯2016决赛]七星填数-next_permutation枚举
- alibaba实习生代码大赛