YBTOJ洛谷P2042:维护数列(平衡树)

文章目录

  • 题目描述
  • 解析
    • 删除区间
    • 插入数列
    • 修改&翻转
    • 区间和&最大子段和
  • 代码

传送门

题目描述

在这里插入图片描述

解析

阴间题…
这不是裸的板子吗?
国赛真的有人能把这题写出来吗…
应该算一道练习作用很强的题了
写完这题,各种平衡树维护区间操作的方法可以说是毕业了吧
题解封装一个split函数的实现技巧值得借鉴
总的来说就一句:不要忘记pushdown和pushup!

删除区间

设删除区间为[l,r]
先把l-1结点splay到根
再把r+1转到l-1的下面
这样[l,r]一定就在r+1的右儿子了
直接删掉即可

插入数列

先在外面用类似线段树的方法建一棵完全平衡的平衡树,然后当单点接到需要的位置(寻找位置与删除类似)

修改&翻转

上标记!
把对应区间转出来,打一个标记即可

区间和&最大子段和

区间和比较无脑
最大子段和做一个经典的小白逛公园就行了

说起来很轻松,本题代码实现在传标记时还是很恶心的
总结一些易错的问题:

  1. 再说一遍:不要漏掉pushup、pushdown!
  2. 注意本题最大子段和可以是负数,所以答案不能和0取max
  3. same操作可以用一个tag记录是否需要下传比laz记录下传的值更加舒服
  4. 反转的时候最大左子段和最大右子段也要交换
  5. 本题卡空间,因此需要动态分配点编号,注意分配之前确定标记都整干净了!
  6. 插入时把一个结点接到父亲上之前要把父亲的标记下传!
  7. 任何操作几乎都要splay一下沿路更新信息

都是辛酸泪…
上代码吧
(实现已经比较简洁了算上debug还是干到了小200行…)

