O(n)RMQ四毛子

2023-11-09

有一种ST表,叫做±1ST表

这种ST表可以在 O ( n ) O(n) O(n) 的时刻内完成建树

其本质就是分块,大块为整除的ST表,小块的差分数组种类不多,完全可以预处理

现在考虑推广到普通的ST表里

我们发现我们真正关心的是数之间的大小关系。但又要使相邻数之间差恰好为±1

考虑什么东西的差为1。树的欧拉环游序点之间的深度差!

现在需要数之间的大小关系,也需要树,那么,我们建笛卡尔树就行了!

总结一下思路

先建出笛卡尔树,求出其欧拉环游序,然后对其欧拉换有序的差进行±1RMQ即可。

#include<bits/stdc++.h>
using namespace std;
//#define int long long
inline int read(){int x=0,f=1;char ch=getchar(); while(ch<'0'||
ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){
x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}return x*f;}
#define Z(x) (x)*(x)
#define pb push_back
//mt19937 rand(time(0));
//mt19937_64 rand(time(0));
//srand(time(0));
#define N 400010
#define M 400010
//#define mo
struct node {
	int val; 
	int dep, dfn, end; 
	int son[2]; 
}T[N], Mn[N][20];
int n, m, i, j, k;
int l, r, top, q, rt; 
int Pos[N], Dif[N], a[N], Log2[N]; 
int b, c; 

int build() {
	int S[N]; 
	for(int i=1; i<=n; ++i) {
		while(top && T[S[top]].val < T[i].val) 
			T[i].son[0]=S[top], --top; 
		if(top) T[S[top]].son[1]=i; 
		S[++top]=i; 
	}
	top=0; 
	return S[1]; 
}

void dfs(int x) {
	a[top]=x; T[x].dfn=top; ++top;  
	for(int i=0; i<=1; ++i)
		if(T[x].son[i]) {
			int y=T[x].son[i]; 
//			printf("%d -> %d\n", x, y); 
			T[y].dep=T[x].dep+1; 
			dfs(y); 
			a[top++]=x; 
		}
	T[x].end=top-1; 
}

node max(node a, node b) {
	if(a.val>b.val) return a; 
	return b; 
}

node min(node a, node b) {
	if(a.val<b.val) return a; 
	return b; 
}

void pre_ST() {
	int i, j, k, l; 
	b=(int)(ceil(log2(top)/2)); c=top/b; 
	Log2[1]=0; 
	for(i=2; i<=top; ++i) Log2[i]=Log2[i>>1]+1; 
	for(i=0; i<c; ++i) {
		Mn[i][0]=T[a[i*b]]; 
		for(j=1; j<b; ++j) 
			Mn[i][0]=max(Mn[i][0], T[a[i*b+j]]); 
	}
	for(k=l=1; l<c; ++k, l<<=1)
		for(i=0, j=l; i+(l<<1)-1<top; ++i, ++j) {
			Mn[i][k]=max(Mn[i][k-1], Mn[j][k-1]); 
		}
}

void pre_small() {
	int i, j, s; 
	for(i=0; i<=c; ++i) {
		for(j=1; j<b && i*b+j<top; ++j)
			if(T[a[i*b+j]].dep<T[a[i*b+j-1]].dep)
				Dif[i]|=(1<<j-1); 
	}
	for(s=0; s<(1<<b-1); ++s) {
		int v=0, mx=0; Pos[s]=0; 
		for(i=1; i<b; ++i) {
			if(s&(1<<i-1)) --v; 
			else ++v; 
			if(v<mx) {
				mx=v; Pos[s]=i; 
			}
		}	
	}
}

node que_small(int l, int r) {
	int p=l/b; 
	int S=((Dif[p]>>(l-p*b))&((1<<r-l)-1)); 
	return T[a[l+Pos[S]]]; 
}

node que_ST(int l, int r) {
	int k=Log2[r-l+1]; 
	return max(Mn[l][k], Mn[r-(1<<k)+1][k]); 
}

int que(int l, int r) {
	if(l>r) return que(r, l); 
	if(l==r) return T[a[l]].val; 
	int pl=l/b, pr=r/b; 
	if(pl==pr) return que_small(l, r).val; 
	node ans=max(que_small(l, pl*b+b-1), que_small(pr*b, r)); 
	if(pl+1<=pr-1) ans=max(ans, que_ST(pl+1, pr-1)); 
	return ans.val; 
}

