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?龙神秘地笑了笑,回答说这是 ...
随机推荐
- 009_go语言中的slices分片
代码演示 package main import "fmt" func main() { s := make([]string, 3) fmt.Println("emp: ...
- Kubernetes Pod OOM 排查日记
一.发现问题 在一次系统上线后,我们发现某几个节点在长时间运行后会出现内存持续飙升的问题,导致的结果就是Kubernetes集群的这个节点会把所在的Pod进行驱逐OOM:如果调度到同样问题的节点上,也 ...
- [机器学习] keras:MNIST手写数字体识别(DeepLearning 的 HelloWord程序)
深度学习界的Hello Word程序:MNIST手写数字体识别 learn from(仍然是李宏毅老师<机器学习>课程):http://speech.ee.ntu.edu.tw/~tlka ...
- 主成分分析PCA数据降维原理及python应用(葡萄酒案例分析)
目录 主成分分析(PCA)——以葡萄酒数据集分类为例 1.认识PCA (1)简介 (2)方法步骤 2.提取主成分 3.主成分方差可视化 4.特征变换 5.数据分类结果 6.完整代码 总结: 1.认识P ...
- Vue 倒计时组件封装
项目中需要用到倒计时的功能,封装了一个组件. 代码解读: 1:created周期中获取传入参数时间的剩余秒数: this.initSecondsLeft() 并绑定间隔事件 intervalEvent ...
- Solon 的 PathVariable 不需注解
相对于 Spring boot 的 path variable :Solon 的不需注解,只要变量名对上即可: //path var demo // @XMapping("e/{p_q}/{ ...
- ssh断连后,保持Linux后台程序连接
#### ssh断连后,如何保持Linux后台程序继续运行?ssh断连后,要想运行在Linux服务器中的程序继续运行,就要用到screen技术.- ##### 新建`session` ```shell ...
- Spark优化之小文件是否需要合并?
我们知道,大部分Spark计算都是在内存中完成的,所以Spark的瓶颈一般来自于集群(standalone, yarn, mesos, k8s)的资源紧张,CPU,网络带宽,内存.Spark的性能,想 ...
- 【AHOI2009】 维护序列 - 线段树
题目描述 老师交给小可可一个维护数列的任务,现在小可可希望你来帮他完成. 有长为N的数列,不妨设为a1,a2,…,aN .有如下三种操作形式: (1)把数列中的一段数全部乘一个值; (2)把数列中的一 ...
- 反制面试官 | 14张原理图 | 再也不怕被问 volatile!
反制面试官 | 14张原理图 | 再也不怕被问 volatile! 悟空 爱学习的程序猿,自主开发了Java学习平台.PMP刷题小程序.目前主修Java.多线程.SpringBoot.SpringCl ...