P9069 [Ynoi Easy Round 2022] 堕天作战 TEST_98 Solution

news2025/1/22 13:37:20

Description

给定长为 n n n 的序列 a = ( a 1 , a 2 , ⋯   , a n ) a=(a_1,a_2,\cdots,a_n) a=(a1,a2,,an),有 m m m 个操作,分以下两种:

  1. modify ⁡ ( l , r , k ) \operatorname{modify}(l,r,k) modify(l,r,k):对于所有满足 i ∈ [ l , r ] ,    a i ≠ k i \in [l,r],\; a_i \neq k i[l,r],ai=k i i i,执行 a i ← a i − k a_i \leftarrow a_i -k aiaik
  2. query ⁡ ( l , r , k ) \operatorname{query}(l,r,k) query(l,r,k):求 ( ∑ i = l r a i )   m o d   2 32 (\sum\limits_{i=l}^r a_i) \bmod 2^{32} (i=lrai)mod232

Limitations

1 ≤ n , m ≤ 5 × 1 0 5 1 \le n,m \le 5\times 10^5 1n,m5×105
1 ≤ l ≤ r ≤ n 1 \le l \le r \le n 1lrn
0 ≤ a i , k ≤ 1 0 9 0 \le a_i,k \le 10^9 0ai,k109
6 s , 512 MB 6\text{s},512\text{MB} 6s,512MB

Solution

考虑上线段树,但 modify ⁡ \operatorname{modify} modify 很棘手,先分讨一下。
下面设线段树当前节点的最小值为 u u u

  • u > k u > k u>k,直接打标记即可。
  • u < k u < k u<k,遍历所有 u ≤ k u \le k uk 的儿子进行修改,由于每个数至多只会修改
    一次,所以这部分耗费 O ( log ⁡ n ) \mathcal{O}(\log n) O(logn) 时间(单次操作)。
  • u = k u = k u=k,最难搞的情况,考虑维护次小值 v v v
  • v > 2 k v > 2k v>2k,次小值无法变成最小值,打上最小值以外的数减去 k k k 的标记即可。
  • v ≤ 2 k v \le 2k v2k,次小值会变成最小值,暴力遍历儿子进行更新,由于 v − k ≤ v 2 v-k \le \frac{v}{2} vk2v,所以 v v v 每次至少减半,这部分耗费 O ( log ⁡ n log ⁡ V ) \mathcal{O}(\log n \log V) O(lognlogV) 时间。( V V V 是值域)

综上述,时间复杂度 O ( n log ⁡ n log ⁡ V ) \mathcal{O}(n\log n \log V) O(nlognlogV)

Code

4.41 KB , 18.88 s , 109.38 MB    (in   total,   C++   20   with   O2) 4.41\text{KB},18.88\text{s},109.38\text{MB}\; \texttt{(in total, C++ 20 with O2)} 4.41KB,18.88s,109.38MB(in total, C++ 20 with O2)

// Problem: P9069 [Ynoi Easy Round 2022] 堕天作战 TEST_98
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P9069
// Memory Limit: 512 MB
// Time Limit: 6000 ms
// 
// Powered by CP Editor (https://cpeditor.org)
#pragma struct pack(1)
#include <bits/stdc++.h>
using namespace std;

using i64 = long long;
using ui64 = unsigned long long;
using i128 = __int128;
using ui128 = unsigned __int128;
using f4 = float;
using f8 = double;
using f16 = long double;

template<class T>
bool chmax(T &a, const T &b){
	if(a < b){ a = b; return true; }
	return false;
}

template<class T>
bool chmin(T &a, const T &b){
	if(a > b){ a = b; return true; }
	return false;
}

const ui64 inf = ULLONG_MAX;
struct Node {
    int l, r, cnt;
    ui64 sum, min, sec, tag, mark;
};
using Tree = vector<Node>;
int ls(int u) { return 2 * u + 1; }
int rs(int u) { return 2 * u + 2; }

inline void pushup(Tree& tr, int u) {
    tr[u].sum = tr[ls(u)].sum + tr[rs(u)].sum;
    if (tr[ls(u)].min == tr[rs(u)].min) {
		tr[u].min = tr[ls(u)].min;
		tr[u].cnt = tr[ls(u)].cnt + tr[rs(u)].cnt;
		tr[u].sec = min(tr[ls(u)].sec, tr[rs(u)].sec);
	}
	else if (tr[ls(u)].min < tr[rs(u)].min) {
	    tr[u].min = tr[ls(u)].min;
	    tr[u].cnt = tr[ls(u)].cnt;
	    tr[u].sec = min(tr[ls(u)].sec, tr[rs(u)].min);
	}
	else {
	    tr[u].min = tr[rs(u)].min;
	    tr[u].cnt = tr[rs(u)].cnt;
	    tr[u].sec = min(tr[ls(u)].min, tr[rs(u)].sec);
	}
}

