洛谷P10716【MX-X1-T4】「KDOI-05」简单的字符串问题(扩展kmp+set+二分+扫描线树状数组)

题目

思路来源

小羊肖恩

题解

羊神这个做法tql,当时只是机械地写,过了之后再想想,才觉得确实是nb

先扩展kmp(Z函数)预处理出来数组,记z[i]为i往后可以和前缀匹配的最大长度

对于每个询问(p,cnt),先离线归位到每个p所在的vector内,

然后扫描线,每个i有一个[i,i+z[i]-1]的区间,在这个区间内以j结尾时,和前缀的lcp固定为[i,j]

对于每个i,先特判只出现一次的情况

对于大于一次的情况,先二分当前最大的可行长度x,

这个需要预处理mn[i][j]表示前i个字符出现j次时的结尾处的最小下标,没有的话就是n+1

对于位置p,出现次数cnt,二分找到满足mn[x][cnt]<=p的最大x,

长度x的前缀出现了cnt次,那么短于长度x的前缀一定出现了cnt次,

然后[i,i+z[i]-1]的扫描线保证了以p结尾的一定是一个当前合法的A串,

这就保证了A在开头和结尾均出现,且总的出现次数不少于cnt次了

一开始暴力双指针预处理mn超时了,后来改成在set上二分就ac了

因为长的串对应的位置一定包含短串,所以可以按zi从大到小释放位置,放进set,在set里二分

代码

