CODEVS-1758-维护数列-NOI2005-splay

最新推荐文章于 2019-01-13 14:48:00 发布
最新推荐文章于 2019-01-13 14:48:00 发布
阅读量1.2k 收藏
点赞数 2
分类专栏： 题解 文章标签： splay
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_21110267/article/details/44873991
版权

描述

要求支持区间插入、区间修改、区间翻转、区间删除、区间求和 和求和最大的子列.


分析

  • 从最开始学完splay做了翻转区间后就想做这个题目, 结果WA了N次后失去调试的信心, 40分收场(这题暴力30分)
  • 快省选了想拿出来再做一下, 因为splay的区间操作这个题算是最全的了, 不做一下的话总担心模版是错的.
  • 然后做了好长时间...终于不耐烦了拿HZWER的改了改, 直到改到所有操作和我的都写的一模一样——除了标记下传. 然后就很开(伤)心地接受了这个现实.
  • 在网上看每个人写的标记下传都不一样, 有的写的简单有的很复杂(比如HZWER). 但一种标记下传对应一种维护方式吧, 看网上很少有写如何处理细节的, 我就写详细点吧.
  • 首先, mls, mrs, maxs分别记录左、右端开始的最大子序列和, 和整段区间的最大子序列和.
  • 规定如果整段区间为负值, mls和mrs可以取0, 而maxs必须取一个存在的数, 这样也是为了维护起来方便吧.
  • 那么mls可以取左子结点的mls、左子结点的sum+根结点的权值+右子结点的mls. mrs同理.
  • maxs可以取左子结点的maxs、右子结点的maxs、左子结点的mrs+根结点权值+右子结点的mls.
  • 这就是mls、mrs最小为0的优势.
  • 在标记下传时下传的实际上是子结点的标记, 而根结点的标记实际上已经在此之前生效了. 这样可以避免在旋转时出错. 旋转里的那个 ch[p][d] 的标记可能没有生效.
  • 遇到带区间修改的题目都要用这种标记下传的方式吧..
  • 又想起原来写过的区间翻转类的题目, 那里的标记下传只需要下传当前节点的即可, 那是因为翻转不影响最大值、不影响求和. 而区间修改会影响.

#include 
   
   
    
    
#include 
    
    
     
     
using namespace std;

const int INF = 0x3f3f3f3f;
const int maxn = 3000000;

int n, m;
int fa[maxn], ch[maxn][2], id[maxn];
int mls[maxn], mrs[maxn], maxs[maxn], s[maxn], rech[maxn], v[maxn], sumv[maxn];
bool mdf[maxn], rev[maxn];

#define lc ch[o][0]
#define rc ch[o][1]

void update(int o) {
	s[o] = s[lc] + s[rc] + 1;
	sumv[o] = sumv[lc] + sumv[rc] + v[o];
	maxs[o] = max(maxs[lc], maxs[rc]);
	maxs[o] = max(maxs[o], mrs[lc] + v[o] + mls[rc]);
	mls[o] = max(mls[lc], sumv[lc] + v[o] + mls[rc]);
	mrs[o] = max(mrs[rc], sumv[rc] + v[o] + mrs[lc]);
}

void pushdown(int o) {
	if(mdf[o]) {
		mdf[o] = rev[o] = 0;
		if(lc) mdf[lc] = 1, v[lc] = v[o], sumv[lc] = v[o]*s[lc];
		if(rc) mdf[rc] = 1, v[rc] = v[o], sumv[rc] = v[o]*s[rc];
		if(v[o] >= 0) {
			if(lc) mls[lc] = mrs[lc] = maxs[lc] = sumv[lc];
			if(rc) mls[rc] = mrs[rc] = maxs[rc] = sumv[rc];
		} else {
			if(lc) mls[lc] = mrs[lc] = 0, maxs[lc] = v[o];
			if(rc) mls[rc] = mrs[rc] = 0, maxs[rc] = v[o];
		}
	}
	if(rev[o]) {
		rev[o]^=1;
		rev[lc]^=1;
		rev[rc]^=1;
		swap(mls[lc], mrs[lc]);
		swap(mls[rc], mrs[rc]);
		swap(ch[lc][0], ch[lc][1]);
		swap(ch[rc][0], ch[rc][1]);
	}
}