inline void addT(Tree& tr, int u, ui64 tag) {
    int len = tr[u].r - tr[u].l + 1;
    tr[u].sum -= tag * len;
    tr[u].min -= tag;
    tr[u].sec -= tag;
    tr[u].tag += tag;
}

inline void addM(Tree& tr, int u, ui64 tag) {
    int len = tr[u].r - tr[u].l + 1;
    tr[u].sum -= tag * (len - tr[u].cnt);
    tr[u].sec -= tag;
    tr[u].mark += tag;
}

inline void pushdown(Tree& tr, int u) {
    if (tr[u].tag) { 
        addT(tr, ls(u), tr[u].tag);
        addT(tr, rs(u), tr[u].tag);
        tr[u].tag = 0;
    }
	if (tr[u].mark) {
		if (tr[ls(u)].min == tr[u].min) {
		    addM(tr, ls(u), tr[u].mark);
		}
		else {
		    addT(tr, ls(u), tr[u].mark);
		}
		
		if (tr[rs(u)].min == tr[u].min) {
		    addM(tr, rs(u), tr[u].mark);
		}
		else {
		    addT(tr, rs(u), tr[u].mark);
		}
		
		tr[u].mark = 0;
	}
}

inline void build(Tree& tr, int u, int l, int r, vector<int>& a) {
    tr[u].l = l;
    tr[u].r = r;
	if (l == r) {
	    tr[u].sum = tr[u].min = a[l];
	    tr[u].cnt = 1;
	    tr[u].sec = inf;
	    return;
	}
	int mid = (l + r) >> 1;
	build(tr, ls(u), l, mid, a);
	build(tr, rs(u), mid + 1, r, a);
	pushup(tr, u);
}

inline void update(Tree& tr, int u, int val) {
    int len = tr[u].r - tr[u].l + 1;
    if (tr[u].min == val && len == tr[u].cnt) {
        return;
    }
	if (tr[u].l == tr[u].r) {
	    tr[u].sum -= val;
	    tr[u].min = tr[u].sum;
	    return;
	}
	if (tr[u].min < val || tr[u].sec <= val * 2) {
		pushdown(tr, u);
		update(tr, ls(u), val);
		update(tr, rs(u), val);
		pushup(tr, u);
		return;
	}
	if (tr[u].min == val) {
	    addM(tr, u, val);
	}
	else {
	    addT(tr, u, val);
	}
}

inline void modify(Tree& tr, int u, int l, int r, int val) {
    if (l <= tr[u].l && tr[u].r <= r) {
        update(tr, u, val);
        return;
    }
    
    int mid = (tr[u].l + tr[u].r) >> 1;
	pushdown(tr, u);
	if (l <= mid) {
	    modify(tr, ls(u), l, r, val);
	}
	if (r > mid) {
	    modify(tr, rs(u), l, r, val);
	}
	pushup(tr, u);
}

inline ui64 query(Tree& tr, int u, int l, int r) {
    if (l <= tr[u].l && tr[u].r <= r) {
        return tr[u].sum;
    }
    
    int mid = (tr[u].l + tr[u].r) >> 1;
    ui64 res = 0;
	pushdown(tr, u);
	if (l <= mid) {
	    res += query(tr, ls(u), l, r);
	}
	if (r > mid) {
	    res += query(tr, rs(u), l, r);
	}
	return res;
}

signed main() {
	ios::sync_with_stdio(0);
	cin.tie(0), cout.tie(0);
	
	int n, m;
	cin >> n >> m;
	
	vector<int> a(n);
	for (int i = 0; i < n; i++) {
	    cin >> a[i];
	}
	
	Tree tree(n << 2);
	build(tree, 0, 0, n - 1, a);
	
	ui64 lst = 0;
	for (int i = 0; i < m; i++) {
	    int op, l, r, x;
	    cin >> op >> l >> r;
	    l ^= lst, r ^= lst;
	    l--, r--;

	    if (op == 1) {
	        cin >> x;
	        x ^= lst;
	        if (x != 0) {
	            modify(tree, 0, l, r, x);
	        }
	    }
	    else {
	        cout << (lst = query(tree, 0, l, r)) << endl;
	        lst &= 1048575;
	    }
	}
	
	return 0;
}

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

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