signed main()
{
//	freopen("in.txt", "r", stdin);
//	freopen("out.txt", "w", stdout);
//	T=read();
//	while(T--) {
//
//	}
	n=read(); q=read(); 
	for(i=1; i<=n; ++i) T[i].val=read(); 
	rt=build(); 
	dfs(rt); 
	pre_ST(); 
	pre_small(); 
	while(q--) {
		l=read(); r=read(); 
		printf("%d\n", que(T[l].dfn, T[r].dfn)); 
	}
	return 0;
}


参考:CSP2021-S 初赛

在这里插入图片描述

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

O(n)RMQ四毛子 的相关文章

  • 纯 CSS 实现高度与宽度成比例的效果

    HTML p class image container img src xxx png p CSS p image container width 100 height 0 padding bottom 60 overflow hidde
  • docker 训练深度学习_Docker:把深度学习装进“口袋”

    熟悉深度学习开发框架的朋友都知道 caffe 的安装麻烦 耗时长 而且很容易因为电脑环境没有配置好导致错误 每次有新同事来 我们都要在他电脑上安装一遍 caffe 经常会因为电脑环境而不断报错 目前TensorFlow受到很多深度学习开发者
  • 阿里云oss视频上传及预览图汇总

    阿里云OSS视频上传的几种方式 1 Web端直传实践简介 https help aliyun com document detail 31923 html spm a2c4g 11186623 6 631 apBNnA 2 JavaScri
  • cacheable 过期设置

    1 增加pom依赖
  • 【Java面试题】IO & NIO

    1 什么是IO流 它是一种数据的流从源头流到目的地 比如文件拷贝 输入流和输出流 输入流从文件中读取数据存储到进程 process 中 输出流从进程中读取数据然后写入到目标文件 2 Java中有几种类型的流 按照单位大小 字符流 字节流 按
  • MCDF实验——Lab3

    Lab3将在Lab2的基础上使用随机约束和环境结构来改进完善实验代码 Lab3中将对generator和initiator之间的数据生成和数据传输的处理进行改进 还将完善何时结束测试 将其主动权交于generator而不再是test组件 在
  • ipython下分析数据

    终端下输入 ipython notebook Jupyter notebook才会有 matplotlib inline import pandas as pd import numpy as np import matplot pyplt
  • C++ Web 编程

    什么是 CGI 公共网关接口 CGI 是一套标准 定义了信息是如何在 Web 服务器和客户端脚本之间进行交换的 CGI 规范目前是由 NCSA 维护的 NCSA 定义 CGI 如下 公共网关接口 CGI 是一种用于外部网关程序与信息服务器
  • Java 优雅判空的方法

    为了避免 NullPointerException 出现 我们经常会进行如下操作 if object null object do 那么每次使用一个不确定的对象 都需要进行一次判空 so 可能你会想 项目中肯定不止你一个人会这样一顿操作 然

