hdu6704 2019CCPC网络选拔赛1003 K-th occurrence 后缀自动机+线段树合并
解题思路:
fail树上用权值线段树合并求right/endpos集合,再用倍增找到待查询串对应节点,然后权值线段树求第k大。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1e5+5;
int n,q;
char s[maxn];
namespace SegTree{
int sum[maxn*100],L[maxn*100],R[maxn*100];
int tot1;
int update(int rt,int l,int r,int pos,int val){
int nrt=++tot1;
L[nrt]=L[rt]; R[nrt]=R[rt]; sum[nrt]=sum[rt]+val;
if(l!=r){
int mid=(l+r)>>1;
if(pos<=mid)L[nrt]=update(L[rt],l,mid,pos,val);
else R[nrt]=update(R[rt],mid+1,r,pos,val);
}
return nrt;
}
int merge(int rt1,int rt2){
if(!rt1 || !rt2)return rt1|rt2;
int nrt=++tot1;
L[nrt]=L[rt1]; R[nrt]=R[rt1]; sum[nrt]=sum[rt1]+sum[rt2];
L[nrt]=merge(L[rt1],L[rt2]);
R[nrt]=merge(R[rt1],R[rt2]);
return nrt;
}
int query(int rt,int l,int r,int k){
if(l==r)return l;
int mid=(l+r)>>1;
if(k<=sum[L[rt]])return query(L[rt],l,mid,k);
return query(R[rt],mid+1,r,k-sum[L[rt]]);
}
}using namespace SegTree;
namespace Suffix_Automaton{
int ch[maxn<<1][26],fa[maxn<<1],len[maxn<<1];
int last,tot;
int rt[maxn],T[maxn<<1];
int Fa[maxn<<1][20];
inline void init(){
last=tot=1;
len[1]=fa[0]=0;
memset(ch[1],0,sizeof(ch[1]));
T[1]=0;
}
inline int newnode(){
++tot;
len[tot]=fa[tot]=0;
memset(ch[tot],0,sizeof(ch[tot]));
T[tot]=0;
return tot;
}
inline void extend(int c,int right){
int p=last,cur=newnode();
len[cur]=len[last]+1;
last=cur;
rt[right]=cur;
T[cur]=update(T[cur],1,n,right,1);
while(p && !ch[p][c]){
ch[p][c]=cur;
p=fa[p];
}
if(!p)fa[cur]=1;
else{
int q=ch[p][c];
if(len[p]+1==len[q])fa[cur]=q;
else{
int clone=newnode();
len[clone]=len[p]+1;
memcpy(ch[clone],ch[q],sizeof(ch[q]));
fa[clone]=fa[q];
fa[q]=fa[cur]=clone;
while(ch[p][c]==q){
ch[p][c]=clone;
p=fa[p];
}
}
}
}
int c[maxn<<1],A[maxn<<1];
inline void init(char *a,int l){
init();
for(int i=1;i<=l;i++)extend(a[i]-'a',i);
for(int i=0;i<=tot;i++)c[i]=0;
for(int i=1;i<=tot;i++)++c[len[i]];
for(int i=1;i<=tot;i++)c[i]+=c[i-1];
for(int i=1;i<=tot;i++)A[--c[len[i]]]=i;
for(int i=tot-1;i>=1;i--)T[fa[A[i]]]=merge(T[fa[A[i]]],T[A[i]]);
for(int i=1;i<=tot;i++)Fa[i][0]=fa[i];
for(int k=1;k<=19;k++)for(int i=1;i<=tot;i++)Fa[i][k]=Fa[Fa[i][k-1]][k-1];
}
inline void solve(int l,int r,int k){
int u=rt[r],length=r-l+1;
if(len[fa[u]]+1>length){
for(int k=19;k>=0;k--)if(len[fa[Fa[u][k]]]+1>length)u=Fa[u][k];
u=fa[u];
}
if(k<=sum[T[u]])printf("%d\n",query(T[u],1,n,k)-length+1);
else printf("-1\n");
}
}using namespace Suffix_Automaton;
int main()
{
//#ifndef ONLINE_JUDGE
// freopen("in.txt","r",stdin);
//#endif
int T;
scanf("%d",&T);
while(T--){
tot1=0;
scanf("%d %d",&n,&q);
scanf("%s",s+1);
init(s,n);
int l,r,k;
while(q--){
scanf("%d %d %d",&l,&r,&k);
solve(l,r,k);
}
}
return 0;
}
hdu6704 2019CCPC网络选拔赛1003 K-th occurrence 后缀自动机+线段树合并的更多相关文章
- hdu6704 2019CCPC网络选拔赛1003 K-th occurrence 后缀数组
题意:给你一个长度为n的字符串,有q个询问,每次询问一个子串s(l,r)第k次出现的位置,若子串出现次数少于k次输出-1. 解题思路:先把SA跑出来,然后对于每次询问可以由l和rank[]找到l在所有 ...
- HDU - 6704 K-th occurrence (后缀数组+主席树/后缀自动机+线段树合并+倍增)
题意:给你一个长度为n的字符串和m组询问,每组询问给出l,r,k,求s[l,r]的第k次出现的左端点. 解法一: 求出后缀数组,按照排名建主席树,对于每组询问二分或倍增找出主席树上所对应的的左右端点, ...
- HDU-6704 K-th occurrence (后缀自动机father树上倍增建权值线段树合并)
layout: post title: HDU-6704 K-th occurrence (后缀自动机father树上倍增建权值线段树合并) author: "luowentaoaa&quo ...
- BZOJ 3065 带插入区间K小值(sag套线段树)
3065: 带插入区间K小值 Time Limit: 60 Sec Memory Limit: 512 MBSubmit: 4696 Solved: 1527[Submit][Status][Di ...
- [BZOJ 3110] [luogu 3332] [ZJOI 2013]k大数查询(权值线段树套线段树)
[BZOJ 3110] [luogu 3332] [ZJOI 2013]k大数查询(权值线段树套线段树) 题面 原题面有点歧义,不过从样例可以看出来真正的意思 有n个位置,每个位置可以看做一个集合. ...
- SPOJ-COT-Count on a tree(树上路径第K小,可持久化线段树)
题意: 求树上A,B两点路径上第K小的数 分析: 同样是可持久化线段树,只是这一次我们用它来维护树上的信息. 我们之前已经知道,可持久化线段树实际上是维护的一个前缀和,而前缀和不一定要出现在一个线性表 ...
- 树上第k小,可持久化线段树+倍增lca
给定一颗树,树的每个结点都有权值, 有q个询问,每个询问是 u v k ,表示u到v路径上第k小的权值是多少. 每个结点所表示的线段树,是父亲结点的线段树添加该结点的权值之后形成的新的线段树 c[ro ...
- 网络赛 I题 Max answer 单调栈+线段树
题目链接:https://nanti.jisuanke.com/t/38228 题意:在给出的序列里面找一个区间,使区间最小值乘以区间和得到的值最大,输出这个最大值. 思路:我们枚举每一个数字,假设是 ...
- ACM-ICPC 2018 南京赛区网络预赛 G Lpl and Energy-saving Lamps(线段树)
题目链接:https://nanti.jisuanke.com/t/30996 中文题目: 在喝茶的过程中,公主,除其他外,问为什么这样一个善良可爱的龙在城堡里被监禁Lpl?龙神秘地笑了笑,回答说这是 ...
随机推荐
- Android 的重要控件 ListView (听说是最难最常用的控件)
这个打字有点慢了,左手受伤了,不过很幸运,左手小拇指没事(这就可以愉快地使用快捷键啦!),虽然有伤,但还是得坚持总结,不只是为自己,还为未来的你们铺路,希望我写的,对你们有帮助. 提前给自己一个祝福: ...
- Python Tricks —— 使用 pywinrm 远程控制 Windows 主机
启用 WinRM 远程服务: winrm quickconfig 很多人学习python,不知道从何学起.很多人学习python,掌握了基本语法过后,不知道在哪里寻找案例上手.很多已经做案例的人,却不 ...
- three.js 着色器材质之变量(二)
上一篇郭先生在例子中用到了着色器变量中的uniform和varying.这篇继续结合例子将一下attribute变量,在使用过程中也发现由于three.js的版本迭代,之前的一些属性和参数已经发生了改 ...
- JS 获取验证码按钮改变案例
HTML代码 <div class="box"> <label for="">手机号</label> <input t ...
- 从0写一个Golang日志处理包
WHY 日志概述 日志几乎是每个实际的软件项目从开发到最后实际运行过程中都必不可少的东西.它对于查看代码运行流程,记录发生的事情等方面都是很重要的. 一个好的日志系统应当能准确地记录需要记录的信息,同 ...
- DAPP开发初探——永存的留言
转载地址 https://blog.csdn.net/qq_33764491/article/details/80570266 前言 最近DAPP的开发貌似很火,学习了区块链的一些知识之后,相信有很多 ...
- menset与fill
menset函数一般只对int型数组进行0.-1的赋值.原因:menset对数组是按字节赋值,对每个字节的赋值是相同的,故int的4个字节全部被赋相同的值,而0正好二进制编码全为0,-1的二进制编码全 ...
- 【Spring注解驱动开发】如何实现方法、构造器位置的自动装配?我这样回答让面试官很满意!
在 冰河技术 微信公众号前面的文章中,我们介绍了如何使用注解来自动装配Spring组件.之前将的都是在来的字段上添加注解,那有没有什么方法可以实现方法.构造器位置的自动装配吗?今天我们就一起来探讨下如 ...
- 9. Technical Overview 技术概述
9.1 Runtime Environment 运行环境 Spring Security 3.0需要Java 5.0运行时环境或更高版本.由于Spring Security的目标是以独立的方式运行,所 ...
- 笔记:Ubuntu安装LAMP环境
一.更换Ubuntu的镜像源 镜像源路径:/etc/apt/sources.list 备份:cp /etc/apt/sources.list /etc/apt/souces.list.bak 使用阿里 ...