相关文章

HTML新春烟花

系列文章 序号目录1HTML满屏跳动的爱心&#xff08;可写字&#xff09;2HTML五彩缤纷的爱心3HTML满屏漂浮爱心4HTML情人节快乐5HTML蓝色爱心射线6HTML跳动的爱心&#xff08;简易版&#xff09;7HTML粒子爱心8HTML蓝色动态爱心9HTML跳动的爱心&#xff08;双心版&#xff09;1…

从结构嵌套的幻梦里:递归与数据构建的精巧和鸣

大家好啊&#xff0c;我是小象٩(๑ω๑)۶ 我的博客&#xff1a;Xiao Xiangζั͡ޓއއ 很高兴见到大家&#xff0c;希望能够和大家一起交流学习&#xff0c;共同进步。 这一节我们来学习递归的相关知识 函数递归 一、什么是递归1.1 递归的思想 二、递归的限制条件三、递归的…

【Linux系统】—— 编译器 gcc/g++ 的使用

【Linux系统】—— 编译器 gcc/g 的使用 1 用 gcc 直接编译2 翻译环境2.1 预处理&#xff08;进行宏替换&#xff09;2.2 编译&#xff08;生成汇编&#xff09;2.3 汇编&#xff08;生成机器可识别代码&#xff09;2.4 链接2.5 记忆小技巧2.6 编译方式2.7 几个问题2.7.1 如何理…

【Unity3D】3D物体摆放、场景优化案例Demo

目录 PlaceManager.cs(放置管理类) Ground.cs(地板类) 和 GroundData.cs(地板数据类) 额外知识点说明 1、MeshFilter和MeshRenderer的Bounds区别 2、Gizmos 绘制一个平行于斜面的立方体 通过网盘分享的文件&#xff1a;PlaceGameDemo2.unitypackage 链接: https://pan.baid…

智能系统的感知和决策

智能系统在感知和决策过程中具备的关键能力表现在智能感知/自主判定上&#xff0c;下面可以从感知的本质、自主判断的含义及其在智能系统中的作用进行深入分析。 1、智能感知&#xff1a;信息获取与理解 智能感知是指智能系统通过传感器或其他数据采集手段获取环境中的信息&…

Spring 中的事件驱动模型

事件驱动的基本了解 事件模式也就是观察者模式&#xff0c;当一个对象改变的时候&#xff0c;所有依赖的对象都会收到一个通知。 Subject&#xff1a;抽象主题 Observer&#xff1a;具体主题 Concrete Subject&#xff1a;抽象观察者&#xff0c;在得到更新通知之后去更新自…

linux-FTP服务配置与应用

也许你对FTP不陌生&#xff0c;但是你是否了解FTP到底是个什么玩意&#xff1f; FTP 是File Transfer Protocol&#xff08;文件传输协议&#xff09;的英文简称&#xff0c;而中文简称为 “文传协议” 用于Internet上的控制文件的双向传输。同时&#xff0c;它也是一个应用程序…

linux-NFS网络共享存储服务配置

1.NFS服务原理 NFS会经常用到&#xff0c;用于在网络上共享存储&#xff0c;这样讲&#xff0c;你对NFS可能不太了解&#xff0c;举一个例子&#xff0c; 加入有三台机器A,B,C&#xff0c;它们需要访问同一个目录&#xff0c;目录中都是图片&#xff0c;传统的做法是把这些 图…

Jenkins 启动

废话 这一阵子感觉空虚&#xff0c;心里空捞捞的&#xff0c;总想找点事情做&#xff0c;即使这是一件微小的事情&#xff0c;空余时间除了骑车、打球&#xff0c;偶尔朋友聚会 … 还能干什么呢&#xff1f; 当独自一人时&#xff0c;究竟可以做点什么&#xff0c;填补这空虚…

消息队列篇--原理篇--Pulsar(Namespace,BookKeeper,类似Kafka甚至更好的消息队列)

Apache Pulusar是一个分布式、多租户、高性能的发布/订阅&#xff08;Pub/Sub&#xff09;消息系统&#xff0c;最初由Yahoo开发并开源。它结合了Kafka和传统消息队列的优点&#xff0c;提供高吞吐量、低延迟、强一致性和可扩展的消息传递能力&#xff0c;适用于大规模分布式系…