代码

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=2e6+100;
const int M=2e6+100;
const ll mod=1ll<<31;
ll read(){ll x=0,f=1;char c=getchar();while(!isdigit(c)){if(c=='-')f=-1;c=getchar();};while(isdigit(c)){x=x*10+c-'0';c=getchar();};return x*f;
}
int n,m,k;
int fre[N],top;
int val[N],siz[N],f[N],lp[N],rp[N],mx[N],sum[N],r,rev[N],tag[N],tr[N][2];
#define ls (tr[x][0])
#define rs (tr[x][1])
void pushup(int x){if(!x) return;siz[x]=siz[ls]+siz[rs]+1;sum[x]=sum[ls]+sum[rs]+val[x];if(!ls&&!rs){lp[x]=rp[x]=max(0,val[x]);mx[x]=val[x];}else if(!ls){lp[x]=max(0,val[x]+lp[rs]);rp[x]=max(rp[rs],sum[rs]+val[x]);mx[x]=max(mx[rs],val[x]+lp[rs]);}else if(!rs){rp[x]=max(0,val[x]+rp[ls]);lp[x]=max(lp[ls],sum[ls]+val[x]);mx[x]=max(mx[ls],val[x]+rp[ls]);}else{lp[x]=max(lp[ls],sum[ls]+val[x]+lp[rs]);rp[x]=max(rp[rs],sum[rs]+val[x]+rp[ls]);mx[x]=max(rp[ls]+val[x]+lp[rs],max(mx[ls],mx[rs]));}
}
void rrr(int x){rev[x]^=1;swap(ls,rs);swap(lp[x],rp[x]);
}
void pushdown(int x){if(tag[x]){tag[x]=rev[x]=0;if(ls){val[ls]=val[x];tag[ls]=1;sum[ls]=siz[ls]*val[x];}if(rs){val[rs]=val[x];tag[rs]=1;sum[rs]=siz[rs]*val[x];	}if(val[x]>=0){if(ls) mx[ls]=lp[ls]=rp[ls]=sum[ls];if(rs) mx[rs]=lp[rs]=rp[rs]=sum[rs];}else{if(ls){lp[ls]=rp[ls]=0;mx[ls]=sum[ls];}if(rs){lp[rs]=rp[rs]=0;mx[rs]=sum[rs];}}}if(rev[x]){rev[x]=0;if(ls) rrr(ls);if(rs) rrr(rs);}
}
void dfs1(int x){if(!x) return;pushdown(x);dfs1(ls);printf("%d ",val[x]);dfs1(rs);pushup(x);
}
void dfs2(int x){if(!x) return;pushdown(x);printf("x=%d val=%d ls=%d rs=%d siz=%d\n",x,val[x],ls,rs,siz[x]);dfs2(ls);dfs2(rs);pushup(x);
}
void debug(int f){printf("------r=%d\n",r);if(f==1) dfs1(r);else dfs2(r);printf("\n\n");
}
int New(int v,int fa){int x=fre[top--];val[x]=mx[x]=sum[x]=v;lp[x]=rp[x]=max(0,v);f[x]=fa;siz[x]=1;tag[x]=rev[x]=0;return x;
}
#define which(x) (tr[f[x]][1]==x)
void rotate(int x){int fa=f[x],gfa=f[fa];pushdown(gfa);pushdown(fa);pushdown(x);int k=which(x),son=tr[x][k^1];f[x]=gfa;if(gfa) tr[gfa][which(fa)]=x;f[fa]=x;tr[x][k^1]=fa;f[son]=fa;tr[fa][k]=son;pushup(fa);pushup(x);pushup(gfa);
}
void splay(int x,int goal){if(x==goal) return;for(int fa;(fa=f[x])!=goal;rotate(x)){if(f[fa]!=goal)which(x)==which(fa)?rotate(fa):rotate(x);}if(!goal) r=x;
}
int a[N];
int build(int l,int r,int fa){if(l>r) return 0;int mid=(l+r)>>1;int x=New(a[mid],fa);tr[x][0]=build(l,mid-1,x);tr[x][1]=build(mid+1,r,x);pushup(x);
//	printf("l=%d r=%d x=%d pl=%d siz=%d\n",l,r,x,mid,siz[x]);return x;
}
int find(int x,int k){pushdown(x);if(siz[ls]>=k) return find(ls,k);else if(siz[ls]+1==k) return x;else return find(rs,k-siz[ls]-1);
}
int split(int l,int rr){int a=find(r,l-1);splay(a,0);int b=find(r,rr+1);splay(b,a);return tr[b][0];
}
void add(int x,int k,int s){pushdown(x);tr[x][k]=s;f[s]=x;pushup(x);
}
char s[150];
void ins(){int pos=read(),tot=read();pos++;for(int i=1;i<=tot;i++) a[i]=read();int now=build(1,tot,0);if(pos==1){int x=find(r,1);splay(x,0);int y=find(r,2);splay(y,x);add(y,0,now);splay(y,0);return;}int x=split(pos,pos);add(x,1,now);pushup(x);splay(now,0);
}
void save(int x){if(!x) return;fre[++top]=x;save(ls);save(rs);
}
void del(){int pos=read(),tot=read();pos++;int x=split(pos,pos+tot-1),fa=f[x];save(x);tr[fa][0]=0;pushup(fa);splay(fa,0);
}
void same(){int pos=read(),tot=read(),v=read();pos++;int x=split(pos,pos+tot-1);val[x]=v;tag[x]=1;sum[x]=siz[x]*v;lp[x]=rp[x]=max(0,sum[x]);splay(x,0);
}
void reverse(){int pos=read(),tot=read();pos++;int x=split(pos,pos+tot-1);rrr(x);//printf("x=%d\n",x);debug(2);splay(x,0);
}
void getsum(){int pos=read(),tot=read();pos++;int x=split(pos,pos+tot-1);printf("%d\n",sum[x]);
}
void maxsum(){printf("%d\n",mx[r]);
}
int main(){for(int i=500050;i>=1;i--) fre[++top]=i;n=read();m=read();for(int i=1;i<=n;i++) a[i]=read();r=build(1,n,0);splay(find(r,1),0);add(r,0,New(-20000,r));splay(find(r,n+1),0);add(r,1,New(-20000,r));//debug(1);for(int i=1;i<=m;i++){//	printf("-------------------------------------\n");scanf(" %s",s+1);if(s[1]=='I') ins();else if(s[1]=='D') del();else if(s[3]=='K') same();else if(s[1]=='R') reverse();else if(s[1]=='G') getsum();else maxsum();//debug(1);}return 0;
}
/*
*/

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

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

相关文章

CAP 2.4版本发布,支持版本隔离特性

