[Ynoi2015]纵使日薄西山
题目大意:
给定一个序列,每次单点修改,然后进行询问。
定义一次操作为,选择一个位置$x$,将这个位置的数和左边、右边两个位置的数(不存在则忽略)各减去1,然后和0取max。
对序列中最大的位置进行一次操作(相同则取最前面的),不断重复,直到所有位置为0为止。
问执行了多少次操作。
询问互相独立(即下一次询问的序列并不是全0)。
解题思路:
在太阳西斜的这个世界里,置身天上之森。等这场战争结束之后,不归之人与望眼欲穿的众人, 人人本着正义之名,长存不灭的过去、逐渐消逝的未来。我回来了,纵使日薄西山,即便看不到未来,此时此刻的光辉,盼君勿忘。————世界上最幸福的女孩
我的愿望是,和珂朵莉一样可爱。
淦不动辣QaQ就窝一个在那瞎分块
一开始的做法是先分块,求出块内的答案,然后瞎分类讨论(自闭了QAQ)。
实际上这个性质挺优美的为啥窝就是发现不了啊
显然一次对某个位置操作到底,和答案是一样的(当你选择这个数进行操作后,其左边和右边就永远不会被操作。而其他数对这个位置无影响)。
所以相当于求所有会被操作的位置上的数的和。
考虑一个单调递增或单调递减的序列,在这上面选的数一定是一个选一个不选(奇偶性相同)。
所以我们对奇数、偶数位置分别用树状数组维护一下,然后用set记录序列的极大、极小值的位置。则可以快速计算出两个极值之间的区间的贡献。
然后对于单点修改,对受到影响的部分区间重新计算贡献即可。
细节挺多的,比如对极小值是否取到的考虑(仅当其左右两个极大值的位置的奇偶性和这个极小值位置相同时,才能取到这个值)。
时间复杂度$O(n\log n)$。
C++ Code:
#include<cstdio>
#include<set>
const int N=1e5+6;
typedef long long LL;
int n,a[N];
struct BIT{
LL b[N];
inline void add(int i,int x){for(;i<N;i+=i&-i)b[i]+=x;}
inline LL ask(int i){LL x=0;for(;i;i^=i&-i)x+=b[i];return x;}
}odd,even;
std::set<int>s;
LL ans=0;
LL get(int x){
if(x<1||x>n)return 0;
auto it=s.find(x);
if(*it==1){
auto nxt=it;++nxt;
if(a[*it]<a[*nxt]){
if(*nxt&1)return odd.ask(*nxt-1)-odd.ask(*it-1);else
return even.ask(*nxt-1)-even.ask(*it-1);
}else return 0;
}
if(*it==n){
auto pre=it;--pre;
if(a[*pre]>=a[*it]){
if(*pre&1)return odd.ask(*it)-odd.ask(*pre-1);else
return even.ask(*it)-even.ask(*pre-1);
}else return a[n];
}
auto pre=it,nxt=it;--pre,++nxt;
if(a[*pre]>=a[*it]&&a[*it]<a[*nxt]){
LL ret=0;
if(*pre&1)ret+=odd.ask(*it-1)-odd.ask(*pre-1);else
ret+=even.ask(*it-1)-even.ask(*pre-1);
if(*nxt&1)ret+=odd.ask(*nxt-1)-odd.ask(*it);else
ret+=even.ask(*nxt-1)-even.ask(*it);
if((*pre&1)==(*it&1)&&(*nxt&1)==(*it&1))ret+=a[*it];
return ret;
}else return 0;
}
int main(){
scanf("%d",&n);
for(int i=1;i<=n;++i)
scanf("%d",a+i),((i&1)?odd:even).add(i,a[i]);
s.insert(-2),s.insert(-1),s.insert(n+2),s.insert(n+3);
s.insert(1);
for(int i=2;i<n;++i)
if(a[i-1]<a[i]&&a[i]>=a[i+1]||a[i-1]>=a[i]&&a[i]<a[i+1])s.insert(i);
s.insert(n);
for(int i:s)
ans+=get(i);
int q;
for(scanf("%d",&q);q--;){
int x,y;
scanf("%d%d",&x,&y);
int nxt=*s.upper_bound(x),pre=*--s.lower_bound(x);
ans-=get(nxt),ans-=get(pre);
nxt=*s.upper_bound(nxt),pre=*--s.lower_bound(pre);
ans-=get(nxt),ans-=get(pre);
if(s.count(x))ans-=get(x),s.erase(x);
s.erase(x-1),s.erase(x+1);
if(x&1)odd.add(x,-a[x]),odd.add(x,y);else
even.add(x,-a[x]),even.add(x,y);
a[x]=y;
for(int i=x-1;i<=x+1;++i){
if(i<1||i>n)continue;
if(i==1||i==n)s.insert(i);else
if(a[i-1]<a[i]&&a[i]>=a[i+1]||a[i-1]>=a[i]&&a[i]<a[i+1])s.insert(i);
}
for(auto l=s.find(pre),r=s.upper_bound(nxt);l!=r;++l)
ans+=get(*l);
printf("%lld\n",ans);
}
return 0;
}
[Ynoi2015]纵使日薄西山的更多相关文章
- 洛谷P5069 [Ynoi2015]纵使日薄西山(树状数组,set)
洛谷题目传送门 一血祭 向dllxl致敬! 算是YNOI中比较清新的吧,毕竟代码只有1.25k. 首先我们对着题意模拟,寻找一些思路. 每次选了一个最大的数后,它和它周围两个数都要减一.这样无论如何, ...
- 【树状数组】【P5069】[Ynoi2015]纵使日薄西山
Description 给定一个长度为 \(n\) 的非负整数序列 \(\{a_n\}\),\(q\) 次操作,每次要么单点修改序列某个值,要么查询整个序列需要操作多少次才能变成全 \(0\). 一次 ...
- [Ynoi2015]即便看不到未来
题目大意: 给定一个序列,每次询问,给出一个区间$[l,r]$. 设将区间内的元素去重后重排的数组为$p$,求$p$中长度为$1\sim 10$的极长值域连续段个数. 长度为$L$的极长值域连续段的定 ...
- [Ynoi2015]盼君勿忘
题目大意: 给定一个序列,每次查询一个区间\([l,r]\)中所有子序列分别去重后的和\(\bmod p\)(每次询问模数不同). 解题思路: 在太阳西斜的这个世界里,置身天上之森.等这场战争结束之后 ...
- [Ynoi2015]我回来了
题目大意: 给定一张无向无权图,每次给定若干个二元组\((x_i,y_i)\),定义点\(u\)满足条件,当且仅当存在\(i\),并满足\(dist(u,x_i)\leqslant y_i\)(\(d ...
- [Ynoi2015]此时此刻的光辉
题目大意: 给定一个序列,每次询问一段区间的数的乘积的约数个数. 解题思路: 在太阳西斜的这个世界里,置身天上之森.等这场战争结束之后,不归之人与望眼欲穿的众人, 人人本着正义之名,长存不灭的过去.逐 ...
- 「Ynoi2015」我回来了
「Ynoi2015」我回来了 这东西已经不是 Ynoi 了,因为太水被嫌弃了. 如何提升自己的数据结构能力?从Ynoi做起 题目链接 其实这个题很小清新的辣,而且不卡常. 由于边权为 \(1\),所以 ...
- 洛谷P5072 [Ynoi2015]盼君勿忘 [莫队]
传送门 辣鸡卡常题目浪费我一下午-- 思路 显然是一道莫队. 假设区间长度为\(len\),\(x\)的出现次数为\(k\),那么\(x\)的贡献就是\(x(2^{len-k}(2^k-1))\),即 ...
- 【题解】Luogu P5071 [Ynoi2015]此时此刻的光辉
众所周知lxl是个毒瘤,Ynoi道道都是神仙题,题面好评 原题传送门 一看这题没有修改操作就知道这是莫队题(我也只会莫队) 我博客里对莫队的简单介绍 一个数N可以分解成\(p_1^{c_1}p_2^{ ...
随机推荐
- 【iOS】网络载入图片缓存与SDWebImage
载入网络图片能够说是网络应用中必备的.假设单纯的去下载图片,而不去做多线程.缓存等技术去优化,载入图片时的效果与用户体验就会非常差. 一.自己实现载入图片的方法 tips: *iOS中全部网络訪问都是 ...
- jsonArray和Java List对象互转,日期处理
List转jsonArray : // 格式化日期 JsonConfig jsonConfig = new JsonConfig(); DSHJsonDateValueProcessor dshJso ...
- luogu2161 [SHOI2009]会场预约
题目大意 随着时间的推移这里有几个任务对应着一段区间.每次要将任务安到时间线上时,要把时间线上已有的与该任务对应区间有交集的区间对应的任务删去.求每次删去的区间个数,以及整个时间线上有几个任务.时间线 ...
- maven工程读取resource下配置文件
maven工程读取resource下配置文件 在maven工程中,我们会将配置文件放到,src/main/resources 下面,例如 我们需要确认resource 下的文件 编译之后存放的位置 ...
- [Contest Hunter#17-C] 舞动的夜晚
[题目链接] http://contest-hunter.org:83/contest/CH%20Round%20%2317/%E8%88%9E%E5%8A%A8%E7%9A%84%E5%A4%9C% ...
- Tool:ProcessOn
ylbtech-Tool:ProcessOn ProcessOn是一个面向垂直专业领域的作图工具和社交网络,成立于2011年6月并于2012年启动.ProcessOn将全球的专家顾问.咨询机构.BPM ...
- php获得本机ipv4地址
if (isset($_ENV["HOSTNAME"])) $MachineName = $_ENV["HOSTNAME"]; else if (isset($ ...
- AOP实现参数的判空问题
不想每次都去判断必传的参数是否为空,写代码太繁琐了,正好最近用了AOP实现权限控制,依葫芦画瓢,现在用它实现参数的判空,至于AOP的原理之类,自己百度了解一下吧 1. NullDisable注解 @D ...
- Android开发利器之ActivityTracker
版权声明:本文为xing_star原创文章,转载请注明出处! 本文同步自http://javaexception.com/archives/113 Android开发利器之ActivityTracke ...
- HDU 4901 DP
我觉得这个DP挺难的...然而这只是lydrainbowcat学长幻灯片上的第一题-- 明天考试要GG. 题意: 给你一个序列,让你选出两个集合S和T.保证S里的数都在T里的数的左边.求一共有多少个集 ...