Python配置MITMPROXY中间人监听配置

1、安装python 环境&#xff0c;此处可以使用conda安装:conda create --name my_new_env python3.12 2、pip安装mitmproxy&#xff1a;pip install mitmproxy&#xff0c;安装后如果使用mitmproxy --version 成功返回结果&#xff0c;说明已经在环境变量路径中&#xff0c;如果…

Java-数据结构-二叉树习题(2)

第一题、平衡二叉树 ① 暴力求解法 &#x1f4da; 思路提示&#xff1a; 该题要求我们判断给定的二叉树是否为"平衡二叉树"。 平衡二叉树指&#xff1a;该树所有节点的左右子树的高度相差不超过 1。 也就是说需要我们会求二叉树的高&#xff0c;并且要对节点内所…

【网络原理】万字详解 HTTP 协议

&#x1f970;&#x1f970;&#x1f970;来都来了&#xff0c;不妨点个关注叭&#xff01; &#x1f449;博客主页&#xff1a;欢迎各位大佬!&#x1f448; 文章目录 1. HTTP 前置知识1.1 HTTP 是什么1.2 HTPP 协议应用场景1.3 HTTP 协议工作过程 2. HTTP 协议格式2.1 fiddler…

基于STM32的智能寝室控制系统设计(论文+源码)

1 .系统整体设计 通过需求分析&#xff0c;本设计基于STM32的智能寝室控制系统整体架构如图2.1所示&#xff0c;整系统利用DHT11温湿度传感器获取室内环境数据&#xff0c;并通过OLED显示&#xff0c;提供用户实时信息&#xff0c;火焰传感器和烟雾传感器用于监测火灾情况&…

日历热力图,月度数据可视化图表(日活跃图、格子图)vue组件

日历热力图&#xff0c;月度数据可视化图表&#xff0c;vue组件 先看效果&#x1f447; 在线体验https://www.guetzjb.cn/calanderViewGraph/ 日历图简单划分为近一年时间&#xff0c;开始时间是 上一年的今天&#xff0c;例如2024/01/01 —— 2025/01/01&#xff0c;跨度刚…

铁电存储器FM25CL64B简介及其驱动编写(基于STM32 hal库)

铁电存储器FM25CL64B简介及其驱动编写&#xff08;基于STM32 hal库&#xff09; 文章目录 铁电存储器FM25CL64B简介及其驱动编写&#xff08;基于STM32 hal库&#xff09;前言一、FM25CL64B简介二、驱动代码1.头文件2.c文件 总结 前言 FM25CL64B是赛普拉斯cypress出品的一款铁…

基于微信小程序的科创微应用平台设计与实现(LW+源码+讲解)

专注于大学生项目实战开发,讲解,毕业答疑辅导&#xff0c;欢迎高校老师/同行前辈交流合作✌。 技术范围&#xff1a;SpringBoot、Vue、SSM、HLMT、小程序、Jsp、PHP、Nodejs、Python、爬虫、数据可视化、安卓app、大数据、物联网、机器学习等设计与开发。 主要内容&#xff1a;…

HarmonyOS Next 最强AI智能辅助编程工具 CodeGenie介绍

随着大模型的兴起&#xff0c;在智能编码领域首先获得了应用。 市面上从Microsoft Copilot到国内阿里通义&#xff0c;字节marscode等&#xff0c;都提供了copilot方式的智能编码工具。HarmonyOS Next作为诞生一年的新事物&#xff0c;由于代码量和文档迭代原因&#xff0c;在智…

WPF2-1在xaml为对象的属性赋值.md

1. AttributeValue方式 1.1. 简单属性赋值1.2. 对象属性赋值 2. 属性标签的方式给属性赋值3. 标签扩展 (Markup Extensions) 3.1. StaticResource3.2. Binding 3.2.1. 普通 Binding3.2.2. ElementName Binding3.2.3. RelativeSource Binding3.2.4. StaticResource Binding (带参…

Appium(四)

一、app页面元素定位 1、通过id定位元素: resrouce-id2、通过ClassName定位&#xff1a;classname3、通过AccessibilityId定位&#xff1a;content-desc4、通过AndroidUiAutomator定位5、通过xpath定位xpath、id、class、accessibility id、android uiautomatorUI AutomatorUI自…