void rotate(int& o, int x)  {
	int p = fa[x], q = fa[p];
	bool d = (ch[p][1] == x), d2 = (ch[q][1] == p);
	if(o == p) o = x; else ch[q][d2] = x;
	fa[x] = q; fa[p] = x; fa[ch[p][d] = ch[x][d^1]] = p;
	ch[x][d^1] = p; update(p); update(x);
}

void splay(int& o, int x) {
	while(o != x) {
		int p = fa[x], q = fa[p];
		if(o != p)
			if(ch[p][0] == x^ch[q][0] == p) rotate(o, x); 
			else rotate(o, p);
		rotate(o, x);
	}
}

void build(int L, int R, int p, int d) {
	if(L > R) return;
	int o = (L+R)>>1;
	ch[p][d] = o; fa[o] = p;
	build(L, o-1, o, 0);
	build(o+1, R, o, 1);
	update(o);
}

int kth(int o, int k) {
	pushdown(o);
	if(s[lc] + 1 == k) return o;
	else if(s[lc] >= k) return kth(lc, k);
	else return kth(rc, k-s[lc]-1);
}


int o, pos, tot, val;

void insert() {
	splay(o, kth(o, pos+1));
	splay(rc, kth(o, pos+2));
	for(int i = n+1; i <= n+tot; i++) scanf("%d", &v[i]);
	build(n+1, n+tot, rc, 0);
	n += tot; //
	update(rc); update(o);
}

void remove() {
	splay(o, kth(o, pos));
	splay(rc, kth(o, pos + tot + 1));
	ch[rc][0] = 0;
	update(rc); update(o);
}

void modify() {
	splay(o, kth(o, pos));
	splay(rc, kth(o, pos + tot + 1));
	int p = ch[rc][0];
	v[p] = val;
	mdf[p] = 1;
	sumv[p] = v[p]*s[p];
	if(val >= 0) mls[p] = mrs[p] = maxs[p] = sumv[p];
	else mls[p] = mrs[p] = 0, maxs[p] = v[p];
	update(rc); update(o);
}

void rever() {
	splay(o, kth(o, pos));
	splay(rc, kth(o, pos + tot + 1));
	if(!mdf[o]) {
		int p = ch[rc][0];
		rev[p]^=1;
		swap(ch[p][0], ch[p][1]);
		swap(mls[p], mrs[p]);
	}
	update(rc); update(o);
}

void query() {
	splay(o, kth(o, pos));
	splay(rc, kth(o, pos + tot + 1));
	printf("%d\n", sumv[ch[rc][0]]);
}

int main()
{
	scanf("%d %d", &n, &m); n += 2;
	v[1] = v[n] = maxs[0] = -INF; //
	for(int i = 2; i < n; i++) scanf("%d", &v[i]);
	build(1, n, 0, 0);
	o = (n+1)>>1;
	
	for(int i = 1; i <= m; i++) {
		char opt[10];
		scanf("%s", opt);
		switch(opt[2]) {
			case 'S': scanf("%d %d", &pos, &tot); insert(); break;
			case 'L': scanf("%d %d", &pos, &tot); remove(); break;
			case 'V': scanf("%d %d", &pos, &tot); rever(); break;
			case 'T': scanf("%d %d", &pos, &tot); query(); break;
			case 'K': scanf("%d %d %d", &pos, &tot, &val); modify(); break;
			case 'X': printf("%d\n", maxs[o]);
		}
	}
	return 0;
}
#include 
     
     
      
      
#include 
      
      
       
       
using namespace std;

const int INF = 0x3f3f3f3f;
const int maxn = 3000000;