随机推荐

  • C语言 - 删除数组中的指定元素(Delete)

    目录 引言 题目说明 举例分析 初始状态 step 1 Step 2 Step 3 Step 4 将最后一个元素赋值为0 Step 5 将i下标向前迁移一位 重复判断 代码实现 总结 引言 这是数组中一个比较简单的题目 今天把它单独拎出来的
  • 自注意力机制超级详解(Self-attention)

    Attention ia all you need论文原文 仅供自我学习使用 Self attention想要解决的问题 目前input都是一个向量 输出可能是一个数值 一个类别 那么如果输入是一排向量呢 且输入的向量数可以改变的话 如何处
  • Transformer(self attention)超详解&示例代码

    目录 说明 数据集介绍 Encoder Input Embedding Positional Encoding Multi Head Attention 第一个Add Norm Feed Forward 第二个Add Norm Decode
  • Git 在 Linux 上和 windows 上的安装

    1 在 Linux 上安装 Centos 上安装 yum y install https repo ius io ius release el7 rpm yum y install git222 git version git versio
  • python如何画多个纵坐标_python-PANDAS绘制多个Y轴

    没有plt的简单解决方案 ax1 df1 plot ax2 ax1 twinx ax2 spines right set position axes 1 0 df2 plot ax ax2 ax3 ax1 twinx ax3 spines
  • CryptoPP的RandomPool算法的使用

    随机数发生器是密码学的一个重要原语 密码学库CryptoPP中提供了一些随机数发生器算法 前面讲述了LC RNG算法的使用 https blog csdn net Lunar Queen article details 81542372 今
  • Python3.6 statsmodels安装

    最近在学习Python 就想着把过程中遇到的比较难解决的问题都记录下来 方便以后查阅 也希望能给有类似问题的小伙伴提供一些帮助 和安装其他模块一样 我最先想到的当然是在命令窗口输入 pip install statsmodels 来安装 没
  • requestFocus() 与 requestFocusFromTouch() 方法的区别

    两个方法都是获取焦点 requestFocus 源码如下 Call this to try to give focus to a specific view or to one of its descendants and give it
  • 在JSP中使用数据库

    需要用到的工具 软件 数据库 MySQL或者SQL Service 以及相对应的jdbc链接器 连接器的对应版本不一样会导致连接失败 eclipse软件 Navivat 没有也行 我用的版本 一 首先连接数据库 在数据库中建立一个名为Boo
  • 史上最全ConstraintLayout使用详解

    今日科技快讯 全球最大半导体代工制造商台积电周四公布了2021年第一季度喜人财报 净利润达49亿美元 同比增长19 然而该公司首席执行官魏哲家警告称 全球芯片短缺可能会持续到2022年 这个问题如今正在严重影响全球工业 作者简介 又要到愉快
  • asp.net 三种基本的参数传递方法(web forms)

    asp net 三种基本的参数传递方法 web forms ASP NET WEB FORMS 给开发者提供了极好的事件驱动开发模式 然而这种简单的应用程序开发模式却给我们带来了一些小问题 举个例子 在传统的ASP应用程序中 你能够通过PO
  • 项目内导入相关py文件时:ModuleNotFoundError: No module named 错误的解决办法

    1 如果需要导入的包和导入包的文件在同一个目录下 直接import即可 例如b py需要用a py的某些代码就可以直接import 2 如果导入的模块是在主程序所在目录的子目录下 可以在子目录中增加一个空白的 init py文件 该文件使得
  • Go面试必会基础题

    文章目录 1 下面代码有什么错误 2 下面代码有什么问题 3 下面代码输出什么 4 下面这段代码输出什么 1 下面代码有什么错误 func main one 0 one 1 参考答案及解析 变量重复声明 不能在单独的声明中重复声明一个变量
  • Instant-ngp 在 window10 上的部署

    一 环境部署 1 Visual Studio 2019 安装 Visual Studio 2019 预览版发行说明 Microsoft Docs 从上述链接下载Community 2019版并安装 界面如下 选择C 选项 右侧复选框默认 点
  • 光立方某一环节编程

    功能 逐一叠加点亮坐标 0 7 Z 其他熄灭 代码 这段代码解释达林顿管和74hc573怎么编程 void CUBE LOG1 uint8 z CLEAR Y 0x80 X 0x01 X 0x00 for z 0 z lt 8 z Z 0x
  • 【成长笔记】html+jqurey全选反选取消--让世界爱上中国造

    一 代码
  • 会话及会话技术、Cookie对象、Session对象 详解

    会话及会话技术 当用户通过浏览器访问Web应用时 通常情况下 服务器需要对用户的状态进行跟踪 例如 用户在网站结算商品时 Web服务器必须根据 请求用户的身份 找到该用户所购买的商品 在Web开发中 服务器跟踪用户信息的技术称为会话技术 下
  • SQL Incorrect syntax near 错误

    案例分析 在我进行Top数据查询的时候用到如下SQL语句 select top TopCount from MD DataCubeDataDetail2 nolock Where 1 1 and seqno gt BeginSeqno an
  • 分布式事务的典型处理方式:2PC、TCC、异步确保和最大努力型

    分布式事务的典型处理方式 2PC TCC 异步确保和最大努力型 标签 分布式事务处理方式2PCTCC异步确保 2017 04 23 16 49 2617人阅读 评论 0 收藏 举报 分类 Java Web 195 目录 1 柔性事务和刚性事
  • O(n)RMQ四毛子

    有一种ST表 叫做 1ST表 这种ST表可以在 O n O n O n 的时刻内完成建树 其本质就是分块 大块为整除的ST表 小块的差分数组种类不多 完全可以预处理 现在考虑推广到普通的ST表里 我们发现我们真正关心的是数之间的大小关系 但