前言自从上次 CAP 2.3 版本发布 以来&#xff0c;已经过去了几个月的时间&#xff0c;这几个月比较忙&#xff0c;所以也没有怎么写博客&#xff0c;趁着2019年到来之际&#xff08;现在应该是2019年开始的时候&#xff09;&#xff0c;CAP也发布了2018年的最后一个大版本 2.4&…

AT2005-[AGC003E]Sequential operations on Sequence【差分,思维】

正题 题目链接:https://www.luogu.com.cn/problem/AT2005 题目大意 开始有一个1∼n1\sim n1∼n依次排列的序列&#xff0c;然后QQQ次&#xff0c;第iii次把序列长度变为aia_iai​&#xff0c;不足的从前往后循环填充。 求最后每个数字的出现次数。 1≤n,q≤105,1≤ai≤10181…

多项式工业基础与全家桶

多项式工业基础与全家桶 开坑待填&#xff0c;放个常数巨大的板子先 别忘了这道题&#xff01;P3338 [ZJOI2014]力 #define Maxn 200005 #define mod 998244353 inline int ksm(int x,int ymod-2) {int ret1;for(;y;y>>1,x1ll*x*x%mod) if(y&1) ret1ll*ret*x%mod;ret…

【周末狂欢赛7】【NOIP模拟赛】七夕祭,齿轮(dfs),天才黑客

文章目录T1题目题解codeT2题目题解codeT3题目题解codeT1 题目 七夕节因牛郎织女的传说而被扣上了「情人节」的帽子。于是TYVJ今年举办了一次线下七夕祭。Vani同学今年成功邀请到了cl同学陪他来共度七夕&#xff0c;于是他们决定去TYVJ七夕祭游玩。 TYVJ七夕祭和11区的夏祭的…

股票买卖 IV

股票买卖 IV 题意&#xff1a; 给定一个长度为 N 的数组&#xff0c;数组中的第 i 个数字表示一个给定股票在第 i 天的价格。 设计一个算法来计算你所能获取的最大利润&#xff0c;你最多可以完成 k 笔交易。 注意&#xff1a;你不能同时参与多笔交易&#xff08;你必须在再…

9.25 模拟

文章目录前言考场总结前言 260pts 100601000 线段树数组要乘4&#xff01;&#xff01;&#xff01; 结果想了一个多小时和暴力没差了。。。 呜呜呜 但总体还可以啦 T3玄学暴力竟然过掉了 告诉我们卡时的重要性 考场 先看题 T1原题啊啊啊&#xff01; 真就签到题了 感觉T2好…

Silence 主题配置代码