int n, m;
int ch[maxn][2], id[maxn];
int mls[maxn], mrs[maxn], maxs[maxn], s[maxn], rech[maxn], v[maxn], sumv[maxn];
bool mdf[maxn], rev[maxn];

#define lc ch[o][0]
#define rc ch[o][1]

void update(int o) {
	s[o] = s[lc] + s[rc] + 1;
	sumv[o] = sumv[lc] + sumv[rc] + v[o];
	maxs[o] = max(maxs[lc], maxs[rc]);
	maxs[o] = max(maxs[o], mrs[lc] + v[o] + mls[rc]);
	mls[o] = max(mls[lc], sumv[lc] + v[o] + mls[rc]);
	mrs[o] = max(mrs[rc], sumv[rc] + v[o] + mrs[lc]);
}

void pushdown(int o) {
	if(mdf[o]) {
		mdf[o] = rev[o] = 0;
		if(lc) mdf[lc] = 1, v[lc] = v[o], sumv[lc] = v[o]*s[lc];
		if(rc) mdf[rc] = 1, v[rc] = v[o], sumv[rc] = v[o]*s[rc];
		if(v[o] >= 0) {
			if(lc) mls[lc] = mrs[lc] = maxs[lc] = sumv[lc];
			if(rc) mls[rc] = mrs[rc] = maxs[rc] = sumv[rc];
		} else {
			if(lc) mls[lc] = mrs[lc] = 0, maxs[lc] = v[o];
			if(rc) mls[rc] = mrs[rc] = 0, maxs[rc] = v[o];
		}
	}
	if(rev[o]) {
		rev[o]^=1;
		rev[lc]^=1;
		rev[rc]^=1;
		swap(mls[lc], mrs[lc]);
		swap(mls[rc], mrs[rc]);
		swap(ch[lc][0], ch[lc][1]);
		swap(ch[rc][0], ch[rc][1]);
	}
}

void rotate(int& o, int d) {
	int p = ch[o][d^1]; ch[o][d^1] = ch[p][d]; ch[p][d] = o;
	update(o); update(p); o = p;
}

int cmp(int o, int k) {
	if(s[lc] + 1 == k) return -1;
	return k < s[lc] + 1 ? 0 : 1;
}

void splay(int& o, int k) {
	pushdown(o);
	int d = cmp(o, k);
	if(d == -1) return;
	if(d == 1) k -= (s[lc] + 1);
	int p = ch[o][d];
	pushdown(p);
	int d2 = cmp(p, k);
	int k2 = (d2 == 0 ? k : k-s[ch[p][0]]-1);
	if(d2 != -1) {
		splay(ch[p][d2], k2);
		if(d == d2) rotate(o, d^1); else rotate(ch[o][d], d);
	}
	rotate(o, d^1);
}

void build(int& o, int L, int R) {
	if(L > R) return;
	o = (L+R)>>1;
	build(lc, L, o-1);
	build(rc, o+1, R);
	update(o);
}

void getInterval(int& o, int pos, int tot) {
	splay(o, pos);
	splay(rc, pos + tot - s[lc]);
}


int o, pos, tot, val;

void insert() {
	getInterval(o, pos+1, 0);
	for(int i = n+1; i <= n+tot; i++) scanf("%d", &v[i]);
	build(ch[rc][0], n+1, n+tot);
	n += tot; //
	update(rc); update(o);
}

void remove() {
	getInterval(o, pos, tot);
	ch[rc][0] = 0;
	update(rc); update(o);
}

void modify() {
	getInterval(o, pos, tot);
	int p = ch[rc][0];
	v[p] = val;
	mdf[p] = 1;
	sumv[p] = v[p]*s[p];
	if(val >= 0) mls[p] = mrs[p] = maxs[p] = sumv[p];
	else mls[p] = mrs[p] = 0, maxs[p] = v[p];
	update(rc); update(o);
}