#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(int i=(a);i<=(b);++i)
#define per(i,a,b) for(int i=(a);i>=(b);--i)
typedef long long ll;
typedef double db;
typedef pair<int,int> P;
#define fi first
#define se second
#define pb push_back
#define dbg(x) cerr<<(#x)<<":"<<x<<" ";
#define dbg2(x) cerr<<(#x)<<":"<<x<<endl;
#define SZ(a) (int)(a.size())
#define sci(a) scanf("%d",&(a))
#define scll(a) scanf("%lld",&(a))
#define pt(a) printf("%d",a);
#define pte(a) printf("%d\n",a)
#define ptlle(a) printf("%lld\n",a)
#define debug(...) fprintf(stderr, __VA_ARGS__)
const int N=2e5+10;
struct BitPre{ // 求前缀和(可改为max等)int n,tr[N];void init(int _n){n=_n;memset(tr,0,(n+1)*sizeof(*tr));}void add(int x,int v){++x;for(int i=x;i<=n;i+=i&-i)tr[i]+=v;}int sum(int x){++x;int ans=0; for(int i=x;i;i-=i&-i)ans+=tr[i];return ans;}
}tr;
int n,m,net[N],ans[N];
char s[N];
P ask[N];
vector<int>q[N],add[N],del[N],mn[N];
void extkmppre(char s[],int len){int i=0,j,pos;net[0]=len;while(i+1<len&&s[i]==s[i+1])i++;net[1]=i,pos=1;rep(i,2,len-1){if(net[i-pos]+i<net[pos]+pos){net[i]=net[i-pos];}else{j=net[pos]+pos-i;if(j<0)j=0;while(i+j<len&&s[j]==s[i+j])j++;net[i]=j,pos=i;}}
}
void init(){set<int>S;vector<vector<int> >in(n+1,vector<int>());rep(i,0,n-1){in[net[i]].pb(i);}per(i,n,1){for(auto &x:in[i]){S.insert(x);}int up=n/i+1;mn[i].resize(up);mn[i][0]=0;mn[i][1]=i-1;//printf("i:%d up:%d\n",i,up);rep(j,2,up-1){if(mn[i][j-1]>n){mn[i][j]=n+1;continue;}auto it=S.upper_bound(mn[i][j-1]);if(it==S.end()){mn[i][j]=n+1;}else{mn[i][j]=(*it)+i-1;}//printf("i:%d j:%d mn:%d\n",i,j,mn[i][j]);}}rep(i,0,n-1){add[i].pb(i);del[i+net[i]-1].pb(i);}
}
bool ok(int len,int cnt,int p){if(!len)return 1;if(1ll*len*cnt>p+1)return 0;return mn[len][cnt]<=p;//printf("len:%d ok:%1d\n",len,o);//"cnt:%d p:%d mn:%d net:%d\n",len,cnt,p,mn[len][cnt],net[p-len+1]);//return o;
}
int main(){sci(n);scanf("%s",s);extkmppre(s,n);// rep(i,0,n-1){// printf("i:%d z[i]:%d\n",i,net[i]);// }init();tr.init(n+1);sci(m);rep(i,1,m){sci(ask[i].fi);sci(ask[i].se);ask[i].fi--;q[ask[i].fi].pb(i);}rep(i,0,n-1){for(auto &x:add[i]){tr.add(x,1);}for(auto &x:q[i]){int p=ask[x].fi,c=ask[x].se;if(c==1){ans[x]=1;continue;}int l=0,r=n;while(l<=r){int mid=(l+r)/2;if(ok(mid,c,i))l=mid+1;else r=mid-1;}//printf("id:%d p:%d c:%d len:%d\n",x,p,c,r);if(r==0)ans[x]=0;else ans[x]=tr.sum(i)-tr.sum(i-r);// rep(j,1,i){// printf("%d ",tr.sum(j)-tr.sum(j-1));// }// puts("");}for(auto &x:del[i]){tr.add(x,-1);}}rep(i,1,m){printf("%d\n",ans[i]);}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/bicheng/44420.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

隧道转发:保护你的数据传输安全

你曾经是否担心过你的数据安全&#xff1f;现代网络威胁不断增加&#xff0c;保护数据传输的安全性变得愈发重要。今天为大家介绍一种强大的网络传输技术——隧道转发&#xff08;Tunneling&#xff09;。这项技术不仅能有效保障数据安全&#xff0c;还能保护用户隐私&#xff…

centOS79中安装nginx12.15

##red## &#x1f534; 大家好&#xff0c;我是雄雄&#xff0c;欢迎关注微信公众号&#xff0c;雄雄的小课堂。 前言 装了这么多&#xff0c;发现Nginx是最简单的&#xff0c;一次性就搞定了。下面我们来看看如何安装 安装Nginx 安装gcc-c编译器 分开运行&#xff1a; yum…

anaconda安装pytorch

&#x1f4da;博客主页&#xff1a;knighthood2001 ✨公众号&#xff1a;认知up吧 &#xff08;目前正在带领大家一起提升认知&#xff0c;感兴趣可以来围观一下&#xff09; &#x1f383;知识星球&#xff1a;【认知up吧|成长|副业】介绍 ❤️如遇文章付费&#xff0c;可先看…

python爬虫入门(三)之HTML网页结构

一、什么是HTML 1、网页的三大技术要素&#xff1a; HTML定义网页的结构和信息&#xff08;骨架血肉&#xff09;CSS定义网页的样式&#xff08;衣服&#xff09;JavaScript定义用户和网页的交互逻辑&#xff08;动作&#xff09; 2、一个最简单的HTML&#xff1a;用<>…

Qt开发 | qss介绍及控件应用 | qss加载方式 | 控件提升 | 鼠标位置与控件位置 | 搜索编辑框 | tab在左文本水平的tabWidget

文章目录 一、qss简介与应用二、QLineEdit qss介绍与使用三、QPushButton qss1.常用qss1.1 基本样式表1.2 背景图片1.3 图片在左文字在右 2.点击按钮弹出菜单以及右侧箭头样式设置3.鼠标悬浮按钮弹出对话框 四、QCheckBox qss妙用&#xff1a;实时打开关闭状态按钮五、QComboBo…

丑数问题,力扣264,坑点

丑数问题&#xff0c;力扣264&#xff0c;坑点 力扣链接 给你一个整数 n &#xff0c;请你找出并返回第 n 个 丑数 。 丑数 就是质因子只包含 2、3 和 5 的正整数。 示例 1&#xff1a; 输入&#xff1a;n 10 输出&#xff1a;12 解释&#xff1a;[1, 2, 3, 4, 5, 6, 8, 9, …

Shell选择结构

文章目录 一、条件判断二、实例 一、条件判断 格式1&#xff1a;test 条件表达式 格式2&#xff1a;[条件表达式] 格式3&#xff1a;[[条件表达式]] -e 是否存在&#xff0c;不管是文件还是目录&#xff0c;只要存在&#xff0c;条件就成立。 -f 是否为普通文件 -d 是否为目录…

《昇思25天学习打卡营第01天|qingyun201003》

打卡 日期 心得 我的主语言并不是Python,以及现在从事的工作也并不是开发&#xff1b;所以对于这个系列的课程&#xff0c;学习起来是较为困难的&#xff0c;所以基于这种情况&#xff0c;该如何进行学习&#xff1f;我的做法是全部交给AI&#xff0c;使用AI一步步解析代码&a…

java.lang.NullPointerException: null cannot be cast to non-null type kotlin.Int

java.lang.NullPointerException: null cannot be cast to non-null type kotlin.Int fun main(args: Array<String>) {var any1: Any?any1 nullval n1 any1 as? Int ?: -2024println(n1)kotlin.runCatching {var any2: Any?any2 nullval n2 any2 as Intprintln(…

python xlsx 导出表格超链接

该Python脚本用于从Excel文件中的第一列提取所有超链接并保存到一个文本文件中。首先&#xff0c;脚本导入必要的库并定义输入和输出文件的路径。然后&#xff0c;它确保输出文件的目录存在。接着&#xff0c;脚本加载Excel文件并选择活动工作表。通过遍历第一列的所有单元格&a…

Internet Download Manager6.42最新下载器互联网冲浪小能手们!

今天我要来种草一个超级棒的宝贝——Internet Download Manager&#xff08;简称 IDM&#xff09;。这个小家伙简直是下载界的“速度与激情”代言人&#xff0c;让我彻底告别了等待的日子。&#x1f389; IDM马丁正版下载如下: https://wm.makeding.com/iclk/?zoneid34275 …

【创作纪念日】我的三周年创作纪念日

光阴荏苒&#xff0c;从我开始在CSDN写作已经整整3年(2024.7.10)&#xff0c;3年前的今天2021.7.10开始记录第一篇文章&#xff0c;到今天的305篇文章&#xff01;没想到&#xff01;积累的力量&#xff01;继续努力&#xff01; 机缘 1、创作初心&#xff0c;对自己工作的经…

299k stars利用Public APIs提升开发效率:探索APILayer提供的开源资源

299k stars利用Public APIs提升开发效率&#xff1a;探索APILayer提供的开源资源 在现代软件开发中&#xff0c;API&#xff08;应用程序接口&#xff09;是实现应用间通信和功能扩展的关键工具。公共API&#xff08;Public APIs&#xff09;则为开发者提供了宝贵的资源&#…

路由器中 RIB 与 FIB 的区别

注&#xff1a; 本文原是前面发的一篇文章的中间注释&#xff0c;但注释跨行太多&#xff0c;影响原文连贯阅读。 RIB&#xff08;Routing Information Base&#xff0c;路由表&#xff09; RIB 存储所有的路由信息&#xff0c;与具体的路由协议无关。所有的路由协议都在这里…

昇思25天学习打卡营第15天|基于 MindSpore 实现 BERT 对话情绪识别

文章目录 昇思MindSpore应用实践1、基于 MindSpore 实现 BERT 对话情绪识别BERT 模型简介数据集数据加载和数据预处理 2、模型训练模型验证 3、模型推理 Reference 昇思MindSpore应用实践 本系列文章主要用于记录昇思25天学习打卡营的学习心得。 1、基于 MindSpore 实现 BERT…

ArduPilot开源代码之AP_OpticalFlow_CXOF

ArduPilot开源代码之AP_OpticalFlow_CXOF 1. 源由2. Library设计3. 重要例程3.1 AP_OpticalFlow_CXOF::init3.2 AP_OpticalFlow_CXOF::update3.3 AP_OpticalFlow_CXOF::detect 4. 总结5. 参考资料 1. 源由 AP_OpticalFlow_CXOF是就是一个光流计&#xff0c;与前面传感模块&…

【网络】为什么SCTP四次握手可以抵御SYN攻击

深入理解SCTP的安全性&#xff1a;从四次握手到抵御SYN攻击 引言 在网络通信的世界中&#xff0c;安全性和可靠性是至关重要的。传统的TCP&#xff08;传输控制协议&#xff09;在建立连接时使用三次握手&#xff0c;但这种机制存在一些安全漏洞&#xff0c;比如SYN攻击。而S…

解决IDEA每次新建项目都需要重新配置maven的问题

每次打开IDEA都要重新配置maven&#xff0c;这是因为在DEA中分为项目设置和全局设置&#xff0c;这个时候我们就需要去到全局中设置maven了。我用的是IntelliJ IDEA 2023.3.4 (Ultimate Edition)&#xff0c;以此为例。 第一步&#xff1a;打开一个空的IDEA&#xff0c;选择左…

数据结构day6链式队列

主程序 #include "fun.h" int main(int argc, const char *argv[]) { que_p Qcreate(); enqueue(Q,10); enqueue(Q,20); enqueue(Q,30); enqueue(Q,40); enqueue(Q,50); show_que(Q); dequeue(Q); show_que(Q); printf(&qu…

stm32按键设置闹钟数进退位不正常?如何解决

&#x1f3c6;本文收录于《CSDN问答解惑-专业版》专栏&#xff0c;主要记录项目实战过程中的Bug之前因后果及提供真实有效的解决方案&#xff0c;希望能够助你一臂之力&#xff0c;帮你早日登顶实现财富自由&#x1f680;&#xff1b;同时&#xff0c;欢迎大家关注&&收…