页面定制 import url(https://fastly.jsdelivr.net/gh/esofar/cnblogs-theme-silence3.0.0-beta3/dist/silence.min.css);media only screen and (max-width:1365px){.github-corner{display:none}}media only screen and (max-width:1259px){#home{width:100%!important}}medi…

AT3857-[AGC020C]Median Sum【背包,bitset】

正题 题目链接:https://www.luogu.com.cn/problem/AT3857 题目大意 给出nnn个数字的一个序列aaa&#xff0c;求它的所有非空子集的和的中位数。 1≤n,ai≤20001\leq n,a_i\leq 20001≤n,ai​≤2000 解题思路 考虑到假设所有数的和为SSS&#xff0c;一个集合的和为xxx&#x…

.NET Core 如何为项目提供高性能解决方案?

本系列&#xff0c;我们将探讨.NET Core 的一些好处&#xff0c;以及它如何为市场提供高性能解决方案&#xff0c;为传统.NET 开发人员和技术人员提供帮助。正文前言随着.NET Core 2.0 在 2016 年首次发布&#xff0c;微软拥有了这个通用、模块化、跨平台开源项目的下一个主要版…

1058. 股票买卖 V

1058. 股票买卖 V 题意&#xff1a; 给定一个长度为 N 的数组&#xff0c;数组中的第 i 个数字表示一个给定股票在第 i 天的价格。 设计一个算法计算出最大利润。在满足以下约束条件下&#xff0c;你可以尽可能地完成更多的交易&#xff08;多次买卖一支股票&#xff09;: …

[2.9训练]【CF909C】Python Indentation,【CF909D】Colorful Points,【CF909E】Coprocessor

文章目录T1&#xff1a;Python Indentation题目题解codeT2&#xff1a;Colorful Points题目题解codeT3&#xff1a;Coprocessor题目题解codeT1&#xff1a;Python Indentation 题目 题目描述 In Python, code blocks don’t have explicit begin/end or curly braces to mark…

模板:2-SAT问题

文章目录前言实现代码所谓2-SAT&#xff0c;就是解决两个SAT的问题 &#xff08;逃&#xff09; 前言 SAT 是适定性&#xff08;Satisfiability&#xff09;问题的简称。一般形式为 k - 适定性问题&#xff0c;简称 k-SAT。而当 k>2 时该问题为 NP 完全的。所以我们只研究 …

CF1481F-AB Tree【构造,背包】

正题 https://www.luogu.com.cn/problem/CF1481F 题目大意 给出nnn个点的一棵树&#xff0c;在每个节点上填a/ba/ba/b&#xff0c;要求恰好有mmm个aaa。要求每个节点到根路径上的字符串种类最少&#xff0c;输出方案。 1≤m≤n≤1051\leq m\leq n\leq 10^51≤m≤n≤105 解题思…

Docker最全教程之使用Tencent Hub来完成CI(十)

本周更新两篇&#xff0c;保证不太监&#xff01;在本系列教程中&#xff0c;笔者希望将必要的知识点围绕理论、流程&#xff08;工作流程&#xff09;、方法、实践来进行讲解&#xff0c;而不是单纯的为讲解知识点而进行讲解。也就是说&#xff0c;笔者希望能够让大家将理论、…

P4169 [Violet]天使玩偶/SJY摆棋子(CDQ分治+欧几里得距离)

P4169 [Violet]天使玩偶/SJY摆棋子(CDQ分治欧几里得距离) 记得上一次欧几里得距离的转化是CF1093G Multidimensional Queries&#xff0c;我们使用了点对在四种方向分别考虑并用 \(\max\) 合并的方法解决&#xff0c;现在使用一种类似的方法。 \(\bigstar\texttt{Trick}\)&…

[2.7]【CF933A】A Twisty Movement【CF926B】Add Points【CF917A】The Monster【CF919E】Congruence Equation

文章目录T1&#xff1a;A Twisty Movement题目题解codeT2&#xff1a;Add Points题目题解codeT3&#xff1a;The Monster题目题解codeT4&#xff1a;Congruence Equation题目题解codeT1&#xff1a;A Twisty Movement 题目 题目 题解 因为aia_iai​1/21/21/2&#xff0c;于…

LIS最长上升子序列

LIS算是比较经典的问题&#xff0c;常用的是O(n^2)的方法 for(int i1;i<n;i){dp[i]1;for(int j1;j<i;j){if(a[j]<a[i])dp[i]max(dp[i],dp[j]1);}mxmax(mx,dp[i]);}我们这里优化成O(nlogn) 我们模拟一个栈stack&#xff0c;每读入一个数&#xff0c;如果这个数大于栈顶…

AT5160-[AGC037C]Numbers on a Circle【贪心,堆】

正题 题目链接:https://www.luogu.com.cn/problem/AT5160 题目大意 给出两个长度为nnn的环序列aaa和bbb&#xff0c;每次你可以让aaa中的一个数变为它和相邻两个的和。 求最少的步数将aaa变为bbb。 1≤n≤105,1≤ai,bi≤1091\leq n\leq 10^5,1\leq a_i,b_i\leq 10^91≤n≤10…

YBTOJ洛谷P3209:平面图判定(2-SAT)

文章目录解析代码传送门解析 关键性质是一个定理&#xff1a;若m>3*n-6&#xff0c;必然不存在合法的平面图 这谁知道啊 不过这题应该往也许图过于稠密时必然无解这方面想 所以我们只需要考虑m、n同阶的情况就行了 这个时候我们直接暴力判断跑2-SAT就行了 代码 #include&…

线性代数 - 矩阵对角化

矩阵对角化 今天听 \(\texttt{m}\color{red}\texttt{yee}\) 嘴的&#xff0c;赶紧来补个学习笔记。 我们有点时候需要计算一个较小矩阵的 \(n\) 次幂&#xff0c;但直接求幂非常不方便&#xff0c;这是会考虑矩阵对角化&#xff0c;将 \(M\) 改写为 \(\mathcal{PDP^{-1}}\)&…