void rever() {
	getInterval(o, pos, tot);
	if(!mdf[o]) {
		int p = ch[rc][0];
		rev[p]^=1;
		swap(ch[p][0], ch[p][1]);
		swap(mls[p], mrs[p]);
	}
	update(rc); update(o);
}

void query() {
	getInterval(o, pos, tot);
	printf("%d\n", sumv[ch[rc][0]]);
}

int main()
{
	scanf("%d %d", &n, &m); n += 2;
	v[1] = v[n] = maxs[0] = -INF; //
	for(int i = 2; i < n; i++) scanf("%d", &v[i]);
	build(o, 1, n);
	
	for(int i = 1; i <= m; i++) {
		char opt[10];
		scanf("%s", opt);
		switch(opt[2]) {
			case 'S': scanf("%d %d", &pos, &tot); insert(); break;
			case 'L': scanf("%d %d", &pos, &tot); remove(); break;
			case 'V': scanf("%d %d", &pos, &tot); rever(); break;
			case 'T': scanf("%d %d", &pos, &tot); query(); break;
			case 'K': scanf("%d %d %d", &pos, &tot, &val); modify(); break;
			case 'X': printf("%d\n", maxs[o]);
		}
	}
	return 0;
}


gohsyi
关注
  • 2
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
[NOI2005] 维护数列（Splay）
hipamp的博客
10-11 375
题意 要求维护一个数列，支持以下操作： 其中，插入总数字个数小于 4e64e64e6，每个数在 [−1000,1000][-1000,1000][−1000,1000] 范围，保证每次操作合法。 分析 这道题是一道 SplaySplaySplay 码农题，细节颇多，让我们一个一个操作分析。 插入操作 我们可以先对要插入的这 tottottot 个数建树，得到树根 trttrttrt，再将 trttrttrt 插入原序列。 删除操作 就是普通的删除，将 r+1r+1r+1 旋到根，将 l−1l-1l−1
【平衡树】【NOI 2005】维护数列
morestep的专栏
07-27 965
339. [NOI2005] 维护数列★★★★☆ 输入文件：seq2005.in 输出文件：seq2005.out 简单对比 时间限制：3 s 内存限制：256 MB 【问题描述】请写一个程序，要求维护一个数列，支持以下 6 种操作：（请注意，格式栏 中的下划线‘ _ ’表示实际输入文件中的空格） 【输入格式】输入文件的第 1 行包含两个数 N 和 M，N 表示初始时数列中数的个数
[NOI2005]维护数列
a6t2007的博客
09-08 100
一道splay综合大板子题。 题面：https://www.lydsy.com/JudgeOnline/problem.php?id=1500 下面是题解： 首先对每个点维护这些量： 1.两个儿子（ch[2]） 2.父节点（fa） 3.当前点权值（vl）和子树权值（sum） 4.修改标记（xg），反转标记（fz）。 5.mx，mxl，mxr（不清楚的建议做小白逛公园）。 ...
【NOI2005】维护数列
Facico的博客
06-28 882
Description　请写一个程序，要求维护一个数列，支持以下6种操作：（请注意，格式栏中的下划线‘ _ ’表示实际输入文件中的空格） 　　1. 插入 INSERT_posi_tot_c1_c2_…_ctot 在当前数列的第posi个数字后插入tot个数字：c1, c2, …, ctot；若在数列首插入，则posi为0 　　2. 删除 DELETE_posi_tot 从当前数列的第posi个数
ACM-NOI-CSP竞赛宝典策略、代码与实战经验共享-acm和noip
最新发布
06-28
### ACM-NOI-CSP竞赛宝典策略、代码与实战经验共享 #### 一、竞赛基础与策略 ##### 1. 竞赛简介 - **ACM/ICPC**：这项竞赛强调团队合作和解决实际问题的能力。参赛队伍通常由三人组成，在限定的五小时内完成8至13...
NOI Linux2.0使用系列视频集-2023-06-14(C).pdf
06-14
【NOI Linux 2.0】是专为信息学奥赛（CSP-J, CSP-S）设计的操作系统，它提供了一个纯净且优化的环境，便于参赛者进行C++编程和算法训练。NOI Linux 2.0的更新旨在提升用户体验，增加新特性，并优化现有功能。以下是...
NOI学习资料整理（依据「NOI大纲」）CSP-J CSP-S NOIP NOI 省选等--2022.03.05(B).pdf
03-05
"NOI学习资料整理（依据「NOI大纲」）CSP-J CSP-S NOIP NOI 省选等" 根据提供的文件信息，我们可以总结出以下知识点： 一、NOI学习资料整理 NOI学习资料整理是基于「NOI大纲」的学习资源整理，涵盖CSP-J、CSP-S、...
信息学竞赛资料学校上课所用，关于CSP-S/NOIP/NOI
01-03
《信息学竞赛全面解析：CSP-S/NOIP/NOI深度探究》 在信息技术日益发展的今天，信息学竞赛已经成为培养青少年计算机科学素养的重要途径。CSP-S（中国计算机学会认证服务提供者-软件能力认证）、NOIP（全国青少年信息...
NOI数学--2021.07.28(E).pdf
08-12
NOI数学竞赛是一场面向中学生的全国性信息学奥林匹克竞赛，其数学部分覆盖了从基础到高级的各类数学知识点，是检验学生数学与编程能力的重要平台。本文档中提到的知识点主要包括数的概念与运算、编码知识、初等数学...
[题解]NOI2005 维护数列
Saramanda的博客
02-14 448
Description Input 输入的第1 行包含两个数N 和M(M ≤20 000)，N 表示初始时数列中数的个数，M表示要进行的操作数目。 第2行包含N个数字，描述初始时的数列。 以下M行，每行一条命令，格式参见问题描述中的表格。 任何时刻数列中最多含有500 000个数，数列中任何一个数字均在[-1 000, 1 000]内。 插入的数字总数不超过4 000
【平衡树】[NOI2005]维修数列
weixin_30412013的博客
01-13 90
Description 请写一个程序，要求维护一个数列，支持以下 6 种操作： 请注意，格式栏 中的下划线‘ _ ’表示实际输入文件中的空格 Input 输入的第1 行包含两个数N 和M(M ≤20 000)，N 表示初始时数列中数的个数，M表示要进行的操作数目。 第2行包含N个数字，描述初始时的数列。 以下M行，每行一条命令，格式参见问题描述中的表格。 任何时刻数列中最多含有5...
【noi2005】维护数列
Evil.livE
08-30 2973
      哎呀，不早了，这道题从9点调到了11：30，总之这道题和之前的【poj3580】非常相似，很欣慰的是这次编写的时候思维非常连贯，一口气写完，debug了2个小时发现一直是读入萎了~~，不知道我的数据是不是noi官网的，总之有些点和题目描述的不一样，比如说GET-MAX操作后面诡异的给我来了几个空格，我就一直在201和215之间徘徊。      依然是裸的splay，发现这次的时间有一些惊险，某个点我要8S才跑得出来，（ms时限有10S），开了release以后神速！每个点都不超过2S，很奇怪啊。
[NOI2005] 维护数列
still_foolish的博客
11-25 1196
[NOI2005] 维护数列时间限制：3 s 内存限制：256 MB 【问题描述】 请写一个程序，要求维护一个数列，支持以下 6 种操作：（请注意，格式栏中的下划线‘_’表示实际输入文件中的空格） 操作编号 输入文件中的格式 说明1.插入 INSERT_posi_tot_c1_c2_…_ctot 在当前数列的第 posi 个数字后插入 tot个数字：c1, c2, …, ctot；
在线评测系统大全-2021-10-5
"该文档是关于2021年10月5日的OJ（Online Judge）系统汇总，主要涵盖了多个在线编程评测平台的信息，包括适合不同水平的用户，特别是针对信奥（如NOI，CSP-J，CSP-S，NOIP）等竞赛的训练。文档内容包括各个OJ的链接...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值